龙格库塔法解二阶微分方程所解方程为 y’’= -y
主程序 M文件
a=0;
b=8*pi;
n=400;
h=(b-a)/n;
xa=3;
ya=0;
T=a:h:b;
X=zeros(1,n+1);
Y=zeros(1,n+1);
X(1)=xa;
Y(1)=ya;
for i=1:n
f1=f(T(i),X(i),Y(i));
g1=g(T(i),X(i),Y(i));
f2=f(T(i)+h/2,X(i)+f1*h/2,Y(i)+g1*h/2);
g2=g(T(i)+h/2,X(i)+f1*...
所解方程为 y’’= -y
主程序 M文件
a=0;
b=8*pi;
n=400;
h=(b-a)/n;
xa=3;
ya=0;
T=a:h:b;
X=zeros(1,n+1);
Y=zeros(1,n+1);
X(1)=xa;
Y(1)=ya;
for i=1:n
f1=f(T(i),X(i),Y(i));
g1=g(T(i),X(i),Y(i));
f2=f(T(i)+h/2,X(i)+f1*h/2,Y(i)+g1*h/2);
g2=g(T(i)+h/2,X(i)+f1*h/2,Y(i)+g1*h/2);
f3=f(T(i)+h/2,X(i)+f2*h/2,Y(i)+g2*h/2);
g3=g(T(i)+h/2,X(i)+f2*h/2,Y(i)+g2*h/2);
f4=f(T(i)+h,X(i)+f3*h,Y(i)+g3*h);
g4=g(T(i)+h,X(i)+f3*h,Y(i)+g3*h);
X(i+1)=X(i)+(f1+2*f2+2*f3+f4)*h/6;
Y(i+1)=Y(i)+(g1+2*g2+2*g3+g4)*h/6;%四阶龙格-库塔方法
end
R1=[T',X',Y']
函数M文件
f.m
function y=f(x,y,z)
y=z;
g.m
function y=g(x,y,z)
y=-y;
本文档为【龙格库塔法解二阶微分方程】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑,
图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。