数值积分方法课件.pptVIP

第五章数值积分方法问题提出计算但是在许多实际问题经常遇到下列情况：（1）原函数存在但不能用初等函数表示；（2）原函数可以用初等函数表示，但结构复杂；（3）被积函数没有表达式，仅仅是一张函数表。解决以上情况的积分问题，最有效的办法为数值积分法。此种方法是利用被积函数在一些离散点处的函数值，而求得满足一定代数精度要求的定积分近似值。

?定积分几何意义：曲边梯形的面积?数值积分的思想：分割、近似、求和取左端点矩形近似

?数值积分公式的一般形式：其中求积节点求积系数仅与求积节点有关求积公式的截断误差或余项：

§5.1插值型求积公式思想用被积函数在区间上的插值多项式近似代替计算设已知函数上的函数值在节点作n次Lagrange插值多项式:

其中余项

则有数值积分公式这是用插值函数代替被积函数导出的定积分近似计算公式，称为插值型数值积分公式。

一、梯形公式n=1时的求积公式

这是用线性插值函数代替被积函数导出的定积分近似计算公式，称为梯形数值积分公式。几何意义用梯形面积近似

截断误差：已知线性插值的截断误差为积分中值定理：连续、不变号

二、Simpson公式将[a,b]二等分，等分节点x=a，x=(a+b)/2，n=2时的求积公式01x=b作为积分节点，构造二次Lagrange插值多2项式L(x)：2

这是用二次插值函数代替被积函数导出的定积分近似计算公式，称为辛普森数值积分公式。几何意义：

Simpson积分公式的截断误差（定理）：积分中值定理：连续、不变号

§5.2复合求积公式复合求积法通常把积分区间等分成若干个子区间，在每个子区间上用低阶的求积公式（如梯形积分公式Simpson积分公式），对所有的子区间求和即得整个区间[a,b]上的积分公式，这种方法称为复合求积法。

5.2.1复化梯形积分将[a,b]分成若干小区间，在每个区间[x,x]上用ii+1梯形积分公式，再将这些小区间上的数值积分累加起来，就得到区间[a,b]上的数值积分。这种方法称为复化梯形积分。★计算公式将[a,b]n等分,h=x-x=(b-a)/n,x=a+ih,i=0,1,2,…,n,i+1ii

记为T(h)或T(f):复化梯形公式n复化梯形公式的几何意义

设复化梯形公式的余项由介值定理余项估计式

5.2.2复化Simpson公式：★计算公式将[a,b]2m等分,m为积分子区间数，记n=2m，n+1为节点总数，h=x-x=(b-a)/n,x=a+ih,i+1iii=0,1,2,…,n,

复化Simpson公式系数首尾为1，奇数点为4，偶数点为2复化Simpson公式的几何意义

设复化Simpson公式的余项由介值定理余项估计式

例：分别利用复化梯形公式、复化Simpson公式计算积分的近似值，要求按复化Simpson公式计算时误差不超过。解：首先来确定步长复化Simpson公式的余项：其中

本题的求法：由归纳法知

解不等式得将区间8等分，分别采用复化Simpson、梯形公式011/81/43/80.9973980.9896880.9767271/25/86/87/810.9588510.9361560.9088580.8771930.841471

复化梯形公式(n=8)复化Simpson公式(n=4)

§5.3数值积分公式的代数精度和Gauss求积公式如果求积公式对一切不高于m次的多项式都恒成立，而对于某个m+1次多项式不能精确成立，则称该求积公式具有m次代数精度。定理求积公式具有次m代数精度的充要条件是为代数精度的判别方法时求积公式精确成立，而不能成为等式。为时求积公式

例2见p73的例5.5

Gauss求积公式一、Gauss积分问题的提法积分公式的一般形式：?前述的求积公式中求积节点是取等距节点，求积系数计算方便，但代数精度要受到限制；?为了提高代数精度，需要适当选择求积节点:①当求积节点个数确定后，不管这些求积节点如何选取，求积公式的代数精度最高能达到多少？②具有最高代数精度的求积公式中求积节点如何选取？

定理的插值型求积公式的代数精度最高不超过2n+1次。形如

这样由方程组的4个方程就能求出4个未知数，得根据定理知三点插值型求积公式的代数精度为5，同理可以去验证三点高斯求积公式二、Gauss求积公式的应用积分见书上的P75.