数值分析课件.pptVIP

前面用复化梯形公式计算此题，满足相同的精度需要将区间 ［0,1］17等分，可见复化抛物形公式的精度的确比复化梯形公式精度高同样也可用  ｜S4m-S2m｜ε 来控制计算的精度. §4.4 龙贝格(Romberg)求积公式 4.4.1 复化梯形公式的逐次分半公式 我们已知的T2n与Tn的关系 结束 于是可以逐次对分形成一个序列{T1,T2,T4,T8,…},此序列收敛 于积分真值 I.当 |T2n-Tn|ε时,取T2n为I的近似值.以上算法称 为复化梯形公式的逐次分半公式. 但由于此序列收敛太慢 ,因 此并不实用.现我们试图将它改造成为收敛快的序列. 如认为 则有 结束 于是有: 记 这样我们从收敛较慢的{Tn}序列推出了收敛较快的{Sn}序列. 可以证明{Sn}序列实际上就是逐次分半的复化抛物形公式序列. 如认为 则有 于是有: 记 这样我们从{Sn}序列又推出了收敛更快的{Cn}序列. 可以证明{Cn}序列实际上就是逐次分半的复化柯特斯公式序列. 结束 如认为 则有 于是有: 记 这样我们从{Cn}序列又推出了收敛更快的{Rn}序列. {Rn}序列也称为龙贝格序列.这样我们从收敛较慢的{Tn}序列只用了一些四则运算,便推出了收敛更快的{Sn}序列, {Cn}序列和{Rn}序列. 这个过程还可继续下去,但已意义不大.我们常将这四个序列排成如下的三角形数表(表4-2) 结束 ﹕ ﹕ ﹕ ﹕ R2 C4 S8 T16 R1 C2 S4 T8 C1 S2 T4 S1 T2 T1 表4-2 该表四个序列都是收敛的. 结束 例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=π相比较已有较好的精度. 结束 §4.5 高斯型求积公式 由前面的讨论已经知道，以a=x0x1…xn=b为节点的N-C求积公式的代数精度一般为n或n+1,这时节点简单地按照闭式等距的方式确定。对一个求积公式而言，如果不固定节点的位置，在节点数目不变的情况下，代数精度能否提高，最多能达到多少?高斯型求积公式讨论的就是最高代数精度的求积公式.先看一个简单的例子，考虑两个节点的求积公式. 4.5.1 最高代数精度的求积公式 结束 * §4.1求积公式 4.1.1 求积公式 结束 对定义在区间［a,b］上的定积分 以上公式多称为牛顿-莱布尼兹公式，F(x)为f(x)的原函数.但有时原函数不能用初等函数表示，有时原函数又十分复杂，难于求出或计算.如被积函数为: 第四章 数值积分 等函数的积分都无法解决，当被积函数为一组数据时，更是无能为力. 为解决定积分的近似计算,从定积分的定义： 这样就避开了求原函数的运算.(4.1)式就叫做求积公式，Ak(k=0,1,…,n)与函数f(x)无关，叫做求积系数,显然要确定一个求积公式，要确定求积结点xk和求积系数Ak，或者说不同的求积结点和求积系数将确定不同的求积公式. 结束 结束 4.1.2 求积公式的余项和代数精度 一般情况下，(4.1)两端并不相等.我们称: (4.2)为求积公式(4.1) 的余项，或截断误差. 为考查一个求积公式的误差，通常用代数精度来表示，如果一个求积公式对于不超过m次的多项式都能够精确成立(R[f]≡0)，而对m+1次以上的多项式不能精确成立，则称该求积公式的代数精度为m. 结束 例如求积公式： 验证当 f(x)=xm，m=0,1,2,3,4 时，是否有R[xm]=0 所以以上求积公式的代数精度为 3. 任何一个求积公式的代数精度至少为零即取f(x)=1时公式应精确成立，这是求积系数应满足的起码条件，可以用它检验一个求积公式的系数的正确性. 4.1.3 矩形求积公式 f(x)= f(a)+ f ’(ξ)(x-a)， ξ在x,a之间,两端积分： 把 f(x)在a处作Taylor展开: 结束 结束 注意到右端第二项积分，设f ′(x)在［a,b］上连续，而x-a在 ［a,b］上不变号(非负)，据积分中值定理有： 于是有左矩形公式： 同理 , f(x)在b点展开,可得右矩形公式: 结束 f(x)在中点(a+b)/2展开,可得中矩形公式 : 不难验证，(4.3)和(4.4)具有零次代数精度，(4.5)具有一次代数精度. 结束 4.1.4 内插求积公式 由插值可知，对任一函数