数值分析801(免费阅读).ppt

湘潭大学数学与计算科学学院 * 稳定性 定义 湘潭大学数学与计算科学学院 * 结论: 湘潭大学数学与计算科学学院 * 计算问题： 隐式计算格式由迭代法去完成. 将上式变形为 记 求 即求隐式方程 的根。 湘潭大学数学与计算科学学院 * 总结 通过对Euler方法的讨论可以看到,微分 方程数值方法的研究应包括以下方面 1.数值计算公式的构造; 2.方法稳定性,收敛性的研究; 3.方法的误差估计; 4.方法的实现等. 湘潭大学数学与计算科学学院 * §3 龙格 - 库塔法 建立高精度的单步递推格式。 单步递推法的基本思想是从 ( xi , yi ) 点出发，以某一斜率沿直线达到 ( xi+1 , yi+1 ) 点。欧拉法及其各种变形所能达到的最高精度为2阶。 ? 考察改进的欧拉法，可以将其改写为： 斜率 一定取K1 K2 的平均值吗？ 步长一定是一个h 吗？ 湘潭大学数学与计算科学学院 * 首先希望能确定系数 ?1、?2、p，使得到的算法格式有2阶精度，即在 的前提假设下，使得 Step 1: 将 K2 在 ( xi , yi ) 点作 Taylor 展开 将改进欧拉法推广为： ) , ( ) , ( ] [ 1 2 1 2 2 1 1 1 phK y ph x f K y x f K K K h y y i i i i i i + + = = + + = + l l Step 2: 将 K2 代入第1式，得到 湘潭大学数学与计算科学学院 * Step 3: 将 yi+1 与 y( xi+1 ) 在 xi 点的泰勒展开作比较 要求 ，则必须有： 这里有 个未知数， 个方程。 3 2 存在无穷多个解。所有满足上式的格式统称为2阶龙格 - 库塔格式。 注意到， 就是改进的欧拉法。 Q: 为获得更高的精度，应该如何进一步推广？ 湘潭大学数学与计算科学学院 * 其中?i ( i = 1, …, m )，?i ( i = 2, …, m ) 和 ?ij ( i = 2, …, m; j = 1, …, i?1 ) 均为待定系数，确定这些系数的步骤与前面相似。 ) ... , ( ... ... ) , ( ) , ( ) , ( ] ... [ 1 1 2 2 1 1 2 32 1 31 3 3 1 21 2 2 1 2 2 1 1 1 - - + + + + + + = + + + = + + = = + + + + = m m m m m m i m i i i i i i m m i i hK hK hK y h x f K hK hK y h x f K hK y h x f K y x f K K K K h y y b b b a b b a b a l l l ? 最常用为四级4阶经典龙格-库塔法 ： 湘潭大学数学与计算科学学院 * 注： ? 龙格-库塔法的主要运算在于计算 Ki 的值，即计算 f 的值。Butcher 于1965年给出了计算量与可达到的最高精度阶数的关系： 7 5 3 可达到的最高精度 6 4 2 每步须算Ki 的个数 ? 由于龙格-库塔法的导出基于泰勒展开，故精度主要受解函数的光滑性影响。对于光滑性不太好的解，最好采用低阶算法而将步长h 取小。 湘潭大学数学与计算科学学院 * §4 收敛性与稳定性 ? 收敛性 定义 　　　 若某算法对于任意固定的 x = xi = x0 + i h，当 h?0 ( 同时 i ? ?) 时有 yi ? y( xi )，则称该算法是收敛的。 例：就初值问题 考察欧拉显式格式的收敛性。 解：该问题的精确解为 欧拉公式为 对任意固定的 x = xi = i h ，有 ? ? 湘潭大学数学与计算科学学院 * ? 稳定性 例：考察初值问题 在区间[0, 0.5]上的解。 分别用欧拉显、隐式格式和改进的欧拉格式计算数值解。 0.0 0.1 0.2 0.3 0.4 0.5 精确解 改进欧拉法 欧拉隐式 欧拉显式 节点 xi 1.0000 ?2.0000 4.0000 ?8.0000 1.6000?101 ?3.2000?101 1.0000 2.5000?10?1 6.2500?10?2 1.5625?10?2 3.9063?