微分方程数值解1.docVIP

1.取步长h 0.1，用改进Euler法解常微分方程初值问题 并将计算结果与准确解相比较。取更小的步长上机计算，结果又如何？ 解： （1）调用文件P641.m，即改进的欧拉法，可以得出计算结果，如下图： 在MATLAB中输入： dz x+y; dz f x,y ; dz x+y; [x,y] P641 f,0,1,1,0.1 ; plot x,y,-.b ; 当取更小的步长时（右图为放大图）： dz f x,y ; dz x+y; [x,y] P641 f,0,1,1,0.1 ; plot x,y,-.b ; hold on; [x,y] P641 f,0,1,1,0.01 ; plot x,y,-k ; legend 改进欧拉法精度0.1,改进欧拉法精度0.01 计算结果与准确解相比较 直接调用文件file.m可得如下图（右图为放大图）： 4.取步长h 0.1，用四阶Runge-Kutta方法解常微分方程初值问题 并和题1的计算结果比较。 解： 调用M文件file3.m,即四阶Runge-Kutta方法，得到的图形如下： 比较：将第一题中的改进欧拉法和第四题的四阶Runge-Kutta方法进行比较，调用文件file4.m：