第六章常微分方程初值问题初步1.pptVIP

第六章 常微分方程初值问题初步 * 6.1 基本理论与Euler方法 6.3用Matlab求解常微分方程的经典例 6.4常微分方程的解析解 6.2 Runge－Kutta方法（龙格－库塔法） 6.5 方程组与高阶问题 6.1 基本理论与Euler方法 function E=euler(f,a,b,ya,N) %f为该问题中的函数f(x,y)。%a,b分别为取值范围的左右端点。 %ya为给定初值y(a)。%N为迭代步数。%h为步长。 %输出值为对应每个节点的近似值。 h=(b-a)/N; T=zeros(1,N+1); Y=zeros(1,N+1); T=a:h:b; Y(1)=ya; for j=1:N Y(j+1)=Y(j)+h*feval(f,T(j),Y(j)); End T=[T Y] 在Euler算法程序的最后改为 Y1=sqrt(1+2*T); T=[T Y Y1]; 可得到解析解和迭代法计算结果比较 function T4=taylor4(df,a,b,ya,N) %df为y的一阶到四阶微商序列=[y y y y]。%y=f(x,y)。%a,b左右端点。%N为迭代步数。%h为步长。%ya为初值。 h=(b-a)/N; T=zeros(1,N+1); Y=zeros(1,N+1); T=a:h:b; Y(1)=ya; for j=1:N D=feval(d