计算方法_数值积分

《计算方法_数值积分》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第5章数值积分若函数f(x)在区间［a,b］上连续且其原函数为F(x),则可用牛顿―莱布尼兹公式，来求定积分。(5―1)求定积分复习函数关系由表格或图形表示,无法求出原函数。定积分计算可能遭遇的三种情况被积函数的原函数不是初等函数被积函数f(x)没有具体的解析表达式被积函数f(x)的原函数F(x)不易找到第5章数值积分从几何上看定积分定积分是曲边梯形的面积图5.1左矩形右矩形(5―2)(5―3)图5.2梯形面积图5.3抛物求积(5―4)(5―5)第5章数值积分近似值§5.1§5.2§5.4牛顿―柯特斯(Ne

2、wton―Cotes)公式复合求积公式龙贝格(Romberg)积分方法5.1牛顿―柯特斯(Newton―Cotes)公式建立数值积分公式最基本的思想是选取一个既简单又有足够精度的函数φ(x),用φ(x)代替被积函数f(x),于是有现用第四章介绍的插值多项式Pn(x)来代替被积函数f(x),即有将积分区间[a,b]n等分，则节点是等距分布的，节点x0,x1,x2,…,xn可表示成xk=x0+kh(k=0,1,…,n)，其中x0=a,xn=b,称为步长。Newton-Cotes公式若Ln(x)为Lagrange

3、插值多项式，则由公式于是令(5.5)公式(5.6)称为等距节点内插求积公式。则有(5.6)求Ak在等距节点前提下，做变换，由，可得而x-xj=(t-j)h(j=0,1,2,…,n)，xk-xj=(k-j)h(j,k=0,1,2,…,n且j≠k)。于是(5.5)式即为记则(5.9)称为牛顿-柯特斯公式。其中Ck(n)叫Cotes系数，Cotes系数与被积函数及积分区间无关。计算柯特斯系数n=1时，有两个Cotes系数n=2时，有三个Cotes系数类似可得，n=3时有四个Cotes系数n=4时，有五个Cotes

4、系数几个常用的牛顿-柯特斯公式n=1时，，此即(5.3)式，为梯形公式。，其中，称为Simpson公式。其中c,d,e为[a,b]的四等分点，称为Cotes公式。n=2时，n=4时，表5―1柯特斯系数柯特斯系数C(n)i仅与n和i有关,与被积函数f(x)无关,且满足(5―15)柯特斯公式对f(x)=1是准确成立的。柯特斯系数的特点例1试分别用梯形公式和辛普森公式计算积分解:利用梯形公式利用抛物线公式原积分的准确值5.1.2误差估计现对牛顿―柯特斯求积公式所产生的误差作一个分析。牛顿―柯特斯求积公式的余项为易

5、知,牛顿―柯特斯求积公式对任何不高于n次的多项式是准确成立的。这是因为f(n+1)(ξ)≡0故Rn(f)≡0(5―10)代数精度一般说来,若某个求积公式对于次数不高于m的多项式都准确成立(即Rn(f)≡0),而对于某一次数为m+1的多项式并不准确成立(即Rn(f)≠0),则称这一求积公式的代数精度为m。牛顿―柯特斯求积公式的代数精度至少为n，若n为偶数，则至少具有n+1次代数精度。通常在基点个数相等的情况下,代数精度愈高,求积公式愈精确。梯形公式、辛普森公式、柯特斯公式分别具有1、3、5次代数精度。例5.1

6、分别利用梯形公式、Simpson公式、Cotes公式计算,n=1,2,3,4,5，并与用牛顿-莱布尼兹公式计算的结果进行比较。解计算结果列于表5-2中。表5-2函数f(x)xx2x3x4x5梯形值0.50.50.50.50.5Simpson值0.50.3333330.250.2083330.1875Cotes值0.50.3333330.250.200.166667准确值0.50.3333330.250.200.166667证由式知,梯形公式的余项为(x-a)(x-b)在区间(a,b)内不变号,f″(ξ)是x

7、的函数且在［a,b］上连续,故根据积分第二中值定理参见有关《数学分析》教材中“一元函数积分学第二中值定理”。知,存在某一η∈(a,b)使定理2(抛物线公式的误差)设f(x)在［a,b］上有连续的四阶导数,则抛物线公式的误差为定理1(梯形公式的误差)设f(x)在区间［a,b］上具有连续的二阶导数,则梯形求积公式的误差为如果在每个子区间上使用梯形公式，就得到复合梯形公式。将积分区间[a,b]N等分后的节点记为xk，xk=a+kh(k=0,1,2,…,N)，在每个子区间[xk,xk+1](k=0,1,2,…，N-

8、1)上应用梯形公式，1.复合梯形公式5.2复合求积公式再求和得：1.复合梯形公式其中xk=a+kh(k=0,1,2,…,N)，1.复合梯形公式复合梯形公式：在每个上用梯形公式：=Tn2.复合Simpson公式如果在每个子区间上使用Simpson公式，就得到复合Simpson公式。将N等分后的每个子区间再对分一次，于是共有2N+1个节点，(k=0,1,2,…,2N)，在每个N等分的子区间[x2k,x2k+2](k