学习资料计算方法习题.ppt

则称f(x, y) 关于y满足利普希茨(Lipschitz)条件， L称为f(x, y) 的利普希茨常数。 定义： 若存在实数L0, 使得 定理 所谓数值解法, 就是寻求解y(x)在一系列离散节点 上的近似值 y1, y2, ?, yi , yi+1 ,?. 相邻两个节点的间距hi=xi+1-xi 称为步长. 今后如不特别说明，总是假定 hi=h(i=1,2,?)为常数, 这时节点为xi=x0+ih(i=0,1,2,?) (等距节点). 初值问题的数值解法一般采取“步进式”，即求解过程顺着节点排列的次序一步一步地向前推进. 描述这类算法，只要给出用已知信息yi, yi-1, yi-2,?计算yi+1的递推公式. 若计算yi+1时只用到前一点的值yi，称为单步法. 若计算yi+1时用到yi+1前面 k 点的值yi, yi-1,?, yi-k+1，称为k步法. §9.2 简单的数值方法 由向前差商公式 一、欧拉法 得 每步计算 只用到 依上述公式逐次计算可得： 几何意义是用一条初始点重合的折线来近似表示曲线。 y P1 (x1, y1) P2 (x2, y2) y=y(x) 0 x0 x1 x2 x3 x 欧拉法的几何意义 由向后差商公式 得 二、后退欧拉法 后退的欧拉公式与欧拉公式有着本质的区别, 后者是关于yi+1的一个直接计算公式，这类公式称作是显式的；前者公式的右端含有未知的yi+1，它实际上是关于yi+1的一个函数方程,这类方程称作是隐式的. 隐式方程通常用迭代法求解，而迭代过程的实质是逐步显式化. 设用欧拉公式 给出迭代初值　 ，用它代入后退欧拉公式式的右端，使之转化为显式，直接计算得 然后再用 代入后退欧拉公式的右端，又有 如此反复进行，得 由于f(x, y)对y满足Lipschitz条件，则 由此可知，只要hL1，迭代法 就收敛到解 三 单步法的局部截断误差与阶 初值问题(1.1),(1.2)的单步法可用一般形式表示为 其中多元函数?与f(x, y )有关，当?含有yi+1时，方法是隐式的，若不含yi+1则为显式方法，所以显式单步法可表示为 ?(x, y, h)称为增量函数，例如对欧拉法有 定义1 设y(x)是初值问题(1.1),(1.2)的准确解, 称 为显式单步法的局部截断误差. Ti+1之所以称为局部的，是假设在xi前各步没有误差.当yi=y(xi)时，计算一步，则有 所以，局部截断误差可理解为用显示单步法计算一步的误差 定义2 若显式单步法的局部截断误差满足 则称显示单步法具有p阶精度或是p阶方法. 欧拉法的局部截断误差为 一阶方法 以上定义对隐式单步法也是适用的.例如，对后退欧拉法，其局部截断误差为 一阶方法 四、 梯形公式 将 区间 积分 忽略高阶项,有 对右端采用梯形公式, 有 注：梯形公式局部截断误差 ， 即梯形公式具有2 阶精度，比欧拉公式有了进步。但注意到该公式是隐式公式. 五、 改进欧拉公式 预测: 先用显式欧拉公式作预测，算出 ) , ( 1 i i i i y x f h y y + = + 校正: 再将 代入隐式梯形公式的右边作校正，得到 1 + i y )] , ( ) , ( [ 2 1 1 1 + + + + + = i i i i i i y x f y x f h y y 改进欧拉公式 这是一种显式格式,它可以表示为嵌套形式。 改进欧拉公式 上页 下页 空议由乙方提供他人呢 空议由乙方提供他人呢 空议由乙方提供他人呢 空议由乙方提供他人呢 空议由乙方提供他人呢 空议由乙方提供他人呢 空议由乙方提供他人呢 空议由乙方提供他人呢 空议由乙方提供他人呢 空议由乙方提供他人呢 空议由乙方提