原创力文档- 1、本文档共16页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
* Range-Kutta公式 捕食者与被捕食者模型 单摆的数学模型 常微分方程边值问题 线性多步法简介 ? ? ? ? ? 《数值分析》 24 初值问题 欧拉公式: yn+1 = yn+ h f(xn , yn) 修改的欧拉公式(Heun方法): k1 = f(xn , yn) , k2 = f( xn+1 , yn+ h k1) (n = 0, 1, 2, ·······, N ) (n = 0, 1, 2, ·······, N ) 2/16 二阶Range-Kutta公式一般形式 yn+1= yn+ h[c1k1+ c2k2] k1=f (xn, yn) k2=f(xn+λ2h, yn+μ21hk1) ( n=0,1,····,N ) c1=0.5, c2 = 0.5, λ2=1, μ21= 1 c1= 0, c2 = 1, λ2=0.5, μ21=0.5 yn+1= yn+ 0.5h[ k1+ k2] k1=f (xn, yn) k2=f(xn+h, yn+hk1) yn+1= yn+ h k2 k1=f (xn, yn) k2=f(xn+.5h, yn+.5hk1) 3/16 三阶Range-Kutta公式一般形式 yn+1= yn+ h[k1+4k2+k3]/6 k1=f(xn,yn), k2=f(xn+0.5h, yn+0.5hk1) k3=f(xn+h, yn – hk1+2hk2) 四阶Range-Kutta公式一般形式 yn+1= yn+ h[k1+2k2+2k3+k4]/6 k1=f(xn,yn), k2=f(xn+0.5h, yn+0.5hk1) k3=f(xn+0.5h, yn+0.5hk2), k4=f(xn+h, yn+hk3) 4/16 例1 数值实验:几种不同求数值解公式的误差比较 n 10 20 30 40 h 0.2 0.1 0.0667 0.05 RK4 6.862e-005 3.747e-006 7.071e-007 2.186e-007 RK3 0.0012 1.529e-004 4.517e-005 1.906e-005 RK2 0.0123 0.0026 0.0011 5.9612e-004 Euler 0.1059 0.0521 0.0342 0.0256 5/16 n=input(input n:= ); f=inline(y-x.*y.^2); h=2/n;x=h:h:2; y0=1;k1=f(0,y0); k2=f(.5*h,y0+0.5*h*k1); k3=f(.5*h,y0+.5*h*k2); k4=f(h,y0+h*k3); y(1)=y0+h*(k1+2*k2+2*k3+k4)/6; for k=1:n-1 xk=x(k);yk=y(k);k1=f(xk,yk); k2=f(xk+h/2,yk+0.5*h*k1); k3=f(xk+h/2,yk+0.5*h*k2); k4=f(xk+h,yk+h*k3); y(k+1)=yk+h*(k1+2*k2+2*k3+k4)/6; end t=-1:h/2:5; yy=1./(t-1+2*exp(-t)); plot(0,1,ro,x,y,ro,t,yy) y1=1./(x-1+2*exp(-x)); error=max(abs(y-y1)) input n:= 20 error = 3.7475e-006 6/16 例2. 捕食者与被捕食者模型 0≤ x ≤ 15 y1(0)=20, y2(0)=20 记 7/16 x(1)=20;y(1)=20; for n=1:length(t)-1 tn=t(n);xn=x(n);yn=y(n); k11=f1(xn,yn);k21=f2(xn,yn); k12=f1(xn+.5*h*k11,yn+.5*h*k21); k22=f2(xn+.5*h*k11,yn+.5*h*k21); k13=f1(xn+.5*h*k12,yn+.5*h*k22); k23=f2(xn+.5*h*k12,yn+.5*h*k22); k14=f1(x
文档评论（0）