数值分析第8章-1.pptVIP

* 8.1 求 积 公 式 8.1.1 求积公式 结束 第八章 数值积分与数值微分 * 对定义在区间［a,b］上的定积分 以上公式多称为牛顿-莱布尼兹公式，F(x)为f(x)的原函数.但有时原函数不能用初等函数表示，有时原函数又十分复杂，难于求出或计算.如被积函数为: 等函数的积分都无法解决，当被积函数为一组数据时，更是无能为力. 为解决定积分的近似计算,从定积分的定义： 结束 * 这样就避开了求原函数的运算.(8.1)式就叫做求积公式，Ak(k=0,1,…,n)与函数f(x)无关，叫做求积系数,显然要确定一个求积公式，要确定求积结点xk和求积系数Ak，或者说不同的求积结点和求积系数将确定不同的求积公式. 8.1.2 求积公式的余项和代数精度 一般情况下，(8.1)两端并不相等.我们称: (8.2)为求积公式(8.1) 的余项，或截断误差. 结束 * 为考查一个求积公式的误差，通常用代数精度来表示，如果一个求积公式对于不超过m次的多项式都能够精确成立(R[f]≡0)，而对m+1次以上的多项式不能精确成立，则称该求积公式的代数精度为m. 例如求积公式： 验证当f(x)=xm，m=0,1,2,3,4时，是否有R[xm]=0 结束 * 所以以上求积公式的代数精度为3. 把f(x)在a处作Taylor展开: f(x)= f(a)+ f ’(ξ)(x-a)， ξ在x,a之间,两端积分： 注意到右端第二项积分，设f ′(x)在［a,b］上连续，而x-a在 ［a,b］上不变号(非负)，据积分中值定理有： 任何一个求积公式的代数精度至少为零即取f(x)=1时公式应精确成立，这是求积系数应满足的起码条件，可以用它检验一个求积公式的系数的正确性. 8.1.3 矩形求积公式 结束 * 于是有左矩形公式(如图) 同理,f(x)在b点展开,可得右矩形公式(如图) : f(x)在中点(a+b)/2展开,可得中矩形公式(如图) : 不难验证，(8.3)和(8.4)具有零次代数精度，(8.5)具有一次代数精度. 结束 * 8.1.4 内插求积公式 由插值可知，对任一函数f(x)(包括表格形式的函数)可用一n次多项式对其插值，即 当Pn (x)为拉格朗日插值多项式时，即 结束 * 其中: 通常将公式(8.6)叫做内插求积公式. §8.2 牛顿-柯特斯公式 为便于上机计算，通常在内插求积公式中取等距节点，即将积分区间[a,b]n等分，即令h=(b-a)/n，且记x0=a,xn=b，则节点为xk=x0+kh(k=0,1,…n)，作变换:t=(x-x0)/h，代入求积系数公式： 结束 * 这种由等距节点的内插求积公式通常叫做牛顿-柯特斯公式,下面介绍几个常用的公式: 取a=x0,b=x1,(即n=1)，代入(8.9)式得 结束 8.2.1 梯形公式 所以梯形公式为 * 从图8.1看到，这是用梯形面积近似代替曲边梯形的面积，对梯形公式的误差估计有如下定理： 定理8.1 设f(x)为二阶连续可微函数，则梯形求积公式的余项为 (证明) 结束 * 其中h=b-a，记成上面形式是为以后复化求积公式余项的一致性.由余项公式立刻可以看出梯形公式的代数精度为1. 例1 利用梯形公式计算 解: 结束 8.2.2 抛物形（辛卜生）公式 * 取a=x0,(a+b)/2=x1,b=x2,(即n=2)，代入(8.9)式得 所以抛物形公式为 其中h=(b-a)/2,上式也可写成: 结束 抛物形公式通常也称为辛普生公式，从图8.2看到，抛物形公式是用抛物线围成的曲边梯形近似代替f(x)围成的曲边梯形. * 定理8.2 设f(x)∈C4[a,b]，则辛普生公式的误差估计为: 直接可以验证抛物形公式代数精度为3(对f(x)为三次以下多项式精确成立). 例2 利用抛物形公式计算 解: * * * * *