求解第一类Fredholm积分方程的多层迭代算法.pdfVIP

2013年8月 计算数学 第35卷第3期 MATHEMATICA Aug．，2013 NUMERICASINICA V01．35，No．3 求解第一类Fredholm积分方程的 多层迭代算法术1) 李繁春 (江西应用技术职业学院基础教学部，江西赣州341000) 杨素华 罗兴钧 彭玉兵 (赣南师范学院数学与计算机科学学院，江西赣州341000) 摘 要 本文先把正则化后的第二类积分方程分解为等价的一对不含积分算子K。瓦、仅含积分算子 咒以及咒+的方程组，再用截断投影方法离散方程组，采用多层迭代算法求解截断后的等价方程 组，并给出了后验参数的选择方法，确保近似解达到最优．与传统全投影方法相比，减少了积分计 算的维数，保持了最优收敛率．最后，算例说明了算法的有效性． 关键词：第一类Fredholm积分方程；Tikhonov正则化；多层迭代算法；后验参数选择 MR(2000)主题分类：65J20，65J10 1．引 言 图像处理以及微分方程中的许多问题可以归结为第一类Fredholm型病态积分方程的求 解．吉洪若夫正则化方法是求解病态积分方程的主要方法【1】．它把病态问题转化为第二类良 态积分方程． ．从计算角度来看，直接利用Tikhonov正则化后的第二类积分方程寻找数值解，计算量巨 大。假设足是积分算子，Tikhonov正则化方法需要生成积分算子瓦+足的离散矩阵表示形式， 而积分算子．IC+．IC的维数是积分算子瓦的两倍，这极大地提高了生成系数矩阵的计算时间，而 且，当我们采用Galerkin方法去离散合成积分算子瓦+瓦时，由于维数的增加，需要花更多的 时间做计算．为解决这一问题，文[2】提出了新算法，将Tikhonov正则化积分方程分成两个等 价的方程组，避免出现积分算子瓦+咒，仅出现积分算子瓦及瓦+【2J．其好处是数值计算减少大 量时间． 本文受文【2]的影响，研究如何快速求解等价的一对方程组．采用的算法是矩阵压缩与分 块迭代算法．与文[2】的本质区别在于采用的算法不一样，文[2]采用配置法离散，压缩系数矩 阵，再采用多层扩充法求解压缩后的一对方程组，本文采用Galerkin离散、压缩系数矩阵，再 采用分块多层迭代算法．受文【2，3】的影响，给出了后验参数的选择方法．从后面可以看出，本 文采用的多层迭代算法，化简后是二步迭代算法，不同于一般的单步多层迭代算法． +2012年9月22日收到． 江西省教育厅科学技术研究资助项目(GJJl0586)． 万方数据 计算数学 2013钽 2．多尺度投影方法 这一节介绍求解用Tikhonov正则化后得到的第二类积分方程的多尺度投影算法．假设 分别用(·，·)与¨|I表示．定义Fredholm积分算子瓦 EE, (删(s)：=五忡∽雄)氓s 这里k(s，t)∈C(E×E)．下面我们考虑第一类积分方程的求解 瓦z=Y， (2．1) 方程 (QJ+庀+瓦)z：=尼。Y6， (2．2) 其中，Ol0是正则化参数，扰动数据Y6∈x满足