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

郑硗势：二维椭圆型偏微分方程反问题 二维椭圆型偏微分方程反问题 摘 要 本文的主要思想来源于CT技术带给我们的重要启示：无须深入物 体内部，只利用物体外部测量的数据即可完成内部图像的重建工作。我 们将拉顿变换及其反变换技术应用于二维椭圆型偏微分方程反问题的求 解，只利用物体在边界上的一些可测值就能快速反演出方程的源函数、 势函数和渗透系数，与传统方法不同的是我们不必知道物体内部的数 据。 ， ／拉顿变换的实质是二元函数(即图像)沿直线的线积分。用一束平 行‘X光束照射物体上区域大小为Q。、图像表达为f(x，y)的某～截面， 入射线与X轴正轴闻夹角为0+丌／2，则经过点(s，曰)或点(x，y)的拉顿 变换可表示为： giJ，口)=彤=』』m，y)8(xcosO+ysinO—J)ax4v —00SoO．0S0万 J州j+fu)exp(iO))du， 甜g(s，p)=町2 ’，～ 一∞ 一∞S∞．0S0石 其中孵表示拉顿变换算子或投影算子。拉顿变换提供了CT中获得投影 数据的基本方法：利用物体边界上的一些可测值而非物体内部的数据获 得f(x，y)的投影值，它们是CT反演的原始数据。 拉顿反变换可表示为： 几∽=st-lg∽耻(嘉)r￡老糕艄 摘要 它提供了进行CT图像重建的方法：在CT投影数据的基础上，利用卷 积反投影方法或滤波反投影方法，可以重建出原始数据的图像。 ， 本文主要讨论二维椭圆型偏微分方程的三个反问题：反源问题、反 势问题和渗透系数反问题。在研究过程中，首先将研究区域限定在 (一1，1)×(一1，1)范围内，增加相应的边界条件，并在边界上测量函数的 一些可测值，然后利用推导的一系列函数关系将这些可测值与函数的拉 顿变换联系起来，最后应用拉顿反变换算法反演出原函数。 1．反源问题： 反源问题讨论方程： △妒(x，J，)=f(x，Y) 研究如何利用边界上的可测值来反演源函数f(x，Y)的分布。通过设定 边界条件妒(一1，Y)=2K(K0)中足的取值，将方程转变为： △y。=／：， (z，Y)∈(一1,1)×(一1,1) ％(一1，Y)=0 I；f，，(1，Y)=0 1坐!!兰!业：o 砂 再利用摄动方法将源函数f(x，Y)的拉顿变换与边界函数联系起来 。f：：厂(‘儿)出：—OO—(1,y-(1,Yo))z(1，y。)一—O妒—(-_1,y—o)+o(￡) 麟 饿 从而可以十分方便地利用这些可测值计算出f(x，y)的拉顿变换，并进 而反演出源函数f(x，Y)。 2．反势问题： 反势问题讨论方程： △妒(x，y)一Z(x，y)妒(x，y)=0 郑晓势：二维椭哩型偏微分方程反问题 研究在V(x，y)摄动的情况下如何利用边界上的可测值来反演势函数 V(x，_y)的分布。通过施加边界条件 19(一l，Y)=3 L：{妒(1，y)=1 Iap(z，±1)／砂=o A●o(x，Y)=s+K(x，Y)妒(J，y)，(x，Y)∈(一1,1)×(一1,1) 妒(一1，Y)=3 p(1，Y)=1 O々o(x，±1)／砂=o 并由摄动关系，推导出 f：：y(‘y。)妒(ty。)dx=旦!!!!!二i；5生!生旦z，(x，y。)l：!。+o(占2) 保持势函数不变，施加另一组边界条件： fy(一1，Y)=2 ， L：{y(1，Y)=0 Ia矿(工，±1)／砂=0 得到关系：