左矩形求积公式

内容提纲（Outline)求积公式的代数精度插值型求积公式

复化求积法第一页第二页，共23页。

为什么要数值积分？在微积分里，按Newton-Leibniz公式求定积分要求被积函数f(x)☞有解析表达式；☞

f(x)的原函数F(x)为初等函数．Whydowedonumericalintegral?第二页第三页，共23页。

问题☎f(x)没有解析表达式，只有数表形式e.g.☎

f(x)有表达式，但原函数不是初等函数e.g.，它们的原函数都不是初等函数.x12345f(x)44.5688.5第三页第四页，共23页。

求定积分就得通过近似计算－数值积分求得积分近似值基本思想是对被积函数进行近似，给出数值积分，同时考虑近似精度。下面首先给出代数精确度的概念第四页第五页，共23页。7.1代数精确度本章讨论的是形如的定积分的数值计算，其中为权函数，要满足5.4节中所提的条件.

第五页第六页，共23页。一般把积分区间n个点{xk}上的函数值f(xk)加权Ak的和作为积分I(f)的近似,

即或记

(2)

第六页第七页，共23页。上式中xk，Ak分别称为求积节点、求积系数.求积系数与被积函数f(x)无关,而与求积节点、求积区间、权函数有关．称公式(2)为n点求积公式，有时也称为一个n点求积公式，为求积公式的误差．用此公式)求积分近似值的计算称为数值积分或数值微分．第七页第八页，共23页。

构造或确定一个求积公式，要讨论解决的问题有

(i)确定求积系数Ak和求积节点n；

(ii)求积公式的误差估计和收敛性．用什么标准来判定两个节点数相同的求积公式的“好”与“差”呢？通常用“代数精确度”的高低作为求积公式“好”与“差”的一个标准．在后面的讨论中我们将看到，节点相同的求积公式，代数精确度越高，求出的积分近似值精确度一般越好．下面给出代数精确度的定义．第八页第九页，共23页。定义１若对任意的，求积公式(2)的误差都满足，则称该求积公式具有n次代数精确度．验证一个求积公式所具有的代数精确度用定义１是极不方便的，为此给出另一个定义．第九页第十页，共23页。定义2若对函数，求积公式(2)精确成立，即而，则称其具有n次代数精确度．因为函数组是的一组基函数，所以两个定义是等价的，但在具体应用时，定义2比定义1要方便的多．第十页第十一页，共23页。例１验证求积公式

具有3次代数精确度．解：当而

有第十一页第十二页，共23页。（1）当

（2）当（3）当第十二页第十三页，共23页。（1）当

故求积公式具有三次代数精确度．第十三页第十四页，共23页。7.2插值型求积公式这一节所讨论的求积公式，都是用在区间[a,b]上对被积函数f(x)作插值所得插值多项式Pn(x)代替被积函数f(x)导出的公式．这一类求积公式的求积节点

xk，就是对f(x)作插值时的插值节点，所以这类求积公式称为插值型求积公式．为简便起见，这节讨论节点分布为等距并且权函数时的插值型求积公式的构造等问题．

第十四页第十五页，共23页。7.2.1Newton-Cotes求积公式一、公式的推导设将积分区间[a,b]n等分，求积节点为，那么，令x=a+th，则t=(x-a)/h，且由可知．由Lagrange插值基函数有而，所以第十五页第十六页，共23页。将n次Lagrange插值多项式Ln(x)代替被积函数f(x)得

记称为Cotes求积系数．它与(3)式中的求积系数Ak相差一个常数b-a即第十六页第十七页，共23页。把Ak代入到(3)式中，得到Newton-Cotes求积公式．例如当n=4,5时，Newton-Cotes公式分别为

n=0,1,2三种情形，在讨论(3)式中的余项R(1,f)后再详细讨论．第十七页第十八页，共23页。二、误差估计求积公式(3)计算出的积分I(f)的近似值In+1(f)的误差多大？若被积函数，记，对n次Lagrange插值余项求积，可得n+1个节点的Newton-Cotes求积公式的误差估计式为

(5)第十八页第十九页，共23页。验证求积公式(3)的代数精确度，不用误差估计的(4)式，而用直接对插值余项求积的形式，即(5)由(5)式，显而易见，当时，因可知，R(1,f)=0，所以我们所n+1点的求积公式(3)至少具有n次的代数精确度．进一步可以证明，当n为偶数时，求积公式(3)的代数精确度可以达到n+1次．第十九页第二十页，共23页。三、几种常见的Newton-Cotes求积公式对n=0,1,2,按公式(3)可以得出下面三种常见的Newton-Cotes求积公式.1.n=0时的矩形求积公式分别以积分区间[a,b]的左、右端点和区间中点，即x=a，b，(a+b)/2为求积节点得到：左矩形求积公式：右矩形求积公式：中矩形求积公式：三个求积公式的误差估计，可将函数f(x)分别在处展开到含f(x)的一阶导数的Taylor公式在区间[a,b]上积分推得．第二十页第二十一页，共23页。2.n=1时的梯形求积公式按Cotes系数公式计算得故求积系数A0,A1为 ，梯形求积公式为记(6)式的几何意义如图7-2所示（见p327）容易验证公式(6)的代数精确度的次数为