代数精度插值求积及复化公式参考资料.pptVIP

? ? ? ? b a a F b F dx x f I ) ( ) ( ) ( 但是，在工程技术领域，在实际使用上述求积分方法时，往 往会遇到下面情况： 1. 函数 f ( x ) 没有具体的解析表达式，只有一些由实验测试 数据形成的表格或 图形。 关于定积分的计算，我们知道，只要求出 f ( x ) 的一个原 函数 F ( x ) ，就可以利用牛顿 — 莱布尼慈（ Newton-Leibniz ）公 式出定积分值： 3. f ( x ) 的结构复杂，求原函数困难，即不定积分难求。 等 3 2 1 , ln 1 , sin , , sin ) ( 2 x x x e x x x f x ? ? ? 2. f ( x ) 的原函数无法用初等函数表示出来，如： 由于以上种种原因，因此有必要研究积分的数值计算方法，进 而建立起上机计算定积分的算法。此外，数值积分也是研究微 分方程和积分方程的数值解法的基础。 数值积分 1.1 构造数值求积公式的基本思想 定积分 I =∫ a b f ( x ) dx 在几何上为 x=a, x=b , y= 0 和 y=f ( x ) 所围成的曲 边梯形的面积。定积分计算之所以困难，是不规则图形的面积。由 积分中值定理，对连续函数 f ( x ) ，在区间 [ a , b ] 内至少存在一点 ? ， 使： ) ( ) ( ) ( ? f a b dx x f I b a ? ? ? ? 也就是说 , 曲边梯形的面积 I 恰好等于 底为 b - a , 高为 f ( ? ) 的规则图形 — 矩形的面 积 ( 图 7-1), f ( ? ) 为曲边梯形的平均高度 , 然 而点 ? 的具体位置一般是不知道的 , 因此难 以准确地求出 f ( ? ) 的值。但是 , 由此可以得 到这样的启发 , 只要能对平均高度 f ( ? ) 提供 一种近似算法 , 便可以相应地得到一种数 值求积公式。 图 7-1 a b ξ ) ( x f y ? ) ( ? f 如用两端点的函数值 f ( a ) 与 f ( b ) 取算术平均值作为平均高度 f ( ? ) 的近似值 , 这样可导出求积公式： 第七章 数值积分与微分 7-3 更一般地在区间 [ a , b ] 上适当选取某些点 x k ( k =0,1,…, n ), 然后 用 f ( x k ) 的加权平均值近似地表示 f ( ? ), 这样得到一般的求积公式： 1) - (7 ) ( ) ( 0 n n k k k b a I x f A dx x f I ? ? ? ? ? ? 其中 , 点 x k 称为求积节点 , 系数 A k 称为求积系数， A k 仅仅与节 点 x k 的选取有关 , 而不依赖于被积函数 f ( x ) 的具体形式。 ( ) ( ( ) ( )) 2 , ( ) ( ) 2 2 b a b a b a I f x dx f a f b a b a b I f x dx b a f ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 梯形公式 取 中矩形公式 另一方面定积分的定义， 0 0 0 ( ) lim ( ) k k n n b k k a Max x k I f x dx f x x ? ? ? ? ? ? ? ? ? ? 其中 ? x k 是 [ a , b ] 的每一个分割小区间的长度 , 它与 f ( x ) 无关 , 去掉 极限 , 由此得到近似计算公式： ? ? ? ? ? ? ? ? ? n k k k n k k k b a x f A x x f dx x f I 0 0 ) ( ) ( ) ( 因此，式（ 7-1 ）可作为一般的求积公式 , 其特点是将积分问 题归结为函数值的计算 , 从而避开了使用牛顿一莱布尼慈公式需 要求原函数的困难 , 适合于函数给出时计算积分 , 也非常便于设计 算法 , 便于上机计算。 求积公式（ 7-1 ）的截断误差为： 0 ( ) ( ) ( ) n b n n k k a k R f R I I f x dx A f x ? ? ? ? ? ? ? ? R n 也称为积分余项 . 1.2 代数精度 定义 1 如果某个求积公式对所有次数不大于 m 的多项式都精确成 立，而至少对一个 m +1 次多项式不精确成，则称该公式具 有 m 次代数精度。 一般来说，代数精度越高，求积公式越好。为了便于应用， 由定义