第一类弱奇异核Fredholm积分方程的数值求解.pdfVIP

第 16卷第 2期 应用泛函分析学报 Vo1．16，No．2 2014年 6月 ACTA ANALYSISFUNCT10NALISAPPLICATA June，2014 DOI：10．3724／SP．J．1160．2014．00167 文章编号：1009—1327(2014)02—0167-05 第一类弱奇异核 Fredholm 积分方程的数值求解 闵 涛，胡 刚，晏立刚 (西安理工大学理学院，西安710054) 摘 要 第一类弱奇异核 Fredholm 积分方程 由于奇异及本质的不适定性，给求解带来很 大难度．本文首先利用克雷斯变换将方程转化，并对转化后的方程进行高斯 一勒让德离散， 得到一离散不适定的线性方程组，结合正则化方法对该类问题进行数值求解．最后给出了 数值模拟，验证 了本文方法的可行性及有效性． 关键词 Fredholm积分方程；正则化；克雷斯变换；高斯一勒让德离散 中图分类号 O241．83 文献标识码 A 1 引言 第一类弱奇异核 Fredholm积分方程由下式给出 厂1 ／ l—Y1一“f(y)dy=g(x) (1) I，一1 其中一1 1， 为常数且满足 0 1，此方程在数学物理反问题 [1-2】，静电学中用测得的半 球形导体上的电势求导体上的电荷密度分布 3【)9]，在几何概率，干涉和体视学 【]等方面有着非常重 要的应用．众所周知，只有极个别的特殊情况能找到其解析解，大量的弱奇异核积分方程需要去寻 找其数值解．而数值求解也面临着两个主要困难，其一是按通常的数值求积公式进行离散是行不通 的，原因是积分在X=Y时是不可积的．其二是第一类Fredholm积分方程本身是不适定的 [5--6】，数 据的微小变化将引起解非常大的变化．受文献 [3]和 6【]的启发，本文首先利用克雷斯变换将方程进 行转化，然后利用高斯 ．勒让德求积公式 []进行离散．最后通过正则化方法进行求解． 2 克雷斯变换 由于积分方程方程 (1)在点Y= =叩，叩∈[--1，1]时，积分不可积，使得对于这类方程不能直接 进行离散，根据文献 [1]1对变量 X和 Y进行变换，变换的具体过程如下：令 { ；， ， 其中 ()= ，()=(一 +t 1 )=2p ，p 为正整数．经验证，当一1 t 1时，变换保证了变量 ∈--【1，1]，Y∈-【1，1]．将变换带入方程 (1) 中得 ． ／I()一叫(s)I一，(叫(s)) (s)ds=9(叫(t)) (2) 给 (2)式两边再同时乘以W )得到 ／f()一w(s)l一，((s)) (s) ()ds=9(()) () 收稿 日期：2013—03—20 资助项目：国家 自然科学基金重大项目；国家自然科学基金 51179151) 作者简介：闵涛 (1963一)，男，陕西西安人，教授，博士，研究方向：科学工程计算与计算机模拟，数学物理反问题等 E-mail：mintaoQxaut．edu．cn；胡刚(1987-)，男，陕西西安人，硕士研窑生，研究方向：科学工程计算与计算机模拟． 168 应用泛 函分 析学 报 第 16卷 设 0(t)=W )， ())，∈(t)=W )9(伽(t))，整理得 -／11 )I一)叫)d8=) 为