第7-3龙格-库塔方法.pptVIP

作业 * §7.3 龙格－库塔方法 一、Taylor展开法 设 截取有限项作为 的近似值,有 Taylor公式 局部截断误差 二、龙格－库塔方法的基本思想 龙格首先提出用n个点处 f 的函数值的线性组合来代替 y 的导数值，然后按泰勒公式展开，确定其中的系数，以提高精度. Taylor公式具有 p阶精度. 当 p=1时,即为Euler公式. Taylor公式需要计算 y 的各阶导数值,当阶数 较高时,求导过程很复杂. 多预报几个点的斜率值，将它们的线性组合作为平均斜率的近似值，有可能构造出精度更高的计算格式. ) , ( 1 1 1 ? ? - = = + + = - - ? i j j ij n i n i i v i i i hk a y h c x f k f v k k 值的个数 所使用的 待定权因子 令 w w 【定理4】（二元函数的泰勒公式） 三、二阶龙格－库塔公式 ———中点公式 ———Heun公式 ———改进的欧拉公式 四、三阶及四阶龙格－库塔公式 三阶龙格—库塔公式有： 三阶Heun公式 四阶龙格—库塔公式有： 例1 P 210 例7.3 为了分析经典R-K公式的计算量和计算精度，将四阶经典R-K公式与一阶显式Euler公式及二阶改进的Euler公式相比较。一般说来，公式的级数越大，计算右端项 f 的次数越多，计算量越大。在同样步长的情况下，Euler方法每步只计算一个函数值，而经典方法要计算4个函数值。四阶R-K法的 计算量差不多是改进的Euler公式的2倍，是显式Euler公式的4倍。 下面的例子中Euler方法用步长 ，二阶改进的Euler法用步长 ，而四阶经典公式用步长 。这样，从 到 三种方法都计算了4个函数制，计算量大体相当。 例2 考虑初值问题 在x=0.5处，三种方法的误差分别是 经典R-K法对多数好条件问题能获得好的效果。 Euler法 改进Euler法 经典R-K法 准确解 h=0.025 h=0.05 h=0.1 0.1 0.096312 0.095123 0 0 0.2 0.183348 0.181193 0 0 0.3 0.262001 0.259085 0 0 0.4 0.333079 0.329563 0 0 0.5 0.397312 0.393337 0 0