求解半正定病态积分方程的多尺度快速Lavrentiev迭代算法.pdfVIP

2012正 赣 南 师范 学 院学 报 No．6 第六期 JournalofGannanNormalUniversity Dec．2Ol2 求解半正定病态积分方程的 多尺度快速 Lavrentiev迭代算法 罗兴钧，谢长发娣，陈维君，李繁春 (赣南师范学院数学与计算机科学学院，江西 赣州 341000) 摘 要：基于多尺度投影方法，构造了求解半正定病态积分方程的截断快速Lawentiev迭代算法，并给出了新 的后验参数选择办法．与传统全投影方法相比，减少了内积计算个数，保持了最优收敛率．最后，算例验证了算法的 有效性． 关键词：病态积分方程；迭代Lavrentiev正则化；截断投影方法；后验参数选择办法 中图分类号：O17 文献标志码：A 文章编号：1004—8332(2012)06—0001—05 1 引言 多数数学物理反问题可以归结为第一类Fredholm型病态积分方程的求解．吉洪若夫正则化方法是求解 病态积分方程的主要方法 J．我们感兴趣的是在有限维空间的近似求解．文献 [2]提出过有限维空间全投影 方法，缺陷是计算量大．文献 [3]采用截断投影方法求解第一类病态半正定积分方程，提出了截断投影方法， 但只给出了先验误差估计，没有给出后验参数选择方法．文献[4]采用截断投影方法求解第一类病态半正定 积分方程，提出了新的正则化参数后验选择准则，证明了近似解达到最优，与传统投影方法相比得到了相同 的最优收敛率，但减少了计算量．然而此算法有缺陷，一是文献[4]中没有证明参数的存在性，二是在精确解 满足某种光滑性的条件下，收敛率是饱和的．要提高收敛率，必须采用其它正则化方法．于是本文将迭代 Lavrentiev正则化方法与截断投影方法相结合，并提出了新的后验参数选择办法，构造了求解半正定病态积 分方程的一种 自适应算法，在精确解满足某种光滑性条件及正则化参数为最佳选择时，与传统全投影方法相 比，既减少了计算量，又保持了收敛率． 假设E是d维欧氏空间 (d 1)的紧集，是Hilbert空间 (E)，用 (·，·)与 lI．I1分别表示Hilbert 空间 的内积及范数．假设 — 是线性半正定紧算子，即对任意 EX，(Ax，) 0．定义Fredholm积分 算子A (Ax)(t)：：I(s，￡)(s)ds=Y(￡)， EE (1．1) 古 这里，KEc(E×E)．以下，考虑第一类积分方程的求解 Ax=Y， (1．2) 即yEX是已知，而 EX是待求的．显然，算子A：x—x是紧算子．因此，方程(1．2)是病态问题．容易得到算 子A的共扼算子A 为 (A )(￡)：=J(t，5)(s)ds=y()， EE． 士 以下，假设YER(A)且(1．2)有一个最小范数解 ，满足 = A ·，lI∞IIsP，JD， 0． (1．3) 在实际问题中，知道的仅是Y的近似值 ，且满足 lIY—Y l【 (1．4) 由于方程(1．2)是不适定的，要使用正则化策略．本文采用迭代 Lavrentiev正则化方法求 的近似解 + 收稿 日期：2012一o6—15 网络出版 日期：2012—12一I1 基金项 目：国家 自然科学基金资助项 目；江西省自然科学基金资助项 目(20114BAB201014)；江西省教育厅科学技术研究资 助项 目(GJJ10586) 作者简介：罗兴钧，男，赣南师范学院数学与计算机科学学院教授，研究方