《计算方法龙贝格高斯求积公式课件.ppt》由会员分享，可在线阅读，更多相关《计算方法龙贝格高斯求积公式课件.ppt（40页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、计算方法龙贝格高斯求积计算方法龙贝格高斯求积公式公式第1页，此课件共40页哦 综合前几节的内容,我们知道梯形公式,Simpson公式,Cotes公式的代数精度分别为1次,3次和5次复合梯形、复合Simpson、复合Cotes公式的收敛阶分别为2阶、4阶和6阶无论从代数精度还是收敛速度,复合梯形公式都是较差的有没有办法改善梯形公式呢?2023/4/102第2页，此课件共40页哦一、复合梯形公式的递推化各节点为复合梯形(Trapz)公式为-(1)-(2)2023/4/103第3页，此课件共40页哦-(3)由(1)(2)两式可42023/4/10(3)式称为递推的梯形公式递推梯形公式加上一个控制精度
2、,即可成为自动选取步长的复化梯形公式优点：梯形法计算简单缺点：收敛慢，为了达到要求的精度，需要二分区间 多次，分点大量增加，计算量很大第4页，此课件共40页哦-(3)2023/4/105第5页，此课件共40页哦则由(1)(2)(3)式,有2023/4/106第6页，此课件共40页哦因此(1)(2)(3)式可化为如下递推公式-(4)上式称为复合变步长梯形求积公式2023/4/107第7页，此课件共40页哦二、外推加速公式由复合梯形公式的余项公式可得由(3)式移项移项合并合并2023/4/108第8页，此课件共40页哦复合Simpson公式表明：由梯形公式前后两次结果的线性组合可构造出精确度较高的
3、辛普森公式2023/4/109第9页，此课件共40页哦-(5)即-(6)当然2023/4/1010第10页，此课件共40页哦11因此由复合Simpson公式的余项可得令-(7)-(8)即当然2023/4/10第11页，此课件共40页哦由复合Cotes公式的余项得令-(9)即当然公式(9)称为龙贝格龙贝格积分公式积分公式2023/4/1012第12页，此课件共40页哦外推加速公式以上整个过程称为Romberg算法将上述结果综合后2023/4/1013第13页，此课件共40页哦其中外推加速公式可简化为-(9)Romberg算法的收敛阶高达m+1的两倍Romberg算法求解步骤Romberg算法的代
4、数精度为m的两倍2023/4/1014第14页，此课件共40页哦上述公式推导说明，T1公式是梯形公式，对于次数不高于1的多项式准确；S1是辛普森公式，对于次数不高于3的多项式准确；C1是柯特斯公式，它对于次数不高于5的多项式准确，每一个公式均由前一公式的适当线性组合得到，精确度都提高2次。因此可以验证，由柯特斯公式C1构造得到的龙贝格公式R1，对次数不高于7次的多项式准确。龙贝格公式计算积分占用内存少占用内存少，精度高精度高2023/4/1015第15页，此课件共40页哦龙贝格公式计算步骤：龙贝格公式计算步骤：其计算过程是将区间逐次分半，加速得到积分近似值，因此称为逐次分半加速法逐次分半加速法
5、2023/4/1016第16页，此课件共40页哦 00.920735510.93979330.946145920.94457350.94608690.946083030.94569090.94608330.94608310.9460831例1 将三个加速公式用于求从表中可以看出三次加速求得R1=0.9460831每位数字都是有效数字172023/4/10第17页，此课件共40页哦例2.182023/4/10解：（1）在区间1,2上用梯形公式得第18页，此课件共40页哦192023/4/10第19页，此课件共40页哦202023/4/10第20页，此课件共40页哦本章作业P97 10212023
6、/4/10第21页，此课件共40页哦 牛顿柯特斯型求积公式是封闭型的（区间区间a,ba,b的的两端点两端点a,ba,b均是求积节点均是求积节点）而且要求求积节点是等距的,受此限制，牛顿柯特斯型求积公式的代数精度只能是n(n为奇数)或n+1(n为偶数)。4.7.1 一般理论一般理论222023/4/10问题：对求积节点也适当的选取,即在求积公式中不仅Ak而且xk也加以选取,是否可以增加自由度，从而可提高求积公式的代数精度？第22页，此课件共40页哦2023/4/1023第23页，此课件共40页哦第24页，此课件共40页哦2023/4/1025第25页，此课件共40页哦2023/4/10定义：26
7、第26页，此课件共40页哦2023/4/10 4.7.2 Gauss求积公式的构造求积公式的构造先从简单情况入手先从简单情况入手27第27页，此课件共40页哦2023/4/1028第28页，此课件共40页哦2023/4/1029第29页，此课件共40页哦2023/4/1030第30页，此课件共40页哦则求积公式为此式对任何不高于3次的多项式f(x)都准确成立2023/4/1031第31页，此课件共40页哦322023/4/10第32页，此课件共40页哦332023/4/10第33页，此课件共40页哦342023/4/10第34页，此课件共40页哦 4.7.3 常用的常用的高斯型求积公式高斯型求积公式352023/4/10由于正交多项式随权函数不同，所以取不同的权函数时，则有各种各样的Gauss型求积公式.第35页，此课件共40页哦（1）高斯勒让德求积公式2023/4/1036第36页，此课件共40页哦372023/4/10（2）高斯切比雪夫求积公式第37页，此课件共40页哦注意：2023/4/1038第38页，此课件共40页哦392023/4/10第39页，此课件共40页哦402023/4/10第40页，此课件共40页哦