用复化梯形公式和外推法求解积分方程-同济大学数学系.pdf

用复化梯形公式和外推法求解积分方程 杨小刚 上海同济大学数学系，上海 200092 摘要：利用数值积分公式对线性积分方程进行离散，舍去误差项，解线性方程组得到积分方程的 近似数值解，然后分析误差，利用Romberg 外推技巧对算法进行加速并提高精度，并在理论上给予证明。 最后探讨了使用Gauss-Legendre 公式的好处和不利的地方，并分析了原因。 关键字：第二类Fredholm 线性积分方程；复化梯形公式；复化Simpson 公式； Romberg 外推技 巧。Gauss-Legendre 公式 0 引言 积分方程分类复杂，各类积分方程的理论的数值解法的差异很大。然而数值积分的算法 已经很成熟，这里借助数值积分的算法对第二类Fredholm 积分方程进行离散进而求出近似 解，由于数值积分的等距节点型算法可以使用 Romberg 外推技巧，这个思想可以应用到积 分方程的求解上，以提高计算精度。在同样多节点的情况下Gauss 型公式具有最高的代数精 度，可以认为采用Gauss 型公式可以得到更好的结果，本文指出在节点较少的情况下的确如 此，当节点过多，由于计算机舍入误差的存在，计算精度会下降。 1 数学模型 第二类Fredholm 线性积分方程的基本格式为： b y(t) k (t, s)y (s)ds f (t) ，a≤t≤b. (1) a 假定f (t) 和k (t, s) 充分光滑，即 f (t) ∈C m [a,b] , k (t, s) ∈C m ([a,b] ×[a,b]) , (2) 其中m ∈N ，m ≥1. 为了保证（1）有唯一解，假定k (t, s) 满足 δ ＝max | k (t, s) |1/(b-a) (3) [a ,b ]2 定理1 如果条件(2)和(3)满足，则积分方程（1）存在唯一解y (t) ∈C m [a,b]. b [1] k (t, s)y (s)ds f (t) m 证明 令T y (t) ＝ ，易见，T 是完备空间C [a,b]到自身的映射， a 由于δ ＝max | k (t, s) |1/(b-a)，于是 [a ,b ]2 b ||Ty 1 ,Ty 2 ||= max | k (t, s)(y 1 (s) y 2 (s))ds | [a ,b ] a ≤δ (b-a)||y 1 , y 2 || 由(3)知：δ (b-a) 1,T 为压缩映射，由压缩映射原理知，T 有唯一不动点，即方程（1）存在 唯一解y (t) ∈C m [a,b]，此证。 2 用数值积分公式离散求解 注意到方程（1）的积分项是一个定积分，可以用数值积分公式进行离散，即用求积公 式离散后代替该定积分。 设a=t 0 t 1 …t n =b , 使用插值型求积公式 b n x (t)dt w x (t ) +E n (x) , x(t) ∈C m [a,b] (4)  j j a j 0 其中w j (j=0,1,…,n)是求积系数，t j (j=0,1,…,n