第2章数值积分优秀课件.ppt

《第2章数值积分优秀课件.ppt》由会员分享，可在线阅读，更多相关《第2章数值积分优秀课件.ppt（62页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第2章数值积分第1页，本讲稿共62页数值微分1.函数f(x)以离散点列给出时，而要求我们给出导数值，2.函数f(x)过于复杂这两种情况都要求我们用数值的方法求函数的导数值微积分中，关于导数的定义如下：自然，而又简单的方法就是，取极限的近似值，即差商第2页，本讲稿共62页向前差商x0 x0+h第3页，本讲稿共62页由Taylor展开因此，有误差第4页，本讲稿共62页向后差商x0-hx0第5页，本讲稿共62页由Taylor展开因此，有误差第6页，本讲稿共62页中心差商x0-hx0 x0+h第7页，本讲稿共62页由Taylor展开因此，有误差第8页，本讲稿共62页由误差表达式，h越小，误差越小，但同

2、时舍入误差增大，所以，有个最佳步长我们可以用事后误差估计的方法来确定设D(h),D(h/2)分别为步长为h,h/2 的差商公式。则时的步长h/2 就是合适的步长第9页，本讲稿共62页f(x)=exp(x)h f(1.15)R(x)h f(1.15)R(x)0.10 3.1630-0.0048 0.05 3.1590-0.00080.09 3.1622-0.0040 0.04 3.1588-0.00060.08 3.1613-0.0031 0.03 3.1583-0.00010.07 3.1607-0.0025 0.02 3.1575-0.00070.06 3.1600-0.0018 0.01

3、3.1550-0.0032例：第10页，本讲稿共62页 插值是建立逼近函数的手段，用以研究原函数的性质。因此，可以用插值函数的导数近似为原函数的导数误差插值型数值微分第11页，本讲稿共62页给定点列 且，求解：例：第12页，本讲稿共62页Taylor展开分析，可以知道，它们都是 称为三点公式第13页，本讲稿共62页数值积分关于积分，有Newton-Leibniz 公式但是，在很多情况下，还是要数值积分：1、函数有离散数据组成2、F(x)求不出3、F(x)非常复杂定义数值积分如下：是离散点上的函数值的线性组合称为积分系数，与f(x)无关，与积分区间和积分点有关第14页，本讲稿共62页为数值积分，

4、为积分，则称数值积分有k阶代数精度是指：两个问题：1、系数ai如何选取，即选取原则2、若节点可以自由选取，取什么点好？代数精度 对任意次数不高于k次的多项式f(x)，数值积分没有误差第15页，本讲稿共62页用插值函数的积分，作为数值积分代数精度由Lagrange插值的误差表达式，有可以看出，至少n 阶代数精度插值型第16页，本讲稿共62页Vandermonde行列式使用尽可能高的代数精度已知求系数所以，要存在唯一，m n，确定一个n 1 阶的方程组待定系数法第17页，本讲稿共62页所以，m=n时存在唯一，且至少n阶代数精度。与节点的选取有关。若数值积分至少n阶代数精度，则系数唯一误差第18页，

5、本讲稿共62页一点数值积分0阶代数精度1阶代数精度例：第19页，本讲稿共62页Newton-Cotes 积分若节点可以自由选取，则，一个自然的办法就是取等距节点。对区间做等距分割。该数值积分称为Newton-Cotes积分第20页，本讲稿共62页设节点步长(b-a)与步长h无关，可以预先求出第21页，本讲稿共62页N 1时梯形公式第22页，本讲稿共62页N 2 时Simpson公式第23页，本讲稿共62页1、梯形公式此处用了积分中值定理误差第24页，本讲稿共62页2、Simpson 公式 注意到，Simpson 公式有3 阶代数精度，因此为了对误差有更精确地估计，我们用3 次多项式估计误差为0

6、第25页，本讲稿共62页一般的有因此，N-C积分，对偶数有n 1 阶代数精度，而奇数为n 阶代数精度第26页，本讲稿共62页复化积分数值积分公式与多项式插值有很大的关系。因此Runge现象的存在，使得我们不能用太多的积分点计算。采用与插值时候类似，我们采用分段、低阶的方法第27页，本讲稿共62页误差做等距节点，复化梯形公式第28页，本讲稿共62页由均值定理知可以看出，复化梯形公式是收敛的。第29页，本讲稿共62页误差做等距节点，复化Simpson公式第30页，本讲稿共62页由均值定理知可以看出，复化Simpson公式是收敛的。第31页，本讲稿共62页 若 一 个 积 分 公 式 的 误 差 满

7、 足 且C 0，则 称该公式是 p 阶收敛的。例：计算解：其中=3.138988494其中=3.141592502运算量基本相同第32页，本讲稿共62页 函数变化有急有缓，为了照顾变化剧烈部分的误差，我们需要加密格点。对于变化缓慢的部分，加密格点会造成计算的浪费。以此我们介绍一种算法，可以自动在变化剧烈的地方加密格点计算，而变化缓慢的地方，则取稀疏的格点。积分的自适应计算第33页，本讲稿共62页先看看事后误差估计以复化梯形公式为例n等分区间2n等分区间近似有：类似，复化Simpson公式第34页，本讲稿共62页自适应计算记为复化一次，2次的Simpson公式控制求第35页，本讲稿共62页是第3

8、6页，本讲稿共62页由前面的事后误差估计式，则，这启发我们，可以用低阶的公式组合后称为一个高阶的公式。类似，Romberg 积分第37页，本讲稿共62页记为以步长为h的某数值积分公式，有第38页，本讲稿共62页有如下的Euler-Maclaurin定理若 为2m阶公式，则Romberg 积分就是不断地用如上定理组合低阶公 式为高阶公式，进而计算积分 Romberg 算法：?T1=)0(0T T8=)3(0T T4=)2(0T T2=)1(0T S1=)0(1T R1=)0(3T S2=)1(1T C1=)0(2T C2=)1(2T S4=)2(1T第39页，本讲稿共62页Lab03 复化积分1

9、.分别编写用复化Simpson积分公式和复化梯形积分公式计算积分的通用程序2.用如上程序计算积分取节点xi,i=0,N，N 为 2k,k=0,1,12，并估计误差3.简单分析你得到的数据第40页，本讲稿共62页重积分的计算 在微积分中，二重积分的计算是用化为累次积分的方法进行的。计算二重数值积分也同样采用累次积分的计算过程。简化起见，我们仅讨论矩形区域上的二重积分。对非矩形区域的积分，大多可以变化为矩形区域上的累次积分。a,b,c,d 为常数，f 在D 上连续。将它变为化累次积分首先来看看复化梯形公式的二重推广第41页，本讲稿共62页做等距节点，x轴，y轴分别有：先计算，将x作为常数，有再将y

10、作为常数，在x方向，计算上式的每一项的积分二重积分的复化梯形公式第42页，本讲稿共62页第43页，本讲稿共62页系数，在积分区域的四个角点为1/4，4个边界为1/2，内部节点为1误差第44页，本讲稿共62页类似前面有：记二重积分的复化Simpson 公式做等距节点，x轴，y轴分别有：m，n为偶数第45页，本讲稿共62页误差第46页，本讲稿共62页Gauss型积分公式 Newton-Cotes 积分公式，可以知道n为偶数时，n 1个点数值积分公式有n 1阶精度。是否有更高的代数精度呢？n个点的数值积分公式，最高可以到多少代数精度？本节会解决这个问题。第47页，本讲稿共62页例：在两点数值积分公式

11、中，如果积分点也作为未知量，则有4个未知量，可以列出4个方程：（以f(x)在-1,1为例）可解出：数值积分公式具有3阶代数精度，比梯形公式1阶代数精度高第48页，本讲稿共62页证明：取易知：也就是说，数值积分公式，对一个2n+2阶的多项式是有误差的，所以，n 1个点的数值积分公式不超过2n 1阶n 个积分点的数值积分公式，最高2n 1阶 定理如何构造最高阶精度的公式？第49页，本讲稿共62页一般性，考虑积分：称为权函数定义两个可积函数的内积为：两个函数正交，就是指这两个函数的内积为0第50页，本讲稿共62页利用Schmidt 正交化过程，变为正交基就可以将多项式基函数第51页，本讲稿共62页以

12、n阶正交多项式的n个零点为积分点的数值积分公式有2n 1阶的代数精度Gauss点Gauss积分，记为Gn(f)证明：若 f 为 2n 1 次多项式，则 为 n 1 次多项式又，仅差一个常数（零点相同）具有一个很好的性质：第52页，本讲稿共62页(2)求出pn(x)的n个零点x1,x2,xn 即为Gauss点.(1)求出区间a,b上权函数为W(x)的正交多项式pn(x).(3)计算积分系数 Gauss型求积公式的构造方法第53页，本讲稿共62页解 按 Schemite 正交化过程作出正交多项式:的2点Gauss公式.求积分例：第54页，本讲稿共62页故两点Gauss公式为 积分系数为P2(x)的

13、两个零点为 第55页，本讲稿共62页 区间-1,1上权函数W(x)=1的Gauss型求积公式,称为Gauss-Legendre求积公式,其Gauss点为Legendre多项式的零点.(1)Gauss-Legendre 求积公式公式的Gauss点和求积系数可在数学用表中查到.几种Gauss 型求积公式由因此,a,b上权函数W(x)=1的Gauss型求积公式为第56页，本讲稿共62页n xkAkn xkAk1 0 260.93246951420.66120938650.23861918610.17132449240.36076157300.46791393462 0.5773502692130.7

14、74596669200.55555555560.888888888970.94910791230.74153118560.405845151400.12948496620.27970539150.38183005050.417959183740.86113631160.33998104360.34785484510.652145154980.96028985650.79666647740.52553240990.18343464250.10122853630.22238103450.31370664590.362683783450.90617984590.538469310100.236926

15、88510.47862867050.5688888889第57页，本讲稿共62页 Gauss 公式的余项：/*设P 为f 的过x0 xn的插值多项式*/*只要P 的阶数不大于2n+1，则下一步等式成立*/插值多项式的余项Q：什么样的插值多项式在 x0 xn 上有 2n+1 阶？第58页，本讲稿共62页A：Hermite 多项式！满足第59页，本讲稿共62页 区间0,)上权函数W(x)=e-x的Gauss型求积公式,称为Gauss-Laguerre 求积公式,其Gauss点为Laguerre多项式的零点.(2)Gauss-Laguerre 求积公式公式的Gauss点和求积系数可在数学用表中查到.

16、由所以,对0,+)上权函数W(x)=1的积分,也可以构造类似的Gauss-Laguerre求积公式:第60页，本讲稿共62页n xkAkn xkAk2 0.58588643763.41421356230.85355339050.146446609450.26356031971.41340305913.59642577107.085810005812.64080084420.52175561050.39866681100.07594244970.00361175870.000023370030.41577455672.29428036026028994508290.71109300990.278

17、51773350.010389256560.22284660411.18893210162.99273632605.77514356919.837467418315.98287398060.45896467930.41700083070.11337338200.01039919750.00026101720.000000898540.32254768961.74576110114.53662029699.39507091230.60315410430.35741869240.03888790850.0005392947第61页，本讲稿共62页(3)Gauss-Hermite 求积公式公式的Ga

18、uss点和求积系数可在数学用表中查到.n xkAkn xkAk2 0.7071067811 0.886226925460.43607741191.33584907042.35060497360.72462959520.15706732030.00453000993 1.224744871300.29540897511.81635900064 0.52464762321.65068012380.80491409000.081312835470.81628788281.67355162872.651961356300.42560725260.05451558280.00097178120.810264617550.95857246462.020182870400.39361932310.01995324210.9453087204区间(-,)上权函数W(x)=的Gauss型求积公式,称为Gauss-Hermite 求积公式,其Gauss点为Hermite多项式的零点.第62页，本讲稿共62页