数值分析--数值积分与数值微分

第5章数值微积分《数值分析》在一元函数的积分学中,我们已经熟知,若函数f(x)在区间［a,b］上连续且其原函数为F(x),则可用牛顿―莱布尼兹公式()()()bafxFbFa第4章数值积分和数值微分§4.1引言第5章数值微积分《数值分析》来求定积分。前面公式虽然在理论上或在解决实际问题中都起了很大的作用,但它并不能完全解决定积分的计算问题。因为定积分的计算常常会碰到以下三种情况:(1)被积函数f(x)的原函数F(x)不易找到。许多很简单的函数,例如等,其原函数都不能用初等函数表示成有限形式。2sin1,,lnxxexx第5章数值微积分《数值分析》(2)被积函数f(x)没有具体的解析表达式。其函数关系由表格或图形表示,无法求出原函数。(3)尽管f(x)的原函数能表示成有限形式但其表达式相当复杂。例如定积分的被积函数的原函数就比较复杂,从数值计算角度来看,计算量太大。41badxx411x第5章数值微积分《数值分析》如图4.1,若用左矩形近似地代替曲边梯形,则得到左矩形公式()()()()()()babafxdxbafafxdxbafb同样可得到右矩形公式第5章数值微积分《数值分析》图4.1第5章数值微积分《数值分析》如图4.2,若用梯形的面积近似地代替曲边梯形的面积,则得到计算定积分的梯形公式()[()()]2babafxdxfafb(1―1)如图4.3,若用抛物线代替曲线f(x),则可得到抛物线公式(或辛普生公式)()[()4()()]2babaabfxdxfaffbb(1―2)第5章数值微积分《数值分析》图4.2图4.3第5章数值微积分《数值分析》.)()3.1(),(d)(,0机械求积公式通常称为的具体形式。被积函数的选取有关，而不依赖仅仅与节点权的权。随节点称为求积系数，亦称伴称为求积节点；式中求积公式一般地xfxAxAxxfAxxfkkkkkknkkba二、代数精度的概念.,1,m次代数精度定义1则称该求积公式具有立成次的多项式等式不准确而对于某一个都准确成立的多项式对于所有次数不超过若一个求积公式mm.,,,)1()0()1(d)(21021011精度尽量高使上述求积公式的代数试确定系数设有求积公式AAAfAfAfAxxf练习.一般方法？第5章数值微积分《数值分析》三、插值型求积公式)()(,,,101010knkknnnfxlxLfffbxxxan项式，就有拉格朗日插值多上已知函数值个互异节点在,d)(d)(d)(0nkkbakbanbafxxlxxLxxf得到(1.5).d)(,d)(0bakknkkkbaxxlwfwxxf其中即得求积公式.称为插值型求积公式第5章数值微积分《数值分析》.d)()!1()(d)()(][0)1(xxxnfxxLxffRbanjjnban它的余项为.d)(0它是插值型求积公式次代数精度至少具有求积公式nfwxxfnkkkba定理1四、求积公式的收敛性和稳定性若中在求积公式,(1.3)定义2baknkkhnxxfxfw,d)()(lim00.(1.3)),(max11是收敛的则称求积公式其中iinixxh第5章数值微积分《数值分析》.]~)([|)~()(|),,,1,0(~)(,)(0knkkknnkkkkkfxfAfIfInkfxfxf则有即有误差设.(1.3),)](~)([|)~()(|),,,0(~)(0,0,0是稳定的则称求积公式就有只要若knkkknnkkxfxfAfIfInkfxf定义3.),,,1,00)3.1(是稳定的则求积公式（中系数若求积公式nAk定理2.)()(~)(||,),,0(~)(,00abAxfxfARnkfxfnkknkkkknkk有时当这是因为第5章数值微积分《数值分析》2牛顿―柯特斯(Newton―Cotes)公式建立数值积分公式最基本的思想是选取一个既简单又有足够精度的函数φ(x),用φ(x)代替被积函数f(x),于是有现用第2章介绍的插值多项式Pn(x)来代替被积函数f(x),即有()()bbaafxdxxdx()()bbnaafxdxPxdx取基点为等距,即a=x0＜x1＜…＜xn=b第5章数值微积分《数值分析》利用拉格朗日插值多项式10,0,1,2,,10,1,2,,kkibahxxknnxxihin()()()nnfxpxRx(2―1)其中00(1)1()()()()()(,)(1)!nnkniikikkinnnxxPxyxxfRxxabn(2―2)第5章数值微积分《数值分析》这里yi=f(xi),对式(2―1)两边积分得00(1)10()()()[]1()()(1)!()bbbnnaaannbkiaikikkibnnaniinifxdxpxdxRxdxxxdxyxxfxdxnayRf第5章数值微积分《数值分析》为牛顿―柯特斯求积公式,Rn(f)为牛顿―柯特斯求积公式的余项。令x=x0+sh,0≤s≤n0(1)101()()()(1)!()nbkiakikkibnnnanbiiaixxadxxxRffxdxnfxdxay(2―3)(2―4)(2―5)我们称第5章数值微积分《数值分析》()00()0000()1(1)()!()!nnniikkinnnikkininnkkibaskadsbacnikskCdsnikskdsinin(2―6)第5章数值微积分《数值分析》称C(n)i为柯特斯求积系数。很显然,当n=1时,可算得1(1)001(1)101(1)212CsdsCsds此时式(2―5)为()[()()]2babafxdxfafb(2―7)第5章数值微积分《数值分析》这是梯形公式。当n=2时,可得2(2)002(2)102(2)2011(1)(2)4614(2)2611(1)46CssdsCssdsCssds于是()[()4()()]62babaabfxdxfaffb(2―8)第5章数值微积分《数值分析》这是抛物线公式。当n=3时,3(3)003(3)103(3)203(3)3011(1)(2)(3)18813(2)(3)6813(1)(3)6811(1)(2)188CsssdsCsssdsCsssdsCsssds第5章数值微积分《数值分析》代入(2―5)式得到求积公式0123()[()3()3()()]8babafxdxfxfxfxfx(2―9)类似地可分别求出n=4,5,…时的柯特斯系数,从而建立相应的求积公式。具体结果见表4―1。从表中可以看出,当n≤7时,柯特斯系数为正;从n≥8开始,柯特斯系数有正有负。因此,当n≥8时,误差有可能传播扩大,牛顿―柯特斯求积公式不宜采用。第5章数值微积分《数值分析》柯特斯系数C(n)i仅与n和i有关,与被积函数f(x)无关,且满足()01nniiC(2―10)事实上,式(2―5)对f(x)=1是准确成立的。第5章数值微积分《数值分析》例1试分别用梯形公式和抛物线公式计算积分解利用梯形公式10.5xdx10.510.5(0.51)20.4267767xdx利用抛物线公式10.510.5(0.540.751)60.43.93403xdx第5章数值微积分《数值分析》原积分的准确值31120.50.520.430964413xdxx第5章数值微积分《数值分析》表4―1第5章数值微积分《数值分析》2.2误差估计现对牛顿―柯特斯求积公式所产生的误差作一个分析。由式(2―4),牛顿―柯特斯求积公式的余项为易知,牛顿―柯特斯求积公式(2―5)对任何不高于n次的多项式是准确成立的。这是因为f(n+1)(ξ)≡0故Rn(f)≡0(1)11()()()(1)!bnnnaRffxdxn第5章数值微积分《数值分析》一般说来,若某个求积公式对于次数不高于m的多项式都准确成立(即Rn(f)≡0),而对于某一次数为m+1的多项式并不准确成立(即Rn(f)0),则称这一求积公式的代数精确度为m。牛顿―柯特斯求积公式的代数精确度至少为n。通常在基点个数相等的情况下,代数精确度愈高,求积公式就愈精确。定理1(梯形公式的误差)设f(x)在区间［a,b］上具有连续的二阶导数,则梯形求积公式的误差为第5章数值微积分《数值分析》由于ω1(x)=(x-a)(x-b)311()()(),(,)121()()()()2babaRffabRffxaxbdx证由式(2―4)知,梯形公式的余项为(2―11)第5章数值微积分《数值分析》在区间(a,b)内不变号,f″(ξ)是x的函数且在［a,b］上连续,故根据积分第二中值定理参见有关《数学分析》教材中“一元函数积分学第二中值定理”。知,存在某一η∈(a,b)使131()()()()21()()12baRffxaxbdxbaf第5章数值微积分《数值分析》定理2(抛物线公式的误差)设f(x)在［a,b］上有连续的四阶导数,则抛物线公式的误差为5(4)22()()(),(,)28801()()()()()3!2babaRffababRffxaxxbdx(2―12)证由式(2―4)知第5章数值微积分《数值分析》§4.3复化求积公式4.3.1复合梯形公式对于定积分,将积分区间［a,b］分成n个相等的子区间［xi,xi+1］,这里步长1,0,1,,1iibahxxinn在每一个子区间［xi,xi+1］上使用梯形公式,则()()()bafxFbFa第5章数值微积分《数值分析》相加后得1311()[()()]()212(,)iixiiixiiihhfxdxfxfxfxx(3―1)110311100()()[()()]()212iinbxaxinniiiiifxdxfxdxhhfxfxf(3―2)若f″(x)在［a,b］上连续,由连续函数的介值定理,存在某一ξ∈(a,b)使得10()()niiffn第5章数值微积分《数值分析》因而33103()()1212()()12niihmhffbaf于是得到复合梯形公式1013()12()[()()2()]2()()()12nbniainhfxdxfxfxfxbaRffn(3―3)其余项为第5章数值微积分《数值分析》例2若用复合梯形公式计算积分问积分区间要等分多少才能保证有五位有效数字解由余项(3―3)式10xedx3()12()12()()()12()()()1()()12nxnbaRffnfxfxfxebafxeeRfn则当0＜x＜1时,有因为又故第5章数值微积分《数值分析》由于原积分的准确值具有一位整数,因此要使近似积分值有五位有效数字,只需取n满足4224110122610lg2lg4lg4lglglg1.82662ennebneebn两边取对数得整理后得到第5章数值微积分《数值分析》4.3.2复化抛物线（辛普森）公式类似复合梯形公式的做法,把区间［a,b］分成n个相等的子区间［x2i,x2i+2］(i=0,1,…,n-1),设每个子区间上的中点为x2i+1(i=0,1,…,n-1),且在每一个子区间［x2i,x2