[2018年最新整理]71-72Newton-Cotes求积公式.pptVIP

华长生制作 传统方法的困境 数值积分的基本思想 数值积分的一般形式 代数精度问题 数值积分的基本思想 The ideas of numerical integration 数值积分的一般形式The General formula 1． Newton-Cotes公式 2．常用的NC公式 3. Newton-Cotes公式的稳定性 常用的NC公式 上式称为Simpson求积公式，也称三点公式或抛物线公式 记为 Simpson公式的余项为 Simpson公式具有3次代数精度 (3). Cotes公式及其余项 Cotes系数为 求积公式为 上式称为Cotes求积公式，也称五点公式 记为 Cotes公式的余项为 Cotes公式具有5次代数精度 一般的有 因此，N-C积分，对偶数有n+1 阶代数精度，而奇数为n 阶代数精度 表中 观察这些公式的代数精度次数，自然会得出结论： 梯形规则简单，有1次代数精度； 再增加一个节点，就是具有3次代数精度的Simpson规则； 而Simpson3-8规则又增加一个节点，精度却没有提高。 所以，人们一般只用到前两个方法。 Cotes系数的性质: 三、Newton-Cotes公式的稳定性(舍入误差) 考察Cotes系数 因此用Newton-Cotes公式计算积分的舍入误差主要由 其值可以精确给定 记 而理论值为 * * 第七章 微积分的数值计算方法 Numerical differentiation And Numerical integration 求函数 f(x) 在区间 [a,b] 上的定积分 是微积分学中的基本问题。 返回章 7.1 基本概念 Basic concepts 对于积分 但是在工程技术和科学研究中,常会见到以下现象: 传统方法的困境 Difficulties The definite integral of a function has no explicit antiderivative or antiderivative is not easy to obtain. So the Newton-Leibniz principle is not suitable to be used. 以上这些现象, Newton-Leibniz很难发挥作用! 只能建立积分的近似计算方法-------- 数值积分正是为解决这样的困难而提出来的， 不仅如此，数值积分也是微分方程数值解法的工具之一。 数值积分----是计算定积分的具有一定精度的近似值的各种计算方法。 从几何上看，就是计算曲边梯形面积的近似值。 由积分中值定理知，在积分区间 内存在一点ξ， 成立 Dxi=xi-xi?1 所以 问题在于点 的具体位置一般是不知道的，因而难以 准确算出 的值. 这样，只要对平均高度 提供一种算法，相应地便 获得一种数值求积方法. 最简单的办法，是用许多小矩形之和近似曲边梯形 的面积，如图7-0所示，这就是----矩形公式： The goal is to approximate the definite integral of f(x) by evaluating f(x) at a finite number of nodes. 图7-0 矩形规则 y x a=x0 x1 x2 xi xi+1 xn-1 xn =b f0 f1 f2 fi fi+1 fn-1 fn f(x) (1) Rectangular rule: 图7-1 梯形规则 x a=x0 x1 x2 xi xi+1 xn-1 xn =b y f0 f1 f2 fi fi+1 fn-1 fn f(x) 如果改用许多小梯形之和近似曲边梯形的面积，如图7-1，就会更精确些，这就是----梯形公式。 (2) Trapezoidal rule: 数值积分的一般形式是： 其中， fi ----是函数f(x)在节点 xi 上的函数值，它可能以列表形式给出，也可以是由函数的解析式计算出的函数值； Ai ----称为节点 xi 上的权系数。 正是由于权系数的构造方法不同，从而决定了数值积分的不同方法。 (3) 最常用的一种方法是利用插值多项式来构造数值求积公式, 具体步骤如下: 不同的 插值方法 有不同的 基函数 插值型求积公式 The quadrature by interpolating polynomial (1)为数值求积公式. Ak为求积系数,且仅与求积节(结)点xk有关. 因此定义代数精度的概念: The degree