数值分析ppt第4章_数值积分与数值微分.ppt

4.1.4 求积公式的收敛性与稳定性 4.2 牛顿—柯特斯公式 由书例2，所求得的两个梯形值T4=0.9445135和T8=0.9456909的精度都很差(与准确值I=0.9460831比较，只有两、三位有效数字)，但如果将它们按上式做线性组合，则新的近似值 却有6位有效数字. 可以直接验证 就是复化辛普森积分公式. Sn的精度为O(h4). 这就是说，用复化梯形法二分前后的两个积分值Tn与T2n ，按上式做线性组合，结果得到了复化辛普森积分公式. 则 同理由辛普森法，用二分前后的两个积分值Sn与S2n，由误差公式即有 不难直接验证就是复化柯特斯积分公式. Cn的精度为O(h6). 则 同理由柯特斯法，用二分前后的两个积分值Cn与C2n，由误差公式即有 这就是复化龙贝格积分公式. Rn的精度为O(h8). 一般我们将这种龙贝格算法做成表格 我们在变步长的过程中运用了三个公式，就能将粗糙的梯形值Tn 逐步加工成精度较高的辛普森值Sn、柯特斯值Cn和龙贝格值Rn. ﹕ ﹕ ﹕ ﹕ R2 C4 S8 T16 R1 C2 S4 T8 C1 S2 T4 S1 T2 T1 见书p135例题3. 例4 利用龙贝格方法计算 3.14159 3.14159 3.14159 3.14094 16 4 3.14159 3.14159 3.14159 3.13899 8 3 3.14212 3.14157 3.13118 4 2 3.13333 3.10000 2 1 3.00000 1 0 R序列 C序列 S序列 T序列 2i i 这一结果与I=π相比较已有较好的精度. 解 计算结果列如下表: 4.4.3 理查森外推加速法 上面讨论说明由梯形公式出发，将区间[a, b]逐次二分可提高求积公式的精度，上述加速过程还可继续下去，其理论依据是梯形公式的余项展开，设 若记Tn=T(h)，当区间[a, b]划分为2n等分时，则有 并且有 可以证明梯形公式余项可展开成级数形式，即 定理4 设f(x)?C∞[a, b]，则有 式中I为积分值，系数?k与h 无关.误差量级为O(h2). 此定理可利用f(x)的泰勒展开推导得到，证略. 定理4表明T(h)≈I是O(h2)阶，若h/2用代替h, 有 用4乘此式，减去上式再除3记为T1(h)，则得 改记为 这里系数βk与h无关，这样构造的T1(h)与积分值I近似的阶为O(h4). 比较T1(h)与Sn可知, 这样构造的序列T1(h),T1(h/2),?.就是辛普森公式序列Sn, S2n,?. 根据 令 则又可进一步从余项展开式中消去h4项, 从而有 这样构造出的{T2(h)}，其实就是柯特斯公式序列，它与积分值I的逼近阶为O(h6). 如此推下去，每加速一次，误差的量级便提高2阶，速度较快，一般地，若记T0(h)=T(h)，则有 误差量级为O(h6)． 误差量级为O(h4)． 如此继续下去，可得 用Ｔm(h)作为I 的近似值, 误差量级为O(h2(m+1)). 经过m(m=1,2,?)次加速后，余项便取下列形式用： 这种处理方法通常称为理查森(Richardson)外推加速方法. 即 又称为逐次分半外推加速求积法，简称外推加速法. 也称为龙贝格求积算法. 以Ｔ0(k)表示二分k次后求得的梯形值, 以Ｔm(k)表示序列{Ｔ0(k)}的m次加速值， 龙贝格求积算法的计算过程如下： (1) 取k=0,h=b-a，求 令1→k(k记区间 [a, b]的二分次数). (2) 求值 ，按梯形递推公式计算Ｔ0(k) . (3) 求计算值，按加速公式逐个求出Ｔ数表的第k行其余各元素Ｔj(k-j) (j=1,2,?,k). (4) 若|Ｔk(0) -Ｔk-1(0)|ε(预先给定的精度)，则终止计算，并取Ｔk(0)≈I; 否则令k+1→k转(2)继续计算. Ｔ数表 ? ? ? ? ? ? ? T4(0) T3(1) T2(2) T1(3) T0(4) 4 T3(0) T2(1) T1(2) T0(3) 3 T2(0) T1(1) T0(2) 2 T1(0) T0(1) 1 T0(0) 0 T4(k-4) T3(k-3) T2(k-2) T1(k-1) T0(k) k 注意计算顺序，第k步子区间长度为h=