CH4数值积分与数值微分—4.1Newton-Cotes公式.pptVIP

华长生制作 第四章 数值积分与 数值微分 张红梅 自动化学院 2010年3月 4.1 Newton-Cotes 公式 用插值多项式 P(x) 代替被积函数或者被微分函数f(x) , 从而导出计算 f(x) 的定积分和微分近似值的公式。 4.1 Newton-Cotes公式 4.2 复合求积法 4.3 Romberg算法 4.4* Gauss求积法 4.5 数值微分 利用等距节点的Lagrange插值 多项式建立插值型求积公式 将求积区间等分为多个子区间,然后在每个子区间用低阶求积公式 利用复合梯形公式,并采用加速技术 用插值多项式代替被求微分的函数 本章要点 对于任意函数的定积分,是否都可以利用Newton-Leibniz 公式求得？ 如果已知 f(x) 的原函数 F(x), 则由 Newton-Leibniz 公式： 可求得式(1)的积分值. (1) f(x)的解析式不存在,只给出了 f(x) 的一些数值； ? f(x)的原函数 F(x) 无法求出,例如 F(x) 不是初等函数； ? f(x)的表达式结构复杂,求原函数较困难. ? 对于积分 解决方案：建立积分近似计算 —— Newton-Leibniz 公式失效的几种情况： 利用插值多项式构造数值求积公式 如何解决上述问题？ 引言 4.1 Newton-Cotes 公式 1-1 插值型求积公式及Cotes系数 使用 Newton-Cotes 公式求解数值积分问题的步骤: 则 f(x) 可表示为 ： (2) 1.将积分区间 [a, b] n 等分, 令 , 得等距节点： 等距节点下使用 Lagrange 插值多项式建立的数值求积公式. Newton-Cotes 公式 : 2.将(2)代入定积分公式 ： (4) (3) (3)式为 插值型求积公式 ? 为 求积系数 ? 为 数值积分余项 ? 利用公式(3)计算 In 的 关键:计算求积系数 Ak 3. 定义近似积分： 如何计算 Ak ？ 4. Ak 的计算 注意: 节点等距 设 由 可知 且有 代入(5)式, 得 (5) (6) 将 Ak 代入(3)式, 称为 Cotes系数. ? 求积区间等分数 注：1.只要已知节点上的函数值, 以及积分区间等分数, 即可求得 积分的近似值； 2.n 取不同的值, 得不同的求积公式。 可知, Cotes 系数与函数 f(x) 以及区间 [a, b] 无关,而只与区间等分数 n 有关。 由 节点下标 得 Newton-Cotes 求积公式: n = 1~8 的Cotes 系数 全正 有正 有负 x0 x1 Newton-Cotes 公式中,n =1,2,4 时的求积公式是最常用、最重要的三个公式。 1. 梯形(trapezoid)公式 Cotes系数为: 求积公式为: 取 即 梯形求积公式/两点公式 则 2. Simpson 公式 Cotes 系数为: 取 则 x0 x1 x2 Simpson求积公式/三点公式/抛物线公式 Simpson求积公式为: 3. Cotes 公式 Cotes系数为: 取 则 Cotes 求积公式为: Cotes 求积公式/五点公式 记为 例：用梯形公式和Simpson公式对幂函数 xm (m=1,2,3,4) 和指数函数 ex 在区间[-2,0]上积分，结果如下： 0.869 6.667 -4 2.667 -2 2 Simpson值 1.135 16 -8 4 -2 2 梯形值 0.865 6.4 -4 2.667 -2 2 精确值 ex x4 x3 x2 x 1 函数f(x) 结论：simpson 公式对更多的函数求积分是准确的。 为了使一个求积公式(近似积分)具有较好的实际计算意义,就希望它对尽可能多的被积函数都准确成立. 但对 m+1 次的多项式却不能准确成立,即只要 则称该求积公式具有 m 次的代数精度. 对不超过 m 次的代数多项式 都准确成立,即 若 的某求积公式 代数精度 第二积分 中值定理 梯形公式具有 1 次代数精度 故 梯形公式的余项为: 由余项公式可知: 1-2 低阶 Newton-Cotes 公式的余项 Simpson 公式的余项为: Simpson公式具有3 次代数精度 Cotes 公式的余项为: Cotes 公式具有5 次代数精度 考察Cotes系数: 因此,用 Newton-Cotes 公式计算积分的舍入误差主要由函数值 的