[理学]数值-8.pptVIP

8.1 求积公式 8.2 牛顿-柯斯特（Newton-Cotes）公式 8.3 复化求积公式 8.4 龙贝格(Romberg)求积公式 8.5 高斯型求积公式 本章小结 当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强调为偶数份。 公式的特点为节点x2k,(k=1,2,…,m-1)作为小区间[x2k, x2k+2]的端点，参与前后两次的辛普生公式的计算，因而系数为2，而奇数节点x2k+1,(k=0,1,…,m-1)因辛普生公式中间点的求积系数为4而保留4，前面的h/3为辛普生公式的公共求积系数。 定理8.4 设函数 f(x)∈C4[a,b]，则 例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}序列。 {Rn}序列也称为龙贝格序列。 这样我们从收敛较慢的{Tn}序列只用了一些四则运算,便推出了收敛更快的{Sn}序列, {Cn}序列和{Rn}序列。我们常将这四个序列排成如下的三角形数表 表8-2 ﹕ ﹕ ﹕ ﹕ R2 C4 S8 T16 R1 C2 S4 T8 C1 S2 T4 S1 T2 T1 该表四 个序列 都是收 敛的。 例5 利用龙贝格方法计算 解:计算结果列如下表: 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=π相比较已有较好的精度。 8.5 高斯型求积公式 由前面的讨论已经知道，以a=x0x1…xn=b为节点的N-C求积公式的代数精度一般为n或n+1,这时节点简单地按照闭式等距的方式确定。对一个求积公式而言，如果不固定节点的位置，在节点数目不变的情况下，代数精度能否提高，最多能达到多少?高斯型求积公式讨论的就是最高代数精度的求积公式。先看一个简单的例子，考虑两 8.5.1 最高代数精度的求积公式 个节点的求积公式。 这里积分区间选为[-1,1] 不失一般性，因为对区间[a,b]，总可用变量替换 将化为[-1,1]上的积分。从而用（8.35）式得到解决。 现在不固定节点x0 , x1的位置为区间端点-1,1，而允许取在(-1,1)内，即选取x0 , x1及A0 ，A1，使(8.35)具有最高的代数精度。为此，分别取 代入(8.35)式，得方程组 取前四个方程构成此方程组是因为只有4个待定数 x0 , x1，A0, A1，即代数精度m=3。 方程组(8.36)有一组解为 即求积公式 具有的代数精度为3。