计算方法-9.5线性多步法.pptVIP

* * 解常微分方程初值问题小结 本章介绍了常微分方程初值问题的基本数值解法。包括单步法和多步法。单步法主要有欧拉法、改进欧拉法和龙格—库塔方法。多步法是阿当姆斯法。它们都是基于把一个连续的定解问题离散化为一个差分方程来求解，是一种步进式的方法。用多步法求常微分方程的数值解可获得较高的精度。 实际应用时，选择合适的算法有一定的难度，既要考虑算法的简易性和计算量，又要考虑截断误差和收敛性、稳定性。 * * 龙格-库塔法较为常用，适用于多步方法中作初值计算和函数f(x,y)较为简单的场合。四阶标准龙格—库塔法精度高，程序简单，易于改变步长，比较稳定，也是一个常用的方法，但计算量较大。当函数f(x,y)较为复杂，可用显式阿当姆斯方法或阿当姆斯预测—校正方法，不仅计算量较小，稳定性也比较好，但不易改变步长。 一般采用龙格—库塔法提供初值y1, y2, y3，然后用阿当姆斯外推公式求得预测值 ，再由阿当姆斯内插值求得校正值yn+1，如此求得的值近似程度好且节省计算量，是一种较好的方法。 * * * * 本节作业 P214 （4） * * * 9.5 线性多步法 单步法计算时只用到前一步的结果，因此只要给定初值，计算就可以进行下去。但是Euler等单步法的精度都较低，龙格-库塔方法虽然可以得到较高的精度，但这类算法为了