陕西师范大学《数值分析》第六章 数值积分与数值微分.pptxVIP

讲授：

用计算机求定积分的方法和估计导数的方法；

重点论述:

插值型求积公式、Newton-Cotes求积公式、Gauss求积公式、复化求积公式、Romberg求积公式及对应的原理、构造、误差分析等。第六章数值积分与数值微分

§6.2基本概念第六章数值积分与微分

问题：定积分计算有公式但不能用公式的定积分要怎样计算呢？怎样把定积分在一定精度下借助计算机计算出来就是本章研究的内容。此外，利用函数在若干个点处的函数值求该函数的导数近似值也是本章介绍的内容。?本章涉及的方法：Newton-Cotes求积公式、Gauss求积公式、复化求积公式、Romberg求积公式等。

1、问题的描述积分离散化处理过程：1）去掉极限号3)减少离散化的误差，做得到离散化形式：

1）求积公式：若存在实数成立称为一个数值求积公式。求积系数求积节点求积余项：求积公式与求积节点个数、求积节点和求积系数有关，显然有很多不同的求积公式！

2）代数精度：则称该求积公式的代数精度为m。这里结论：一个求积公式的代数精度越大，则该求积公式越好。

确定代数精度的方法：分别代入余项公式：则对应的代数精度为m-1。

例6.1确定求积公式的代数精度。解：故本题求积公式代数精度为1。

例6.2确定求积公式的参数A，B，C，使它具有尽可能高的代数精度，并指出相应的代数精度。解：要先求出具体的求积公式，再判断代数精度。

故所求的求积公式为为确定其代数精度，再取故所求的求积公式具有二阶代数精度。

§6.3插值型求积公式第六章数值积分与微分

借助多项式插值函数来构造的求积公式称为插值型求积公式。

一般选用不同的插值公式就可以得到不同的插值型求积公式。1）基本思想利用被积函数的插值函数代替被积函数来构造求积公式。

2）构造原理考虑函数y=f(x)在n个求积节点上的函数值构造的一个n-1次Lagrange插值多项式

若舍去等式右端的积分项，得到

插值型求积公式：插值型求积公式的求积余项插值型求积公式的代数精度至少为n-1。

例1.已知求这3点为求积节点在[0,1]上的插值型求积公式。解：法1：3个点的3个插值多项式基函数为

代入具体数值计算有:

故所求的插值型求积公式为法2：过这3个点的2次插值多项式为故所求的插值型求积公式为

3）分析（1）结论定理1：插值型求积公式的代数精度至少为n-1。推论：插值型求积公式的求积系数之和为b-a，即：

（2）Newton-Cotes求积公式可以得到[a,b]上的n个等距节点Cotes系数

通常称为n点的Newton-Cotes公式。也称为等距节点求积公式。可以事先算出cotes系数，如n=2的两个系数为所以2点的Newton-Cotes公式为

常用的Newton-Cotes公式梯形公式Simpson公式或抛物线公式

（3）n点Newton-Cotes公式的代数精度定理n点Newton-Cotes公式的代数精度至少为n-1,节点个数n为奇数时，对应的Newton-Cotes求积公式的代数精度至少为n。

(4)梯形公式与Simpson公式的余项广义积分中值定理I、梯形公式余项

抛物线公式的余项（4）Newton-Cotes公式的数值稳定性则有n8时，不稳定。

§6.4Gauss求积公式第六章数值积分与微分

1）基本思想在插值型求积公式中利用求积节点的选取构造具有最高代数精度的求积公式。问题：n点的插值型求积公式的代数精度至少是n-1,那么是否还能提高其代数精度呢？若能，其代数精度最大能是多少？

2）构造原理考虑一般的求积公式这里ρ(x)是权函数。权函数ρ(x)要非负权函数ρ(x)=1时就是前面的求积公式。定理:求积公式(*)的代数精度最大为2n-1。

证明的任意一组求积节点。代入求积公式，有余项由定义得求积公式的代数精度不大于2n-1.?

为证明本定理，设f(x)是任意一个2n-1次多项式，由多项式除法有式中q(x),r(x)都是次数小于n的多项式因为求积公式(*)是具有n个节点的插值型求积公式，故其代数精度至少为n-1，故有

可以用q(x)的这种任意性，选择求积节点正交多项式理论可知上面的节点是存在唯一的，它们就是在[a,b]上关于权的n次正交多项式的零点。至此得出定理结论。

结论与概念:定理：n点插值型求积公式的代数精度至少是n-1，至多为2n-1。定义:若n点的求积公式具有2n-1次代数精度，