Abel变换数值反演的积分算子方法.pdfVIP

V01．29(2009) 数学杂志 No．3 J．ofMath．(PRC) Abel变换数值反演的积分算子方法。 吕小红，吴传生 (武汉理工大学理学院，湖北武汉，430070) 摘要：本文研究了Abel变换的数值反演问题．利用Abel变换的理论反演公式与数值求导 的积分算子法相结合的方法．对反演公式中奇异积分合理处理。获得Abel变换数值反演的一种 算法，并进行了理论分析与数值实验．结果表明该算法具有计算简单、数值稳定等优点． 关键词：Able变换；积分算子l数值反演；奇异积分 MR(2000)主题分类号；65J20；65D25中图分类号：0241．2；0241．4 文献标识码，A 文章编号：0255—7797(2009)03—0383-04 1 引言 Abel变换或Abel积分方程在力学、光谱学、地震学、等离子物理及散射理论等众多领 域有着广泛的应用‘川．但Abel积分方程的求解问题是一个典型的不适定问题，即数据测量 中微小的误差会导致解的巨大误差r引，因此必须采用特殊的方法处理．前解决的方法主要 有两类：一是直接对Abel积分方程应用正则化的方法嘲I二是寻找一种稳定的数值微分方 法然后通过反演公式获得原方程稳定的近似解[4’5]． 我们考虑下述类型的Abel积分方程的求解问题： Tu一鬲妥IF竺≤刮t一厂(s)，0≤s≤1 u’(1) 一而Jo百轳一八’≮5≮1 其中，，(o)=0，口∈(o，1)为常数，甜(￡)为待求函数． 在一定条件下，问题(1)有如下的理论反演公式[5] (2) 础)一T-1f一志J。尚ds，0≤￡≤1 分公式获得问题(1)稳定的近似解‘4。，但其计算量太大．文献[5]提出一种解决奇异积分的 方法，但它在端点及附近点的逼近效果不太好．为解决奇异积分和端点附近的缺陷，对式 (2)右端项作分部积分可得 此时式(3)中的积分不再有奇异性，本文首先采用积分算予方法求得／(o)，／，(s)的稳定近 似，然后用简单的复化梯形积分公式对式(3)进行数值积分，从而得到“(￡)的稳定的近似． 数值实验结果表明该方法具有计算简单、精度高、数值稳定等优点． *收稿日期：2008—08-12接收日期：2009-03-16 基金项目：中国博士后基金资助项目(No．20080431004) 作者简介：吕小红(1969一)。男，湖北黄陂，硕士生。讲师。从事反问题研究． E-mail：Iv】ciaohon92008@sina．com 万方数据 384 数 学 杂 志 Vol_29 2求导的积分算子方法 在实际问题中，方程(1)右端通常为带有扰动数据的^(s)，满足条件 I，一^峙≤艿， (4) 其中艿为已知误差水平．下面我们给出由扰动数据^(s)求／(o)、／，(s)的稳定近似的积分 算子方法． J(x)dx=1，且满足 假定，(s)∈c[口，63，J(z)∈cr[一1，1]为非负的偶函数，I J(z)一雨———I一(1一z2)7． j一。(1一tz刈z 定义积分算子 (历似s)=半f。∥∽弛+缸)出 的近似[6]． 们考虑将厂(5)延拓至[口一Jll，6+J11]， F(s)={