刘伟201120370122.doc

下载文档

2
0
约1.73千字
约 6页
2016-10-21 发布于重庆
举报
版权申诉
保障服务

刘伟201120370122.doc

1、本文档共6页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

刘伟201120370122

用改进的欧拉公式，求解常微分方程的初值问题：用四阶龙格—库塔法求解微分方程的初值问题改进的欧拉方法%改进的欧拉方法 clear all,close all,clc format long h=0.01; a=0.1; b=0.4; N=(b-a)/h; N=round(N); y=zeros(1,N+1); y(1)=1.0000; x=a:h:b; for n=1:N t=y(n)+h*y(n).^2; y(n+1)=y(n)+h*(y(n).^2+t.^2)/2; end Y=-1./(x-11/10); %精确解 Error=abs(Y-y); disp([' x 数值解精确解误差']) for k=1:N+1 %disp([num2str(y(k)) ' ' num2str(Y(k))]) fprintf('%.4f %.4f %.4f %.4d\n',x(k),y(k),Y(k),Error(k)); if mod(k,5)==0 disp([' x 数值解精确解误差']) end end scatter(x,y,'o'),hold on scatter(x,Y,'r.'); grid on; axis([0.05,0.45,0.95,1.45]); title('数值解与精确解对比');%作出数据点和拟合曲线的图形 xlabel('x'); ylabel('y'); legend('数值解','精确解','Location','SouthEast');%图形标注 gtext('109040669 杨永富')改进的欧拉方法四阶龙格—库塔法： %四阶龙格—库塔法 clear all close all clc format long h=0.01; a=2.0; b=2.6; N=(b-a)/h; N=round(N); y=zeros(1,N+1); y(1) = 1; x=a:h:b; for n=1:N K1=x(n)/y(n); K2=(x(n)+0.5*h)/(y(n)+0.5*h*K1); K3=(x(n)+0.5*h)/(y(n)+0.5*h*K2); K4=(x(n)+h)/(y(n)+h*K3); y(n+1)=y(n)+h*(K1+2*K2+2*K3+K4)/6; end format short Y=sqrt(x.*x-3); Error=abs(Y-y); disp([' x 数值解精确解误差']) for k=1:N+1 %disp([num2str(y(k)) ' ' num2str(Y(k))]) fprintf('%.4f %.4f %.4f %.4d\n',x(k),y(k),Y(k),Error(k)); if mod(k,5)==0 disp([' x 数值解精确解误差']) end end scatter(x,y,'o'),hold on scatter(x,Y,'r.'); grid on; axis([1.9,2.7,0.9,2]); title('数值解与精确解对比(四阶龙格—库塔法)');%作出数据点和拟合曲线的图形 legend('数值解','精确解','Location','SouthEast');%图形标注 gtext('109040669 杨永富') xlabel('x'); ylabel('y');