论文部分内容阅读
波导计算在集成光学的研究中具有重要作用。针对更多实际应用的周期非均匀的波导结构，本文发展了一些基于DtN映射的M算子方法和逆基本解算子方法(IFS)。与传统的算法相比，这些算法引入DtN映射，将微分方程边值问题转化为算子初值问题，并且充分利用了光波导结构的几何结构特性，引入只与光波导结构有关的中间量，构建算子递推格式，对周期结构可大量降低计算量。在计算算子M时，用三点四阶差分格式近似计算二阶导数的方法的产生系数矩阵为块状三对角矩阵，较Chebyshev配置法[1]大大提高了计算效率，且在边界面处充分利用原方程的信息，减少因近似计算带来的误差。用Riccati法[2]求解逆基本解算子的方法，更适用于介质折射率光滑性不好的情况。这些算法主要是针对非齐次Helmholtz方程[3]，也适用于齐次的情形。算例表明，这些算法均是高效的算法。