精选优质文档-倾情为你奉上常微分方程组求解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=
Copyright © 2018-2021 Wenke99.com All rights reserved
工信部备案号:浙ICP备20026746号-2
公安局备案号:浙公网安备33038302330469号
本站为C2C交文档易平台,即用户上传的文档直接卖给下载用户,本站只是网络服务中间平台,所有原创文档下载所得归上传人所有,若您发现上传作品侵犯了您的权利,请立刻联系网站客服并提供证据,平台将在3个工作日内予以改正。