《常微分方程数值解的若干Matlab函数文件实用教案》由会员分享，可在线阅读，更多相关《常微分方程数值解的若干Matlab函数文件实用教案（7页珍藏版）》请在装配图网上搜索。
1、 function x,y=euler(f,tspan,y0,n) % 解初值问题： =f(x,y),y(a)=y0。 % 使用 n 步的 Euler 法，步长 h=(b-a)/h。 a=tspan(1);b=tspan(2);h=(b-a)/h; x=(a+h:h:b); y(1)=y0+h*feval(f,a,y0); for i=2:n y(i)=y(i-1)+h*fevol(f,x(i-1),y(i-1); end x=a x; y=y0 y;y 8.7常微分方程数值(shz)解的若干Matlab函数文件8.7.1 Euler 方法(fngf)的 Matlab 函数文件第1页/共6页第
2、一页，共7页。8.7.2 经典(jngdin)Runge-Kutta 法的 Matlab 函数文件function x,y =rk4(f,tspan,y0,n) % 解初值问题：y =f(x,y),y(a)=y0。 % 使用(shyng) n 步 4 阶 R-K 法。 a=tspan(1);b=tspan(2);h=(b-a)/n; k1=h*feval (f,a,y0); k2=h*feval (f,a+h/2,y0+k1/2); k3=h*feval (f,a+h/2,y0+k2/2); k4=h*feval (f,a+h,y0+k3); y(1)=y0+k1/6+k2/3+k3/3+k4
3、/6; 第2页/共6页第二页，共7页。for i=1:n-1 k1=h*feval(f,x(i),y(i); k2=h*feval(f,x(i)+h/2,y(i)+k1/2); k3=h*feval(f,x(i)+h/2,y(i)+k2/2); k4=h*feval(f,x(i)+h/2,y(i)+k3/2); y(i+1)=y(i)+k1/6+k2/3+k3/3+k4/6 ； end x=a x; y=y0 y;第3页/共6页第三页，共7页。function x,y =ad3 (f,tspan,y0,n) a=tspan(1);b=tspan(2); h=(b-a)/n;hh=h/12; x
4、=(a+h:h:b); %使用中点方法(fngf)提供开始值。 z0=feval (f,a,y0); k1=h*z0;k2=h*feval (f,a+h/2,y0+k1/2) ; y(1)=y0+k2; z(1)=feval (f,x(1)+h/2,y(1)+k1/2); 8.7.3 三阶 Adams方法的 Matlab 函数(hnsh)文件第4页/共6页第四页，共7页。k1=h*z(1);k2=h*feval (f,x(1)+h/2,y(1)+k1/2) ; y(2)=y(1)+k2 ; % 继续使用(shyng)三阶 Adams 方法。 z(2)=feval (f,x(2),y(2); y
5、(3)=y(2)+hh*(23*z(2)-16*z(1)+5*feval (f,a,y0) ; for i=3:n-1 z(i) =feval (f,x(2),y(2) ; y(i+1)=y(i)+hh(23*z(i)-16*z(i-1)+5z(i-2) : end x=a x ; y=y0 y; 第5页/共6页第五页，共7页。感谢您的观看(gunkn)！第6页/共6页第六页，共7页。NoImage内容(nirng)总结function x,y=euler(f,tspan,y0,n)。% 解初值问题：y =f(x,y),y(a)=y0。第3页/共6页。y(2)=y(1)+k2。z(2)=feval (f,x(2),y(2)。z(i) =feval (f,x(2),y(2)。第5页/共6页。感谢您的观看(gunkn)第七页，共7页。
常微分方程数值解的若干Matlab函数文件实用教案
