matlabeuler法解常微分方程.pdfVIP

禁止复制 北京石油化工学院 Onlyunited Euler 法解常微分方程 Euler 法解常微分方程算法： Step 1 分别取积分上限、积分下限、步长 Step 2 计算n  n  h 判断n  b 是否成立，成立转到 Step 3，否则继续进行 Step 4 Step 3 计算y n1  y n  hf (x n , y n ) Step 4 y n1  y n  hf (x n , y n ) Euler 法解常微分方程算程序： function euler2(fun,y0,A,h) %fun--y %y0初值 %Ax取值范围 %ax左区间端点值 %bx右区间端点值 %h给定步长 x=min(A); b=max(A); y=y0; while xb-h b=y; y=y+h*feval(fun,x,b) x=x+h; end  2 y  y  xy 例：用Euler 法计算下列初值问题（取步长 h=0.2） (0  x  0.6) y (0)  1  输入：fun=inline(-y-x*y^2) euler2(fun,1,[0 0.6],0.2) 得到： y = 0.8000 y = 0.6144 y = 禁止复制 北京石油化工学院 Onlyunited 0.4613 指导教师： 年 月 日 改进 Euelr 法解常微分方程 改进 Euler 法解常微分方程算法： Step 1 分别取积分上限、积分下限、步长 Step 2 取一个以 h 为步长，a，b 分别为左右端点的矩阵 Step 3 （1）做显性Euler 预测y i1  y i  hf (x n , y n ) h （2）将y 带入y  y  [f (x , y )  f (x , y )] i1 i1 i i i i1 i1 2 Step 4 计算n  n  h 判断n  b 是否成立，成立返回 Step 3，否则继续进行 Step 5 h Step 5 y  y  [f (x , y )  f (x , y )] i1 i i i i1 i1 2 改进 Euler 法解常微分方程算程序： function gaijineuler2(fun,y0,A,h) %fun--y %y0初值 %Ax 取值范围 %ax 左区间端点值 %bx 右区间端点值 %h给定步长 a=min(A); b=max(A); x=a:h:b; y(1)=y0; for i=1:length(x)-1 w1=feval(fun,x(i),y(i)); y(i+1)=y(i)+h*w1; w2=feval(fun,x(i+1),y(i+1)); y(i+1)=y(i)+h*(w1+w2)/2; end x=x