資源描述:
《Matlab解微分方程組》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、求解微分方程可以按如下程序得到數(shù)值解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;輸出求解微分方程組可以按如下程序得到數(shù)值解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)注意函數(shù)dsolve只可用于求解線性常微分方程(組)