有限元近似误差对EIT正逆问题影响的定量研究及算子分解方法的应用

来源 :第四军医大学 | 被引量 : 2次 | 上传用户:bmw335
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
电阻抗断层成像(Electrical Impedance Tomography, EIT)在数学上实际是一类椭圆型偏微分方程正反问题。已知模型的几何结构、激励信号的源参数、阻抗分布求解模型的电势分布称之为正问题(forward problem);根据边界电压的测量值求阻抗分布或变化,被称为逆问题(inverse problem)。正问题的研究主要包括数学物理理论描述、数学模型建立和大规模方程组的求解等。逆问题在数学上是不适定的、病态的,即很小的误差会对解产生很大的影响,这导致了图像重建问题对于正问题的误差包括测量误差、模型误差高度敏感。近几十年,测量误差的研究颇为广泛,而模型误差的研究报道不多。通常,模型误差主要来源于求解域的有限元近似、未知的接触阻抗及未知的被测组织的具体形状等。这些模型误差对于EIT问题是至关重要的,决定着重建效率和图像质量。其中有限元的近似误差是数学与工程理论基础研究的一般问题,且在EIT领域的研究报道很少,因此本文的研究立足点为有限元近似误差对EIT正逆问题的定量影响。有限元近似误差主要取决于网格剖分密度和单元插值函数阶数两个方面。本文在大量的数值试验的基础上,具体从以下几个方面做初步讨论:(1) EIT正问题有限元解的误差分析通过不同的剖分规模、不同阶次的插值函数,在二维正方形区域选择特殊的椭圆型偏微分方程例子,比较三角形与四边形剖分的数值结果。结果表明,计算精度随着单元数目和插值函数阶次的增加而提高;随着插值函数阶次的提高,计算精度增加的倍数也逐次提高;并且单元数目愈多、插值函数阶数愈高,计算时耗愈大。推导二维均匀圆域模型、三层同心圆域模型的解析解公式,应用解析方法对三层同心圆的数值解进行讨论,比较FEM数值解与解析解,定量地分析FEM数值解的计算精度、收敛阶、相对误差等数值指标。结果表明,对于某一特定的高阶单元如二阶单元、三次单元,增加剖分规模对于其计算精度的提升作用有限,不同的剖分规模误差只在小数点后三位有变化;对于某一特定的剖分规模,应用高阶单元对提高精度作用不大,如二次、三次、四次单元间的误差只在小数点后两位有变化。(2)有限元剖分对二维、三维图像重建的影响通过仿真实验,给出了三种评价图像的数值指标:重建图像质量函数D、结构相似度SSIM和误差总和TE,运用非迭代的线性近似算法—NOSER和迭代算法—基于TV正则化算法,分别定量地讨论了剖分规模对重构图像的影响。数值试验表明,在NOSER算法中,随着二维圆域的网格剖分不断加密,三种数值指标显示图像质量越来越高;通过基于极大后验估计的误差模型的误差补偿后,图像质量有所改善;加入不同信噪比的噪声后,虽然噪声对重建图像造成了诸多伪影、斑迹;但是纵向比较来看;剖分不断加密,图像质量依然越来越高。当网格剖分规模达到一定程度时,图像质量提高幅度有限,故对于线性近似的非迭代算法来说,选择中分辨网格即可。采用迭代算法,低分辨网格的模型随着迭代次数增加图像质量变高,误差总和TE减小;而高分辨网格的模型在迭代过程中收敛速率变慢或者不收敛,误差增大,重建过程困难。故对于基于全变差正则化的迭代重构算法选择低分辨网格和合适的迭代次数就可以重建高质量的图像。针对三维全电极模型,随着单元数目、节点数目的增加,边界电压误差减少;用于计算正问题的网格从低分辨到高分辨分为七个等级,用于重建的网格均采用中分辨率网格。数值试验表明,从低分辨网格到高分辨网格,图像质量越来越高,伪影减少,边界越发清晰。(3)算子分解方法在EIT中的应用算子分解方法是图像重建方法中的一种直接方法。本文得到了算子分解方法的主要数学模型和理论构架,以及算子分解方法的正则化理论。通过仿真实验,对于多目标模型算子分解方法能得到较好的重建结果。综上,本文从数值分析的角度,在大量的数值模拟仿真试验的基础上,定量地考察了有限元剖分、不同阶数的插值函数对正问题结果的影响;结合均匀场域的解析解公式对数值解进行了误差分析;结合不同重建算法定量地讨论了有限元剖分规模对图像重建的影响;讨论了求解EIT逆问题三大类方法之一直接法中的算子分解方法,为今后更为全面的EIT算法的定量评估和误差分析打下了前期的理论基础。
其他文献
长输管道沿线地质地貌错综复杂,自然条件恶劣,容易受外界载荷(如地质滑坡、坍塌、疏松土体沉降、湿陷性黄土、地震、水毁等)影响导致大跨度悬空,给管道的安全运行带来严重隐
随着科技的发展,人们需要不断开发高性能的新材料来满足日新月异的各种高新技术需要,尤其是仿生材料越来越引起材料学家们的重视,仿生材料工程也是开发新材料最重要的手段之
本论文研究工作主要在考察真空带式干燥设备结构特点的基础上,对干燥过程进行数学模拟,并以香蕉为研究物料,分析了干燥参数对干燥特性的影响。 首先介绍了真空带式连续干燥设
在压气机叶栅内气流三维流动过程中,二次流的产生和发展不仅直接影响叶栅中的能量损失和压气机的效率,也影响了压气机的稳定性。翼刀技术以直接影响流道内气流的横向流动和叶展
一些抗生素对某些特定微生物的生长繁殖有抑制作用,抗生素对微生物的作用效果好坏反映了抗生素效价。现在国内外使用最多的测定抗生素效价的方法是管蝶法和药敏纸片法,这两种
炼油厂常压塔装置是石油加工的基本装置也是重要装置,其安全运行决定着整个炼油厂的正常作业与安全。腐蚀是常压装置的最普遍的问题,常压系统的腐蚀影响着炼化装置的正常运行及
引起脑损伤的原因有多种,缺血、缺氧、中毒、感染以及神经性疾病等因素都有可能造成脑细胞的死亡(destruction)或退化(degeneration),引起脑损伤。其中,由中风和心脏停搏引起的
在环境污染和能源问题日趋严峻的背景下,空气制冷循环作为能满足大气生态环境要求的最佳制冷循环己经引起了人们的关注。而目前我国居民用空调使用的大都是以R-22制冷剂为主的
液压传动与控制技术在各个行业中的推广和渗透对液压系统本身提出了更高的要求。通过对实际应用中出现的一些问题进行调查研究,设计出了一套具有压力滤波和稳流功能的液压元件
随着现代高新技术的发展,在化工、航天、制冷、电子、生物工程等工业领域广泛地存在着微细通道内流动与换热现象,同时,由于微细通道流动与换热出现了一些不同于常规大直径管道的