计算方法龙贝格高斯求积公式

计算方法龙贝格高斯求积公式第一页，共四十页，2022年，8月28日

综合前几节的内容,我们知道梯形公式,Simpson公式,Cotes公式的代数精度分别为1次,3次和5次复合梯形、复合Simpson、复合Cotes公式的收敛阶分别为2阶、4阶和6阶无论从代数精度还是收敛速度,复合梯形公式都是较差的有没有办法改善梯形公式呢?2023/2/212第二页，共四十页，2022年，8月28日一、复合梯形公式的递推化各节点为复合梯形(Trapz)公式为--------(1)--------(2)2023/2/213第三页，共四十页，2022年，8月28日--------(3)由(1)(2)两式可42023/2/21(3)式称为递推的梯形公式递推梯形公式加上一个控制精度,即可成为自动选取步长的复化梯形公式优点：梯形法计算简单缺点：收敛慢，为了达到要求的精度，需要二分区间多次，分点大量增加，计算量很大第四页，共四十页，2022年，8月28日--------(3)2023/2/215第五页，共四十页，2022年，8月28日则由(1)(2)(3)式,有2023/2/216第六页，共四十页，2022年，8月28日因此(1)(2)(3)式可化为如下递推公式-------(4)上式称为复合变步长梯形求积公式2023/2/217第七页，共四十页，2022年，8月28日二、外推加速公式由复合梯形公式的余项公式可得由(3)式移项合并2023/2/218第八页，共四十页，2022年，8月28日复合Simpson公式表明：由梯形公式前后两次结果的线性组合可构造出精确度较高的辛普森公式2023/2/219第九页，共四十页，2022年，8月28日--------(5)即--------(6)当然2023/2/2110第十页，共四十页，2022年，8月28日11因此由复合Simpson公式的余项可得令--------(7)--------(8)即当然2023/2/21第十一页，共四十页，2022年，8月28日由复合Cotes公式的余项得令--------(9)即当然公式(9)称为龙贝格积分公式2023/2/2112第十二页，共四十页，2022年，8月28日外推加速公式以上整个过程称为Romberg算法将上述结果综合后2023/2/2113第十三页，共四十页，2022年，8月28日其中外推加速公式可简化为--------(9)Romberg算法的收敛阶高达m+1的两倍Romberg算法求解步骤Romberg算法的代数精度为m的两倍2023/2/2114第十四页，共四十页，2022年，8月28日上述公式推导说明，T1公式是梯形公式，对于次数不高于1的多项式准确；S1是辛普森公式，对于次数不高于3的多项式准确；C1是柯特斯公式，它对于次数不高于5的多项式准确，每一个公式均由前一公式的适当线性组合得到，精确度都提高2次。因此可以验证，由柯特斯公式C1构造得到的龙贝格公式R1，对次数不高于7次的多项式准确。龙贝格公式计算积分占用内存少，精度高2023/2/2115第十五页，共四十页，2022年，8月28日龙贝格公式计算步骤：其计算过程是将区间逐次分半，加速得到积分近似值，因此称为逐次分半加速法2023/2/2116第十六页，共四十页，2022年，8月28日

00.920735510.93979330.946145920.94457350.94608690.946083030.94569090.94608330.94608310.9460831例1

将三个加速公式用于求从表中可以看出三次加速求得R1=0.9460831每位数字都是有效数字172023/2/21第十七页，共四十页，2022年，8月28日例2.182023/2/21解：（1）在区间[1,2]上用梯形公式得第十八页，共四十页，2022年，8月28日192023/2/21第十九页，共四十页，2022年，8月28日202023/2/21第二十页，共四十页，2022年，8月28日本章作业P9710212023/2/21第二十一页，共四十页，2022年，8月28日

牛顿—柯特斯型求积公式是封闭型的（区间[a,b]的两端点a,b均是求积节点）而且要求求积节点是等距的,受此限制，牛顿—柯特斯型求积公式的代数精度只能是n(n为奇数)或n+1(n为偶数)。

4.7.1一般理论4.7高斯求积公式222023/2/21问题：对求积节点也适当的选取,即在求积公式中不仅Ak而且xk也加以选取,是否可以增加自由度，从而可提高求积公式的代数精度？

第二十二页，共四十页，2022年，8月28日2023/2/2123第二十三页，共四十页，2022年，8月28日第二十四页，共四十页，2022年，8月28日2023/2/2125第二十五页，共四十页，2022年，8月28日2023/2/21定义：26第二十六页，共四十页，2022年，8月28日2023/2/21

4.7.2Gauss求积公式的构造先从简单情况入手27第二十七页，共四十页，2022年，8月28日2023/2/21①②28第二十八页，共四十页，2022年，8月28日2023/2/2129第二十九页，共四十页，2022年，8月28日2023/2/2130第三十页，共四十页，2022年，8月28日则求积公式为此式对任何不高于3次的多项式f(x)都准确成立2023/2/2131第三十一页，共四十页，2022年，8月28日322023/2/21第三十二页，共四十页，2022年，8月28日332023/2/21第三十三页，共四十页，2022年，8月28日342023/2/21第三十四页，共四十页，2022年，8月28日

4.7.3常用的高斯型求积公式352023/2/21由于正交多项式随权函数不同，所以取不同的权函数时，则有各种各样的Gauss型求积公式.第三十五页，共四十页，2022年，8月28日（1）高斯－勒让德求积公式2023/2/2136第三十六页，共四十页，2022年，8月28日372023/2/2