推广的Abe1积分方程的高精度算法.pdfVIP

第 6期 徐征顺 ，等 ：推广的Abel积分方程的高精度算法 推广的Abel积分方程的高精度算法 徐征顺 (四川I大学数学学院，成都 610064) 摘 要：本文提 出了求推广的Abel积分方程数值解的高精度算法，该算法节省运行时间，精 度高且数值解有后验估计．数值算倒表 明算法计算结果与理论分析吻合． 关键词 ：推广的Abel积分方程 ；数值计算 ；高精度算法；误差估计 中图分类号 ：O241．83 文献标识码：A 文章编号 ：0490—6756(2014)06—1125—11 Efficientalgorithm forgeneralizedAbelintegralequations X己，Zheng—Shun (CollegeofMathematics，SichuanUniversity，Chengdu610064，China) Abstract：Anefficientalgorithm isproposedforsolving thegeneralizedAbelintegralequations．theory andnumericalexamplesshow thatthepresentmethodispowerful，savingstoragecapacityandhashigh accuracy．M oreover，aposteriorierrorestimateisderived． Keywords：GeneralizedAbelintegralequations；Numericalsolution；Highaccuracy；Errorestimate 引 言 考虑如下推广的Abel积分方程 ： )出 一 厂( 0a 1，o≤ z≤ l,fl l (1) 其 中K(x，)、厂(-z)是已知函数 ， (f)是待求函数 ，且K(x，)充分光滑．方程 (1)在理论和工程领域均有 广泛应用 ，诸如半导体 、扩散理论 、冶金、流体力学、人 口动态学、地震学、热传导等许多实际问题 中都需要 求方程 (1)的数值解． 虽然通常的积分方法不能有效的求一般 的奇异积分方程 ，但求各种特殊类型的奇异积分方程数值解 的方法却很多，而针对弱奇异积分方程的求解方法更多[1_8]，如刘亚平和吕涛口提出的关于方程 I 鲁 ()出一f(x)，oa1，o≤≤1的高精度算法．方程』： ()d一，(z)，0a 1，0≤ ．z≤ 1，0 1也可化为上述类型．值得注意的是 ，方程 (1)必须做变换才能应用文献-I3]中 提供 的方法 ，这正是本文的主要工作．本文提出了求推广的Abel积分方程数值解的高精度算法 ，该算法 节省运行时间，精度高且数值解有后验估计．数值算例表明，算法计算结果与理论分析吻合． 2 预备知识 对端点带奇异性的积分方程 收稿 日期 ：2013—11-14 作者简介 ：徐征顺 (1981一)，男，山东曹县人，硕士，主要研究方 向为微分方程数值解．E—mail：xuzhengshundfrj@126．com 1126 四川大学学报 (自然科学版) 第 51卷 ((GG))一一IIGG((zz))dx=：II ((1一一LLzz))g((zz))dx，，卢卢一一11，，口口一一1，，gg(())ECN-Eo，，l1]，， J0 J0 Lyness和 NinhamE推导出了如下的渐近展开式 ： 耋G叫2HG()一 +。 等l ～l +。o、)， 薯GHcG=茎 +。厶 + O(hN)，