大连理工大学 计算机科学计算 第五章1.pdfVIP

第 5 章 插值函数的应用 5.1 基于插值公式的数值积分 5.1.1 数值求积公式及其代数精度 5.1.2 复化求积公式 5.1.3 数值微分公式 5.2 Gauss型求积公式 5.2.1 基于Hermite插值的Gauss型求积公式 5.2.2 常见的Gauss型求积公式和数值稳定性 5.3 外推加速原理和Romber算法 5.3.1 逐次分半算法 5.2.2 外推加速公式和Romber算法 5.1.1 数值求积公式及其代数精度 由 Newton-Leibniz公式，连续函数 f (x)在 [a,b] 上的定积分 b f x ( ) dx ( )F a ( ) ∫a F b − 其中 F(x) 是f (x) 的原函数。但是大多数实际问题，N-L公式已经 无能为力。 常常遇到的困难是： F (x)不能用初等函数表示，即f (x)找不到的原函数； 1 sin x 2 1 f (x ) , ( ) , −x f (x ) , ln x f x x f (x ) e , 1−k 2 sin2 x f (x)没有解析表达式，用表格方式给出时； 大多数的无穷积分，除特殊的无穷积分外。 虽然找到f (x) 的原函数，但是它比被积函数复杂的多 2 x dx 1 2 x+ 1 + x 4 + 1 arcsin x 2 +c 4 4 ln 1−x (1+x ) 2 2 ∫ 4 2 1−x 2 2 1+x 上述的积分就只能利用数值积分公式进行近似计算。 设f (x) 是定义在[a,b]上的可积函数，考虑带权积分 b I f ( ) ∫a (ρ)xf (x ) dx （5-1 ） 其中权函数 r(x) 在[a, b]上非负可积，且至多有有限个零点。 本节只讨论r(x) ≡1的情形。 所谓数值求积就是用 n （5-2 ） f ( x ) I (f n ) ∑Ak k k 0 近似计算 I(f ) 的值。 求 公式（5-2 ）称为数值求积公式， 数值求积公式 积 系 n 数 I (f ) A n ∑f ( x ) k k k 0 A ( 0, 1, , )k