计算方法 教学课件 作者 江爱民 主编 张有正 副主编 nn04.pptVIP

上一页 下一页 返回首页 第4章 数值积分 4.1 数值积分的基本思路 4.2 牛顿-柯特斯公式 4.3 梯形求积公式 4.4 辛卜生公式及龙贝格求积法 4.5 高斯求积公式 4.1 数值积分的基本思路 4.1.1 数值积分的必要性 计算定积分 仅有牛顿—莱布尼兹公式 是不够的，这主要有以下几方面原因： 4.1.1 数值积分的基本思想 方法一 利用函数的级数展开，转化为简简单函数的积分问题，对于这一种解法对函数本身的要求还是比较高的，即函数能够按某一种方式展开，比如展开成幂级数或三角级数等，然后用级数的部分和代替函数进行积分。 对定积分，我们常用以下方法进行 方法二 而多项式分的积分则是我们熟悉的。下面我们主要讨论这一方法。 4.2 牛顿柯特斯公式 思路 利用插值多项式 则积分易算。 ? 在[a, b]上取 a ? x0 x1 … xn ? b，做 f 的 n 次插值多项式 ，即得到 Ak（x 由 决定， 与 无关。 节点 f x 误差 ?k ? 当节点等距分布时： 令 Cotes系数 注：Cotes 系数仅取决于 n 和 k，可查表得到。与 f x 及区间[a, b]均无关。 即 n 2: n 1: 梯形公式 辛卜生公式 到这里，积分问题似乎已经解决了，只要用高阶插值多项式替代f x 进行积分，但是高次插值有Runge 现象，故采用分段低次插值——分段低次合成的 Newton-Cotes 复合求积公式。 4.3 梯形求积公式 4.3.1 基本公式 1阶牛顿—柯特斯公式即为梯形求积公式 其误差为: h b-a 注:从余项中看出，误差的大小取决于h的大小，当比较大h 时精度是不高的，而误差太大就会失去数值计算的实际意义. 2.复化的梯形公式 1. 梯形公式 在每个 上用梯形公式： Tn /*中值定理*/ 将[a,b]分成n等分， 4.3.2 变步长梯形公式 4.4 辛卜生公式及龙贝格求积法 4.4.1 辛卜生基本求积公式 2阶的牛顿-柯特斯公式 即为辛卜生公式 余项： 将[a,b]n等分 4.4.2 复化辛卜生公式 4 4 4 4 4 Sn n级辛卜生公式 4.4.3 龙贝格求积法 在复化的梯形公式中，将区间[a,b]n等分,h b-a /n ,其余项: 再将区间对分，得到： 令： 由 来计算 I 效果是否好些？ 考察 Sn 即： 余项： 比较： 一般有： Romberg 公式 ? Romberg 算法： ? ? ? ? ? T1 0 0 T ? T8 3 0 T ? T4 2 0 T ? T2 1 0 T ? S1 0 1 T ? R1 0 3 T ? S2 1 1 T ? C1 0 2 T ? C2 1 2 T ? S4 2 1 T 四阶牛顿-柯特斯公式 以上所述的方法称为龙贝格求积法，可以描述如下 此列是用梯形公式计算的 上一页 下一页 返回首页