数值分析第7章数值积分与数值微分.pdf

第七章 数值积分与数值微分 §1 数值积分基本概念 §2 插值型求积公式 §3 复化求积算法 §4 Romberg求积法 §5 Gauss型求积公式 §6 数值微分 2009-09-26 1 §7.1 数值积分基本概念 由积分学基本定理Newton-Leibenize公式有： b ∫ f (x)dx = F (b) − F (a) a 但在应用中常会碰到如下情况： ①f (x) 的原函数无法用初等函数给出 ②f (x)用表格形式给出 ③虽然f (x) 的原函数能用初等函数表示, 但表达式过于复杂。 2009-09-26 2 一.求积公式 b n I = ∫a f ( x ) dx ≈ ∑ A k f ( x k ) = I n k = 0 b n I = f (x)dx = ∑A f (x ) + R[f ] = I + R[f ] ∫ k k n a k =0 x (0k,1=,2, Ln) 其中R [f ]称为求积公式的余项， k A (0k,1=,2 , L n ) 称为求积节点，k 称为求积系数。 A k 仅与求积节点 x k 的选取有关，而不依赖与被积 函数f (x) 的具体形式。 2009-09-26 3 二.求积公式的代数精确度 衡量一个求积公式好坏的标准。 b n 定义：如果求积公式 f (x)dx ≈ ∑A f (x ) ∫ k k a k =0 对于一切不高于m次的代数多项式准确成立， 而对于某个m+1次多项式并不准确成立， 则称上述求积公式具有m 次代数精确度，或称 为具有m次代数精度。 2009-09-26 4 求积公式的代数精确度 （续） 如果要构造具有m次代数精度的求积公式， 只要令它对于f (x ) =1, x , x 2 , L, x m 都能准确成立即可。 b