原创力文档- 1、本文档共47页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
例1观察向前欧拉、向后欧拉算法计算情况。与精确解进行比较。误差有多大？解：1）解析解：y=x+e-xy=dsolve(Dy=-y+x+1,y(0)=1,x)1、欧拉方法第29页,共47页，星期六，2024年，5月2）向前欧拉法：yn+1=yn+h(-yn+xn+1)=(1-h)yn+hxn+h3）向后欧拉法：yn+1=yn+h(-yn+1+xn+1+1)转化yn+1=(yn+hxn+1+h)/(1+h)y’=f(x,y)=-y+x+1;1、欧拉方法第30页,共47页，星期六，2024年，5月x1(1)=0;y1(1)=1;y2(1)=1;h=0.1;%(died.m)fork=1:10x1(k+1)=x1(k)+h;y1(k+1)=(1-h)*y1(k)+h*x1(k)+h;y2(k+1)=(y2(k)+h*x1(k+1)+h)/(1+h);endx1,y1,y2,%（y1——向前欧拉解，y2——向后欧拉解）x=0:0.1:1;y=x+exp(-x)%（解析解）plot(x,y,x1,y1,k:,x1,y2,r--)1、欧拉方法第31页,共47页，星期六，2024年，5月x精确解向前欧拉向后欧拉01110.11.004811.00910.21.01871.011.02640.31.04081.0291.05130.41.07031.05611.08300.51.10651.09051.12090.61.14881.13141.16450.71.19661.17831.21320.81.24931.23051.26650.91.30661.28741.324111.36791.34871.3855（1）步长h=0.1的数值解比较表计算结果第32页,共47页，星期六，2024年，5月（2）步长h=0.01的数值解比较表x精确解向前欧拉向后欧拉01110.11.00481.00441.00530.21.01871.01791.01950.31.04081.03971.04190.41.07031.06901.07170.51.10651.10501.10800.61.14881.14721.15040.71.19661.19481.19830.81.24931.24751.25110.91.30661.30471.308411.36791.36601.3697结论：显然迭代步长h的选取对精度有影响。第33页,共47页，星期六，2024年，5月图形显示有什么方法可以使精度提高？向后欧拉法返回第34页,共47页，星期六，2024年，5月对方程y’=f(x,y),两边由xi到xi+1积分，并利用梯形公式，有：2、使用数值积分即梯形法：第35页,共47页，星期六，2024年，5月梯形公式改进欧拉公式yn+hf(xn,yn)返回第36页,共47页，星期六，2024年，5月以例1为例，用改进欧拉公式编程计算，再与精确解的比较。yn+1=yn+(h/2)*[(-yn+xn+1)+(-yn+1+xn+1+1)]=yn+(h/2)*[(-yn+xn+1)-(yn+h*(-yn+xn+1))+xn+1+1]=yn+(h/2)*[(1-h)*xn+xn+1+2-h+(h-2)*yn]died1.m使用改进欧拉公式的例第37页,共47页，星期六，2024年，5月x精确解向前欧拉向后欧拉改进欧拉011110.11.004811.00911.0050.21.01871.011.02641.0190.31.04081.0291.05131.04120.41.07031.05611.08301.07080.51.10651.09051.12091.10710.61.14881.13141.16451.14940.71.19661.1783
文档评论（0）