二维椭圆型偏微分方程反势问题.pdfVIP

二维椭圆型偏微分方程的反势问题 郑晓势 张玉海 亓桂明 （山东省计算中心 250014 ） 摘要 将 Radon 变换及其反投影变换原理应用于二维椭圆型偏微分方程反势问题的求解，从 另一个角度解决了小扰动情况下椭圆型偏微分方程的反势问题。 关键词 椭圆型偏微分方程 拉顿（Radon）变换 反势问题 反源问题 中图分类号 O175.25 1 引言 椭圆型偏微分方程在许多物理研究领域 （如量子力学、电磁学、流体动力学和弹性力学 等）都有重要应用。快速而准确地求解其反问题，能够帮助我们了解上述物理现象的过程和 本质，从而更好地研究和控制它们。在文献[1]中系统地讨论了二维椭圆型偏微分方程的反 源问题（Invserse Source Problem ），本文将把拉顿（Radon ）变换及其反投影变换原理应用 于二维椭圆型偏微分方程的另一类反问题：反势问题（Invserse Potential Problem ），以解决 小扰动情况下的反势问题。本文讨论如下二维椭圆型偏微分方程： ∆ϕ −Vϕ 0 （1） 其中 V 称为势，即物理学中常说的势函数，其位置和大小分布的估测有利于更好地了解和 控制物体。反势问题可陈述如下：在ϕ(x , y ) 已知或部分已知、V 发生微小扰动的情况下， 求势函数 V(x,y) 的分布。 由于势函数在许多研究领域都有很好的应用，因此人们在反势问题上做了大量的理论和 实际工作。过去对反势问题的讨论大都局限于某个具体的物理现象和方程，如引力探测和地 震探测中重力势的反问题[2] 、心电图仪中电势的反问题[3] 、旋转抛物线和无限抛物柱面中 的反势问题[4] 等，从不同侧面解决了势函数的反问题，其中以迭代法[5] 和最小二乘法[6] 为 多，这些方法一般都需要事先知道区域内部某物理量的具体分布，且对数据精度要求较高。 但在许多实际问题中，可能根本无法测量区域内部某物理量的具体数值，而且测量值都存在 不同程度的噪声，从而在不同程度上限制了我们准确而快捷地求解势函数。 拉顿变换及其反投影变换是图像处理中的一种重要研 究方法，许多图像重建便是有效地利用了这种方法。它不必 知道图像内部的具体细节，仅利用图像在边界上的一些投影 值，即可很好地反演出原图像。它成功地应用于层析 X 射 线摄影[7] 中。本文将利用它研究二维椭圆型偏微分方程，不 必知道区域内部物理量的具体数值，仅通过物体在边界上的 一些可测值来反演势函数 V(x,y) 。 2 预备知识—拉顿变换及其反变换 设函数 f (x ,y ) ，其中(x , y ) ∈R 2 , 则f (x ,y ) 的拉顿 变换（如图 1）为 1 +∞ rf (s,θ) ∫−∞f (s cosθ−t sinθ,s sinθ+t cosθ)dt (2) 此积分变换显然是f (x ,y ) 在方向θ上的一个投影。 ˆ ˆ 令f (x ,y ) 的Fourier变换为f (ω ,ω ) ，而r (s,θ) 关于s 的Fourier变换为r (ω,θ) 。 1 2 f f [9] 则有 +∞ +∞ ˆ −j ωs rf (ω,θ