4.1-4.2Newton-Cotes求积公式.ppt

华长生制作 传统方法的困境 数值积分的基本思想 数值积分的一般形式 代数精度问题 数值积分的基本思想 数值积分的一般形式 1． Newton-Cotes公式 2．常用的NC公式 3. Newton-Cotes公式的稳定性 4. 复化求积法 常用的NC公式 * * 第四章 微积分的数值计算方法 求函数 f(x) 在区间 [a,b] 上的定积分 是微积分学中的基本问题。 4.1 基本概念 对于积分 但是在工程技术和科学研究中,常会见到以下现象: 传统方法的困境 以上这些现象, Newton-Leibniz很难发挥作用! 只能建立积分的近似计算方法-------- 数值积分正是为解决这样的困难而提出来的， 不仅如此，数值积分也是微分方程数值解法的工具之一。 数值积分----是计算定积分的具有一定精度的近似值的各种计算方法。 从几何上看，就是计算曲边梯形面积的近似值。 最简单的办法，是用许多小矩形之和近似曲边梯形 的面积，如图4-0所示，这就是----矩形公式： 图4-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) 图4-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) 如果改用许多小梯形之和近似曲边梯形的面积，如图4-1，就会更精确些，这就是----梯形公式。 (2) 数值积分的一般形式是： 其中， fi ----是函数f(x)在节点 xi 上的函数值，它可能以列表形式给出，也可以是由函数的解析式计算出的函数值； Ai ----称为节点 xi 上的权系数。 正是由于权系数的构造方法不同，从而决定了数值积分的不同方法。 (3) 最常用的一种方法是利用插值多项式来构造数值求积公式, 具体步骤如下: 插值型求积公式 (1)为数值求积公式. Ak为求积系数, 且仅与求积节(结)点xk有关. 因此定义代数精度的概念: 定义1. 若求积公式 则称该求积公式具有m次的代数精度. 代数精度也称 代数精确度 例1. 试确定下面积分公式中的参数使其代数精确度尽量高. 解: 因此 所以该积分公式具有3次代数精确度 4.2 Newton-Cotes求积公式 1、Newton-Cotes数值求积公式 Newton-Cotes公式是指等距节点下使用Lagrange插值 多项式建立的数值求积公式 各节点为 其中 而 因此对于定积分 有 令 即有 n阶Newton-Cotes求积公式 Newton-Cotes公式的余项(误差) 注意是等距节点 所以Newton-Cotes公式化为 2、低阶Newton-Cotes公式及其余项 在Newton-Cotes公式中,n=1,2,4时的公式是最常用也 最重要三个公式,称为低阶公式 (1). 梯形(trapezia)公式及其余项 Cotes系数为 求积公式为 上式称为梯形求积公式,也称两点公式，记为 梯形公式的余项为 第二积分 中值定理 梯形(trapezia)公式具有1次代数精度 故 (2). Simpson公式及其余项 Cotes系数为 求积公式为 上式称为Simpson求积公式，也称三点公式或抛物线公式 记为 Simpson公式的余项为 Simpson公式具有3次代数精度 (3). Cotes公式及其余项 Cotes系数为 求积公式为 上式称为Cotes求积公式，也称五点公式 记为 Cotes公式的余项为 Cotes公式具有5次代数精度 表中 观察这些公式的代数精度次数，自然会得出结论： 梯形规则简单，有1次代数精度； 再增加一个节点，就是具有3次代数精度的Simpson规则； 而Simpson3-8规则又增加一个节点，精度却没有提高。 所以，人们一般只用到前两个方法。 * * * *