介质导电率成像的反演问题及算法实现

来源 :东南大学 | 被引量 : 0次 | 上传用户:oppoyy
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
电阻抗成像(EIT,Electrical Impedance Tomography)是当今生物医学工程领域的重大研究课题之一,数学上对应于一类经典的椭圆型方程的系数反问题.传统的EIT 技术由于仅利用了介质外部的直接测量数据,得到的图像的分辨率较低,很难满足应用问题的需要.近年来发展起来的核磁共振成像技术(MREIT,Magnetic Resonance Electrical Impedance Tomography)是改进成像精度的一个新的研究方向,它利用间接测量到的介质内部的磁场信息来反演介质内部的导电率,由此有望得到稳定的高精度的成像结果。数学上而言,它是一类由椭圆型方程边值问题解的内部信息来反演方程的待定系数的新的反问题.由于利用了解在内部的有关数据,而不仅是边界数据,因此问题的不适定性有所减弱. 本文讨论该反问题在一个具有明确医学背景的模型上的算法-调和B<,z>算法的实现问题.该模型可以用下面的带积分边值条件的椭圆型方程定解问题来描述:在分析了算法基本思想和引起算法不适定性的因素的基础上,本文重点考虑了利用正则化方法处理该算法中扰动输入数据的求导不适定性的问题,在此基础上提出了B<,z>反演方法对扰动数据的数值实现.通过计算机模拟实验验证,取得了令人满意的结果。 本文的工作包含了以下四个方面的内容.首先,与目前已有的模型问题相比,我们讨论了一个更为现实的模型,利用的是边界电极上的总的输入电流的强度,数学上这是一个带有积分边界条件的定解问题.其次,我们把此问题转化为一个标准的椭圆型方程混合边值问题,自编了有限元程序得到了有效的正问题的数值解,实现了反演数据的模拟,同时也为B<,z>反演方法中方程的迭代求解提供了一个必要的工具.再次,我们考虑了利用扰动输入数据B<,z>时,算法的稳定实现问题.调和B<,z>反演方法中的关键一步是要对输入数据进行Laplace运算.众所周知要对扰动数据求数值微分是个不适定的问题,对此我们应用基于Tikhonove正则化方法的五次样条函数的数值微分方法,提出了一个稳定的数值求导方法,大大减弱了迭代算法中的不适定性,并在适当的假设下得到了误差估计.最后,我们给出了所提出反演方法的数值实现,尤其是对扰动数据的实现效果,验证了本文的理论结果。
其他文献
本文以我国目前中小型水库坝体的使用现状为切入点,从设计与实际形成效果的联系等几个方面进行分析,总结出一些切实可行的设计方法,总结了坝体渗流分析的各种方法,并综述了水力学
期刊
期刊
历史教育与社会主义核心价值观教育二者有其内在的联系,本文从史学与社会主义核心价值观教育的关系、史学在高校群团组织开展核心价值观教育的优势、高校群团组织利用史学功用
数学教学中的一题多解是一种备受推崇的教学手段,但随之而来的问题是教师对一题多解背后的目的、数学价值等把握不准,导致灌输式罗列解题方法,即存在目的不明确、重点未突出、能
期刊
请下载后查看,本文暂不支持在线获取查看简介。 Please download to view, this article does not support online access to view profile.
我国城市化的快速发展在城市中最直观的表现之一就是城市建筑的高层化。高层建筑的建设己经越来越成为我国城市发展过程中的一个普遍现象。随之而来的高层建筑与其他建筑间的
期刊
二维矩形切割问题广泛存在于各种工业部门中,如生产钢铁制品、纸张、木材、皮革、玻璃等行业。在所有这些应用中,通常都是先生产出少数几种标准规格的大件产品,然后再切割成用户
期刊
本文考虑如下带Hardy。项的半线性椭圆问题的球对称解及非球对称解的存在性.这里Ω={x|x∈R,n≥3,a