计算方法公选课课件(第六章).pptVIP

第六章 数值积分 一、问题的提出 按照定义，定积分 是利用牛顿-莱布尼兹（Newton-Leibnitz ）公式 （6-1） 来计算的。其中F(x)是被积函数f(x)的原函数，F’(x)=f(x)。可是实际中经常遇到的一种情况是被积函数f(x)的原函数F(x)不存在（例如椭圆积分），因此无法利用（6-1）式计算 ；另一种情况是被积函数f(x)的原函数F(x)虽然存在，但很复杂，因此利用（6-1）式计算 很麻烦；还有一种情况是被积函数f(x)不是以数学公式的形式给出，而是以表格的形式给出，因此公式（6-1）就 更无法利用了。为此，我们讨论计算 近似值的数值计算方法--数值积分。 二、基本求积公式 1、梯形求积公式 用a,b两点作一次插值多项式 用 代替f(x)在区间[a,b]上作积分，得 （6-2） 从图6.1看到，这是用梯形面积近似代替曲边梯形的面积，所以式（6-2）叫做梯形求积公式。 图6.1 根据定理2.1有 其中 是依赖于x的函数，两边积分得 假定 在[a,b]上连续，而(x-a)(x-b)在[a,b]上恒小于零，利用积分中值定理，在[a,b]上存在一点 ，使 因此，梯形求积公式的误差估计为 （6-3） 2、抛物线求积公式 把区间[a,b]二等分，过a,b和 三点作二次插值多项式 用 代替f(x)在区间[a,b]上作积分，得 （6-4） 图6.2 从图6.2看到，这是用抛物线 所围成的曲边梯形面积近似代替f(x)所围成的曲边梯形的面积，所以式（6-4）叫做抛物线求积公式，有时也叫做辛浦生（Simpson）公式。 利用和梯形求积公式的误差估计相同的方法，若假定 在[a,b]上连续，则得抛物线求积公式的误差估计为 （6-5） 3、牛顿-柯特斯（ Newton-Cotes）公式 把区间[a,b]n等分，其分点为 过这n+1个节点，构造一个n次插值多项式 其中 。用 代替f(x)在区间[a,b]上作积分，得 （6-6） ? 公式（6-6）叫做牛顿-柯特斯（ Newton-Cotes）公式。其中的系数 叫做牛顿-柯特斯系数。下面的表6.1给出了n等于1至8的牛顿-柯特斯系数 n 1 2 3 4 5 6 7 8 4、例题 例1 试用梯形求积公式、抛物线求积公式和n=4的牛顿-柯特斯公式计算定积分 解：（1）梯形求积公式 （2）抛物线求积公式 （3） 的牛顿-柯特斯公式 该定积分的准确值为 三、复化求积公式 1、复化梯形求积公式 把区间[a,b]n等分，其分点为 对每个小区间[xk,xk+1]用梯形求积公式，则得 （6-7） Tn称为复化梯形求积公式，下标n表示将区间n等分。若在每个小区间上直接用梯形求积公式的误差估计（6-3），则 由于f”(x)在[a,b]上连续，利用连续函数的性质知在[a,b]中存在一点 ，使 这样就得到了复化梯形求积公式的误差估计 （6-8） 2、复化抛物线求积公式 因为抛物线求积公式用到了区间的中点，所以在构造复化抛物线求积公式时必须把区间[a,b]等分