原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
6．1
一、题目：
用欧拉公式和四阶龙格-库塔法分别求解下列初值问题；
二、算法分析：
Euler公式为：
四阶龙格-库塔法为：
这里我们设h=0.1
三、程序设计及运行结果：
（1）f11.m
global y;%定义y为数组
global x;
a=0;b=1;h=0.1;alpha=1;f=inline(0.9*y/(1+2*x));
x=a:h:b;
n=(b-a)/h;
y(1)=alpha;
disp(euler公式求得的解为：)
for i=2:n+1
y(i)=y(i-1)+h*feval(f,x(i-1),y(i-1));
sprintf(x=%3.1f,y=%9.7f,x(i),y(i))
end
disp(四阶龙格-库塔法求得的解为：)
for i=2:n+1
k1=feval(f,x(i-1),y(i-1));
k2=feval(f,x(i-1)+h/2,y(i-1)+k1*h/2);
k3=feval(f,x(i-1)+h/2,y(i-1)+k2*h/2);
k4=feval(f,x(i-1)+h,y(i-1)+k3*h);
y(i)=y(i-1)+h/6*(k1+2*k2+2*k3+k4);
sprintf(x=%3.1f,y=%9.7f,x(i),y(i))
end
运行结果为：euler公式求得的解为：
x=0.1,y=1.0900000
x=0.2,y=1.1717500
x=0.3,y=1.2470768
x=0.4,y=1.3172249
x=0.5,y=1.3830861
x=0.6,y=1.4453250
x=0.7,y=1.5044519
x=0.8,y=1.5608688
x=0.9,y=1.6148989
x=1.0,y=1.6668064
四阶龙格-库塔法求得的解为：
x=0.1,y=1.0855051
x=0.2,y=1.1634777
x=0.3,y=1.2355334
x=0.4,y=1.3027862
x=0.5,y=1.3660420
x=0.6,y=1.4259056
x=0.7,y=1.4828446
x=0.8,y=1.5372290
x=0.9,y=1.5893579
x=1.0,y=1.6394763
x
euler公式求得的解
四阶龙格-库塔法求得的解
x=0.1
y=1.0900000
y=1.0855051
x=0.2
y=1.1717500
y=1.1634777
x=0.3
y=1.2470768
y=1.2355334
x=0.4
y=1.3172249
y=1.3027862
x=0.5
y=1.3830861
y=1.3660420
x=0.6
y=1.4453250
y=1.4259056
x=0.7
y=1.5044519
y=1.4828446
x=0.8
y=1.5608688
y=1.5372290
x=0.9
y=1.6148989
y=1.5893579
x=1.0
y=1.6668064
y=1.6394763
（2）f12.m
global y;%定义y为数组
global x;
a=0;b=1;h=0.1;alpha=2;f=inline(-x*y/(1+x^2));
x=a:h:b;
n=(b-a)/h;
y(1)=alpha;
disp(euler公式求得的解为：)
for i=2:n+1
y(i)=y(i-1)+h*feval(f,x(i-1),y(i-1));
sprintf(x=%3.1f,y=%9.7f,x(i),y(i))
end
disp(四阶龙格-库塔法求得的解为：)
for i=2:n+1
k1=feval(f,x(i-1),y(i-1));
k2=feval(f,x(i-1)+h/2,y(i-1)+k1*h/2);
k3=feval(f,x(i-1)+h/2,y(i-1)+k2*h/2);
k4=feval(f,x(i-1)+h,y(i-1)+k3*h);
y(i)=y(i-1)+h/6*(k1+2*k2+2*k3+k4);
sprintf(x=%3.1f,y=%9.7f,x(i),y(i))
end
运行结果为：
euler公式求得的解为：
x=0.1,y=2.0000000
x=0.2,y=1.9801980
x=0.3,y=1.9421173
x=0.4,y=1.8886645
x=0.5,y=1.8235382
x=0.6,y=1.7505966
x=0.7,y=1.6733644
x=0.8,y=1.5947500
x=0.9,y=1.5169573
x=1.
文档评论（0）