歐拉法的应用.docVIP

1.取h=0.2用Euler格式、隐式Euler格式和改进的Euler格式解微分方程（0≤x≤1）,并比较精度（解析解y=）。 1.运用欧拉法，其主程序段（M函数）如下。 function [ x,y ] = naler( dyfun,xspan,y0,h) x=xspan(1):h:xspan(2); y(1)=y0; for n=1:length(x)-1; y(n+1)=y(n)+h*feval(dyfun,x(n),y(n)); end x=x;y=y; 运行过程及中间结果： clear; dyfun=inline(y-2*x/y); [x,y]=naler(dyfun,[0,1],1,0.2); p=[x,y] p = 0 1.0000 0.2000 1.2000 0.4000 1.3733 0.6000 1.5315 0.8000 1.6811 1.0000 1.8269 2.运用隐式的欧拉法，其主程序段（M函数）如下。 function [ x,y ] = yseuler( dyfun,xspan,y0,h) x=xspan(1):h:xspan(2); y(1)=y0; for n=1:length(x)-1; y(n+1)=iter(dyfun,x(n+1),y(n),h); end x=x;y=y function y=iter(dyfun,x,y,h) y0=y;e=1e-4;k=1e+4;y=y+h*feval(dyfun,x,y);y1=y+2*e;k=1; while abs(y-y1)e y1=y; y=y0+h*feval(dyfun,x,y); k=k+1; if kk error(μü′ú·￠é￠); end end [x,y]=yseuler(dyfun,[0,1],1,0.2);[x,y] y = 1.0000 1.1641 1.3014 1.4146 1.5019 1.5561 ans = 0 1.0000 0.2000 1.1641 0.4000 1.3014 0.6000 1.4146 0.8000 1.5019 1.0000 1.5561 3.运用改进的欧拉法，其主程序段（M函数）如下。 function [ x,y ] =zseuler( dyfun,xspan,y0,h) x=xspan(1):h:xspan(2);y(1)=y0; for n=1:length(x)-1; k1=feval(dyfun,x(n),y(n)); y(n+1)=y(n)+h*k1; k2=feval(dyfun,x(n+1),y(n+1)); y(n+1)=y(n)+h*(k1+k2)/2; end x=x;y=y; [x,y]=zseuler(dyfun,[0,1],1,0.2);[x,y] ans = 0 1.0000 0.2000 1.1867 0.4000 1.3483 0.6000 1.4937 0.8000 1.6279 1.0000 1.7542 2 3 4 num=[5 100];den=[1 8 32 80 100]; sys=tf(num,den);step(sys); axis=([0 10 0 1.2]); grid on; 5 无干扰时 有干扰时