第第2章章 数值积分与数值微分数值积分与数值微分牛顿牛顿-科特斯科特斯(Newton-Cotes)公式公式及其复合求积公式及其复合求积公式2021/6/161牛顿牛顿- -科特斯公式科特斯公式q 等距节点的插值型求积公式称为等距节点的插值型求积公式称为牛顿牛顿- -科特斯公式科特斯公式：取等距节点：取等距节点：xi = a + i h， ，i = 1, 2, , n令令 x = a + t h 得：得：q 插值型求积公式插值型求积公式其中其中2021/6/162牛顿牛顿- -科特斯公式（续）科特斯公式（续）注：注：Cotes 系数系数仅取决于仅取决于 n 和和 i，可通过查表得到。可通过查表得到。与被积函数与被积函数 f (x) 及积分区间及积分区间 a, b 均无关。均无关。科特斯科特斯(Cotes)系数系数q 牛顿牛顿- -科特斯公式科特斯公式：2021/6/163几个常见公式几个常见公式n = 1:代数精度代数精度 = 1梯形求积公式梯形求积公式n = 2:代数精度代数精度 = 3抛物线求积公式抛物线求积公式Simpson求积公式求积公式n = 4:科特斯科特斯(Cotes)求积公式求积公式TSC2021/6/164科特斯系数表科特斯系数表2021/6/165系数特点和稳定性系数特点和稳定性q 科特斯系数具有以下特点：科特斯系数具有以下特点：(1) (2) (3) 当当 n 8 时，出现负数，时，出现负数，稳定性得不到保证稳定性得不到保证。而且。而且当当 n 较大时，由于较大时，由于Runge现象，现象，收敛性也无法保证收敛性也无法保证。故故一般一般不采用不采用高阶的牛顿高阶的牛顿- -科特斯求积公式科特斯求积公式。q 当当 n 7 时，牛顿时，牛顿- -科特斯公式是稳定的。科特斯公式是稳定的。2021/6/166牛顿牛顿- -科特斯公式的代数精度科特斯公式的代数精度定理定理当当 n 为偶数时，牛顿科特斯公式至少有为偶数时，牛顿科特斯公式至少有 n+1 阶阶代数精度。代数精度。证证：只要证明当只要证明当 n 为偶数时，公式对为偶数时，公式对f (x)xn+1精确成立。精确成立。由插值型求积公式的误差公式得由插值型求积公式的误差公式得作变量代换作变量代换 x = a + t h，并将，并将 xi = a + i h 代入得代入得再作变量代换再作变量代换 t = n - s，得，得又又n 偶数偶数2021/6/167余项余项q 梯形公式的余项梯形公式的余项中值定理中值定理q Simpson公式的余项公式的余项三次三次Hermite插值插值2021/6/168余项的一般形式余项的一般形式定理定理(1) 若若 n 为偶数，为偶数， f (x) Cn+2a, b ，则存在，则存在 (a, b) 使得使得设设 ，则有，则有(2) 若若 n 为奇数，为奇数， f (x) Cn+1a, b ，则存在，则存在 (a, b) 使得使得2021/6/169举例（一）举例（一）q 例：例：分别用梯形公式和分别用梯形公式和simpson公式计算积分公式计算积分 解：解：a0, b1, f (x) = e -x ，由由 simpson 公式可公式可得得 由由梯形公式可梯形公式可得得 与精确值与精确值 0.6321 相比相比得误差分别为得误差分别为 0.0518 和和 0.0002。2021/6/1610复合求积公式复合求积公式q 提高积分计算精度的常用两种方法提高积分计算精度的常用两种方法 用用 复合公式复合公式 用用 非等距节点非等距节点q 复合求积公式：复合求积公式：将积分区间分割成多个小区间，然将积分区间分割成多个小区间，然后在每个后在每个小区间小区间上使用低次牛顿科特斯求积公式。上使用低次牛顿科特斯求积公式。q 将将a, b 分成分成 n 等分等分 xi , xi+1 ，其中节点，其中节点(i = 0, 1, , n)2021/6/1611复合梯形公式复合梯形公式q 复合梯形公式复合梯形公式:Tnq 余项：余项：， (a, b)2021/6/1612复合复合simpson公式公式q 复合复合simpson公式公式:Snq 余项：余项：， (a, b)444442021/6/1613复合科特斯公式复合科特斯公式q 复合复合cotes公式公式:Cnq 余项：余项：， (a, b)2021/6/1614举例（二）举例（二）解：解：q 例：例：设设 ，利用下表中的数据分别用复合梯，利用下表中的数据分别用复合梯形公式和复合形公式和复合simpson公式计算积分公式计算积分 xi01/82/83/84/85/86/87/81.0f (xi )10.9970.9900.9770.9540.9360.9090.8770.8412021/6/1615h 很小时的误差很小时的误差 i (xi, xi+1 )(h 0)定积分定义定积分定义即即同理同理2021/6/1616 收敛速度与误差估计收敛速度与误差估计 若若一一个个积积分分公公式式的的误误差差满满足足 且且C 0，则称该公式是则称该公式是 p 阶收敛阶收敛的。的。例：例：计算计算解：解：其中其中= 3.138988494其中其中= 3.141592502运算量基运算量基本相同本相同2021/6/1617Q: 给定精度给定精度 ，如何取，如何取 n ?例如：要求例如：要求 ，如何判断，如何判断 n = ？?上例中若要求上例中若要求 ，则，则即：取即：取 n = 4092021/6/1618通常采取将区间通常采取将区间不断对分不断对分的方法，即取的方法，即取 n = 2k上例中上例中2k 409 k = 9 时，时，T512 = 3.14159202注意到区间再次对分时注意到区间再次对分时可用来判断迭代可用来判断迭代 是否停止。是否停止。Q: 给定精度给定精度 ，如何取，如何取 n ?2021/6/16192.3 龙贝格算法龙贝格算法n梯形法的递推化梯形法的递推化n龙贝格算法龙贝格算法n理查森外推加速法理查森外推加速法2021/6/16201 梯形法的递推化梯形法的递推化方法思路方法思路 : :复化求积方法可提高求积精度，实际计算复化求积方法可提高求积精度，实际计算时可以将步长逐次分半。时可以将步长逐次分半。在每个子区间在每个子区间 x xk k,x,xk+1k+1 经过二分只增加了一经过二分只增加了一个分点个分点x xk+1/2k+1/2=1/2(x=1/2(xk k+x+xk+1k+1),),用复化梯形公式用复化梯形公式求得该子区间上的积分值为求得该子区间上的积分值为2021/6/1621注意，这里h=(a+b)/n代表二分前的步长。将每个子区间上的积分值相加得从而可导出下列递推公式从而可导出下列递推公式1 梯形法的递推化梯形法的递推化2021/6/1622 龙贝格算法龙贝格算法龙贝格积分法是在计算梯形和序列的龙贝格积分法是在计算梯形和序列的基础上应用了线性外推的加速方法，基础上应用了线性外推的加速方法，由此构成的一种具有超线性收敛的自由此构成的一种具有超线性收敛的自动积分法动积分法 2021/6/1623基本思想基本思想根据复化梯形公式的余项表达式可知2021/6/1624将上式移项整理，可得将上式移项整理，可得可以做这样的补偿可以做这样的补偿基本思想基本思想2021/6/1625同理同理由此得到由此得到同理同理基本思想基本思想2021/6/1626由此法，可得如下三角形数表由此法，可得如下三角形数表梯形梯形辛卜生辛卜生柯特斯柯特斯龙贝格龙贝格T0T3T2T1S0 S2S1 C0 C1 D0基本思想基本思想2021/6/1627样条插值积分样条插值积分q 用三次样条插值函数用三次样条插值函数 S(x) 近似被积函数近似被积函数 f (x) ，从而得到样条插值积分公式。，从而得到样条插值积分公式。(i = 0, 1, , n) 将将a, b 分分 n 等分等分 ， ， 设设 S(xi)mi ，则，则 S (x) 在在 xi , xi+1 上为满足以下条件上为满足以下条件的三次多项式：的三次多项式：，由三次由三次 Hermite 插值多项式公式插值多项式公式（P.46）可得可得2021/6/1628样条插值积分（续）样条插值积分（续）于是有于是有 由于由于 S (x) 在在 xi , xi+1 上为三次多项式，所以上为三次多项式，所以simpson公式精确成立，即公式精确成立，即于是得积分公式于是得积分公式Un2021/6/1629样条插值积分（续）样条插值积分（续）UnTn余项：余项：只增加计算两端点的导数，计算精度即由只增加计算两端点的导数，计算精度即由O(h2) 提高到提高到O(h4) 2021/6/1630 结束语结束语若有不当之处，请指正，谢谢！若有不当之处，请指正，谢谢！