常微分方程初值问题数值解法研究.ppt

例2 用改进欧拉法解初值问题 区间为?0,1?,取步长h=0.1 解: 改进欧拉法的具体形式 本题的精确解为 ,计算结果略。 例题 * §9. 3 龙格-库塔方法 1 龙格-库塔(Runge-Kutta)法的基本思想 Euler公式可改写成 则yi+1的表达式y(xi+1)与的Taylor展开式的前两项完全相同,即局部截断误差为 。 改进的Euler公式又可改写成 * 上述两组公式在形式上有一个共同点:都是用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）法的基本思想。 龙格-库塔法的基本思想 * 4 四阶龙格—库塔法 如果需要再提高精度，用类似上述的处理方法，只需在区间 上用四个点处的斜率加权平均作为平均斜率k*的近似值，构成一系列四阶龙格—库塔公式。具有四阶精度，即局部截断误差是 。 由于推导复杂，这里从略，只介绍最常用的一种四阶经典龙格—库塔公式。 四阶龙格-库塔法 * 四阶龙格-库塔法 * （20） 本章介绍了常微分方程初值问题的基本数值解法。包括单步法和多步法。单步法主要有欧拉法、改进欧拉法和龙格—库塔方法。它们都是基于把一个连续的定解问题离散化为一个差分方程来求解，是一种步进式的方法。 实际应用时，选择合适的算法有一定的难度，既要考虑算法的简易性和计算量，又要考虑截断误差和收敛性、稳定性。 本章小结 * 第9章 常微分方程初值问题数值解法 §9.1 引言 §9.2 简单的数值方法 §9.1　引言 包含自变量、未知函数及未知函数的导数或微分的方程称为微分方程。在微分方程中, 自变量的个数只有一个, 称为常微分方程。自变量的个数为两个或两个以上的微分方程叫偏微分方程。微分方程中出现的未知函数最高阶导数的阶数称为微分方程的阶数。如果未知函数y及其各阶导数 都是一次的,则称它是线性的,否则称为非线性的。 * 在高等数学中，对于常微分方程的求解，给出了一些典型方程求解析解的基本方法，如可分离变量法、常系数齐次线性方程的解法、常系数非齐次线性方程的解法等。但能求解的常微分方程仍然是有限的，大多数的常微分方程是不可能给出解析解。 譬如 这个一阶微分方程就不能用初等函数及其积分来表达它的解。 引言 * 再如，方程 的解 ,虽然有表可查,但对于表上没有给出 的值,仍需插值方法来计算 引言 * 从实际问题当中归纳出来的微分方程，通常主要依靠数值解法来解决。本章主要讨论一阶常微分方程初值问题 ( 1 ) 在区间a ≤ x ≤ b上的数值解法。 可以证明,如果函数在带形区域 R=a≤x≤b, -∞＜y＜∞}内连续，且关于y满足李普希兹(Lipschitz)条件，即存在常数L(它与x,y无关)使 引言 * 引言 * 对R内任意两个 都成立,则方程( 1 )的解 在?a, b?上存在且唯一。 §9. 2 简单的数值方法 对常微分方程初值问题(1)式的数值解法，就是要算出精确解y(x)在区间?a,b?上的一系列离散节点 处的函数值 的近似值 。相邻两个节点的间距 称为步长，步长可以相等，也可以不等。本章总是假定h为定数，称为定步长，这时节点可表示为 * 数值方法的基本思想 * 数值解法需要把连续性的问题加以离散化，从而求出离散节点的数值解。 对常微分方程数值解法的基本出发点就是离散化。其数值解法有两个基本特点，它们都采用