基本信息
文件名称:计算方法龙贝格高斯求积公式.ppt
文件大小:2.3 MB
总页数:40 页
更新时间:2025-07-09
总字数:约2.12千字
文档摘要

*①②*第28页,共40页,星期日,2025年,2月5日计算方法龙贝格高斯求积公式第1页,共40页,星期日,2025年,2月5日综合前几节的内容,我们知道梯形公式,Simpson公式,Cotes公式的代数精度分别为1次,3次和5次复合梯形、复合Simpson、复合Cotes公式的收敛阶分别为2阶、4阶和6阶无论从代数精度还是收敛速度,复合梯形公式都是较差的有没有办法改善梯形公式呢?**第2页,共40页,星期日,2025年,2月5日一、复合梯形公式的递推化各节点为复合梯形(Trapz)公式为--------(1)--------(2)**第3页,共40页,星期日,2025年,2月5日--------(3)由(1)(2)两式可**(3)式称为递推的梯形公式递推梯形公式加上一个控制精度,即可成为自动选取步长的复化梯形公式优点:梯形法计算简单缺点:收敛慢,为了达到要求的精度,需要二分区间多次,分点大量增加,计算量很大第4页,共40页,星期日,2025年,2月5日--------(3)**第5页,共40页,星期日,2025年,2月5日则由(1)(2)(3)式,有**第6页,共40页,星期日,2025年,2月5日因此(1)(2)(3)式可化为如下递推公式-------(4)上式称为复合变步长梯形求积公式**第7页,共40页,星期日,2025年,2月5日二、外推加速公式由复合梯形公式的余项公式可得由(3)式移项合并**第8页,共40页,星期日,2025年,2月5日复合Simpson公式表明:由梯形公式前后两次结果的线性组合可构造出精确度较高的辛普森公式**第9页,共40页,星期日,2025年,2月5日--------(5)即--------(6)当然**第10页,共40页,星期日,2025年,2月5日*因此由复合Simpson公式的余项可得令--------(7)--------(8)即当然*第11页,共40页,星期日,2025年,2月5日由复合Cotes公式的余项得令--------(9)即当然公式(9)称为龙贝格积分公式**第12页,共40页,星期日,2025年,2月5日外推加速公式以上整个过程称为Romberg算法将上述结果综合后**第13页,共40页,星期日,2025年,2月5日其中外推加速公式可简化为--------(9)Romberg算法的收敛阶高达m+1的两倍Romberg算法求解步骤Romberg算法的代数精度为m的两倍**第14页,共40页,星期日,2025年,2月5日上述公式推导说明,T1公式是梯形公式,对于次数不高于1的多项式准确;S1是辛普森公式,对于次数不高于3的多项式准确;C1是柯特斯公式,它对于次数不高于5的多项式准确,每一个公式均由前一公式的适当线性组合得到,精确度都提高2次。因此可以验证,由柯特斯公式C1构造得到的龙贝格公式R1,对次数不高于7次的多项式准确。龙贝格公式计算积分占用内存少,精度高**第15页,共40页,星期日,2025年,2月5日龙贝格公式计算步骤:其计算过程是将区间逐次分半,加速得到积分近似值,因此称为逐次分半加速法**第16页,共40页,星期日,2025年,2月5日00.920735510.93979330.946145920.94457350.94608690.946083030.94569090.94608330.94608310.9460831例1将三个加速公式用于求从表中可以看出三次加速求得R1=0.9460831每位数字都是有效数字**第17页,共40页,星期日,2025年,2月5日例2.**解:(1)在区间[1,2]上用梯形公式得第18页,共40页,星期日,2025年,2月5日**第19页,共40页,星期日,2025年,2月5日**第20页,共40页,星期日,2025年,2月5日本章作业P9710**第21页,共40页,星期日,2025年,2月5日牛顿—柯特斯型求积公式是封闭型的(区间[a,b]的两端点a,b均是求积节点)而且要求求积节点是等距的,受此限制,牛顿—柯特斯型求积公式的代数精度只能是n(n为奇数)或n+1(n为偶数)。4.7.1一般理论4.7高斯求积公式**问题:对求积节点也适当的选取,即在求积公式中不仅Ak而且xk也加以选取,是否可以增加自由度,从而可提高求积公