数值积分与数值微分.pptx
《数值积分与数值微分.pptx》由会员分享,可在线阅读,更多相关《数值积分与数值微分.pptx(152页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、进行计算,但在工程计算和科学研究中,经常会遇进行计算,但在工程计算和科学研究中,经常会遇到被积函数到被积函数f(x)的下列一些情况:的下列一些情况:的原函数的原函数对定积分对定积分的被积函数的被积函数已知,在高等数学中可用牛顿已知,在高等数学中可用牛顿莱布尼兹公式莱布尼兹公式4.1 数值积分概论数值积分概论 实际问题当中常常要计算积分,有些数值方法,实际问题当中常常要计算积分,有些数值方法,如微分方程和积分方程的求解,也都和积分计算相如微分方程和积分方程的求解,也都和积分计算相联系联系.4.1.1 数值求积的基本思想数值求积的基本思想第1页/共152页(4)f(x)本身没有解析表达式,其函数关
2、系由表格本身没有解析表达式,其函数关系由表格或图形给出,列如为实验或测量数据或图形给出,列如为实验或测量数据.(2)f(x)的原函数不能用初等函数形式表示,例如的原函数不能用初等函数形式表示,例如(3)f(x)的原函数虽然可用初等函数形式表示,但的原函数虽然可用初等函数形式表示,但其原函数表示形式相当复杂,例如其原函数表示形式相当复杂,例如(1)f(x)复杂,求原函数困难,列如复杂,求原函数困难,列如第2页/共152页 以上的以上的 4种情况都不能用牛顿种情况都不能用牛顿莱布尼兹公莱布尼兹公式方便地计算该函数的定积分,满足不了实际需式方便地计算该函数的定积分,满足不了实际需要,因此,有必要研究
3、定积分的数值计算问题;要,因此,有必要研究定积分的数值计算问题;另外,对一些函数的求导问题,其求导、微分也另外,对一些函数的求导问题,其求导、微分也相当复杂,也有必要研究求导、微分的数值计算相当复杂,也有必要研究求导、微分的数值计算问题。本章主要介绍数值求积分和数值求微分的问题。本章主要介绍数值求积分和数值求微分的方法。方法。第3页/共152页 由积分中值定理由积分中值定理,对连续函数对连续函数f(x),在区间在区间a,b内至少存在一点内至少存在一点,使,使只要对平均高只要对平均高度度 f()提供一提供一种种近似算法近似算法,便便可相应地获得可相应地获得一种一种数值求积数值求积方法方法.即所谓
4、即所谓矩矩形公式形公式.第4页/共152页 例如例如,用区间用区间a,b两端点的函数值两端点的函数值 f(a)与与f(b)的的算术平均值作为算术平均值作为f()的近似值的近似值,可导出可导出求积公式求积公式这便是人们所熟知的这便是人们所熟知的梯形公式梯形公式.第5页/共152页 如果改用区间如果改用区间a,b的中点的中点 c=(a+b)/2 处的函数值处的函数值f(c)近似代替近似代替f(),则又可导出所谓则又可导出所谓(中中)矩形公式矩形公式第6页/共152页 一般地一般地,在区间在区间a,b上适当选取点上适当选取点xk(k=0,1,n),然后用然后用 f(xk)的的加权平均值加权平均值作为
5、作为f()的近似值的近似值,可得到可得到更为更为一般的求积公式一般的求积公式 其中:点其中:点xk叫叫求积节点求积节点,系数系数Ak叫叫求积系数求积系数.Ak仅与节仅与节点点xk的选取有关的选取有关,而与被积函数而与被积函数 f(x)无关无关.求积公式的求积公式的截断误差截断误差为为 R(f)又称为又称为求积余项求积余项.这类数值积分方法通常称为这类数值积分方法通常称为机械求积机械求积,其特点是,其特点是将积分求值问题归结为函数值的计算,这就避开了牛将积分求值问题归结为函数值的计算,这就避开了牛-莱公式寻求原函数的困难莱公式寻求原函数的困难.第7页/共152页4.1.2 代数精度的概念代数精度
6、的概念 定义定义1 如果求积公式如果求积公式(1)对所有次数不超过对所有次数不超过m的多项式都精确成立;的多项式都精确成立;(2)至少对一个至少对一个m+1次多项式不精确成立,次多项式不精确成立,则称则称该公式具有该公式具有m次代数精度次代数精度(或或代数精确度代数精确度).数值求积方法的近似方法,为要保证精度,我数值求积方法的近似方法,为要保证精度,我们自然希望求积公式能对们自然希望求积公式能对“尽可能多尽可能多”的函数准确的函数准确地成立,这就提出了所谓代数精度的概念地成立,这就提出了所谓代数精度的概念.第8页/共152页 一般来说,代数精度越高,求积公式越好。一般来说,代数精度越高,求积
7、公式越好。结论结论 一个求积公式具有一个求积公式具有m次代数精度的次代数精度的充要条充要条件件是该求积公式是该求积公式:(1)对对xk(k=0,1,m)精确成立;精确成立;(2)对对xm+1不精确成立不精确成立.故一般地,要验证一个求积公式具有故一般地,要验证一个求积公式具有m次代数次代数精度,只要令对于精度,只要令对于 f(x)=1,x,xm求积公式精确成立求积公式精确成立等式就行等式就行.第9页/共152页 即对于求积公式即对于求积公式 给定给定n+1个互异的求积节点个互异的求积节点 x0,x1,xn-1,xn,令求积公式对令求积公式对 f(x)=1,x,xn 精确成立精确成立,即得即得求
8、解该方程组即可确定求积系数求解该方程组即可确定求积系数Ak,所得到的求积公所得到的求积公式式至少具有至少具有n 次代数精度次代数精度.第10页/共152页 解解 当当 f(x)=1时时,此时公式精确成立。此时公式精确成立。例例1 验证梯形公式验证梯形公式具有一次代数精度。具有一次代数精度。当当 f(x)=x时,时,公式也精确成立公式也精确成立.当当 f(x)=x2 时,时,公式对公式对x2不精确成立不精确成立.故由定理故由定理1知知,梯形公式的代数精度为梯形公式的代数精度为1次次.第11页/共152页 例例2 确定求积公式中的待定系数,使其代数精确定求积公式中的待定系数,使其代数精度尽量高,并
9、指明求积公式所具有的代数精度度尽量高,并指明求积公式所具有的代数精度.解解 令令 f(x)=1,x,x2 代入公式两端并令其相等,得代入公式两端并令其相等,得 解得解得第12页/共152页得得求积公式求积公式为为令令 f(x)=x3,得,得令令 f(x)=x4,得,得故故求积公式求积公式具有具有3次代数精度次代数精度.第13页/共152页 如果我们事先选定求积节点如果我们事先选定求积节点xk,譬如,以区间,譬如,以区间a,b的等距分点作为节点,这时取的等距分点作为节点,这时取m=n求解方程组求解方程组即可确定求积系数即可确定求积系数Ak,而使求积公式至少具有,而使求积公式至少具有 n次次代数精
10、度代数精度.本章第本章第2节介绍这样一类求积公式,梯形节介绍这样一类求积公式,梯形公式是其中的一个特例公式是其中的一个特例.如为了构造出上面的求积公式,原则上是一个如为了构造出上面的求积公式,原则上是一个确定参数确定参数xk和和Ak的代数问题的代数问题.方程组方程组(1.4)实际上是一实际上是一2n+2个参数的非线性方程组,此方程组当个参数的非线性方程组,此方程组当n1时求时求解非常困难,但当解非常困难,但当n=0及及n=1的情形还是可以通过求的情形还是可以通过求解方程组得到相应的求积公式解方程组得到相应的求积公式.下面对下面对n=0讨论求积讨论求积公式的建立及代数精确度公式的建立及代数精确度
11、.第14页/共152页 此时求积公式为此时求积公式为其中,其中,x0及及A0为待定参数为待定参数.根据代数精确度定义可令根据代数精确度定义可令f(x)=1,x,由方程组知,由方程组知.得得得到的求积公式就是得到的求积公式就是(1.2)式的中矩形公式式的中矩形公式.再令再令f(x)=x2,代入,代入(1.4)式的第三式式的第三式第15页/共152页说明说明(1.2)式对式对f(x)=x2不精确成立,故它的代数精确不精确成立,故它的代数精确度为度为1.方程组方程组(1.4)是根据是根据(1.3)式的求积公式得到的,式的求积公式得到的,按照代数精确度的定义,如果求积公式中除了按照代数精确度的定义,如
12、果求积公式中除了f(xi)还有还有(x)在某些节点上的值,也同样可得到相应的在某些节点上的值,也同样可得到相应的求积公式求积公式.第16页/共152页 例例1 给定形如下面的求积公式,试确定系数给定形如下面的求积公式,试确定系数A0,A1,B0,,使公式具有尽可能高的代数精确度,使公式具有尽可能高的代数精确度.解解 根据题意可令根据题意可令f(x)=1,x,x2分别代入求积公式分别代入求积公式使它精确成立:使它精确成立:第17页/共152页解得解得当当f(x)=x3时,上式右端为时,上式右端为1/3,而左端是,而左端是于是有求积公式于是有求积公式故积分公式对故积分公式对f(x)=x3不精确成立
13、,其代数精确度为不精确成立,其代数精确度为2.第18页/共152页4.1.3 插值型的求积公式插值型的求积公式设给定一组节点设给定一组节点且已知且已知f(x)在这些节点上的函数值在这些节点上的函数值 f(xk),则可求得则可求得f(x)的拉格朗日插值多项式的拉格朗日插值多项式(因为因为Ln(x)的原函数易求的原函数易求)其中其中lk(x)为插值基函数为插值基函数,取取由上式确定系数的公式称为由上式确定系数的公式称为插值型求积公式插值型求积公式.即即则则 f(x)Ln(x)第19页/共152页插值型求积公式积分法几何表示第20页/共152页由插值余项定理由插值余项定理,其求积余项为其求积余项为其
14、中其中=(x)如果求积公式如果求积公式(1.5)是插值型的,按照插值余是插值型的,按照插值余项式子,对于次数不超过项式子,对于次数不超过n的多项式的多项式f(x),其余项,其余项 R(f)等于零,因而等于零,因而这时求积公式至少具有这时求积公式至少具有n次代数次代数精度精度.第21页/共152页 反之,如果求积公式至少具有反之,如果求积公式至少具有n次代数精度,次代数精度,则它必定是插值型的则它必定是插值型的.事实上,这时求积公式对于事实上,这时求积公式对于插值基函数插值基函数 lk(x)应准确成立,即有应准确成立,即有注意到注意到lk(xj)=kj,上式右端实际上即等于,上式右端实际上即等于
15、Ak,因,因而下面式子成立而下面式子成立.第22页/共152页 定理定理1 具有具有n+1个节点的数值求积公式个节点的数值求积公式(1.5)是是插值型求积公式插值型求积公式的的充要条件充要条件为为:该公式该公式至少具有至少具有n次代数精度次代数精度.综上所述,我们有结论为综上所述,我们有结论为 这时令这时令f(x)=1代入又有结论为代入又有结论为 结论结论 对插值型求积公式的系数必有对插值型求积公式的系数必有第23页/共152页4.1.4 求积公式的余项求积公式的余项若求积公式若求积公式(1.3)的代数精确度为的代数精确度为m,则由求积公,则由求积公式余项的表达式式余项的表达式(1.7)可以证
16、明余项形如可以证明余项形如其中其中K为不依赖于为不依赖于f(x)的待定参数的待定参数.这个结果表明当这个结果表明当f(x)是次数小于等于是次数小于等于m的多项式时,由于的多项式时,由于f(m+1)(x)=0,故此时故此时R f=0,即求积公式,即求积公式(1.3)精确成立精确成立.而当而当f(x)=xm+1时,时,f(m+1)(x)=(m+1)!,(1.8)式左端式左端 Rf 0,故可求得,故可求得第24页/共152页代入余项公式代入余项公式(1.8)式可以得到更细致的余项表达式式可以得到更细致的余项表达式.例如梯形公式例如梯形公式(1.1)的代数精确度为的代数精确度为1,可以证明,可以证明它
17、的余项表达式为它的余项表达式为其中其中于是得到梯形公式于是得到梯形公式(1.1)的余项为的余项为第25页/共152页对中矩形公式对中矩形公式(1.2),其代数精确度为,其代数精确度为1,可以证,可以证明它的余项表达式为明它的余项表达式为其中其中于是得到中矩形公式于是得到中矩形公式(1.2)的余项为的余项为第26页/共152页 例例2 求例求例1中求积公式中求积公式的余项的余项.解解 由于此求积公式的代数精确度为由于此求积公式的代数精确度为2,故余项,故余项表达式为表达式为R=K(),令,令f(x)=x3,得,得()=3!,于是有,于是有故得故得第27页/共152页其中其中h=max(xi-xi
18、-1),则称求积公式,则称求积公式(1.3)是是收敛的收敛的.4.1.5 求积公式的收敛性与稳定性求积公式的收敛性与稳定性 定义定义2 在求积公式在求积公式(1.3)中,若中,若 在求积公式在求积公式(1.3)中,由于计算中,由于计算f(xk)可能产生误可能产生误差差k,实际得到,实际得到 ,即,即 .记记如果对任给小正数如果对任给小正数 0,只要误差,只要误差|k|充分小就有充分小就有它表明求积公式它表明求积公式(1.3)计算是计算是稳定的稳定的,由此给出,由此给出第28页/共152页 定义定义3 对任给小正数对任给小正数 0,若存在,若存在 0,只要,只要 就有就有(1.12)式成立,则式
19、成立,则称求积公式称求积公式(1.3)是是稳定的稳定的.证明证明 对任给对任给 0,若取若取=/(b-a),对所有对所有k都有都有故求积公式故求积公式(1.3)是稳定的是稳定的.证毕证毕.定理定理2 若求积公式若求积公式(1.3)中所有系数中所有系数Ak0,则此,则此求积公式是稳定的求积公式是稳定的.则有则有定理定理2表明只要求积系数表明只要求积系数Ak0,就能保证计算的,就能保证计算的稳定性稳定性.第29页/共152页 4.2 牛顿牛顿柯特斯公式柯特斯公式 为便于上机计算,通常在内插求积公式中我们为便于上机计算,通常在内插求积公式中我们通常取等距节点,即将积分区间通常取等距节点,即将积分区间
20、a,b划分划分n等分,等分,即令步长即令步长h=(b-a)/n,且记,且记x0=a,xn=b,则节点记为,则节点记为xk=x0+kh(k=0,1,n),然后作变换,然后作变换:t=(x-x0)/h,代入代入求积系数公式,将会简化计算求积系数公式,将会简化计算.第30页/共152页4.2.1 柯特斯系数与辛普森公式柯特斯系数与辛普森公式设将积分区间设将积分区间a,b划分成划分成 n等分等分,步长步长h=求积节点取为求积节点取为xk=a+kh(k=0,1,n),由此构造插值型由此构造插值型求积公式求积公式,则其求积系数为则其求积系数为引入变换引入变换 x=a+th,则有则有(k=0,1,n)(k=
21、0,1,n)第31页/共152页其中其中记记于是得求积公式于是得求积公式称为称为n 阶牛顿阶牛顿-柯特斯柯特斯(Newton-Cotes)公式公式.显然显然,柯特斯系数与被积函数柯特斯系数与被积函数 f(x)和积分区间和积分区间a,b无关无关,且为容易计算的多项式积分且为容易计算的多项式积分.称为称为柯特斯系数柯特斯系数.第32页/共152页n11/21/221/64/61/631/83/83/81/847/9032/90 12/90 32/907/90519/28875/28850/28850/28875/28819/288641/840216/84027/840272/84027/8402
22、16/84041/840常用的柯特斯系数表常用的柯特斯系数表第33页/共152页 当当n=1时,时,柯特斯系数柯特斯系数为为这时的这时的牛顿牛顿-柯特斯公式柯特斯公式为一阶求积公式,就是我们为一阶求积公式,就是我们所熟悉的所熟悉的梯形公式梯形公式,即,即第34页/共152页 当当n=2时,时,柯特斯系数柯特斯系数为为相应的相应的牛顿牛顿-柯特斯公式柯特斯公式为二阶求积公式,就是为二阶求积公式,就是辛普辛普森森(simpson)公式公式(又称为又称为抛物形求积公式抛物形求积公式),即,即第35页/共152页式中式中(k=0,1,2,3,4),h=(b-a)/4.n=4 时的时的牛顿牛顿-柯特斯公
23、式柯特斯公式就特别称为就特别称为柯特斯公柯特斯公式式.其形式是其形式是 在在柯特斯系数表柯特斯系数表中中(见书见书p104)看到看到n 7时,时,柯特柯特斯系数斯系数出现负值,于是有出现负值,于是有第36页/共152页特别地,假定特别地,假定则有则有这表明在这表明在b-a1时,初始误差将会引起计算结果误差时,初始误差将会引起计算结果误差增大,即计算不稳定,故增大,即计算不稳定,故n 7的的牛顿牛顿-柯特斯公式柯特斯公式是是不用的不用的.第37页/共152页4.2.2 偶阶求积公式的代数精度偶阶求积公式的代数精度 作为插值型求积公式,作为插值型求积公式,n阶阶牛顿牛顿-柯特斯公式柯特斯公式至至少
24、具有少具有n次代数精度次代数精度(推论推论1).实际的代数精度能否进实际的代数精度能否进一步提高呢?一步提高呢?先看先看辛普森公式辛普森公式,它是二阶,它是二阶牛顿牛顿-柯特斯公式柯特斯公式,因此至少具有二次代数精度因此至少具有二次代数精度.进一步用进一步用f(x)=x3进行检进行检验,按验,按辛普森公式辛普森公式计算得计算得第38页/共152页另一方面,直接求积得另一方面,直接求积得这时有这时有S=I,即,即辛普森公式辛普森公式对不超过三次的多项式均对不超过三次的多项式均能精确成立,又容易验证它对能精确成立,又容易验证它对f(x)=x4通常是不精确通常是不精确的的(如取如取a=0,b=1进行
25、验证有,进行验证有,S=5/24I=1/5),因此,因此,辛普森公式辛普森公式实际上实际上具有三次代数精度具有三次代数精度.一般地,我们可以证明下述论断:一般地,我们可以证明下述论断:第39页/共152页 定理定理3 n 阶牛顿阶牛顿-柯特斯公式的代数精度至少为柯特斯公式的代数精度至少为 证明证明 由定理由定理1已知,无论已知,无论n为奇数或偶数,插值为奇数或偶数,插值型求积公式都至少具有型求积公式都至少具有n次代数精度次代数精度.因此我们证明因此我们证明n为偶数的情形,即对为偶数的情形,即对n+1次多项式余项为零次多项式余项为零.令令n=2k,设设为任一为任一n+1次多项式,其最高次系数为次
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数值 积分 微分
限制150内