第四章数值积分与数值微分.pptVIP

?收敛速度与误差估计【定义4.2】（p阶收敛）：若一个积分公式的误差满足，且，则称该公式是p阶收敛的。~~~**第60页，共106页，星期日，2025年，2月5日【例】利用数据表01/83/81/25/83/47/811/422.265492.460002.876403.200003.506853.764703.938464计算积分解：这个问题有明显的答案**第61页，共106页，星期日，2025年，2月5日取n=8用复化梯形公式=3.138988494取n=4用复化Simpson公式=3.141592502运算量基本相同**第62页，共106页，星期日，2025年，2月5日复化梯形公式的误差估计给定精度，如何取?【例】要求，如何判断n=？1、误差先验估计式记则**第63页，共106页，星期日，2025年，2月5日?上例中若要求，则即：取n=409通常采取将区间不断对分的方法，即取n=2k上例中2k?409?k=9时，T512=32、误差后验估计式**第64页，共106页，星期日，2025年，2月5日复化Simpson公式的误差估计1、误差先验估计式2、误差后验估计式**第65页，共106页，星期日，2025年，2月5日复化Cotes公式的误差估计1、误差先验估计式2、误差后验估计式**第66页，共106页，星期日，2025年，2月5日§4.3Romberg算法一、梯形法的递推化(变步长)二、龙贝格求积公式三、理查德森外推法**第67页，共106页，星期日，2025年，2月5日1、首先将区间n等分：2、再将区间2n等分，即步长减半：一、梯形法的递推化(变步长)基本思想：将积分区间逐次分半，比较前后两次的近似值终止法则：前后两次近似值的误差小于已知精度**第68页，共106页，星期日，2025年，2月5日上述条件满足，程序终止；否则，继续分半计算。3、终止条件：由复化梯形公式的余项知:变化不大时由此得到近似关系式:误差控制条件:*第69页，共106页，星期日，2025年，2月5日P85页例：根据如下函数值表，利用复化梯形公式计算积分的近似值，要求误差不超过。01/81/43/810.9973980.9896880.9767271/25/86/87/810.9588510.9361560.9088580.8771930.841471**第70页，共106页，星期日，2025年，2月5日然后将区间二等分，利用递推公式求出*递推公式进一步二分积分区间，类似可求出如此不断二分并利用递推公式，可得下表中的结果解：先在整个区间上用梯形公式第71页，共106页，星期日，2025年，2月5日k表示二分次数，区间数1234567890.93979330.94451350.94569090.94598500.94605960.94607690.94608150.94608270.9460830由表中可以看出，对分8次和对分7次之间的差因而是满足精度要求的解。误差控制条件**第72页，共106页，星期日，2025年，2月5日二、龙贝格求积公式1.龙贝格积分思想由上节分析知，用复化梯形公式计算积分值,的误差大约为：令由复化梯形公式知:误差补偿**第73页，共106页，星期日，2025年，2月5日令**第74页，共106页，星期日，2025年，2月5日梯形公式的加速方法：利用复化梯形公式前后两次的积分近似值和，按照上式作出的线性组合是具有更高精度的积分值。龙贝格积分公式正是由此产生！上述公式说明：**第75页，共106页，星期日，2025年，2月5日**第28页，共106页，星期日，2025年，2月5日二、截断误差与代数精度1、截断误差*