数值积分公式

数值积分公式第一页，共三十二页，编辑于2023年，星期六内容提纲

数值积分的必要性求积公式及其代数精度插值型求积公式Newton-Cotes公式及数值稳定性复化求积公式及误差估计第二页，共三十二页，编辑于2023年，星期六

数值积分的必要性本章主要讨论如下形式的一元函数积分在微积分里，按Newton-Leibniz公式求定积分要求被积函数f(x)☞有解析表达式；☞

f(x)的原函数F(x)为初等函数．第三页，共三十二页，编辑于2023年，星期六实际问题1.

f(x)的原函数F(x)不能用初等函数表示例如函数:考虑一个实际问题:建筑上用的一种铝制波纹瓦是用一种机器将一块平整的铝板压制而成的.第四页，共三十二页，编辑于2023年，星期六

假若要求波纹瓦长4英尺,每个波纹的高度(从中心线)为1英寸,且每个波纹以近似2π英寸为一个周期.

求制做一块波纹瓦所需铝板的长度L.

这个问题就是要求由函数f(x)=sinx给定的曲线,从x=0到x=48英寸间的弧长L.由微积分学我们知道,所求的弧长可表示为:上述积分称为第二类椭圆积分,它不能用普通方法来计算.第五页，共三十二页，编辑于2023年，星期六2.

有些被积函数其原函数虽然可以用初等函数表示成有限形式,但表达式相当复杂,计算极不方便.例如函数并不复杂,但它的原函数却十分复杂:第六页，共三十二页，编辑于2023年，星期六

3.f(x)没有解析表达式，只有数表形式:

x12345f(x)44.5688.5

这些都说明,通过原函数来计算积分有它的局限性,因而,研究关于积分的数值方法具有很重要的实际意义.第七页，共三十二页，编辑于2023年，星期六求积公式及其代数精度

求积公式的概念积分值

在几何上可解释为由x=a,x=b,y=0和y=f(x)

所围成的曲边梯形的面积.积分计算之所以有困难，就是因为这个曲边梯形有一条边y=f(x)是曲的.

第八页，共三十二页，编辑于2023年，星期六

依据积分中值定理,对于连续函数f(x)

,在[a,b]内存在一点ξ,使得

称f(ξ)为区间[a,b]的平均高度.问题在于点ξ的具体位置一般是不知道的.这样,只要对平均高度f(ξ)提供一种算法,相应地便获得一种数值求积方法.第九页，共三十二页，编辑于2023年，星期六如果简单地选取区间[a,b]的一个端点或区间中点的高度作为平均高度,这样建立的求积公式分别是:左矩形公式:

I(f)≈(b-a)f(a)右矩形公式:

I(f)≈(b-a)f(b)中矩形公式:

I(f)≈(b-a)f[(a+b)/2]第十页，共三十二页，编辑于2023年，星期六此外,众所周知的梯形公式:

I(f)≈(b-a)[f(a)+f(b)]/2和Simpson公式:

I(f)≈(b-a)[f(a)+4f((a+b)/2)+f(b)]/6则分别可以看作用a,b,c=(a+b)/2,

三点高度的加权平均值

[f(a)+f(b)]/2和[f(a)+4f(c)+f(b)]/6作为平均高度f(ξ)的近似值.第十一页，共三十二页，编辑于2023年，星期六更一般地,取区间[a,b]内n+1个点

{xi},(i=0,1,2,…n)处的高度{f(xi)}(i=0,1,…,n)通过加权平均的方法近似地得出平均高度f(ξ),这类求积方法称为机械求积:第十二页，共三十二页，编辑于2023年，星期六或写成:数值积分公式求积系数求积节点(1)第十三页，共三十二页，编辑于2023年，星期六记称(2)为数值求积公式,(3)为求积公式余项(误差).构造或确定一个求积公式，要讨论解决的问题有(i)

确定求积系数Ak和求积节点xk

；(ii)

求积公式的误差估计和收敛性为了构造形如式(2)的求积公式,需要提供一种判定求积方法精度高低准则第十四页，共三十二页，编辑于2023年，星期六求积公式的代数精度定义1

称求积公式(2)具有m次代数精度,如果它满足如下两个条件:

(i)对所有次数≤m次的多项式,有

(ii)存在m+1次多项式,使得定义1中的条件(i),(ii)等价于:第十五页，共三十二页，编辑于2023年，星期六插值型求积公式

在积分区间[a,b]上取n+1个节点xi，i=0,1,2,…,n,作f(x)的n次代数插值多项式（拉格朗日插值公式）:则有

为插值余项于是有第十六页，共三十二页，编辑于2023年，星期六取称(4)式为插值型求积公式,求积系数Ak由(5)式确定.(4)(5)Ak由节点决定，与f(x)

无关。第十七页，共三十二页，编辑于2023年，星期六第十八页，共三十二页，编辑于2023年，星期六推论1

求积系数满足:误差定理1

形如的求积公式至少有n

次代数精度

该公式为插值型（即：）第十九页，共三十二页，编辑于2023年，星期六取节点为等距分布：令Cotes系数注：Cotes系数仅取决于n

和k，可查表得到。与f(x)及区间[a,b]均无关。Newton-Cotes公式由此构造的插值型求积公式称为Newton-Cotes公式,此时求积系数第二十页，共三十二页，编辑于2023年，星期六记则求积公式(4)变为(6)(7)(8)称(8)式为n阶闭型Newton-Cotes求积公式.第二十一页，共三十二页，编辑于2023年，星期六Newton-Cotes公式的误差为:与x有关注意:由(6)式确定的Cotes系数只与j和n有关,与f(x)和积分区间[a,b]无关,且满足:(9)第二十二页，共三十二页，编辑于2023年，星期六定理2当阶数n为偶数时,Newton-Cotes公式(8)至少具有n+1次代数精度.证明

只需验证当n为偶数时,Newton-Cotes公式对f(x)=xn+1的余项为零.由于f(x)=xn+1,所以f(n+1)(x)=(n+1)!.由式(9)得引进变换t=u+n/2,因为n为偶数,故n/2为整数,于是有据此可断定R(f)=0,因为上述被积函数是个奇函数.第二十三页，共三十二页，编辑于2023年，星期六Newton-Cotes公式的数值稳定性

现在讨论舍入误差对计算结果产生的影响.设用公式

近似计算积分时,其中计算函数值f(xj)有误差εj

(j=0,1,2,…,n).设计算Cj(n)没有误差,中间计算过程中的舍入误差也不考虑,则在式(10)的计算中,由εj引起的误差为(10)第二十四页，共三十二页，编辑于2023年，星期六如果Cj(n)都是正数,并设故en是有界的,即由εj引起的误差受到控制,不超过ε的(b-a)倍,保证了数值计算的稳定性.而当n>7时,Cj(n)将出现负数,保证数值稳定性.因此高阶公式不宜采用,有实用价值的仅仅是几种低阶的求积公式.将随n增大,因而不能则有第二十五页，共三十二页，编辑于2023年，星期六n=1:TrapezoidalRule/*令x=a+th,h=ba,用中值定理*/代数精度=1n=2:Simpson’sRule代数精度=3n=3:Simpson’s3/8-Rule,代数精度=3,n=4:CotesRule,代数精度=5,n

为偶数阶的Newton-Cotes

公式至少有n+1

次代数精度。第二十六页，共三十二页，编辑于2023年，星期六复化型求积公式高次插值有Runge

现象,高阶Newton-Cotes公式会出现数值不稳定,低阶Newton-Cotes公式有时又不能满足精度要求.解决这个矛盾的办法是将积分区间[a,b]分成若干小区间,在每个小区间上用低阶求积公式计算,然后将它们加起来,这就是复化求积方法.第二十七页，共三十二页，编辑于2023年，星期六

复化梯形公式：在每个上用梯形公式：=

Tn/*中值定理*/第二十八页，共三十二页，编辑于2023年，星期六

复化梯形公式积分法第二十九页，共三十二页，编辑于2023年，星期六收敛性由上述的误差估计式可知，当f(x)C2[a,b]

时,只要h0时数列Tn(f)