二章节数值微分和数值积分.ppt

数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS第二章 数值微分和数值积分数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS数值微分1. 函数f(x)以离散点列给出时，而要求我们给出导数值，2. 函数f(x)过于复杂这两种情况都要求我们用数值的方法求函数的导数值微积分中，关于导数的定义如下：自然，而又简单的方法就是，取极限的近似值，即差商数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS由Taylor展开因此，有误差向前差商数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS由Taylor展开因此，有误差向后差商数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS由Taylor展开因此，有误差中心差商数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS由误差表达式，h越小，误差越小，但同时舍入误差增大，所以，有个最佳步长我们可以用事后误差估计的方法来确定设D(h),D(h/2)分别为步长为h,h/2的差商公式。

则时的步长h/2就是合适的步长数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICSf(x)=exp(x)hf’(1.15)R(x)hf’(1.15)R(x)0.103.1630-0.00480.053.1590-0.00080.093.1622-0.00400.043.1588-0.00060.083.1613-0.00310.033.1583-0.00010.073.1607-0.00250.023.1575-0.00070.063.1600-0.00180.013.1550-0.0032例：数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS 插值是建立逼近函数的手段，用以研究原函数的性质因此，可以用插值函数的导数近似为原函数的导数误差插值型数值微分数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS给定点列且，求解：例：数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICSTaylor展开分析，可以知道，它们都是称为三点公式三点公式数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS数值积分关于积分，有Newton-Leibniz公式但是，在很多情况下，还是要数值积分：1、函数有离散数据组成2、F(x)求不出3、F(x)非常复杂定义数值积分如下：是离散点上的函数值的线性组合称为积分系数，与f(x)无关，与积分区间和积分点有关数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS代数精度为数值积分，为积分，则称数值积分有k阶代数精度是指：两个问题：1、系数ai如何选取，即选取原则2、若节点可以自由选取，取什么点好？定义 对任意次数不高于k次的多项式f(x)，数值积分没有误差数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS用插值函数的积分，作为数值积分代数精度代数精度由Lagrange插值的误差表达式，，有可以看出，至少n阶代数精度插值型数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICSVandermonde行列式使用尽可能高的代数精度已知求系数所以，要存在唯一，m＝＝n，确定一个n＋＋1阶的方程组是否有更好的方法使得代数精度为至少为n+1阶？数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS所以，m=n时存在唯一，且至少n阶代数精度。

与节点的选取有关若数值积分至少n阶代数精度，则系数唯一误差数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS一点数值积分0阶代数精度1阶代数精度例：数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICSNewton-Cote’s Newton-Cote’s 积分积分若节点可以自由选取，则，一个自然的办法就是取等距节点对区间做等距分割该数值积分称为Newton-Cote’s积分数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS设节点步长(b-a)与步长h无关，可以预先求出数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICSN＝1时梯形公式数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICSN＝2时Simpson公式数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS1、梯形公式此处用了积分中值定理误差数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS2、Simpson公式 注意到，Simpson公式有3阶代数精度，因此为了对误差有更精确地估计，我们用3次多项式估计误差为0数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS一般的有因此，N-C积分，对偶数有n＋1阶代数精度，而奇数为n阶代数精度数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS复化积分数值积分公式与多项式插值有很大的关系。

因此Runge现象的存在，使得我们不能用太多的积分点计算采用与插值时候类似，我们采用分段、低阶的方法数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS误差误差做等距节点，复化梯形公式数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS由均值定理知可以看出，复化梯形公式是收敛的数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS误差误差做等距节点，复化Simpson公式数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS由均值定理知可以看出，复化Simpson公式是收敛的数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS　　　　 若若一一个个积积分分公公式式的的误误差差满满足足 且且C   0，，则则称该公式是称该公式是 p 阶收敛阶收敛的。

的~~~例：例：计算计算解：解：其中其中= 3.138988494其中其中= 3.141592502运算量基运算量基本相同本相同数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS3、积分的自适应计算 函数变化有急有缓，为了照顾变化剧烈部分的误差，我们需要加密格点对于变化缓慢的部分，加密格点会造成计算的浪费以此我们介绍一种算法，可以自动在变化剧烈的地方加密格点计算，而变化缓慢的地方，则取稀疏的格点数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS①先看看事后误差估计 （不同的误差表达式，事后误差估计式是不同的）以复化梯形公式为例n等分区间2n等分区间近似有：类似，复化Simpsom公式数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS②自适应计算记为复化一次，2次的Simpson公式控制求数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS是数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS4、Romberg积分由前面的事后误差估计式，则，这启发我们，可以用低阶的公式组合后称为一个高阶的公式。

类似，数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS记为以步长为h的某数值积分公式，有数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS有如下的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数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS重积分的计算 在微积分中，二重积分的计算是用化为累次积分的方法进行的。

计算二重数值积分也同样采用累次积分的计算过程简化起见，我们仅讨论矩形区域上的二重积分对非矩形区域的积分，大多可以变化为矩形区域上的累次积分a,b,c,d为常数，f在D上连续将它变为化累次积分首先来看看复化梯形公式的二重推广数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS做等距节点，x轴，y轴分别有：先计算，将x作为常数，有再将y作为常数，在x方向，计算上式的每一项的积分二重积分的复化梯形公式二重积分的复化梯形公式数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS系数，在积分区域的四个角点为1/4，4个边界为1/2，内部节点为1误差数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS类似前面有：记二重积分的复化二重积分的复化Simpson公式公式做等距节点，x轴，y轴分别有：m，n为偶数数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS误差数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICSGauss型积分公式 Newton-Cote’s积分公式，可以知道n为偶数时，n＋＋1个点数值积分公式有n＋＋1阶精度。

是否有更高的代数精度呢？n个点的数值积分公式，最高可以到多少代数精度？本节会解决这个问题数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS例：在两点数值积分公式中，如果积分点也作为未知量，则有4个未知量可以列出4个方程： （以f(x)在[-1,1]为例）可解出：可以看出，数值积分公式具有3阶代数精度，比梯形公式1阶代数精度高数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICSn个积分点的数值积分公式，最高2n－－1阶证明：取易知：也就是说，数值积分公式，对一个2n+2阶的多项式是有误差的，所以，n＋1个点的数值积分公式不超过2n＋1阶如何构造最高阶精度的公式？定理数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS一般性，考虑积分：称为权函数定义两个可积函数的内积为：两个函数正交，就是指这两个函数的内积为0数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS利用Schmidt正交化过程，变为正交基就可以将多项式基函数数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS以n阶正交多项式的n个零点为积分点的数值积分公式有2n－1阶的代数精度Gauss点Gauss积分，记为Gn(f)证明：若f为2n－－1次多项式，则为n－－1次多项式又，仅差一个常数（零点相同）具有一个很好的性质：数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS(2)求出pn(x)的n n个零点x1 , x2 , … xn 即为Gsuss点. (1)求出区间[a,b]上权函数为W(x)的正交多项式pn(x) .(3)计算积分系数 Gauss型求积公式的构造方法数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS解解 按 Schemite 正交化过程作出正交多项式: 的2点Gauss公式.求积分例：数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS故两点Gauss公式为 积分系数为P2(x)的两个零点为 数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS 区间[-1,1]上权函数W(x)=1的Gauss型求积公式,称为Gauss-Legendre求积公式求积公式,其Gauss点为Legendre多项式的零点. (1) Gauss-Legendre求积公式求积公式公式的Gauss点和求积系数可在数学用表中查到 .几种几种Gauss型求积公式型求积公式由因此,[a,b]上权函数W(x)=1的Gauss型求积公式为数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICSnxkAknxkAk1026±0.9324695142±0.6612093865±0.23861918610.17132449240.36076157300.46791393462±0.577350269213±0.774596669200.55555555560.88888888897±0.9491079123±0.7415311856±0.405845151400.12948496620.27970539150.38183005050.41795918374±0.8611363116±0.33998104360.34785484510.65214515498±0.9602898565±0.7966664774±0.5255324099±0.18343464250.10122853630.22238103450.31370664590.36268378345±0.9061798459±0.538469310100.23692688510.47862867050.5688888889数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS  Gauss 公式的余项：公式的余项：/* 设设P为为f 的过的过x0 … xn的插值多项式的插值多项式 *//*只要只要P 的阶数不大于的阶数不大于2n+1，，则下一步则下一步等式成立等式成立*/插值多项式插值多项式的余项的余项Q：：什么样的什么样的插值多项式插值多项式在在 x0 … xn 上上有有 2n+1 阶？阶？数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICSA：：Hermite 多项式！多项式！ 满足满足数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS 区间[0,)上权函数W(x)=e-x的Gauss型求积公式,称为Gauss-Laguerre求积公式求积公式,其Gauss点为Laguerre多项式的零点. (2) Gauss-Laguerre求积公式求积公式公式的Gauss点和求积系数可在数学用表中查到 .由所以,对[0, + )上权函数W(x)=1的积分,也可以构造类似的Gauss-Laguerre求积公式:数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICSnxkAknxkAk20.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.0005392947数 学 系University of Science and Technology of ChinaDEPARTMENT OF MATHEMATICS (3) Gauss-Hermite求积公式求积公式公式的Gauss点和求积系数可在数学用表中查到 .nxkAknxkAk2±0.70710678110.88622692546±0.4360774119±1.3358490704±2.35060497360.72462959520.15706732030.00453000993±1.224744871300.29540897511.81635900064±0.5246476232±1.65068012380.80491409000.08131283547±0.8162878828±1.6735516287±2.651961356300.42560725260.05451558280.00097178120.81026461755±0.9585724646±2.020182870400.39361932310.01995324210.9453087204区间(- ,)上权函数W(x)= 的Gauss型求积公式,称为Gauss-Hermite求积公式求积公式, 其Gauss点为Hermite多项式的零点. 。