Matlab解微分方程组

《Matlab解微分方程组》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、求解微分方程可以按如下程序得到数值解clear;clc;[t,x]=ode45('xprim1',[01],1);plot(t,x,'o-');xlabel('timet0=0,tt=1');ylabel('xvaluesx(0)=1');functionxprim=xprim1(t,x)xprim=-x.^2;输出求解微分方程组可以按如下程序得到数值解clear;clc;[t,x]=ode45('xprim1',[020],[30;20]);plot(t,x);xlabel('timet0=0,tt=20');ylabel('xvaluesx1(0)=30,x2(0)=20')

2、;functionxprim=xprim1(t,x)xprim=[x(1)-0.1x(1)x(2)+0.01t;-x(2)+0.02x(1)x(2)+0.04t];输出求的通解输入命令dsolve('Du=1+u^2','t')输出i-itan(C3+t)求微分方程的特解输入命令y=dsolve('D2y+4Dy+29y=0','y(0)=0,Dy(0)=15','x')输出y=(3sin(5x))/exp(2x)求微分方程组通解输入命令[x,y,z]=dsolve('Dx=2x-3y+3z','Dy=4x-5y+3z','Dz=4x-4y+2z','t');x=simplify

3、(x)y=simplify(y)z=simplify(z)注意函数dsolve只可用于求解线性常微分方程（组）