常微分方程组求解matlab程序代码(共2页).docx

精选优质文档-倾情为你奉上常微分方程组求解matlab程序代码clc,clear,close allCT=1;1;1;h=0.25;k,X,Y,wucha,P=RK4z(dydx,0.1,60,CT,h),plot(X,Y(:,1),g-,X,Y(:,2),b*-,X,Y(:,3),mp-)xlabel(轴it x); ylabel(轴it y)grid onlegend(方程解z1的曲线,方程解z2的曲线, 方程解z3的曲线)function dY=dydx(X,Y)dY(1)=Y(2);dY(2)=Y(3);dY(3)=Y(3)*X(-1)-3*X(-2)*Y(2)+2*X(-3)*Y(1)+9*X3*sin(X);dY=dY(1);dY(2);dY(3);endfunction k,X,Y,wucha,P=RK4z(dydx,a,b,CT,h)n=fix(b-a)/h);X=zeros(n+1,1); Y=