3数值积分资料.pptVIP

第二章 数值积分 /*Numerical Integration*/ 近似计算 但是在许多实际问题经常遇到下列情况： （1）原函数存在但不能用初等函数表示； （2）原函数可以用初等函数表示，但结构复杂； （3）被积函数没有表达式，仅仅是一张数表。 ?几何意义： 曲边梯形的面积 - * f(x)没有解析表达式，只有数表形式 e.g. f(x)有表达式，但原函数不是初等函数 e.g.　　　　　，　　　　　　　　 它们的原函数都不是初等函数. x 1 2 3 4 5 f(x) 4 4.5 6 8 8.5 * 取左端点矩形近似 ?数值积分的思想： 分割、近似、求和 取右端点矩形近似 复化型求积公式 - * ?数值积分公式的一般形式： 其中 求积节点 求积系数 仅与求积节点有关 求积公式的截断误差或余项： * 代数精度的判别方法 ?求积公式的代数精度（/*Algebraic Precision */） 如果求积公式 对一切不高于m次的多项式都恒成立，而对于某个m+1次多项式不能精确成立，则称该求积公式具有m次代数精度。 求积公式 具有次m代数精度的充要条件是 为 时求积公式精确成立，而 为 时求积公式不能成为等式。 * 求积系数的特征： ?求积公式的收敛性和稳定性 若 则称求积公式（*）是收敛的。 设 有舍入误差 ，实际计算的求积公式为： * 两者的误差为 其中 求积系数全为正时，公式是稳定的 * §2 Newton—Cotes公式 一、插值型求积公式/*Integration Formula of Interpolation Type*/ 思想 用被积函数 在区间 上的插值多项式近似代替计算 作n次Lagrange插值多项式: 设已知函数 在节点 上的函数值 * 其中 插值型求积公式： 余项 * 形如 的求积公式至少 有n次代数精度的充要条件是它是插值型求积公式。 * 二、 Newton—Cotes求积公式 Newton—Cotes公式是插值型求积公式的特殊形式： 求积节点 取等距分布： * 步长 Cotes系数 Newton—Cotes公式： * n=1时的求积公式 梯形公式/*Trapezoidal Formula */ 1次代数精度 用梯形面积近似 - * n=2时的求积公式 3次代数精度 Simpson公式 用抛物形面积近似 - * n=4时的求积公式 Cotes公式 5次代数精度 近似等于曲边梯形的面积 - * 例1：分别利用梯形公式、 Simpson公式、 Cotes公式 计算积分 的近似值。 解： * 注： ? Newton—Cotes公式中 时不能使用； ? 当Newton—Cotes公式中 时不能满 足高精度要求。 对于Newton—Cotes求积公式 当n为奇数时至少具有n次代数精度；当n为偶数时至少具有n+1次代数精度。 * 第二积分中值定理 三、 前述三种求积公式的余项 ?梯形公式 设 连续 * ? Simpson公式 构造次数不超过3次的多项式 ，满足： 其中 设 连续 * 第二积分中值定理 * （补充：Newton—Cotes求积公式的误差估计） （1）当n为偶数时，如果 ，则 其中 （2）当n为奇数时，如果 ，则 其中 * ? Cotes公式 设 连续 Newton—Cotes求积方法的缺陷： 从余项公式可以看出，要提高求积公式的代数精 度，必须增加节点个数，而节点个数的增加，会导致 （1）插值多项式出现Runge现象； （2）Newton—Cotes数值稳定性不能保证。（n7） *