数值计算方法第2版 第6章 数值积分和数值微分.ppt

6.5 数值微分 当f(x)在[a,b]有连续的6阶导数时，复化柯特斯公式的误差 6.3 变步长求积和龙贝格算法 复化求积公式能提高精度，但要给出步长，步长太大精度低，步长太小，计算量大。实际计算用变步长计算，在步长逐次二分过程中，反复利用复化求积公式进行计算，直到所求积分值满足精度要求为止。 将积分区间等分成n个子区间，则有n+1个分点 对子区间再增加一个新节点 ，区间增加1倍，有 对子区间 运用梯形公式，有 6.3.1 变步长梯形法则 比较 取 （ 允许截断误差ε ） 在步长逐次二分的过程中，校验上式，取满足精度的 。 若将区间再分半，为 则有 6.3.2 龙贝格（Romberg)求积法 梯形法的加速 梯形法计算简单，精度较低，收敛慢，当把区间分成n等份，用复化公式计算积分的近似值为 ，截断误差为 当 时，T2n即为所求的近似值。 是T2n 的修正项，它与T2n 之和比T2n 、Tn更接近与真值，即它是一种补偿。 取 设f″(x)在 [a,b] 连续且变化不大时，有f″(ξn )≈ f″(ξ2n ) ,可得近似式 验后误差估计式 下面说明 将Tn ,T2n的表达式 代入，有 2 辛卜生法的加速 当把区间分成n等份，用复化辛卜生公式计算积分的近似值为 ，截断误差为 若将区间再分半，为 则有 设 连续且变化不大时，有 ,可得近似式 具有5次代数精度。 3 龙贝格公式 （柯特斯法的加速） 当把区间分成 n 等份，用复化柯特斯公式计算积分的近似值为 ，截断误差为 若将区间再分半，为 则有 设 连续且变化不大时，有 ,可得近似式 具有7次代数精度。 龙贝格积分法可以按下面表的顺序进行： 当对角线上最后两个相邻项满足 时,可停止计算并取 作为所求积分的近似值。 n ? ? ? ? ? ? ? ? 1 1/2 1/2 ? ? ? ? ? ? ? 2 1/6 4/6 1/6 ? ? ? ? ? ? 3 1/8 3/8 3/8 1/8 ? ? ? ? ? 4 7/90 16/45 2/15 16/45 7/90 ? ? ? ? 5 19/288 25/96 25/144 25/144 25/96 19/288 ? ? ? 6 41/840 9/35 9/280 34/105 9/280 9/35 41/840 ? ? 7 751/17280 3577/17280 1323/17280 2989/17280 2989/17280 1323/17280 3577/17280 751/17280 ? 8 989/ 28350 5888/ 28350 -928/ 28350 10496/ 28350 -4540/ 28350 10496/ 28350 -928/ 28350 5888/ 28350 989/ 28350 柯特斯求积系数表： 例如：n=1时，有 n=2时，有 柯特斯系数的性质 (2) 系数有对称性。 (3) 当n≥8时开始出现负值的柯特斯系数。 (1) 取f(x)≡1，则 f(n+1)(x)≡0, Rn(f)≡0，于是 梯形公式 当n=1时, 有 相当于用直线P(x)代替f(x)计算积分。 3 常用的低阶牛顿-柯特斯公式 抛物线(辛卜生)公式 牛顿－柯特斯求积公式 当n=2时有 相当于用过两个端点和中点的二次 抛物线P(x)代替f(x)计算积分。 辛卜生公式的几何意义 柯特斯公式 牛顿－柯特斯求积公式 当n=4时有 6.2.2 牛顿－柯特斯公式的代数精度 当 f(x)是 1,x,x2,…, xm 时，准确成立，但当f(x)= xm+1时，不准确成立，则称求积公式的代数精确度（简称代数精度）为m。 复习 定义 求积公式 （Ai与f(x)无关）