MATLAB使用自定义的欧拉法求解常微分方程组

House

7 n( J; W' V( }5 D! F" b: X) o& J
. d6 E7 c/ r- o* a5 R+ B7 OMATLAB使用自定义的欧拉法求解常微分方程组
7 P: D1 }6 s- k5 i0 F& N

% D( T( |5 k  j" z
) n+ a5 T6 E: s# ]3 m0 H( T: wfunction F=f(t,Y)
9 p* I" Z" q- y: U/ t5 j5 o' I% 定义待求的微分方程组
x=Y(1);
y=Y(2);
f1=3*y;
; v* M; K4 i, Y* {$ If2=(1-x^2)*y-x;
F=[f1;f2];
end
# t6 t. Z7 u9 Q
3 h& y! e( t4 {) `/ g
%% 定义计算的步长, 设置变量的初始值
" J; A4 S4 M' t2 Yclear;clc;close all
Delta=0.001;      % 定义步长
$ q+ m% |) e! p0 r. u$ x7 \t=0: Delta :20;     % 定义自变量 t
n=length(t);
Y(:,1)=[2;0];     % 定义 x y 的初始值

%% 自定义欧拉法, 求解微分方程组
for k=1:n-1
    Y(:,k+1)=Y(:,k)+Delta*f(t(k),Y(:,k));
end
- n% U( {* L0 t5 rx=Y(1,: );
8 Q4 Y+ R3 _) w8 @$ xy=Y(2,: );
4 h0 h7 x7 g9 l/ K/ J! \
%% 绘制 x y 的求解结果
figure
set(gcf,'units','normalized','position',[0.15 0.2 0.7 0.6]);  % 设置 figure 窗口的位置和尺寸
! I7 ]' z0 W" I3 [( Dsubplot(1,2,1)
# H1 ^3 T, `/ k  A( a) I3 J/ ?plot(t,x,'b')
/ U9 c" z: h. k& U' xxlabel('t')
ylabel('x')
( }2 M5 R; g- R0 s9 K, b- s! n
$ Y3 Z5 C3 n: h; G) K! ssubplot(1,2,2)
plot(t,y,'r')
xlabel('t')
* B8 g& ^$ z  d$ Kylabel('y')
* K8 S% z- {4 D8 b! y+ k2 c

wu68aq

谢谢楼主分享的代码。