第四章数值积分和数值微分Chapter4NumericalIntegration1.pptVIP

第四章 数值积分与数值微分 /*Chapter4 Numerical Integration and numerical differentiation */ * 近似计算 机械公式 矩形公式 梯形公式 4.1. 数值积分概论 1. 数值积分的基本思想 解：逐次检查公式是否精确成立 代入 P0 = 1： = 代入 P1 = x ： = 代入 P2 = x2 ： ? 代数精度 = 1 例 求代数精确度 梯形公式 如果某个求积公式对于次数不超过m的多项式均能准确成立，但对于m+1次多项式就不准确成立，则称该求积公式具有m次代数精度。 2. 代数精度的概念 思路 利用插值多项式 则积分易算。 ? 在[a, b]上取 a ? x0 x1 … xn ? b，做 f 的 n 次插值多项式 ，即得到 Ak 由 决定， 与 无关。 节点 f (x) 插值型积分公式 /*interpolatory quadrature*/ 误差 3. 插值型的求积公式 定理1：形如 的求积公式至少有 n 次代数精度 ? 该公式为插值型（即： ） ? 当节点等距分布时： 令 Cotes系数 注：Cotes 系数仅取决于 n 和 i，可查表得到。与 f (x) 及区间[a, b]均无关。 4.2 Newton-Cotes牛顿-柯特斯公式 牛顿-柯特斯公式 Cotes系数表 可查下表 n = 1: n = 2: Simpson公式 代数精度 = 3 定理 当阶n为偶数时，newton-cotes公式至少具有n+1次代数精确度. n = 4: Cotes公式 求积公式的余项 （2） (-1)n+k=(-1)n-k 特点： Cotes系数特点： 1.5 Newton—Cotes公式的数值稳定性和收敛性 1、数值稳定性 若某个求积公式的舍入 设 的近似值为 ， 公式是稳定的；否则，若影响较大，则称为不稳定的。 误差，即f(xk)的误差对数值积分的结果影响较小，则称该数值求积 以下推导说明 Newton—Cotes公式的数值稳定性 2、收敛性 4.3 复合求积公式 高次插值不稳定，故采用分段低次插值 ? 分段低次合成的 Newton-Cotes 复合求积公式。 ? 复合梯形公式： 在每个 上用梯形公式： ? 复合 Simpson 公式： = Sn 考察复合梯形公式，区间二分前后： 4.4 龙贝格求积公式/*Romberg Integration*/ 理查森（Richardson）外推法 理查森外推加速方法 *