3数值积分与数值微分(12)

1、2020/8/1,1,数值积分与数值微分,引言,一、 问题的提出,既使函数是以解析的形式给出，但由于其表达式比较复杂，我们在高数中所学的方法很难甚至无法计算出其积分或微分的准确值，如,（1）,（2）,在工程技术和科学研究中，很多情况下变量间的函数关系是以数表的形式给出的，无法运用我们学过的方法计算其积分和导数！,2020/8/1,2,二、 数值微积分概述,利用函数在一些点上的函数值，计算出该函数的积分或微分满足一定精度要求的近似值！,（1）,（2）,数值微积分方法是其它数值方法，如微分方程数值解法等的必备基础。,2020/8/1,3,一、 数值积分公式及其代数精度,数值积分,1、数值积分公式,

2、根据积分中值定理,可以通过在区间a,b内选择 的近似值得到积分的近似计算公式:,(1),(2),(3),(4),矩形公式, 梯形公式,2020/8/1,4,数值求积公式的一般形式为, 求积节点,其中, 求积系数（与f(x)无关）,（1）求积公式由求积节点和求积系数唯一确定；,（1）,注：,（2）求积系数和求积区间、求积节点有关。,2020/8/1,5,为求积公式的截断误差或求积余项。,称,定义:,注: 求积余项刻画了求积公式的计算精度!,2、代数精度,定义:,若求积公式（1）对任何次数不高于m的多项式都准确成立，而对于m+1次的多项式不能准确成立，则称求积公式（1）具有m阶代数精度。,2020

3、/8/1,6,二、 插值型求积公式,设函数f(x)在区间a,b上由定义，已知节点,和相应的函数值,从而可以确定插值函数,1、定义及一般形式,Lagrange插值多项式,以pn(x)作为f(x)的近似，可得, 插值型求积公式,其中,(2),注：求积系数仅与求积区间和求积节点有关!,2020/8/1,7,插值型求积公式的余项,令 ，构造插值型求积公式,2、Newton-Cotes公式,将求积区间a,bn等分，步长为 ，,其中,Newton-Cotes公式,Newton-Cotes系数,2020/8/1,8,事实上，,注: Newton-Cotes系数只与n，k有关,而与f(x)和求积区间a,b无关

4、。并且,2020/8/1,9,几种低阶Newton-Cotes公式的系数,梯形公式(1阶),辛普生(Simpson)公式(3阶),n=4时，称为Cotes公式。(5阶),2020/8/1,10,三、 复合求积公式,考虑到多项式插值的Runge现象，通常不用高阶的插值型积分公式！,为了提高计算精度，我们把积分区间分成若干个子区间，每个子区间上的积分使用少结点的Newton-Cotes求积公式计算，然后再把结果相加，这就是复合求积的思想，所得到的公式就是复合求积公式。,2020/8/1,11,1、复合求积公式,将求积区间a,bn等分，,再相加得到复合梯形公式：,(1) 复合梯形公式,则复合梯形公式

5、的余项为,算法收敛，且数值稳定!,问题9: 实现复合梯形求 积公式，画出流程图！,2020/8/1,12,(2) 复合辛普生(Simpson)公式,再相加得到复合辛普生公式：,将求积区间a,b2n等分，,复合Simpson公式的余项有表达式,算法收敛，且数值稳定!,2020/8/1,13,2、区间逐次分半的思想以复合梯形公式为例,给出精度要求后，一般很难确定把求积区间多少等分，就可以利用复化梯形公式得到所需的积分值。通常等分区间太多则计算量增大，等分区间少，则达不到精度要求！为克服此困难，考虑如下方法:,将当前的每个小求积区间 二等分，从而得到2n个小求积区间，区间长度为 ，再利用复化梯形公式

6、来计算积分值，记为T1,把原来的积分值记为T0,则有,因此，可以由| T1-T0 |作为精度控制条件!,若满足精度要求则停止， T1即为所求，否则重复上述过程。,2020/8/1,14,1、插值型求积算法适用于有限区间a,b上连续函数的积分问题； 2、复合求积公式是数值稳定的，龙贝格积分算法是插值型求积算法中最有效的算法。 3、算法的评价指标：代数精度、求积余项，收敛阶，数值稳定性和函数值的计算次数！ 4、对于振荡积分，即,插值型求积算法不是很理想！此时，先将f(x)用样条函数逼近，再分部积分。,插值型算法总结,2020/8/1,15,数值微分,插值函数数值微分法,思想方法：,以插值多项式近似代替函数，以插值多项式在节点上的导数值近似代替函数在节点上的导数值。,推导公式：,用此方法求微商，可以先求出插值多项式，然后各点上的微商就可以同时求出。,下面给出几个常用的数值微分公式。,设f(x)有插值多项式为Pn(x)，则,2020/8/1,16,2、二阶微商的三点公式,一阶微商的三点公式(留作练习!),1、一阶微商的两点公式,2020/8/1,17,3、一阶微商的五点公式,注：,(1)三点公式与五点公式中有中点项。由于中点的导数值的表达式中不含有中点的函数值项，且函数值项的系数不大，因此选取节点的方法