基于积分方程几类不适定问题正则化方法和数值解

《基于积分方程几类不适定问题正则化方法和数值解》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、东南大学学位论文独创性声明本人声明所呈交的学位论文是我个人在导师指导下进行的研究工作及取得的研究成果。尽我所知，除了文中特别加以标注和致谢的地方外，论文中不包含其他人已经发表或撰写过的研究成果，也不包含为获得东南大学或其它教育机构的学位或证书而使用过的材料。与我一同工作的同志对本研究所做的任何贡献均已在论文中作了明确的说明并表示了谢意。研究生签名：_j邋日期：勿12．舡《东南大学学位论文使用授权声明东南大学、中国科学技术信息研究所、国家图书馆有权保留本人所送交学位论文的复印件和电子文档，可以采用影印，缩印或其他复制手段保存论文。本人电子文档的内容和纸质论文的内容相一致。除在保密期内的保密论

2、文外，允许论文被查阅和借阅，可以公布(包括刊登)论文的全部或部分内容。论文的公布(包括刊登)授权东南大学研究生院办理。研究生签名鳞弦f厶Co。Ⅸ摘要本文考虑了三个典型不适定问题的基于积分方程的正则化求解方法，分别是：一元函数高阶导数的数值计算问题、核磁共振成像中调和展算法在输入数据有误差时反演结果的稳定件和误差估计、有界区域卜．的逆时热传导问题．这三个不适定问题的共同特点就是问题的解都满足一个积分方程．具体而言，第一个问题和第三个问题都满足一个第一类的线十牛积分方程，而第二个问题满足一个第二类的非线性积分方程，但积分项中由B构成的核函数的数值计算是不适定的．本文由五个部分组成．第一章，介绍

3、不适定问题和儿种常用币则化方法及其诈则化参数选取的基本理论以及已有的结论．概述了本文研究的三类不适定问题国内外已有的研究现状，并在此基础卜-阐明了本文研究的目的、主要内容和创新点．第二章，研究一元函数的二阶导数具有有限个跳跃间断点时，二阶导数的数值计算问题．这是对计算充分光滑函数的二阶导数已有工作的本质推广．借助十二阶导数算子的基本解，将求一元函数二阶导数的问题转化为一个第一类线件积分方程的求解问题．对此不适定的问题，考虑剑二阶导数是分片连续的，本文采用全变分_F则化方法求其币则化解，给出了该方法的收敛性分析，并利用Bregman距离得到了先验条件下收敛速度的估计．另外我们还证明了若一元函

4、数的二阶导数有脉冲奇件，则利用本章所给方法计算的二阶导数的诈则化解的L2范数会爆破．这一结论为确定非光滑函数的不可导点提供了理论依据．最后通过数值算例说明了该方法的有效件和数值稳定件．第三章，考虑核磁共振成像(MREIT)技术中的调和B算法在输入数据有扰动时的稳定性和误差估计．调和成算法的基本思想是基于散度型椭网方程并利用介质内部磁场的一个分量B：的信息来重建导电率．该算法借助于二维拉普拉斯方程的基本解，得卅导电率满足一个第二类非线t牛积分方程(该积分方程中核函数需要计算晚的Laplace运算)，进而用迭代算法求其近似解．调和玩算法虽然在最近十年来得到了快速的发展，而且从数值模拟和模型试验

5、结果来看，可以得到物体高分辨率、高精确度的导电率，但其严格的理论分析却非常少．已有的结果是在精确的输入数据下，该算法的误差估计及其收敛十牛的严格证明．但在输入数据有误差时迭代算法的误差估计还未考虑，这个问题是和算法的收敛件紧密相连的．在本章中我们假设：1)由扰动输入数据计算其L印lace运算的误差是可计算的(该误差依赖_J二采取的正则化方法)；2)待求的导电率仃+满足

6、

7、V盯+怯。f矗1冬E，其中E是充分小的常数．在这样的假发条件下，我们讨论了输入数据有误差时，调和晚算法的误差估计．第四章，研究有界区域卜．的逆时热传导问题，给m了曲种正则化方法：1)转化为第一类nedholm积分方程的方法

8、．首先由热传导过程的特征函数系统构造了一个积分算子，将逆时热传导问题转化为一个第一类的线件nedholnl积分方程求解问题．对此问题，我们考虑了初始温度不连续情形的数值求解，这个问题的研究还在初始阶段．和第二章一样利用全变分币则化方法进行求解，给m了该方法的收敛性分析，并借助于Bregman距离得到了收敛速度的估计．和第二章相比，这里的积分核函数具有更高的光滑件，因而求解更为困难．文中给m的数值结果验证了本文关于收敛速度的估计．2)多参数终值数据拟合法．该方法是将扰动的终值数据按热传导过程的有限项特征函数系展开，这里展开式的项数和有限项展开的误差看成两个诈则化参数．利用近似展开的系数，我们

9、可以显式构造}H逆时问题的正则化解．在初值有界的先验假设下，我们得到了诈则化解关于所有时刻收敛率的一致估计．最后通过数值算例说明了该方法的有效性和数值稳定性．此方法的优点在于：首先，任意时刻的温度分布具有统一的显式表达式；其次，此方法所有的不适定性都集中在终值数据的拟合卜，本质卜．仅是对终值数据进行诈则化，因此计算量大为降低；最后，对求解区域的维数、热传导问题边值条件的类型没有严格的限制．最后，在第五章我们对全文作了简单