数值积分和数值微分yjs0000市公开课一等奖百校联赛特等奖课件.pptx

科大研究生学位课程数值分析数值分析 计算定积分有微积分基本公式但很多函数找不到原函数，如等。而实际上，有很多函数只知一些离散点函数值，并无表示式，这就需要利用已知条件求出近似值。第第5章章 数值积分与数值微分数值积分与数值微分 x12345f(x)44.5688.51第1页科大研究生学位课程数值分析数值分析数值积分数值积分/Numerical Integration/定义数值积分以下：是离散点上函数值线性组合是离散点上函数值线性组合称为求积系数,与f(x)无关,与积分区间和求积节点求积节点相关称为数值积分公式称为数值积分公式数值积分问题可分解为下述三个问题：1、求积公式详细结构问题；、求积公式详细结构问题；(包含包含xi选取和选取和Ai结构结构)3、准确性程度衡量标准问题。、准确性程度衡量标准问题。2、余项预计问题、余项预计问题(亦即误差预计问题亦即误差预计问题)；求积公式误差求积公式误差 RfI*fIf2第2页科大研究生学位课程数值分析数值分析1、处理第一个问题；节点、处理第一个问题；节点xi 和和系数系数Ai怎样选取，即选取标准怎样选取，即选取标准两个目标：两个目标：1、余项预计问题；求积公式误差、余项预计问题；求积公式误差 RfI*fIf尽可能小。尽可能小。2、求积公式代数精度尽可能高。、求积公式代数精度尽可能高。2、处理第二个问题；依赖插值多项式余项预计公式。、处理第二个问题；依赖插值多项式余项预计公式。3、对于第三个问题；引进代数精度概念、对于第三个问题；引进代数精度概念3第3页科大研究生学位课程数值分析数值分析 定义定义5.1 若求积公式对(x)=xj(j=0,1,2,m)都准确成立,但对(x)=xm+1不准确成立,即则称此公式含有含有m次代数精度次代数精度.可见,若公式含有m次代数精度,则公式对全部次数不超出m多项式都准确成立.注意：注意：、求积公式误差是计算精度度量标志，而代数精度、求积公式误差是计算精度度量标志，而代数精度是求积公式优良性能标志。是求积公式优良性能标志。2、求积公式误差小，不代表代数精度高。代数精度高，、求积公式误差小，不代表代数精度高。代数精度高，也不代表求积公式误差小。它们没有必定联络。也不代表求积公式误差小。它们没有必定联络。4第4页科大研究生学位课程数值分析数值分析例例 1 确定形如求积公式，使其代数精度尽可能高。数值求积公式为 解解 令公式对(x)=1,x,x2都准确成立,则 A0+A1+A2=3 A1+3A2=4.5 A1+9A2=9解之得:A0=0,A1=9/4,A2=3/4.例例2 试确定参数A0,A1,A2,使求积公式含有尽可能高代数精度,并问代数精度是多少?5第5页科大研究生学位课程数值分析数值分析 解得:A0=A2=1/3,A1=4/3.求积公式为 当(x)=x3时,左=0,右=0,公式也准确成立.解解 令公式对(x)=1,x,x2 都准确成立,则 A0+A1+A2=2 -A0+A2=0 A0+A2=2/3当(x)=x4时,左=2/5,右=2/3,公式不准确成立.所以,此公式代数精度为3.6第6页科大研究生学位课程数值分析数值分析例例3 试确定参数A0,A1和x0,x1,使求积公式含有尽可能高代数精度,并问代数精度是多少?解解 令公式对(x)=1,x,x2,x3都准确成立,则 A0+A1=2 A0 x0+A1x1=0 A0 x02+A1x12=2/3 A0 x03+A1x13=0解得:求积公式为求积公式代数精度为3。7第7页科大研究生学位课程数值分析数值分析5.1 插值型求积公式插值型求积公式思思绪绪利用利用插值多项式插值多项式 ，则积分易算。，则积分易算。在在a,b上取上取 a x0 x1 xn b，做，做 f n 次插值多次插值多项式项式 ，即得到，即得到Ak由由 决定，决定，与与 无关。无关。节点节点 f(x)插值型积分公式插值型积分公式误差误差8第8页科大研究生学位课程数值分析数值分析例：例：对于对于a,b上上1次插值，有次插值，有考查其代数精度。考查其代数精度。f(x)abf(a)f(b)梯形公式梯形公式/*trapezoidal rule*/解解：逐次检验公式是否准确成立逐次检验公式是否准确成立代入代入 P0=1：=代入代入 P1=x：=代入代入 P2=x2：代数精度代数精度=1定理定理：形如形如 求积公式最少有求积公式最少有 n 次代数精度次代数精度 该该公式公式为为插值型插值型（即：（即：）9第9页科大研究生学位课程数值分析数值分析 为了简化计算,取等距节点xk=a+kh,(k=0,1,2,n,则有 令 则有称为Newton-Cotes公式公式.Ck(n)称为Cotes系数.(5.6)它不但与它不但与函数函数f(x)无关，而且与无关，而且与积分区间积分区间a,b无关。无关。10第10页科大研究生学位课程数值分析数值分析 例例1 设(x)C2a,b,求n=1时Newton-Cotes公式并预计误差.解解 计算Cotes系数于是有5.2 几个惯用求积公式几个惯用求积公式从几何上看:用梯形面积近似曲边梯形面积。所以公式=T也称为梯形公式梯形公式,记为T.11第11页科大研究生学位课程数值分析数值分析称之为Simpson公式公式或抛物线公式抛物线公式，记为S.结构三次多项式H3(x),使满足 H3(a)=(a),H3(b)=(b),于是有轻易证实Simpson公式对不高于三次三次多项式准确成立准确成立,即这时插值误差为=S.例例2.设(x)C4a,b,求n=2时Newton-Cotes公式并预计误差.解解 计算Cotes系数12第12页科大研究生学位课程数值分析数值分析于是有13第13页科大研究生学位课程数值分析数值分析 因为结构Newton-Cotes公式需要Cotes系数,将其列表以下:14第14页科大研究生学位课程数值分析数值分析 （3）牛顿求积公式：代数精度代数精度=315第15页科大研究生学位课程数值分析数值分析例例3 求n=4Newton-Cotes公式及误差.解解 查表可得于是有称之为Cotes公式，公式，记为C。其误差为其中,xk=a+kh,k=0,1,2,3,4,h=(b-a)/4.代数精度代数精度=516第16页科大研究生学位课程数值分析数值分析 普通地,Newton-Cotes公式截断误差为 例例4 用梯形公式、Simpson公式和Cotes公式求积分 近似值。解解 IT=1/2*(4+2)=3IS=1/6*(4+12.8+2)=3.13333IC=1/90*(28+14)=3.1421217第17页科大研究生学位课程数值分析数值分析5.3 复化求积公式复化求积公式高次插值有高次插值有Runge 现象现象，故采取分段低次插值，故采取分段低次插值 分段低次合成分段低次合成 Newton-Cotes 复化复化求积公式。求积公式。复化梯形公式：复化梯形公式：在每个在每个 上用梯形公式：上用梯形公式：=Tn18第18页科大研究生学位课程数值分析数值分析可见,复化梯形公式是收敛。而且,要使|RTn|,只要假如记M2=复化梯形公式误差为，则有 若在每个小区间上积分采取Simpson公式,则可得到复化复化Simpson公式公式:19第19页科大研究生学位课程数值分析数值分析 复化复化 Simpson 公式：公式：44444=Sn误差为复化Simpson公式也是收敛。20第20页科大研究生学位课程数值分析数值分析假如记M4=，则有复化Simpson公式也是收敛,而且,要使|RSn|,只要21第21页科大研究生学位课程数值分析数值分析 例例 已知函数分别用复化梯形公式、复化Simpson公式计算积分 解解数据表xk(xk)xk(xk)xk(xk)013/80.97672673/40.90885171/80.99739781/20.95885117/80.87719261/40.98961585/80.936155610.8414710近似值。I准确到数点后7位值是0.9460831。22第22页科大研究生学位课程数值分析数值分析 例例 利用复化梯形公式和复化Simpson公式分别计算上例中定积分,若使精度=10-6,问各需取n为多少?解解 因为(x)=,所以有于是有对复化梯形公式,若使|RTn|10-6,只要故应取n=167.对复化Simpson公式,若使|RSn|10-6,只要 故只需取n=3.实际上,S3=0.9460838.23第23页科大研究生学位课程数值分析数值分析变步长求积方法变步长求积方法 实际积分计算问题，极难依据误差实际积分计算问题，极难依据误差|Rf|0时时,总有总有Rf-0.这说明这说明,只需只需 h 充分小充分小,必可满足误差要求。必可满足误差要求。所以为计算积分，通常采取所以为计算积分，通常采取逐步缩小步长逐步缩小步长方法。方法。利用两种步长计算积分时，为了降低计算函数f(x)次数，通常取 h*=h/2.比如应用复化梯形求积公式时，注意当前步长为h时,有 即先任取步长即先任取步长h 进行计算，然后取较小步长进行计算，然后取较小步长h*进行计算进行计算，假如两次计算结果相差较大，则取更小步长进行计算，如此假如两次计算结果相差较大，则取更小步长进行计算，如此下去，直到相邻两次计算结果相差不大为止，取最小步长算下去，直到相邻两次计算结果相差不大为止，取最小步长算出结果作为积分值。出结果作为积分值。这种方法称为变步长积分法。这种方法称为变步长积分法。24第24页科大研究生学位课程数值分析数值分析可见步长减半时 这表明算出T(h)后，为算T(h/2)，只需计算新增节点xi-1/2=a+(i-1/2)h(i=1,n)处函数值f(xi-1/2)，将它们和乘新步长h/2，再加上T(h)二分之一。利用T(h)和T(h*)还可近似误差预计,称之事后误差预计事后误差预计.25第25页科大研究生学位课程数值分析数值分析对于复化梯形公式n等分区间2n等分区间近似有：引入龙贝格求积方法。26第26页科大研究生学位课程数值分析数值分析由此得因为 首先，若|T2n-Tn|3,则有近似误差|I*-T2n|.5.4 Romberg求积公式求积公式所以有 其次，(4T2n-Tn)/3应比Tn和T2n近似程度更好.事实上,有27第27页科大研究生学位课程数值分析数值分析其中,xk=a+kh,k=0,1,2,n,h=(b-a)/n,而且有于是有所以有逐次分半复化梯形公式递推公式：而且，要使 ,只要=Sn28第28页科大研究生学位课程数值分析数值分析复化Simpson公式能加工成更高精度公式吗?由复化Simpson公式误差预计式有：29第29页科大研究生学位课程数值分析数值分析所以有由此得 首先，若|S2n-Sn|15,则有近似误差|I*-S2n|.其次,(16S2n-Sn)/15应比Sn和S2n近似程度更好.事实上,有(16S2n-Sn)/15=Cn类似地,因为30第30页科大研究生学位课程数值分析数值分析所以有由此得 首先，若|C2n-Cn|63,则有近似误差|I-C2n|.其次,(64C2n-Cn)/63应比Cn和C2n近似程度更好.记(64C2n-Cn)/63=Rn,称为Romberg求积公式求积公式.31第31页科大研究生学位课程数值分析数值分析 用Tm(k)(m=1,2,3,4)分别表示把区间2k等分复化梯形公式,复化Simpson公式,复化Cotes公式和Romberg求积公式.而且，要使|I*-Tm(k)|,只要|Tm(k)-Tm(k-1)|(4m-1)(m=1,2,3,4).则有若对Romberg求积公式作组合也有 32第32页科大研究生学位课程数值分析数值分析普通有：普通有：Romberg 序列序列 Romberg 算法：算法：?T1=)0(1T T8=)3(1T T4=)2(1T T2=)1(1T S1=)0(2T R1=)0(4T S2=)1(2T C1=)0(3T C2=)1(3T S4=)2(2T33第33页科大研究生学位课程数值分析数值分析 实际计算可按下表次序进行 k区间等分数n=2k梯形公式T1(k)Simpson公式T2(k)Cotes公式T3(k)Romberg公式T4(k)01234124816T1(0)T1(1)T1(2)T1(3)T1(4)T2(0)T2(1)T2(2)T2(3)T3(0)T3(1)T3(2)T4(0)T4(1)例例 利用Romberg积分公式计算积分 34第34页科大研究生学位课程数值分析数值分析 解解 按递推公式计算,结果以下可见,因为|T1(4)-T1(3)|=0.0019531,应有|I-T1(4)|0.000651033.kn=2kT1(k)T2(k)T3(k)T4(k)012341248163.00000003.10000003.13117653.13898853.14094163.13333333.14156873.14159253.14159263.14211773.14159413.14159263.14158583.1415926因为|T2(3)-T2(2)|=0.0000001,应有|I-T2(3)|0.00000000666.因为|T3(2)-T3(1)|=0.0000015,应有|I-T3(2)|0.0000000238.因为|T4(1)-T4(0)|=0.0000068,应有|I-T4(1)|=2n=2时就极难求解时就极难求解.故普通不经过解非线性方程求故普通不经过解非线性方程求 ，而从分析高斯点特征来结构高斯求积公式而从分析高斯点特征来结构高斯求积公式.此方法称为用待定系数法结构高斯求积公式此方法称为用待定系数法结构高斯求积公式.利用正交多项式结构高斯求积公式利用正交多项式结构高斯求积公式.39第39页科大研究生学位课程数值分析数值分析 定理定理是高斯点充分必要条件是以这些节点为零点多项式是高斯点充分必要条件是以这些节点为零点多项式与任何次数不超出与任何次数不超出 多项式多项式 带权带权 正交，正交，即插值型求积公式插值型求积公式(5.1)(5.1)节点节点为了结构高斯求积公式节点为了结构高斯求积公式节点,有下述结论有下述结论.40第40页科大研究生学位课程数值分析数值分析(2)求出pn(x)n个零点x0,x1,xn 即为Gsuss点.(1)求出区间a,b上权函数为正交多项式pn(x).(3)再计算积分系数 Gauss型求积公式结构方法型求积公式结构方法一是采取施密特正交化方法.另是借用现成正交多项式函数组.怎样求正交多项式n(x).41第41页科大研究生学位课程数值分析数值分析解解 按 Schemite 正交化过程作出正交多项式:2点Gauss公式.求积分例：=x42第42页科大研究生学位课程数值分析数值分析故两点Gauss公式为 积分系数为P2(x)两个零点为 一些现成正交多项式组有 43第43页科大研究生学位课程数值分析数值分析 区间-1,1上权函数(x)=1Gauss型求积公式,称为Gauss-Legendre求积公式求积公式,其Gauss点为Legendre多项式零点.几个几个Gauss型求积公式型求积公式 (1)Gauss-Legendre求积公式求积公式 公式Gauss点和求积系数可在数学用表中查到.nxkAknxkAk10260.93246951420.66120938650.23861918610.17132449240.36076157300.467913934620.5773502692130.774596669200.55555555560.888888888970.94910791230.74153118560.405845151400.12948496620.27970539150.38183005050.417959183740.86113631160.33998104360.34785484510.652145154980.96028985650.79666647740.52553240990.18343464250.10122853630.22238103450.31370664590.362683783450.90617984590.538469310100.23692688510.47862867050.568888888944第44页科大研究生学位课程数值分析数值分析例 用3点Gauss公式计算积分 解解 查表得x0=-0.7745966692,x1=0,x2=0.7745966692,A0=A2=0.5555555556,A1=0.8888888889,所以有 Gauss-Legendre求积公式余项为 误差为 实际上,I*=2sin1=1.68294197,误差为|R|=6.15810-5.用Simpson公式,则有I*1.69353487,误差为|R|=1.0610-2.45第45页科大研究生学位课程数值分析数值分析因为所以,a,b上权函数(x)=1Gauss型求积公式为 例 用3点Gauss公式计算积分结果远比Simpson公式结果准确.解解 这里Gauss点和积分系数与上例相同,所以 求积误差可表示为46第46页科大研究生学位课程数值分析数值分析 区间0,)上权函数(x)=e-xGauss型求积公式,称为Gauss-Laguerre求积公式求积公式,其Gauss点为Laguerre多项式零点.(2)Gauss-Laguerre求积公式求积公式 公式Gauss点和求积系数可在数学用表中查到.nxkAknxkAk20.58588643763.41421356230.85355339050.146446609450.26356031971.41340305913.59642577107.085810005812.64080084420.52175561050.39866681100.07594244970.00361175870.000023370030.41577455672.29428036026028994508290.71109300990.27851773350.010389256560.22284660411.18893210162.99273632605.77514356919.837467418315.98287398060.45896467930.41700083070.11337338200.01039919750.00026101720.000000898540.32254768961.74576110114.53662029699.39507091230.60315410430.35741869240.03888790850.000539294747第47页科大研究生学位课程数值分析数值分析Gauss-Laguerre求积公式为 求积公式误差为 因为 所以,对0,+)上权函数(x)=1积分,也能够结构类似Gauss-Laguerre求积公式:48第48页科大研究生学位课程数值分析数值分析 区间(-,)上权函数(x)=Gauss型求积公式，(3)Gauss-Hermite求积公式求积公式 公式Gauss点和求积系数可在数学用表中查到.nxkAknxkAk20.70710678110.886226925460.43607741191.33584907042.35060497360.72462959520.15706732030.004530009931.224744871300.29540897511.816359000640.52464762321.65068012380.80491409000.081312835470.81628788281.67355162872.651961356300.42560725260.05451558280.00097178120.81026461755l0.9585724646l2.028704l00.39361932310.01995324210.9453087204称为Gauss-Hermite求积公式求积公式,其Gauss点为Hermite多项式零点.Gauss-Hermite求积公式为49第49页科大研究生学位课程数值分析数值分析数值微分数值微分就是用离散方法近似地求出函数在某点导数值就是用离散方法近似地求出函数在某点导数值.按照按照Taylor展开原理可得展开原理可得其中其中h步长。步长。5.6 数值微分数值微分 后一个数值微分方法称为中点方法，它其实是前两种后一个数值微分方法称为中点方法，它其实是前两种方法算术平均方法算术平均.但它误差阶却由但它误差阶却由 提升到提升到 50第50页科大研究生学位课程数值分析数值分析 设Ln(x)是(x)以ax0 x1xnb为节点n次Lagrange插值多项式,则取当(x)Cn+1+ka,b时,有5.6.2 插值型数值微分插值型数值微分尤其,当k=1时有假如仅限定在节点xi处求导,则有 51第51页科大研究生学位课程数值分析数值分析如取n=1线性插值L1(x)=(x-x0)(x1)-(x-x1)(x0)/h,(其中h=x1-x0)可得数值微分二点公式:如取n=2等距节点(x2-x1=x1-x0=h)抛物线插值:L2(x)=(x-x1)(x-x2)(x0)-2(x-x0)(x-x2)(x1)+(x-x0)(x-x1)(x2)/2h2 则有 L2(x)=(2x-x1-x2)(x0)-2(2x-x0-x2)(x1)+(2x-x0-x1)(x2)/2h2 L2(x)=2(x0)-4(x1)+2(x2)/2h2 L2(x)=(x0)-2(x1)+(x2)/h2 52第52页科大研究生学位课程数值分析数值分析可得数值微分三点公式:53第53页科大研究生学位课程数值分析数值分析练习练习 题题第第116页页 习题习题55.1-5.4,5.7-5.105.13-5.14 54第54页