常微分方程数值解讲解.pptVIP

　　若将(4.2)展开式写成 则 称为局部截断误差主项. 以上定义对隐式单步法(4.1)也是适用的.例如，对隐式欧拉法(3.3)其局部截断误差为 这里p=1是1阶方法，局部截断误差主项为 　　同样对梯形法(3.5)有 所以梯形方法(3.5)是二阶的. 其局部截断误差主项为 * 分别将 y(xn+1), y(xn-1) 在 xn 点处用泰勒公式展开： 二步欧拉法的局部截断误差： 二步欧拉法： 2 阶精度 * 梯形公式： y(xn+1) 在 xn 点处展开： 2 阶精度 * 各种欧拉法的比较 方法 精度 评述 显式欧拉法 1 最简单，精度低 隐式欧拉法 1 不便计算，稳定性好 二步欧拉法 2 需要两步初值，且第 2 个初值只能由其它方法给出，可能对后面的递推精度有影响 梯形公式 欧 拉 法 2 精度有所提高，但为隐式，需要迭代求解，计算量大 * 提高计算精度的方法 增加展开项数 增加y(xn+1)函数的泰勒展开式的阶数，即逐次取其高次项，这样就可以减少截断误差。其误差系列为 O(h2), O(h3), O(h4), … 缩小步长值h 由于差分间隔减小，即h值变得很小，可使计算精度提高。在理论上步长值h越小越好，而且h趋近于0最好。但是，h越小会使计算次数增加，浪费计算时间，而且还会使舍入误差增大 * 10.4.2单步法的收敛性