光纤的模式MATLAB图模拟.docxVIP

光纤模式图 matlab 模拟 ************** 摘要 : 光纤通信是现代化通信的支柱 , 在光纤通信中 , 光纤是最重要的部件之一。本文利 用电磁波动理论推导了光在光纤中的传输模式的本征方程 , 并使用 Matlab 软件绘出不同条件下的模式图 . 关键词 : 光纤模式 ; 电磁波动 ;Matlab 一、引言 对光纤中光的传播理论的研究 , 可以有多种方法 , 比如射线法 , 标量近似分析法等 , 但为了更广泛地描述光纤波导中光的传播 , 更详细地研究光纤的传输特性 , 就必须运用波动光学理论对光纤进行分析 . 本文从麦克斯韦方程的求解出发推导光纤的传播模式本征方 程并利用 Matlab 模拟其模式图 . 要对光在光纤中的传播特性有详细的理解 , 必须依靠麦克斯韦方程22 , 结合问题中的边2 界条件 , 求解电磁矢量场 . 求解的方法一般是 W a k0 n2 :1 、先求出亥姆霍兹方程组以及电磁场纵向 分量 Ez和Hz的具体形式 .2 、把 Ez 和Hz有具体形式代入麦克斯韦方程以求取其他电磁场横向 分量 E 、 r 、 H 、、利用界面上电磁场 E 和 H 切向连续条件，求取模式本征方程 [1] . E 二、波动方程 由麦克斯韦方程组 , 我们知道 , 光纤中电磁场的波动方程可以写成 : (1) 式中参量ε表示介质的介电常数 , μ表示介质的磁导率 . 对于在圆柱形光纤中传播的电磁波 . 电场和磁场具有如下形式的函数关系： (2) 式中β为光纤中导波沿 z 轴方向的传播常数 , 其值由纤芯———包层界面处的电磁场边界条件决定 . 不同的β值对应于不同阶的导波模式 , 它们的场分布也不同 . 将式 (2) 代入波动方程式 (1) 中, 可得到矢量亥姆霍兹方程 , 即 (3) 在柱坐标系中 , 只有沿 z轴方向的单位矢量与场点位置无关，所以，在柱坐标系中，只有Ez和Hz才满足标量亥姆霍兹方程 , 可得 (4) 解得方程 (4), 可得到电场和磁场的纵向分量 Ez和Hz. E ( x, y) E (r , )e j ( z t ) (5) H (x, y) H (r , )ej ( z t ) (6) 式中 (r ≤a) 表示在纤芯内部 , 称为归一化横向传播常数 . 其大小随纤芯内 场的不同模式而变 .(ra) 表示在纤芯的外部区域 , 称为归一化横向衰减 常数 , 其大小及符号反映了包层中场的状态 . 有了场分量 Ez和Hz的表达式 , 再利用麦氏方 程组即可求出场的2其它四2个分2量 E 、 Er 、 H 2 、Hr. 2 2 H 0 E n k0 E 0 H n k0 三、本征方程及模式图 本征方程 光纤中传播模式及传输特性都是由它自身的本征方程确定的 . 在光纤的基本参量 n1,n2,a ,k0 已知的条件下 , U,W仅与传播常数β有关 . 用所导出的各个区域中电磁场的表达式 , 再利用电磁场切向分量在纤芯 - 包层界面上 (r=a) 连续的条件 , 就可以救出模式本征方程 , 也称特征方程 . (7) 各类模式 根据 Jm(u) 的振荡特性 , 对于一特定的 m值 , 本征方程存在着 n个根 . 当m=0时，若 , 相应于 , 只有磁场纵向分量 . (8) 当m=0时，若 , 相应于 , 只有磁场纵向分量 . (9) 当m0时，混合模式 HE模和 EH模 (10) 其中 : Matlab 模拟流程图及模式图 为了分析导波模的传输特性 , 就需要得知各模式传播常数β随光纤归一化频率 V的变化情况 . 这可通过对本征方程 (10) 求解而得出 . 其解可写为 (11) 式中 方程 (11) 是超越方程 , 在截止和远离截止的情况下 , 可以将它简化成简单的形式求解 , 得出各种矢量模式的截止频率 Vc和远离截止时的 u值, 从而进行传输特性的分析 . 而在一般情况下 ( 不局限于截止和远离截止两种状态 ) , 就需使用计算机对特征方程 (7) 求数值解 . 其计算流程图如图 1所示 . 图2分别给出了依照此流程绘制出的 TE模、HE11模及 EH11模的模式 图 . 图1 计算β / k0 -V 曲线的程序流程图 (a) (b) 图 2 β/ K0 - V 曲线模式 四、结束语 (c) 光纤中的传播理论已为人们充分了解 , 本文直接由麦克斯韦方程组出发 , 精确地求出 图 电磁场各个分量 , 根据光纤芯———包层界面处电磁场的边界条件 , 在计算机上通过数 值求解 , 得到传播常数及光纤中的场分布模式 , 因而对光纤中传导模的描述完整 , 结果明 A0,B0 Ez (a)EH 0 11 模 确 [4] . 0, H z 五、参考文献：