[理学]第4章_数值积分和数值微分.ppt

4.1.4 求积公式的余项 4.1.5 求积公式的收敛性与稳定性 4.2 牛顿—柯特斯公式 即 又称为逐次分半外推加速求积法，简称外推加速法. 也称为龙贝格求积算法. 以Ｔ0(k)表示二分k次后求得的梯形值, 以Ｔm(k)表示序列{Ｔ0(k)}的m次加速值， 龙贝格求积算法的计算过程如下： (1) 取k=0,h=b-a，求 令1→k(k记区间 [a, b]的二分次数). (2) 求值 ，按梯形递推公式计算Ｔ0(k) . (3) 求加速值，按加速公式逐个求出Ｔ数表的第k行其余各元素Ｔj(k-j) (j=1,2,?,k). (4) 若|Ｔk(0) -Ｔk-1(0)|ε(预先给定的精度)，则终止计算，并取Ｔk(0)≈I; 否则令k+1→k转(2)继续计算. Ｔ数表 ? ? ? ? ? ? ? T4(0) T3(1) T2(2) T1(3) T0(4) 4 T3(0) T2(1) T1(2) T0(3) 3 T2(0) T1(1) T0(2) 2 T1(0) T0(1) 1 T0(0) 0 T4(k) T3(k) T2(k) T1(k) T0(k) k 注意计算顺序，第k步子区间长度为h=(b-a)/2k. 可以证明，如果f(x)充分光滑，那么T数表每一列的元素及对角线元素均收敛到所求的积分值I，即 对于f(x)不充分光滑的函数也可以用龙贝格算法计算，只是收敛慢一些，这时也可以直接使用复合辛普森公式计算.课本例6. 4.6 高斯求积公式 由前面的讨论已经知道，以a=x0x1?xn=b为节点的N-C求积公式的代数精度一般为n或n+1，这时节点简单地按照闭区间等距的方式确定。 对一个求积公式而言，如果不固定节点的位置，在节点数目不变的情况下，代数精度能否提高，最多能达到多少? 这里高斯型求积公式讨论的就是最高代数精度的求积公式. 4.6.1 一般理论 前面给出的机械求积公式 含有2n+2个待定参数xk,Ak(k=0,1,?,n). 当xk为等距节点时得到的插值求积公式其代数精度至少为n次，如果适当选取xk(k=0,1,?,n)，有可能使求积公式具有2n+1次代数精度，这类求积公式称为高斯(Gauss)求积公式. 例8 对于求积公式 试确定节点 及系数，使其具有尽可能高的代数精度。 这里ρ(x)为权函数，类似的求积公式为 在这个求积公式里Ak(k=0,1,?,n)为不依赖于f(x)的求积系数， xk(k=0,1,?,n)为求积节点，可适当选取xk及Ak(k=0,1,?,n)使(6.1)式具有2n+1次代数精度. 定义4 如果求积公式(6.1)具有2n+1次代数精度，则称其节点xk(k=0,1,?,n)称为高斯点，相应公式(6.1)称为高斯(Gauss)求积公式. 根据定义要使 (6.1)具有2n+1次代数精度，只要取f(x)=xm，对m=0,1,?,2n+1，精确成立，则得 当给定权函数ρ(x)，求出右端积分，则可由(6.2)式解得xk及Ak(k=0,1,?,n). 求解非线性方程组(6.2)较复杂，通常n1就很难求解，故一般不通过解方程(6.2)求解xk及Ak(k=0,1,?,n)，而从分析高斯点的特性来构造高斯求积公式. 定理5 对插值型求积公式 与一切次数不超过n的多项式p(x)带权ρ(x)正交, 即 节点xk (k=0,1,2,?,n) 是高斯点的充要条件是这些节点为零点的多项式 证明 (必要性) 设p(x)∈Hn, 则p(x)ωn+1(x)∈H2n+1, 因此，如果xk (k=0,1,2,?,n)是高斯点，则求积公式(5.1)对于p(x)ωn+1(x)精确成立，即有 因ωn+1(xk)=0 (k=0,1,2,?,n)，故结论成立. (充分性) 对任意f(x)∈H2n+1, 用ωn+1(x)去除 f(x), 则可表示成 其中，q(x)为商式，r(x)为余式，q(x),r(x)均为不超过n次的多项式，于是有 由已知得 于是有 由于(6.1)是插值型求积公式，故对r(x)∈Hn精确成立 再注意到在节点处ωn+1(xk)=0 (k=0,1,2,?,n)，知道有等式r(xk)= f(xk) (k=0,1,2,?,n)，从而有 可见积分公式(6.1)对一切次数不超过2n+1的多项式均精确成立. 因此, xk(k=0,1,2,?,n)为高斯点. 证毕.