湖南大学_科学计算_数值积分案例.ppt

求积公式的收敛性与稳定性 1.收敛性定义： 2稳定性定义 （定义7.3） 3.结论：若Ak0,则求积公式稳定。（证明略） 误差估计 定理7.3：n 为偶数阶的Newton-Cotes 公式至少有 n+1 次代数精度。 证明：只证对f(x)=xn+1,R(f)=0. 由插值余项定理知 等式两边积分得 由于f(x)?C2[a,b] ，且(x-a)(x-b)在[a,b]上非正 (不变 号),故根据积分中值定理知,至少存在一点??[a,b], 使 抛物线求积公式的代数精度为3,为此构造三次多项式P3(x)，满足 P3(a)=f(a), 则 等式两边从a到b积分得 由于P3(x)是三次多项式，故抛物线求积公式对它准确成立，即 这样 由于f(x)?C4[a,b] ，且 在[a,b]上非 正(不变号),故根据积分中值定理知,至少存在一点 ??[a,b],使 解 (1)将f(x)在a处展开，得 两边在[a, b]上积分，得 由于x-a在[a, b]上不变号，故有??[a,b], 使 从而得 将f(x)在(a+b)/2处展开，得 两边积分，得 由于 在[a, b]上不变号，故有??[a,b], 使 练习 P269 练习1,2,6 讨论题 1. ：若Ak0，则 2.导下列矩形求积公式 3计算Cotes系数(N≤6) 7.3 复合求积公式 /* Composite Quadrature */ 从Newton-Cotes求积公式的余项可知，被积函数所用的插 值多项式的次数越高，相应的求积公式的精度也越高，而高次 插值多项式稳定性很差，从而导致高次插值求积公式很不稳定。 因此在实际应用中，为了既能提高结果的精度，又使算法简便 且易在电子计算机上实现，往往采用复合求积的方法。所谓复 合求积，就是先将积分区间分成几个小区间，并在每个小区间 上用低阶牛顿—柯特斯公式计算积分的近似值，然后对这些近 似值求和，从而得到所求积分的近似值。由此得到的一些具有 更大实用价值的数值求积公式，统称为复合求积公式。 例如，先将区间[ a，b ] n 等分，记分点为 其中 称为步长，然后在每个小区间[ xi-1，xi ] 上应用梯形公式即 余项为： 就可导出复合梯形公式 若将所得积分近似值记成Tn ，并注意到x0=a，xn=b，则上式即为 余项为 仿上，可得复合辛普生公式 也可以像教材上介绍的将区间[a,b]分为2n等分，得到辛普生公式 例7.4.利用梯形公式和辛普生公式求 和复合柯特斯公式 例7.4。解： 若 f (x) 在积分区间 [ a，b ] 上分别具有二阶、四阶和六阶连续导数,则复合求积公式的余项分别为 定理 证明 只对复合梯形公式证明余项公式。 由于 在 上连续，故由定理知，对每个小区间上积 分 使用梯形公式时，所得近似值的误差为 故 即 由上式再根据已知条件容易判定复合梯形公式是收敛的。 而事实上只要 就可得到收敛性。 因为 由介值定理知，在 中必有点 ，使 故余项公式成立。 复化Simpson公式证明余项公式,见p345公式7.26 7.4 变步长求积公式 /* Composite Quadrature */ 复合求积公式的截断误差随着n的增大而减小，但是对于一 个给定的积分，选定了某种积分方法之后，如何选择适当的n， 使得计算结果达到预先选定的精度？由误差的要求可以定出n ， 但事先计算时很难估计：应边算边估计进而加密。 接下来以复合梯形求积公式为例，探讨如何避免事先误差 估计中的高阶导数的计算，以及事后误差估计中步长的选择问 题。 研究对象： 上的定积分 。 取2n等分计算 n等分计算 若n充分大， 在 上连续且 则 根据上式，当 时， 具体计算过程： 每次检验 当