《六章数值微分和数值积分》由会员分享，可在线阅读，更多相关《六章数值微分和数值积分（58页珍藏版）》请在金锄头文库上搜索。
1、泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College第六章 数值微分和数值积分泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College数值微分1. 函数f(x)以离散点列给出时，而要求我们给出导数值，2. 函数f(x)过于复杂这两种情况都要求我们用数值的方法求函数的导数值微积分中，关于导数的定义如下：自然，而又简单的方法就是，取极限的近似值，即差商泰山学院 信息科学技术系Department of
2、 Information Science and Technology, Taishan College由Taylor展开因此，有误差向前差商泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College由Taylor展开因此，有误差向后差商泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College由Taylor展开因此，有误差中心差商泰山学院 信息科学技术系Department of Information
3、 Science and Technology, Taishan College由误差表达式，h越小，误差越小，但同时舍入误差增大，所以，有个最佳步长我们可以用事后误差估计的方法来确定设D(h),D(h/2)分别为步长为h,h/2的差商公式。则时的步长h/2就是合适的步长泰山学院 信息科学技术系Department of Information Science and Technology, Taishan Collegef(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
4、.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例：泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College 插值是建立逼近函数的手段，用以研究原函数的性质。因此，可以用插值函数的导数近似为原函数的导数误差插值型数值微分泰山学院 信息科学技术系Department of Information Science a
5、nd Technology, Taishan College给定点列且，求解：例：泰山学院 信息科学技术系Department of Information Science and Technology, Taishan CollegeTaylor展开分析，可以知道，它们都是称为三点公式三点公式泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College数值积分关于积分，有Newton-Leibniz公式但是，在很多情况下，还是要数值积分：1、函数有离散数据组成2、F(x)求不出3、F(x)非常复杂定
6、义数值积分如下：是离散点上的函数值的线性组合称为积分系数，与f(x)无关，与积分区间和积分点有关泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College代数精度为数值积分，为积分，则称数值积分有k阶代数精度是指：两个问题：1、系数ai如何选取，即选取原则2、若节点可以自由选取，取什么点好？定义 对任意次数不高于k次的多项式f(x)，数值积分没有误差泰山学院 信息科学技术系Department of Information Science and Technology, Taishan Colleg
7、e用插值函数的积分，作为数值积分代数精度代数精度由Lagrange插值的误差表达式，有可以看出，至少n阶代数精度插值型泰山学院 信息科学技术系Department of Information Science and Technology, Taishan CollegeVandermonde行列式使用尽可能高的代数精度已知求系数所以，要存在唯一，mn，确定一个n1阶的方程组是否有更好的方法使得代数精度为至少为n+1阶？泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College所以，m=n时存在唯一
8、，且至少n阶代数精度。与节点的选取有关。若数值积分至少n阶代数精度，则系数唯一误差泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College一点数值积分0阶代数精度1阶代数精度例：泰山学院 信息科学技术系Department of Information Science and Technology, Taishan CollegeNewton-Cotes Newton-Cotes 积分积分若节点可以自由选取，则，一个自然的办法就是取等距节点。对区间做等距分割。该数值积分称为Newton-Cotes
9、积分泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College设节点步长(b-a)与步长h无关，可以预先求出泰山学院 信息科学技术系Department of Information Science and Technology, Taishan CollegeN1时梯形公式泰山学院 信息科学技术系Department of Information Science and Technology, Taishan CollegeN2时Simpson公式泰山学院 信息科学技术系Department of
10、 Information Science and Technology, Taishan College1、梯形公式此处用了积分中值定理误差泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College2、Simpson公式 注意到，Simpson公式有3阶代数精度，因此为了对误差有更精确地估计，我们用3次多项式估计误差为0泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College一般的有因此，N-C积分，对
11、偶数有n1阶代数精度，而奇数为n阶代数精度泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College复化积分数值积分公式与多项式插值有很大的关系。因此Runge现象的存在，使得我们不能用太多的积分点计算。采用与插值时候类似，我们采用分段、低阶的方法泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College误差误差做等距节点，复化梯形公式泰山学院 信息科学技术系Department of Informati
12、on Science and Technology, Taishan College由均值定理知可以看出，复化梯形公式是收敛的。泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College误差误差做等距节点，复化Simpson公式泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College由均值定理知可以看出，复化Simpson公式是收敛的。泰山学院 信息科学技术系Department of Informat
13、ion Science and Technology, Taishan College 若若一一个个积积分分公公式式的的误误差差满满足足 且且C 0，则则称该公式是称该公式是 p 阶收敛阶收敛的。的。例：例：计算计算解：解：其中其中= 3.138988494其中其中= 3.141592502运算量基运算量基本相同本相同泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College3、积分的自适应计算 函数变化有急有缓，为了照顾变化剧烈部分的误差，我们需要加密格点。对于变化缓慢的部分，加密格点会造成计算的
14、浪费。以此我们介绍一种算法，可以自动在变化剧烈的地方加密格点计算，而变化缓慢的地方，则取稀疏的格点。泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College先看看事后误差估计 （不同的误差表达式，事后误差估计式是不同的）以复化梯形公式为例n等分区间2n等分区间近似有：类似，复化Simpsom公式泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College自适应计算记为复化一次，2次的Simpson公式控制
15、求泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College是泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College4、Romberg积分由前面的事后误差估计式，则，这启发我们，可以用低阶的公式组合后称为一个高阶的公式。类似，泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College记为以步长为h的某数值积
16、分公式，有泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College有如下的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泰山学院 信息科学技术系Department of Info
17、rmation Science and Technology, Taishan College重积分的计算 在微积分中，二重积分的计算是用化为累次积分的方法进行的。计算二重数值积分也同样采用累次积分的计算过程。简化起见，我们仅讨论矩形区域上的二重积分。对非矩形区域的积分，大多可以变化为矩形区域上的累次积分。a,b,c,d为常数，f在D上连续。将它变为化累次积分首先来看看复化梯形公式的二重推广泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College做等距节点，x轴，y轴分别有：先计算，将x作为常数，
18、有再将y作为常数，在x方向，计算上式的每一项的积分二重积分的复化梯形公式二重积分的复化梯形公式泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College系数，在积分区域的四个角点为1/4，4个边界为1/2，内部节点为1误差泰山学院 信息科学技术系Department of Information Science and Technology, Ta
19、ishan College类似前面有：记二重积分的复化二重积分的复化Simpson公式公式做等距节点，x轴，y轴分别有：m，n为偶数泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College误差泰山学院 信息科学技术系Department of Information Science and Technology, Taishan CollegeGauss型积分公式 Newton-Cotes积分公式，可以知道n为偶数时，n1个点数值积分公式有n1阶精度。是否有更高的代数精度呢？n个点的数值积分公式，
20、最高可以到多少代数精度？本节会解决这个问题。泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College例：在两点数值积分公式中，如果积分点也作为未知量，则有4个未知量可以列出4个方程： （以f(x)在-1,1为例）可解出：可以看出，数值积分公式具有3阶代数精度，比梯形公式1阶代数精度高泰山学院 信息科学技术系Department of Information Science and Technology, Taishan Collegen个积分点的数值积分公式，最高2n1阶证明：取易知：也就是说，数
21、值积分公式，对一个2n+2阶的多项式是有误差的，所以，n1个点的数值积分公式不超过2n1阶如何构造最高阶精度的公式？定理泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College一般性，考虑积分：称为权函数定义两个可积函数的内积为：两个函数正交，就是指这两个函数的内积为0泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College利用Schmidt正交化过程，变为正交基就可以将多项式基函数泰山学院 信息科学
22、技术系Department of Information Science and Technology, Taishan College以n阶正交多项式的n个零点为积分点的数值积分公式有2n1阶的代数精度Gauss点Gauss积分，记为Gn(f)证明：若f为2n1次多项式，则为n1次多项式又，仅差一个常数（零点相同）具有一个很好的性质：泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College(2)求出pn(x)的n n个零点x1 , x2 , xn 即为Gsuss点. (1)求出区间a,b上权函数
23、为W(x)的正交多项式pn(x) .(3)计算积分系数 Gauss型求积公式的构造方法泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College解解 按 Schemite 正交化过程作出正交多项式: 的2点Gauss公式.求积分例：泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College故两点Gauss公式为 积分系数为P2(x)的两个零点为 泰山学院 信息科学技术系Department of Info
24、rmation Science and Technology, Taishan College 区间-1,1上权函数W(x)=1的Gauss型求积公式,称为Gauss-Legendre求积公式求积公式,其Gauss点为Legendre多项式的零点. (1) Gauss-Legendre求积公式求积公式公式的Gauss点和求积系数可在数学用表中查到 .几种几种Gauss型求积公式型求积公式由因此,a,b上权函数W(x)=1的Gauss型求积公式为泰山学院 信息科学技术系Department of Information Science and Technology, Taishan Colleg
25、enxkAknxkAk10260.93246951420.66120938650.23861918610.17132449240.36076157300.467913934620.5773502692130.774596669200.55555555560.888888888970.94910791230.74153118560.405845151400.12948496620.27970539150.38183005050.417959183740.86113631160.33998104360.34785484510.652145154980.96028985650.79666647740
26、.52553240990.18343464250.10122853630.22238103450.31370664590.362683783450.90617984590.538469310100.23692688510.47862867050.5688888889泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College Gauss 公式的余项：公式的余项：/* 设设P为为f 的过的过x0 xn的插值多项式的插值多项式 */*只要只要P 的阶数不大于的阶数不大于2n+1，则下一步则下一步等式成立等
27、式成立*/插值多项式插值多项式的余项的余项Q：什么样的什么样的插值多项式插值多项式在在 x0 xn 上上有有 2n+1 阶？阶？泰山学院 信息科学技术系Department of Information Science and Technology, Taishan CollegeA：Hermite 多项式！多项式！ 满足满足泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College 区间0,)上权函数W(x)=e-x的Gauss型求积公式,称为Gauss-Laguerre求积公式求积公式,其Gau
28、ss点为Laguerre多项式的零点. (2) Gauss-Laguerre求积公式求积公式公式的Gauss点和求积系数可在数学用表中查到 .由所以,对0, + )上权函数W(x)=1的积分,也可以构造类似的Gauss-Laguerre求积公式:泰山学院 信息科学技术系Department of Information Science and Technology, Taishan CollegenxkAknxkAk20.58588643763.41421356230.85355339050.146446609450.26356031971.41340305913.59642577107.08
29、5810005812.64080084420.52175561050.39866681100.07594244970.00361175870.000023370030.41577455672.29428036026028994508290.71109300990.27851773350.010389256560.22284660411.18893210162.99273632605.77514356919.837467418315.98287398060.45896467930.41700083070.11337338200.01039919750.00026101720.0000008985
30、40.32254768961.74576110114.53662029699.39507091230.60315410430.35741869240.03888790850.0005392947泰山学院 信息科学技术系Department of Information Science and Technology, Taishan College (3) Gauss-Hermite求积公式求积公式公式的Gauss点和求积系数可在数学用表中查到 .nxkAknxkAk20.70710678110.886226925460.43607741191.33584907042.35060497360.7
31、2462959520.15706732030.004530009931.224744871300.29540897511.816359000640.52464762321.65068012380.80491409000.081312835470.81628788281.67355162872.651961356300.42560725260.05451558280.00097178120.810264617550.95857246462.020182870400.39361932310.01995324210.9453087204区间(- ,)上权函数W(x)= 的Gauss型求积公式,称为Gauss-Hermite求积公式求积公式, 其Gauss点为Hermite多项式的零点.
