光纤的模式MATLAB图模拟

1、光纤模式图matlab模拟*摘要：光纤通信是现代化通信的支柱,在光纤通信中,光纤是最重要的部件之一。本文利用电磁波动理论推导了光在光纤中的传输模式的本征方程，并使用Matlab 软件绘出不同条件下的模式图.关键词：光纤模式；电磁波动；Matlab一、引言对光纤中光的传播理论的研究，可以有多种方法，比如射线法，标量近似分析法等，但 为了更广泛地描述光纤波导中光的传播，更详细地研究光纤的传输特性，就必须运用波动 光学理论对光纤进行分析.本文从麦克斯韦方程的求解出发推导光纤的传播模式本征方 程并利用Matlab模拟其模式图.要对光在光纤中的传播特性有详细的理解，必须依靠有图斯韦点程，结合问题中的边

2、界条件，求解电磁矢量场.求解的方法一般是：1、先求出亥W霍酬0n组以及电磁场纵向 分量Ez和HZ的具体形式.2、把Ez和HZ有具体形式代入麦克斯韦方程以求取其他电磁场横向 分量E、Er、H、利用界面上电磁场E和H切向连续条件，求取模式本征方程1.二、波动方程由麦克斯韦方程组，我们知道，光纤中电磁场的波动方程可以写成：式中参量e表示介质的介电常数，仙表示介质的磁导率.对于在圆柱形光纤中传播的电磁波.电场和磁场具有如下形式的函数关系：(2)式中B为光纤中导波沿z轴方向的传播常数，-值由纤芯包层界面处的电磁场边界条件决定.不同的B值对应于不同阶的导波模式，它们的场分布也不同.将式(2)代入波动方程式

3、(1)中，可得到矢量亥姆霍兹方程，即在柱坐标系中，只有沿z轴方向的单位矢量与场点位置无关，所以，在柱坐标系中，只有 Ez和H才满足标量亥姆霍兹方程，可得(4)解得方程(4),可得到电场和磁场的纵向分量Ezf口Hz.E(x,y) E(r, )e j( z % H (x, y) H (r, )ej( z t)式中(r & a)表示在纤芯内部，称为归一化横向传播常数.其大小随纤芯内场的不同模式而变.(ra)表示在纤芯的外部区域，称为归一化横向衰减常数，其大小及符号反映了包层中场的状态.有了场分量Ez和Hz表达式，再利用麦氏方 程组即可求出场咱它伴来范 E、E、H 2HHr. n2k2H o 三、本征

4、方程及模式图本征方程光纤中传播模式及传输特性都是由它自身的本征方程确定的.在光纤的基本参量n1,n2,a ,k0已知的条件下,U,W仅与传播常数B有关.用所导出的各个区域中电磁场的表 达式，再利用电磁场切向分量在纤芯-包层界面上(r=a)连续的条件，就可以救出模式本征 方程，也称特征方程.各类模式根据Jm(u)的振荡特性,对于一特定的m直，本征方程存在着n个根.当m=0寸，若，相应于,只有磁场纵向分量.当m=0寸，若，相应于，只有磁场纵向分量.当m0寸，混合模式HE真和EH奠其中：Matlab模拟流程图及模式图为了分析导波模的传输特性，就需要得知各模式传播常数B随光纤归一化频率V勺变化情况.这

5、可通过对本征方程(10)求解而得出.其解可写为式中方程(11)是超越方程，在截止和远离截止的情况下，可以将它简化成简单的形式求解， 得出各种矢量模式的截止频率 Vcffi远离截止时的u值，从而进行传输特性的分析.而在一 般情况下(不局限于截止和远离截止两种状态)，就需使用计算机对特征方程(7)求数值解. 其计算流程图如图1所示.图2分别给出了依照此流程绘制出的TE真、HEi模及EHi模的模式 图.图1计算0/ k0 -V曲线的程序流程图(a)(b)图2 (3/ K0 - V曲线模式四、结束语 (c)光纤中的传播理论已为人们充分了解，馋直接由麦克斯韦方程组出发，精确地求出 电磁场各个分量，根据光

6、纤芯包层界面处电磁场的边界条件，在计算机上通过数值学解，得到传播常数及心得反勺智布模小 ,脯户,纤1M导模的描述完整，结果明 五、参考文献：1余守宪等.导波光学物理基础M.北京：北方交通大学出版社,2002.2陈军等.光学电磁理论叫.北京：科学出版社,2005.3陈抗生等.微波与光波队B教0 M.赤江:薪朴吠电出版社,2000.4薛苏云等.阶跃折射率光纤的电磁场模式研究J.河海大学常州分校学报.(4):16-20.附录一 ：Matlab实现TE模程序format longclearn1 = n2= a= k=(aA2)*(n1A2-n2A2);u1=(aA2)*(n1A2);w1=(aA2)*

7、(n2A2);delta 1 =(n 1 A2+n2A2)/(2*n 1A2);delta2=(n 1 A2-n2A2)/(2*n 1A2);i=1;n=n2;for V=0:6kO2=VA2/k;k01=sqrt(k02);for BeiTa=n:n1%n nn=n+ U2=u 1 *kO2-(aA2)*kO2*(BeiTaA2);U=sqrt(U2);W2=(aA2)*kO2*(BeiTaA2)-w1 *k02;W=sqrt(W2);if(U=0 | W=0) break;%disp(sss)elsez1=n1A2*besselj(1 ,U)/(U*besselj(O,U);z2=n2A2

8、*besselk(1,W)/(W*besselk(0,W);z3=1/(UA2)+delta1*(1/WA2-z2)-sqrt(delta2A2*(1/WA2-z2)A2+(BeiTa/n1)A2*(V/(U*W )A4);if(abs(z1+z2) x(i 尸 V;y(i)=BeiTa;i=i+1;z1； z2;%disp(nnn) n=BeiTa; break;endendendendplot(x,y);axis(0 );附录二:Matlab实现HEi模程序format long clearn1 = n2=a=k=(aA2)*(n1A2-n2A2);u1=(aA2)*(n1A2);w1=(

9、aA2)*(n2A2);delta1=(n1A2+n2 A2)/(2*n1A2);delta2=(n1A2-n2A2)/(2*n1A2);i=1;n=n2;for V=0:1k02=VA2/k;k01=sqrt(k02);for BeiTa=n:n1 %n nn=n+U2=u1*k02-(aA2)*k02*(BeiTaA2);U=sqrt(U2)W2=(aA2)*k02*(BeiTaA2)-w1*k02;W=sqrt(W2)if(U=0 | W=0) break;%disp(sss)elsez1=besselj(0,U)/(U*besselj(1,U);z2=besselk(0,W)/(W*b

10、esselk(1,W);z3=1/(UA2)+delta1*(1/WA2+z2)+sqrt(delta2A2*(1/WA2+z2)A2+(BeiTa/n1)A2*(V/(U*W )A4);if(abs(z1-z3) x(i)=V;y(i)=BeiTa;i=i+1;z1;z2;%disp(nnn) n=BeiTa; break;end end endendfor V=1:6k02=VA2/k;k01=sqrt(k02);for BeiTa=n:n1 %nU2=u1*k02-(aA2)*k02*(BeiTaA2);U=sqrt(U2);W2=(aA2)*k02*(BeiTaA2)-w1*k02;W

11、=sqrt(W2);if(U=0 | W=0) break;elsez1=besselj(0,U)/(U*besselj(1,U);z2=besselk(0,W)/(W*besselk(1,W);z3=1/(UA2)+delta1*(1/WA2+z2)+sqrt(delta2A2*(1/WA2+z2)A2+(BeiTa/n1)A2*(V/(U*W )A4);if(abs(z1-z3) x(i)=V;y(i)=BeiTa;i=i+1;z1;z2; n=BeiTa; break; end end endendplot(x,y);axis(0 );附录三：Matlab实现EH1模程序format l

12、ong clearn1 = n2= a= k=(aA2)*(n1A2-n2A2); u1=(aA2)*(n1A2);w1=(aA2)*(n2A2);delta1=(n1A2+n2 A2)/(2*n1A2);delta2=(n1A2-n2A2)/(2*n1A2);i=1;n=n2;for V=0:1k02=VA2/k;k01=sqrt(k02);for BeiTa=n:n1%nnn=n+U2=u1*k02-(aA2)*k02*(BeiTaA2);U=sqrt(U2)W2=(aA2)*k02*(BeiTaA2)-w1*k02;W=sqrt(W2)if(U=0 | W=0)break;%disp(s

13、ss)elsez1=besselj(0,U)/(U*besselj(1,U);z2=besselk(0,W)/(W*besselk(1,W);z3=1/(UA2)+delta1*(1/WA2+z2)+sqrt(delta2A2*(1/WA2+z2)A2+(BeiTa/n1)A2*(V/(U*W )A4);if(abs(z1-z3)x(i)=V;y(i)=BeiTa;i=i+1;z1;z2;%disp(nnn)n=BeiTa;break;endendendendfor V=1:6k02=VA2/k;k01=sqrt(k02);for BeiTa=n:n1%nU2=u1*k02-(aA2)*k02*(BeiTaA2);U=sqrt(U2);W2=(aA2)*k02*(BeiTaA2)-w1*k02;W=sqrt(W2);if(U=0 | W=0)break;elsez1=besselj(0,U)/(U*bess