一类变系数亥姆霍兹算子的特征模计算及应用研呢究-计算数学专业毕业论文.docx

摘 摘 要 以物理中的光波导为背景，数学上，本文分别对无界变系数亥姆霍兹方程 的模式求解问题及反散射问题的数值算法两类问题进行了研究。第一类问题， 对于缓变光波导，本文推导得到了使用完美匹配层(perfectly matched layer，简 称PML)边界的TM模式下传播常数满足的三角矩阵近似的色散方程；进一步 地，对一般的光波导，本文得到了使用PML边界的TE和TM模式下全矩阵 近似的高精度的色散方程。另外，本文给出了泄漏模和Berenger模的解析的渐 近解公式，当求解这两种模在大模(蚓较大)时的高精度解时，它们可以作为使 用牛顿迭代法高精度求解色散方程时的初始值。对于另一类问题，反介质的重 建问题，本文给出了较好的数值反演算法。 具体地说，对于无界的光波导，本文使用PML边界条件将无界问题有界 化。由于讨论的是波导非均匀的情况，所以使用了微分转移矩阵来推导芯层的 转移矩阵。本文首先推导了TM模式的微分转移矩阵，并对该微分转移矩阵进 行积分，对积分后的形式进行矩阵指数运算。但是因为指数运算难以给出简单 明确的表达式，所以对波导缓变的情况，本文使用上、下三角矩阵近似微分转 移矩阵，最终推导得到了PML边界下TM模式的两个三角矩阵近似的色散关 系，并对波导的大模情况，本文给出了泄漏模和Berenger模的解析的渐近解公 式。因为在波导模式的分析中，相对于小模，高精度求解大模要困难得多，所 以使用牛顿法求解泄漏模和Berenger模的高精度解时，它们的渐近解可以作 为迭代初始值。在另一个研究工作中，为了进一步提高特征值的精度，本文不 对微分转移矩阵进行任何近似，直接对其积分，并对积分后的形式进行矩阵变 换，然后近似地进行矩阵指数的计算，最终结合PML边界条件，推导得到了 特征模式(传播模、泄漏模、Berenger模)满足的全矩阵近似的高精度的解析的 色散关系。这个色散关系对一般折射率连续变化且可导的波导都适用，不在局 限于缓变波导的情况。而且，在一定的条件下，可以精确地计算矩阵指数，此 时上面得到的高精度的色散关系变为精确的色散关系。 对于波导中的反介质重建问题，基于已有的理论方法，本文给出了一种有 效的反演数值算法。首先对于每一个固定的入射频率，使用数值搜索和牛顿迭 万方数据 摘 摘 要 代，得到了一系列高精度的特征值。其次使用多频信息，结合理论表达式的特 点，建立了规模相对小的线性系统。最后通过分析奇异值的分布情况，主要使 用奇异值分解，迭代正则化等技术对线性系统进行正则化反演求解。数值实验 部分对复杂程度不同的杂质(含有低频到高频傅里叶模式的杂质)进行数值验 证，结果表明本文提出了一个较好的反演数值算法。 关键词：光波导，亥姆霍兹方程，完美匹配层，微分转移矩阵，色散关系，泄 漏模，Berenger模，渐近解，反散射问题，线性系统，反演计算，正则化方法 万方数据 AbstractBased Abstract Based the waveguide in physics，in the thesis，modal solution of Helmholtz equation with variable coeffcients in unbounded domain iS firstly considered．then the numerical algorithm of some inverse scattering problem is studied．For the first kind of problem，by using the perfectly matched layer(PML)to terminate the unbounded domain，two asymptotic dispersion relations for TM case proposed，for the waveguides with slowly varied refractive index in the transverse direction。A more general situation，when the refractive jndex profile in the iS continuous and derivable function，iS also considered．Just as the above situation，the thesis gives the dispersion relation with higher precision for both TM and TE cases．Moreover，the asymptotic solutions for Berenger modes and lea