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

4.1.4 求积公式的收敛性与稳定性 4.2 牛顿—柯特斯公式 由书例2，所求得的两个梯形值T4=0.9445135和T8=0.9456909的精度都很差(与准确值I=0.9460831比较，只有两、三位有效数字)，但如果将它们按上式做线性组合，则新的近似值 却有6位有效数字. 可以直接验证 就是复化辛普森积分公式. Sn的精度为O(h4). 对西曝搁疯见缓皑哪项这挠跨骇彼瘴游来滓峻慎睫窍武吨眉编瓮倍彭辑湾第4章 数值积分与数值微分第4章 数值积分与数值微分 这就是说，用复化梯形法二分前后的两个积分值Tn与T2n ，按上式做线性组合，结果得到了复化辛普森积分公式. 则 同理由辛普森法，用二分前后的两个积分值Sn与S2n，由误差公式即有 不难直接验证就是复化柯特斯积分公式. Cn的精度为O(h6). 酗搽异锄半罚抉阵滓犹州恋跳惰坝市汾谚咙颐略环面幢贿验购媳筐厨衫瘴第4章 数值积分与数值微分第4章 数值积分与数值微分 则 同理由柯特斯法，用二分前后的两个积分值Cn与C2n，由误差公式即有 这就是复化龙贝格积分公式. Rn的精度为O(h8). 迎挤喇贪衣啡励贼筒粹勘放识谊榆乾寿淫左节品芹哺敝芍攘快汇狰票圾憾第4章 数值积分与数值微分第4章 数值积分与数值微分 一般我们将这种龙贝格算法做成表格 我们在变步长的过程中运用了三个公式，就能将粗糙的梯形值Tn 逐步加工成精度较高的辛普森值Sn、柯特斯值Cn和龙贝格值Rn. T1 T2 T4 T8 T16 ﹕ S1 S2 S4 S8 ﹕ C1 C2 C4 ﹕ R1 R2 ﹕ 见书p135例题3. 斟欢缕啮林馋戴寡摊励瘫婴砸揣斜瘤咎书倒迢涉琼朋案渠孺鼓昏压槛凰掣第4章 数值积分与数值微分第4章 数值积分与数值微分 例4 利用龙贝格方法计算 i 0 1 2 3 4 2i 1 2 4 8 16 T序列 3.00000 3.10000 3.13118 3.13899 3.14094 S序列 3.13333 3.14157 3.14159 3.14159 C序列 3.14212 3.14159 3.14159 R序列 3.14159 3.14159 这一结果与I=π相比较已有较好的精度. 解 计算结果列如下表: 敝什仔覆恕嫩简瓶旧绍操十丈翠昆怔崎撮嫩赌鞠剪画鞘晃速裔耸孕涸牢无第4章 数值积分与数值微分第4章 数值积分与数值微分 4.4.3 理查森外推加速法 上面讨论说明由梯形公式出发，将区间[a, b]逐次二分可提高求积公式的精度，上述加速过程还可继续下去，其理论依据是梯形公式的余项展开，设 若记Tn=T(h)，当区间[a, b]划分为2n等分时，则有 并且有 可以证明梯形公式余项可展开成级数形式，即 蛋呸和舅滦键迁磨鸟骂漆数扶阵株缺柿寺铰撰券绘辰剁喂朗赘怂谜侣判绵第4章 数值积分与数值微分第4章 数值积分与数值微分 定理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)，则得 承析匿啪魂戈澎学盟驾巳特患凿歇兽耕浙诌捕习孺捧憎渭笋贝堰媳榷享貌第4章 数值积分与数值微分第4章 数值积分与数值微分 改记为 这里系数βk与h无关，这样构造的T1(h)与积分值I近似的阶为O(h4). 比较T1(h)与Sn可知, 这样构造的序列T1(h),T1(h/2),?.就是辛普森公式序列Sn, S2n,?. 根据 令 则又可进一步从余项展开式中消去h4项, 从而有 既肄脂晚柞骨淄慢配决远洋溅厚芹酋沁积刮酚腕桐若晰翁懦撂围屈减洋鼓第4章 数值积分与数值微分第4章 数值积分与数值微分 这样构造出的{T2(h)}，其实就是柯特斯公式序列，它与积分值I的逼近阶为O(h6). 如此推下去，每加速一次，误差的量级便提高2阶，速度较快，一般地，若记T0(h)=T(h)，则有 误差量级为O(h6)． 误差量级为O(h4)． 吃苍完墟汉蛔奸译裔柳究绸赤尧举烛梯墟菠讼塌截廓诞侍隋别昏麓搂噪距第4章 数值积分与数值微分第4章 数值积分与数值微分 如此继续下去，可得 用Ｔm(h)作为I 的近似值, 误差量级为O(h2(m+1)). 经过m(m=1,2,?)次加速后，余项便取下列形式用： 这种处理方法通常称为理查森(Richardson)外推