计算方法ch7讲解.pptVIP

R阶龙格和库塔（ Runge-Kutta)法的一般形式： 练习 取步长 ，求解初值问题 用经典的四阶龙格—库塔法求 的值。 第7章 常微分方程数值解法 7.1 欧拉法和改进的欧拉法 7.2 龙格-库塔法 7.3 线性单步法的收敛性和稳定性 引言 7.4 线性多步法 7.5 一阶常微分方程组和高阶方程 7.3 线性单步法的收敛性和稳定性 任何一种数值方法，必须能提供比较可靠的计算结果，收敛性和稳定性从不同的角度描述了数值方法的可靠性。由于多步法的收敛性和稳定性比较复杂，本节只讨论单步法的收敛性和稳定性。 2. 稳定性 1. 收敛性 1. 收敛性 O(hp+1), 类似的也可证明其他方法的收敛性 7.3 线性单步法的收敛性和稳定性 任何一种数值方法，必须能提供比较可靠的计算结果，收敛性和稳定性从不同的角度描述了数值方法的可靠性。由于多步法的收敛性和稳定性比较复杂，本节只讨论单步法的收敛性和稳定性。 2. 稳定性 1. 收敛性 检验稳定性的方法： 计算机求解时，每步都有可能有舍入误差。 定义：若在计算过程中某部引入的舍入误差，在以后的计算 （或传播）中被压缩、衰减，则称)是数值稳定的；若在传播中 误差被放大，则称是数值不稳定的； 使用检验方程 ，其中 为常数。 2. 稳定性 稳定性不仅与算法有关，而且与方程中的函数f(x,y)也有关，讨论比较复杂。为简单起见，通常只针对模型方程讨论。 若一个数值方法对模型方程式稳定的，并不能保证该方法对任何方程都稳定，但某个方法若对模型方程都不稳定，也就很难用于其他方程的求解。 检验方程 试验方程 模型方程 使用检验方程 ，其中 为常数。 2. 稳定性 该模型方程的初值问题的解是 所以 显然 （0,1） 　　　　　　，任何单步法用于它都不稳定 使用检验方程 ，其中 小于零。 定义7-3 向前欧拉公式的稳定性 向后欧拉公式的稳定性 四阶R-K公式的稳定性 o 向前欧拉公式的稳定性 o 向后欧拉公式的稳定性 若一个数值方法的绝对稳定区域包括了复平面λh的左半个平面，则称该方法是Ａ－稳定的。这时计算的数值稳定性比较好。 显然，向前欧拉法不是Ａ－稳定的。向后欧拉法是Ａ－稳定的。 只要选取适当的h,就能满足稳定条件。比如:λ是负实数时，取h=-2/ λ就可以满足条件。 四阶R-K公式的稳定性 显然 绝对稳定区域为 因此，四阶龙格－库塔法是条件稳定的。 2. 稳定性 稳定性的好坏依次为： (3)经典 Runge-Kutta方法， (1)改进的Euler方法， (2)隐式Euler方法， (4)显式Euler和预估－校正Euler方法。 第7章 常微分方程数值解法 7.1 欧拉法和改进的欧拉法 7.2 龙格-库塔法 7.3 线性单步法的收敛性和稳定性 引言 7.4 线性多步法 7.5 一阶常微分方程组和高阶方程 7.4 线性多步法 7.4.1 四阶阿达姆斯(Adams)外插公式 7.4.2 四阶阿达姆斯(Adams)内插公式 多步法的概念 7.4.3 初始出发值的计算 7.4.4 阿达姆斯预测-校正公式 与其他算法的实例比较 单步法----计算yn+1时只使用yn的值。 多步法----计算yn+1时使用前面的k个yi 值，即由 yn-k+1 , yn-k+2 , …,yn-1, yn计算yn+1。(k=1,2,……) 线性多步法----计算yn+1的公式由yn-k+1 , yn-k+2 , …,yn-1, yn 的线性组合表达。 多步法的概念 7.4.1 四阶阿达姆斯(Adams)外插公式 设想用yn-3 ,yn-2 ,yn-1 ,yn 的值计算yn+1。为方便，讨论由 出发计算y(x+h), 由初值问题的方程y′=f(x,y(x))两边从x到x+h积分, 可得到等价的积分方程 设想运用数值积分方法，取 x-3h, x-2h, x- h, x 为插值基点 做 f(s,y(s)) 的三次拉格朗日插值，用它近似计算上式的积分 。 这样得到的数值积分公式是f(s,y(s))在4个插值基点处的函数值的线性组合。 由于 f(x-ih,y(x-ih))= y′(x-ih)，所得到的计算 y(x+h) 的近似公式形为： 为达到四阶精度，希望确定参数b0 ,b1 ,b2 ,b3使满足 运用在 h=0 处的泰勒公式得 求得 代入前面的公式得 令x=xn 并记 为使误差等于O(h5)，令h, h2,h3,h4 的系数为0， 得方程组：