论文部分内容阅读
本文主要工作是应用局部间断有限元方法(LDG)求解偏微分方程及利用基于偏微分方程的正则化算法研究生物医学成像中的去噪问题。主要分成五部分。第一部分主要研究了具有复杂几何结构的量子定向耦合器上的量子传输现象的数值求解问题。这个现象可以由一个静态Schrodinger方程来解释,我们运用最小耗散LDG方法来求解该方程。为了保证数值格式的稳定性,本文在边界数值流通量上添加惩罚项。此外,量子传输本身的物理性质决定了其频率的变化主要在y方向上,因而针对该问题的求解,我们除了应用基于多项式基函数的LDG方法外,还应用了基于指数基函数的LDG方法来减少计算量。数值算例验证了该算法的有效性。第二部分针对含波动算子的非线性Schrodinger (NLSW)方程给出了守恒的LDG格式。NLSW方程的一个重要性质是能量守恒,而非守恒的数值格式很容易导致数值解爆破,因此我们在空间上使用LDG方法进行离散,而在时间上使用Crank-Nicholson格式进行离散,从而得到一个全离散的守恒数值格式并在理论上给出了能量守恒的严格证明。同时,我们还进一步给出了线性情形下半离散LDG方法最优误差估计的证明。数值结果证明了守恒数值格式在长时间数值模拟中的优越性。第三部分采用高阶保正LDG格式求解含有爆破解的抛物方程。对于该类方程,如果初始条件及源项都是正的,则根据极大值原理,方程的解也是正的,而不加保正限制器的高阶格式会导致错误的爆破时间和爆破区间。由于该问题是Dirichlet边界条件,为了得到最优收敛阶,我们需要在边界数值流通量上添加惩罚项。数值结果表明,保正LDG格式能够准确捕捉到爆破时间和爆破区间。第四部分主要研究基于偏微分方程正则化算法的去噪方法,主要应用于乳腺超声波弹性成像问题。传统的临床设备只能得到较为准确的轴向(平行于波)位移测算,而对于侧向(垂直于波)的位移测算质量是较低的。但是在多种临床弹性成像应用中(如模型重建及温度成像)同时获取轴向和侧向的精确超声波测算位移却是很重要的。因此,在该问题的研究中,我们主要使用在临床设备上获取的传统超声回声数据来提升侧向散斑追踪的精确度。该算法已被用于计算机模拟数据、组织模拟仿真及人体数据测试。第五部分研究如何减少生物医学领域中磁共振成像中的噪音。由于在相位对比磁共振成像中,时域解析的3D(time-resolved3D)相位对比磁共振成像或磁共振造影(PC-MRA/MRI)只用来判定基本的信息,如:临床环境下动脉瘤(aneurysm)的尺寸和流率。而计算误差和潜在的缺陷都会对磁共振的结果产生不利的影响。但是要表征内动脉瘤血液动力学(intra-aneurismal hemodynamics)中的复杂流体,则需要更精确的血管成像。因此这里主要采用基于偏微分方程正则化的去噪算法减少噪音。计算机模拟数据验证了该算法的性能。