计算方法习题.ppt

第9章 常微分方程初值问题数值解法;其中f (x,y)是已知函数，(1.2)是定解条件也称为 初值条件。;则称f(x, y) 关于y满足利普希茨(Lipschitz)条件， L称为f(x, y) 的利普希茨常数。;定理; 所谓数值解法, 就是寻求解y(x)在一系列离散节点; 初值问题的数值解法一般采取“步进式”，即求解过程顺着节点排列的次序一步一步地向前推进. 描述这类算法，只要给出用已知信息yi, yi-1, yi-2,?计算yi+1的递推公式.;§9.2 简单的数值方法;由向前差商公式;每步计算; 几何意义是用一条初始点重合的折线来近似??示曲线。 y P1 (x1, y1) P2 (x2, y2) y=y(x) 0 x0 x1 x2 x3 x ;由向后差商公式; ; 设用欧拉公式;由于f(x, y)对y满足Lipschitz条件，则;三 单步法的局部截断误差与阶; 定义1 设y(x)是初值问题(1.1),(1.2)的准确解, 称; 定义2 若显式单步法的局部截断误差满足; 以上定义对隐式单步法也是适用的.例如，对后退欧拉法，其局部截断误差为;四、 梯形公式;;五、 改进欧拉公式 ;;例：用改进欧拉公式求解初值问题; 于是有 由y0=1计算得;由微分中值定理，有; 欧拉公式可改写成 ;将改进的欧拉公式推广为; 改进的Euler 公式推广为二阶Runge-Kutta公式带来这样的启示： 若在[xi, xi+1]上多取几个点的斜率值，然后将它们的线性组合作为平均斜率的近似值，则有可能构造出具有更高精度的计算公式。 --------Runge-Kutta方法的基本思想。;显式N级龙格-库塔公式