Matlab解微分方程组Matlab解微分方程组
求解微分方程
2,xx',, ,x(0)1,,
可以按如下程序得到数值解
clear;clc;
[t,x]=ode45('xprim1',[0 1],1); plot(t,x,'o-');
xlabel('time t0=0,tt=1'); ylabel('x values x(0)=1'); function xprim=xprim1(t,x) xprim=-x.^2;
输出
求解微分方程组
xxxxt'0.10.01,,,,1112,xxxxt'0.020.04,,,...
Matlab解微分方程组
求解微分方程
2,xx',, ,x(0)1,,
可以按如下程序得到数值解
clear;clc;
[t,x]=ode45('xprim1',[0 1],1); plot(t,x,'o-');
xlabel('time t0=0,tt=1'); ylabel('x values x(0)=1'); function xprim=xprim1(t,x) xprim=-x.^2;
输出
求解微分方程组
xxxxt'0.10.01,,,,1112,xxxxt'0.020.04,,,,,2212 ,x(0)30,1,
,x(0)20,,2
可以按如下程序得到数值解
clear;clc;
[t,x]=ode45('xprim1',[0 20],[30;20]); plot(t,x);
xlabel('time t0=0,tt=20'); ylabel('x values x1(0)=30,x2(0)=20'); function xprim=xprim1(t,x) xprim=[x(1)-0.1*x(1)*x(2)+0.01*t;
-x(2)+0.02*x(1)*x(2)+0.04*t];
输出
du2求的通解 ,,1udt
输入命令
dsolve('Du=1+u^2','t')
输出
i
-i
tan(C3 + t)
求微分方程的特解
2,dydy,,,4290y,2 dxdx,
,yy(0)0,'(0)15,,,
输入命令
y=dsolve('D2y+4*Dy+29*y=0','y(0)=0,Dy(0)=15','x')
输出
y =(3*sin(5*x))/exp(2*x)
dx,,,,233xyz,dt,dy,,,,453xyz求微分方程组通解 ,dt,
dz,,,,442xyz,dt,
输入命令
[x,y,z]=dsolve('Dx=2*x-3*y+3*z','Dy=4*x-5*y+3*z','Dz=4*x-4*y+2*z','t'
);
x=simplify(x)
y=simplify(y)
z=simplify(z)
注意函数dsolve只可用于求解线性常微分方程(组)
本文档为【Matlab解微分方程组】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑,
图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。