第3章_常微分方程的差分方法.pptVIP

稳定性在微分方程的数值解法中是一个非常重要的问题。因为微分方程初值问题的数值方法是用差分格式进行计算的，而在差分方程的求解过程中，存在着各种计算误差，这些计算误差如舍入误差等引起的扰动，在传播过程中，可能会大量积累，对计算结果的准确性将产生影响。这就涉及到算法稳定性问题。 算法稳定性 例 3 对初值问题 证明用梯形公式求得的近似解为 并证明当步长h?0时,yn收敛于精确解 证明: 解初值问题的梯形公式为 ∵ ∴ 整理成显式 反复迭代,得到 ∵ ∴ 由于 ，有 ∴ 证毕 作业题！ 1、列出下列初值问题的欧拉格式： （1） 0=x=0.4, y(0)=1, 取h=0.2 （2） 1=x=1.2 y(1)=1, 取h=0.1 2、用梯形方法做 1=x=2 y(1)=2, 取h=0.1， 另外用改进欧拉法计算。 基本思想: Euler公式可改写成 则yi+1的表达式y(xi+1)与的Taylor展开式的前两项完全相同,即局部截断误差为 改进的Euler公式又可改写成 龙格-库塔（Runge-Kutta）法 上述两组公式在形式上有一个共同点:都是用f(x,y)在某些点上值的线性组合得出y(xi+1)的近似值yi+1,而且增加计算的次数f(x,y)的次数,可提高截断误差的阶。如欧拉公式:每步计算一次f(x,y)的值,为一阶方法。改进欧拉公式需计算两次f(x,y)的值，它是二阶方法。它的局部截断误差为 。 于是可考虑用函数f(x,y)在若干点上的函数值的线性组合来构造近似公式，构造时要求近似公式在(xi,yi)处的Taylor展开式与解y(x)在xi处的Taylor展开式的前面几项重合，从而使近似公式达到所需要的阶数。既避免求偏导,又提高了计算方法精度的阶数。或者说,在 这一步内多预报几个点的斜率值，然后将其加权平均作为平均斜率，则可构造出更高精度的计算格式，这就是龙格—库塔（Runge-Kutta）法的基本思想。 二阶龙格—库塔法 在 上取两点xi和 ,以该两点处的斜率值k1和k2的加权平均(或称为线性组合)来求取平均斜率k*的近似值K，即 式中:k1为xi点处的切线斜率值， k2为 点处的切线斜率值,比照改进的欧拉法,将 视为 ，即可得 对常微分方程初值问题(3.1)式的解 y=y(x),根据微分中值定理，存在点 ，使得 式中 K可看作是y=y(x)在区间 上的平均斜率。所以可得计算公式为： （3.14） 将y(xi)在x=xi处进行二阶Taylor展开： （3.15） 也即 将 在x=xi处进行一阶Taylor展开： 将以上结果代入（3.14）得： （3.16） 对式(3.15)和(3.16)进行比较系数后可知,只要 （3.17） 成立,格式(3.14)的局部截断误差就等于 有2阶 精度 式(3.17)中具有三个未知量,但只有两个方程,因而有无穷多解。若取 ,则p=1，这是无穷多解中的一个解，将以上所解的值代入式(3.14)并改写可得 不难发现，上面的格式就是改进的欧拉格式。凡满足条件式（3.17）有一簇形如上式的计算格式，这些格式统称为二阶龙格—库塔格式。因此改进的欧拉格式是众多的二阶龙格—库塔法中的一种特殊格式。 若取 ,则 ，此时二阶龙格-库塔 法的计算公式为 此计算公式称为变形的二阶龙格—库塔法。式中 为区间 的中点。 三阶龙格-库塔法 为了进一步提高精度，设除 外再增加一点 并用三个点 , , 的斜率k1,k2,k3加权平均 得出平均斜率k*的近似值，这时计算格式具有形式: （3.18） 为了预报点 的斜率值k3,在区间 内有两 个斜率值k1和k2可以用,可将k1,k2加权平均得出 上的平均斜率,从而得到 的预报值 于是可得 运用Taylor展开方法选择参数 ,可以使格式(3.18)的局部截断误差为 ,即具有三阶精度，这类格式统称为三阶龙格—库塔方法。下列是其中的一种，称为库塔（Kutta）公式。 （3.19） 四阶龙