电阻抗成像的数值模拟与分析.pdf

DoctoralDissertation ShandongUniversity 初始猜测解，这样，均匀分布下的Jacobi矩阵可预先计算好并存储起来，实际计算时只 进行一步迭代，从而有助于快速求解．但这样的方法也限制了其适用范围和求解精度． 本文作者在袁益让教授的精心指导下，在大量的数值试验的基础上，选择Levenberg- Marquardt方法与偏微分方程数值解法相结合，对电阻抗成像进行数值模拟研究．L—M 方法是一种信赖域方法，也是一种正则化的迭代算法，因而可以用来稳定地求解非线性 model)和 反问题(【18，19】)．本文主要考虑两种电阻抗成像模型：连续模型(continuum 电极模型(electrode model)．分别应用有限元、体积元和块中心有限体积方法在规则和 不规则区域，对二维和三维EIT问题进行求解．其中体积元和块中心有限体积方法是 首次用于求解电阻抗成像问题．对于连续模型，我们用精确解来进行数值模拟，验证模 型的正确性和算法的可靠性及可行性，然后应用到电极模型上去，并在电极模型下给出 了适用于以上三种数值方法(对称格式)的一类电流模式，将计算量最大的Jacobi矩阵 的计算彻底简化． 全文共分三章． 第一章介绍基于有限元方法的电阻抗成像的数值模拟和分析．§1．2—1．5研究了基 于有限元的二维EIT问题的数值模拟的方法和技术．采用了两种不同的网格：三角形 网格和四边形网，给出了伪单元刚度矩阵的计算公式，进行分析比较．在矩形区域上，采 用连续模型，对边界条件施以不同幅度的扰动，分别用不同空间步长的数值解和精确 解进行对比，验证了连续模型的正确性和算法的可靠性及可行性．数值实验表明，当步 长缩小时，反问题的数值解能逼近精确解；数值解关于Neumann边界条件的稳定性优 于Diriehlet边界条件，四边形网格优于三角形网格．§1．5在电极模型上的应用，研究两 种网格上EIT问题的各种数值表现．特别比较研究了固定发射电极模式和轮换发射电 极模式下的数值表现，结果表明轮换发射电极模式优于固定电极模式．在§1．6-1．9研究 了三维EIT问题的数值模拟的方法和技术．针对基于有限元方法的迭代算法计算量大 的特点，给出了基于伪单元刚度矩阵不变性的计算技巧，可大大减少有限元方法所需的 数值积分的计算．提出了一类特殊的电流施加模式，可充分利用每一步迭代中已有的 计算结果来简单计算目标函数Jacobi矩阵，而不再需要另外求解任何大规模线性方程 组．最后给出三维不规则区域上对电阻抗成像的实际数值模拟的结果．§1．乒1．5的主要 Mathematica and 内容已投稿((Acta Scientia》，§1．6-1．9的主要内容已投稿((Computers Mathematicswith Applications)}，并已投修改稿． 第一章的创新之处有：(1)首次采用精确解对EIT问题进行数值模拟，尤其是三 维EIT问题．而前人所做的都是二维的，而且是用有限元方法求得的近似解来代替边 界测量(边界条件)来进行数值模拟．(2)首次提出伪单元刚度矩阵方法，简化了有限元 方法中数值积分的计算．(3)本文提出的Jacobi矩阵的计算方法适合多次迭代，并使每 次迭代中求解正问题的次数降到最低． 第二章介绍基于有限体积元方法的电阻抗成像的数值模拟和分析．本章§2．2—2．5 为二维区域上的数值模拟和分析．将体积元方法用于二维电阻抗成像的数值模拟，给出 DoctoralDissertation ShandongUniversity 了Jacobi矩阵的计算公式．并在正方形区域上，用精确解进行数值模拟，结果验证了连 续模型的正确性和算法的可靠性和可行性．在2．5应用于电极模型，给出了在扇形区域 上图像重建的仿真结果．本章§2．争2．9为三维区域上的数值模拟和分析．对于三维椭圆 方程Neumann边值问题，提出了四面体单元上的一类对称体积元格式，并证明了格式 的半正定性及解的存在性；引入单元形状矩阵的概念，简化了系数矩阵的计