找回密码
 注册
关于网站域名变更的通知
查看: 641|回复: 1
打印 上一主题 下一主题

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

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-3-10 10:36 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
本帖最后由 House 于 2020-3-10 17:43 编辑
- \% J/ ]( S& s& L
0 X9 g7 A  `4 }& [3 |MATLAB使用自定义的欧拉法求解常微分方程组' {8 a  S6 Y# g* B
! o& D: I5 C# z" W$ v: ?1 G4 y% b
0 w8 r8 B9 B2 W% i3 v

  q/ f$ n2 K6 U' Z1 Bfunction F=f(t,Y)
+ J; ~3 n# L1 ~- h& v% 定义待求的微分方程组  t- z) u4 v' K" g5 X9 R
x=Y(1);
" f4 z$ b( i0 ?: v# p* Z3 H8 my=Y(2);
0 \0 r# ?( ?2 [) |: g$ g* H( y. Mf1=3*y;
2 V- h. G( ?9 af2=(1-x^2)*y-x;
- }9 ~7 r. T2 I8 WF=[f1;f2];! V# A) J3 ?6 I4 V) _
end: d8 C  X% P6 Q& `* M
( r2 C1 m$ F4 g7 r9 n3 m
4 W. x' W) F: H% o
%% 定义计算的步长, 设置变量的初始值3 ~, \& P& q  B5 L8 m2 R
clear;clc;close all! H+ ~( O. C3 C5 ]# Q
Delta=0.001;      % 定义步长
' ]) g: H  E0 V0 i. H( Kt=0: Delta :20;     % 定义自变量 t) M) |% u1 V3 ]7 t7 ^
n=length(t);
- j3 ]3 @! d; J2 j9 WY(:,1)=[2;0];     % 定义 x y 的初始值$ ^, B: ?  J4 C9 p1 c9 L
" F/ F1 L( ]/ Z
%% 自定义欧拉法, 求解微分方程组# F  |& L2 e, Y2 j- Z4 B/ A! E
for k=1:n-1
. O% `) V8 J$ Y    Y(:,k+1)=Y(:,k)+Delta*f(t(k),Y(:,k));
/ E& b% _4 r. D; C- qend! Y) x. r& N! b8 C$ }+ _
x=Y(1,: );9 q% V, _( j. V0 N# G
y=Y(2,: );) E' d1 w5 ^- T; N8 N
' r8 w3 l# e9 m' b9 e
%% 绘制 x y 的求解结果
2 v" R7 b. S4 r' R0 I  F) sfigure9 g" Y" g: c" x+ \6 @" x3 _4 A
set(gcf,'units','normalized','position',[0.15 0.2 0.7 0.6]);  % 设置 figure 窗口的位置和尺寸
; `4 ?+ r) Y2 A3 s  Qsubplot(1,2,1)1 L) E, q$ K5 `/ A" x" Y, U1 M
plot(t,x,'b')( ^6 ?: U3 ~  X
xlabel('t')
8 w. b$ C( J  F' k3 y* a* K+ z+ Mylabel('x'), S, I8 q/ g5 ~1 w0 {, _; }

* x. ]" P3 {1 V5 F+ zsubplot(1,2,2)+ x# t- R9 a: f) A2 R
plot(t,y,'r')$ d) _& S: f9 S, F8 A% _& c
xlabel('t')- P4 d! @2 o! f' ^
ylabel('y')) x1 W5 }4 S; f! O: l
% Z% C1 V: S4 q) U8 s- R$ R+ }# v  e
. Y1 m0 R* K2 K, B: ~

2 r( ~+ `6 K4 a( a

该用户从未签到

2#
发表于 2020-3-10 17:43 | 只看该作者
谢谢楼主分享的代码。
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-10-6 01:12 , Processed in 0.109375 second(s), 26 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表