数值分析第8章-2.pptVIP

* 结束 §8.3 复化求积公式 * 从求积公式的余项的讨论中我们看到，被积函数所用的插值多项式次数越高，对函数光滑性的要求也越高.另一方面，插值节点的增多(n的增大)，在使用牛顿-柯特斯公式时将导致求积系数出现负数(当n≥8时,牛顿.柯特斯求积系数会出现负数)因而在实际应用中往往采用将积分区间划分成若干个小区间，在各小区间上采用低次的求积公式(梯形公式或抛物形公式)，然后再利用积分的可加性，把各区间上的积分加起来，便得到新的求积公式，这就是复化求积公式的基本思想.为叙述方便，我们仅讨论各小区间均采用同一低次的求积公式的复化求积公式对各小区间也可分别采用不同的求积公式，也可推出新的求积公式，读者可按实际问题的具体情况讨论 结束 8.3.1. 复化梯形公式 用n+1个分点将区间[a，b]n等分。每个区间长 在[xk，xk+1]上用梯形公式，则 * Tn叫做复化梯形求积公式，下标n表示将积分区间等分的份数. 结束 从公式的特点可以看出，内节点xk(k=1,2,…n-1)作为小区间的端点参与前、后两个小区间的计算，因而系数为2，端点a与b只参与一次计算，系数为1. * 如果在Tn的基础上，将各小区间对分，这时节点数为2n+1，分段数为2n.记新的分点的函数值的和为σn，则T2n应为原内节点与新增节点函数值的和的两倍，加上两端点a,b的函数值之和再乘上新区间长度的一半，即 结束 定理8.3 设 f (x)∈C2[a,b]，复化梯形公式的截断误差 (证明) 这一复化梯形求积公式的余项在形式上与(8.13)式相同，不同的是，这里的h=(b-a)/n，而(8.13)式中的h=b-a. * 从这一公式可以看出，将区间对分后，原复化梯形公式的值Tn作为一个整体保留.只需计算出新分点的函数值，便可得出对分后的积分值，不需重复计算原节点的函数值，从而减少了计算量. 利用复化求积公式的余项，我们可以估计出在满足精度的要求下，应将积分区间等分多少份，即n取多少.这种误差估计方法称为事前误差估计.如例8.3 例3 利用复化梯形公式计算 使其误差限为10-4，应将区间［0,1］几等分? 结束 解: 因为被积函数 * 取n=17可满足要求. 另一方法是利用公式前后两次计算结果的差来估计误差的，即用｜T2n-Tn｜＜ε,这是因为 结束 * 当f ″(x)在[a,b]上连续，并且假定当n充分大时有f″(ξ)≈f″(η)，则 结束 * 这种误差估计方法通常叫做事后误差估计，在计算机上用来控制计算精度常用这一方法，有的也把这种方法叫做步长的自动选取或逐次对分的方法. 因此当｜T2n-Tn｜＜ε时，可认为 结束 * 将积分区间[a,b]2m等分，n=2m，节点为xk=a+kh(k=0,1,2,…,2m)，h=(b-a)/2m.在每两个小区间［x2k,x2k+2］(k=0,1,2,…,m-1)上用 抛物形公式，则有: 8.3.2 复化抛物形公式 S2m叫做复化抛物形求积公式，下标2m表示积分区间等分的份数，2m强调为偶数份. 结束 * 定理8.4 设函数f(x)∈C4[a,b]，则 公式的特点为节点x2k,(k=1,2,…,m-1)作为小区间[x2k, x2k+2]的端点，参与前后两次的辛普生公式的计算，因而系数为2，而奇数节点x2k+1,(k=0,1,…,m-1)因辛普生公式中间点的求积系数为4而保留4，前面的h/3为辛普生公式的公共求积系数. 例4 利用复化抛物形公式计算 使其误差限为10-4，应将区间［0,1］几等分? 解:利用例3的结果 因此只需将区间［0,1］二等分，即取m=1(n=2). 结束 * 前面用复化梯形公式计算此题，满足相同的精度需要将区间 ［0,1］17等分，可见复化抛物形公式的精度的确比复化梯形公式精度高同样也可用 ｜S4m-S2m｜ε 来控制计算的精度. 结束 * §8.4 龙贝格(Romberg)求积公式 我们已知的T2n与Tn的关系 8.4.1 复化梯形公式的逐次分半公式 于是可以逐次对分形成一个序列{T1,T2,T4,T8,…},此序列收敛于积分真值I.当|T2n-Tn|ε时,取T2n为I的近似值.以上算法称为复化梯形公式的逐次分半公式.但由于此序列收敛太慢,因此并不实用.现我们试图将它改造成为收敛快的序列. 结束 * 如认为 则有 于是有: 记 这样我们从收敛较慢的{Tn}序列推出了收敛较快的{Sn}序列. 可以证明{Sn}序列实际上就是逐次分半的复化抛物形公式序列. 如认为 则有 于是有: 记 这样我们从{Sn}序列又推出了收敛更快的{Cn}序列. 可以证明{Cn}序列实际上就是逐次分半的复化柯特斯公式序列. 结束 * 如认为 则有 于是有: 记 这样我们从{Cn}序列又推出了收敛更快的{Rn}