介质波导的数值模拟.docVIP

介质波导的数值模拟 摘要：我们从麦氏方程出发，采用标量法对波导中电磁波满足的微分方程进行数值求解，并对其结果给出形象的图形表示。对于求解微分方程，我们采用有限差分法，将求解微分方程转化为等价的求解矩阵特征值和相应特征向量的问题，利用MATLAB在处理矩阵运算方面的强大功能，求解出电场在介质波导截面的分布。此外，我们对不同结构、不同折射率分布的波导进行了模拟分析，作为扩展给出了金属波导和光纤的近似数值模拟方法。本文深入浅出的介绍了这种介质波导数值模拟方法，学过电磁场理论并对MATLAB有一定了解的同学只要参考这篇文章便能很快学会这种方法。 第一部分 原理 电磁波方程的导出 在无源简单介质中，麦克斯韦方程写成： 从两个旋度方程消去E或H得到 其中 设，取截面的二维矢量波方程为 其中 即得 如果我仅仅考虑电场的标量方程，则电场大小E(Ex或Ey)满足： 令 得 ----------------------------------------------------------（*） 这便是我们进行数值求解的微分方程。 有限差分法 有限差分法的基本思想是利用网格线将定解区域（场域）离散化为网格离散节点的集合，然后，基于差分原理的应用，以离散点上函数的差商来近似替代该点的偏导数，这样，待求的偏微分方程定解问题可转化为一组相应的差分方程组的问题。根据差分方程组（代数方程组），解出各个离散点上的待求函数值，即为所求定解问题的离散解，再应用插值方法便可从离散解得到定解问题在中个场域上的近似解。对于包括电磁场在内的物理场，应用有限差分法进行数值计算的步骤如下： 采用一定的网格划分方式离散化场域； 基于差分原理的应用，对场域内偏微分方程以及定解条件进行差分离散化处理（构造差分格式）； 对所建立的差分格式（即与偏微分方程定解对应的代数方程组），进行计算机求解。 有关这方面的问题可参考偏微分方程数值解的书，下面结合我们的问题，介绍应用有限差分法求解方程（*）。 问题的求解 我们考虑的介质波导结构如图1所示方型波导生长在SiO2 衬底上，芯层折射率大于包层折射率（如图中所示n1n2）。 显然，如果芯层折射率比包层折射率大的多，电磁波将被限制在芯层中传播，在包层介质n2中，电磁波已经很弱，因此我们将包层介质与空气及衬底边界的电场设为零，在这样的假设下，求解介质波导截面电场分布就转化成下面的微分方程求解问题： 由于波导形状规则，我们很容易将其作网格划分，网格线交点（节点）处的电场大小就是我们要求解的电场离散解（如图2所示）。每一个节点的电场大小都是未知数（除了边界点外），我们求解的是（Nx-2）*（Ny-2）个未知数E（i,j）（i=2,3,4…Nx-1; j=2,3,4…Ny-1），下面我们从偏微分方程（*）中“提取”信息,构造求离散解所需的方程组。 我们分析偏微分方程（*）。首先将偏导数差分化，考虑函数f（x）,取小量△x=h,则 同样的二阶微分, 同样对E(x,y) 的二阶偏导有, _ 代入偏微分方程（*）得： 如图2所示，我们采用E（i,j）表示节点处的电场，则有： 如果记图中A点的电场为为E（i,j），则上式给出了节点B、C、D、E处电场和节点A处电场的关系，即所谓的五点差分格式。对所有节点列出这种关系式，并将其写成矩阵的形式，得到： 其中X是由各节点电场E（1，1），E（1，2）…组成的Nx*Ny 个元素的列向量，A是Nx*Ny行、Nx*Ny列的矩阵，其每一行对应一个节点的五点差分格式方程。作为例子，我们给出如图3所示网格（节点处折射率均为n）的矩阵方程： 从简单的例子中可以看出矩阵A是个数字分布有规律的对称的庞大的稀疏矩阵，在我们的程序中给出了一种巧妙的方法来构造这个矩阵。 至此，我们的问题就转化为求解矩阵A的特征值以及相应的特征向量，从电磁波理论上讲，这里的一个特征向量对应一种电磁场在波导中的模式。而求解求矩阵A的特征值以及相应的特征向量，在?MATLAB中，只要用一个语句就行了（用eigs()函数或者eig()）。求出特征向量后，进行归一化，就得到离散节点处电场相对强弱的分布，再用MATLAB中的作图函数直观地显示出来。 值得注意的一点是我们似乎没有对芯层折射率进行考虑，其实不然，在稀疏矩阵A中已经包含每个节点对应的折射率n（i,j）的信息。正是基于这一点，我们可以考虑多种介质复杂结构的波导的数值模拟，在程序中使用者可以自由地设计波导结构和各个部分的折射率。当然，由于这种模拟方法的局限性，不能对超出矩形和其组合的波导，如三角形，圆形进行数值模拟。自由是相对的。 金属波导和光纤的模拟 我们这里不是从金属波导或光纤中电磁波