光纤通信-阶跃光纤模式色散曲线计算光纤材料色散计算、-光纤光栅特性分析

《光纤通信》大作业阶跃光纤模式色散曲线计算光纤材料色散计算光纤光栅特性分析指导老师:专业:姓名：学号：一、阶跃光纤模式色散曲线计算1.原理分析模式色散是指多模式传输时同一波长分量的各种传导模式的相位常数不同，群速度不同，引起到达终端的脉冲展宽的现象。对光纤中光的传播理论的研究,可以有多种方法,比如射线法,标量近似分析法等,但为了更广泛地描述光纤波导中光的传播,更详细地研究光纤的传输特性,就必须运用波动光学理论对光纤进行分析。要对光在光纤中的传播特性有详细的理解,必须依靠麦克斯韦方程,结合问题中的边界条件,求解电磁矢量场.求解的方法一般是:1、先求出亥姆霍兹方程组以及电磁场纵向分量Ez和Hz的具体形式。2、2、把Ez和Hz有具体形式代入麦克斯韦方程以求取其他电磁场横向分量E、Er、H、Hr。3、3、利用界面上电磁场E和H切向连续条件，求取模式本征方程。对于时谐场光波，在均匀、无源介质中，同样满足矢量的亥姆霍兹方程（式1-1,1-2）022EkE（1-1）022HkH（1-2）对于在圆柱形光纤中传播的电磁波.电场和磁场具有如下形式的函数关系（式1-3,1-4）：)(),(),(tzjerEyxE（1-3）)(),(),(tzjerHyxH（1-4）代入，得到亥姆霍兹方程（式1-5,1-6）。02202zzEnkE(1-5)02202zzHnkH(1-6)经过推导，可以得到关于r的贝塞尔方程或修正的贝塞尔方程（式1-7）。（1-7）最后可解得电场和磁场的纵向分量Ez和Hz。);()()();()()(areraWKWKAareraUJUJAEimmmimmmz0)(])[()()(2222120222rRmrnkrrRrrrRr光纤中传播模式及传输特性都是由它自身的本征方程确定的。在光纤的基本参量n1,n2,a,k0已知的条件下,U,W仅与传播常数β有关.用所导出的各个区域中电磁场的表达式,再利用电磁场切向分量在纤芯-包层界面上(r=a)连续的条件,就可以救出模式本征方程,也称特征方程（式1-8）。（1-8）当m=0时，若A=0,B0,相应于Ez=0,Hz0,只有磁场纵向分量。TE模：(1-9)当m0时，混合模式HE模和EH模。(1-10)为了分析导波模的传输特性,就需要得知各模式传播常数β随光纤归一化频率V的变化情况.这可通过对本征方程求解而得出。2.相关代码TE：formatlongclearn1=1.45n2=1.447a=4.5k=(a^2)*(n1^2-n2^2);u1=(a^2)*(n1^2);w1=(a^2)*(n2^2);delta1=(n1^2+n2^2)/(2*n1^2););()()();()()(areraWKWKBareraUJUJBHimmmimmmz222204202221)11()()()(])()()()(][)()()()([WUkmUWVkmWWKWKnUUJUJnWWKWKUUJUJmmmmmmmm0)()(1)()(10000WKWKWUJUJU)]11)(11(44[)1(21)1(21222222122221212221212122WUWnnUmKnnKnnKnnJdelta2=(n1^2-n2^2)/(2*n1^2);i=1;n=n2;forV=0:0.01:6k02=V^2/k;k01=sqrt(k02);forBeiTa=n:0.00001:n1%nnn=n+0.00001U2=u1*k02-(a^2)*k02*(BeiTa^2);U=sqrt(U2);W2=(a^2)*k02*(BeiTa^2)-w1*k02;W=sqrt(W2);if(U==0||W==0)break;%disp('sss')elsez1=n1^2*besselj(1,U)/(U*besselj(0,U));z2=n2^2*besselk(1,W)/(W*besselk(0,W));z3=1/(U^2)+delta1*(1/W^2-z2)-sqrt(delta2^2*(1/W^2-z2)^2+(BeiTa/n1)^2*(V/(U*W))^4);if(abs(z1+z2)0.01)x(i)=V;y(i)=BeiTa;i=i+1;z1;z2;%disp('nnn')n=BeiTa;break;endendendendplot(x,y);axis([06.51.4471.45]);HE11：formatlongclearn1=1.45n2=1.447a=4.5k=(a^2)*(n1^2-n2^2);u1=(a^2)*(n1^2);w1=(a^2)*(n2^2);delta1=(n1^2+n2^2)/(2*n1^2);delta2=(n1^2-n2^2)/(2*n1^2);i=1;n=n2;forV=0:0.001:1k02=V^2/k;k01=sqrt(k02);forBeiTa=n:0.00001:n1%nnn=n+0.000001U2=u1*k02-(a^2)*k02*(BeiTa^2);U=sqrt(U2)W2=(a^2)*k02*(BeiTa^2)-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*besselk(1,W));z3=1/(U^2)+delta1*(1/W^2+z2)+sqrt(delta2^2*(1/W^2+z2)^2+(BeiTa/n1)^2*(V/(U*W))^4);if(abs(z1-z3)0.01)x(i)=V;y(i)=BeiTa;i=i+1;z1;z2;%disp('nnn')n=BeiTa;break;endendendendforV=1:0.01:6k02=V^2/k;k01=sqrt(k02);forBeiTa=n:0.00001:n1%nU2=u1*k02-(a^2)*k02*(BeiTa^2);U=sqrt(U2);W2=(a^2)*k02*(BeiTa^2)-w1*k02;W=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/(U^2)+delta1*(1/W^2+z2)+sqrt(delta2^2*(1/W^2+z2)^2+(BeiTa/n1)^2*(V/(U*W))^4);if(abs(z1-z3)0.01)x(i)=V;y(i)=BeiTa;i=i+1;z1;z2;n=BeiTa;break;endendendendplot(x,y);axis([06.51.4471.45]);3.仿真结果以下两幅图是分别是TE模和HE11模的b~V关系曲线。理论上讲，HE11模的归一化截止频率应为0，对所有的V数都存在，然而仿真过程中有近似计算，因此出现了一些偏差。从图中可以看出TE模的Vc=2.4,与理论计算数值相当。tehe11二、光纤材料色散计算1.原理分析材料色散是由材料自身特性造成。光波不同，折射率不同，光传输速度不同。对于单模光纤，由于只有一个模式在光纤中传输，所以不存在输入光脉冲的模式混合色散。但是，由于实际光源不是纯单色光，单模光纤的折射率又随耦合进光纤的不同波长的光而改变，所以不同波长光的时间延迟不同，从而使输出光脉冲展宽。这种色散取决于材料折射率的波长特性和光源的谱线宽度。上图形象地说明了由于光纤材料色散引起的输出脉冲展宽。所有发射光源都是在一定波长范围内发射的非单色光，当各种波长的光进入纤芯后，由于波长与折射率有关，所以在光纤波导中的光以不同的群速度Vg(λm)在纤芯内以基模传输，波长短的波(频率高)速度慢，波长长的波(频率低)速度快，所以它们到达光纤末端的时间也不同，导致输出脉冲展宽。通过解波动方程和利用边界条件确定传播常数，求得色散表达式。单模光纤还利用了基模的特点和折射率差很小的事实进行简化，从而获得近似结果和有用的经验公式。光纤越长，因材料色散引起输出脉冲展宽越大。因此材料色散用单位长度的展宽来表示(式2-1)。(2-1)Dm(λ)是材料色散系数，由材料折射率的二阶导数(式2-2)给出。)(mDL(2-2)材料的折射率可用Sellmerier公式(式2-3)近似表示：(2-3)对石英光纤，其中：B1=0.696166300B2=0.407942600B3=0.897479400C1=4.67914826×10−3µm2C2=1.35120631×10−2µm2C3=97.9340025µm2则可通过所给参数，求出材料色散系数与波长的对应关系，画出曲线。2.相关代码symsl%定义波长变量lB1=0.696166300B2=0.407942600B3=0.897479400C1=4.67914826*10^3C2=1.35120631*10^2C3=97.9340025c=299792458%定义石英光纤各参数值y=1+(B1*l^2)/(l^2-C1)+(B2*l^2)/(l^2-C2)+(B3*l^2)/(l^2-C3)%中间变量yn=sqrt(y)%n=根号y定义折射率n与波长l的函数关系dn1=diff(n,l)%函数n对波长l求导数dn2=diff(n,l)%再次求导，得到n对l的二阶导数ddd=l*dn2/c%ddd为材料色散系数l1=linspace(1300,1600,300);%定义波长变量的取值范围dddl1=subs(ddd,l,l1)%在波长定义域内的材料色散系数222)(dndcDm32232222122121)(CBCBCBn函数值plot(l1,dddl1)%画图3.仿真结果扩大波长λ的取值范围，能得到，在1.27微米波长处，材料色散系数为零。三、光纤光栅特性分析1.原理分析光纤光栅是利用光纤中的光敏性而制成的。所谓光敏性，是指强激光（在10ns~40ns脉冲内产生几百毫焦耳的能量）辐照掺杂光纤时，光纤的折射率将随光强的空间分布发生相应的变化，变化的大小与光强成线性关系。如用特定波长的激光干涉条纹（全息照相）从侧面辐照掺锗光纤，就会使其内部折射率呈现周期性变化，就像一个布拉格光栅，称为光纤光栅。光纤布拉格光栅是一小段光纤，一般几毫米长，其纤芯折射率经两束相互干涉的紫外光（峰值波长为240nm）照射后产生周期性地调制，干涉条纹周期由两光束之间的夹角决定，大多数光纤的纤芯对于紫外光来说是光敏的，这就意味着将纤芯直接曝光于紫外光下将导致纤芯折射率永久性变化。这种光纤布拉格光栅的基本特性就是以共振波长为中心的一个窄带光学滤波器。该共振波长称为布拉格波长。啁啾光纤光栅的栅格周期不是常数，而是沿轴向变化。不同的栅格周期对应不同的布拉格反射波长，不同波长的入射光在啁啾光纤光栅的不同位置反射,如下图所示。布拉格反射波长随光纤光栅的位置而变化，在某点Z处所对应的布拉格反射波长如下（式3-1）)()(2)(zznzeffB(3-1)啁啾光纤光栅的折射率表示如下(式3-2)：其中是布拉格周期，Φ(z)表示折射率的相位，通常用来描述光栅的啁啾量。)]}(2cos[1{)(zznzneffeff(3-2)啁啾光纤光栅的折射率随z变化曲线如图所示:对于一定长度的均匀光纤Bragg光栅，其反射谱中主峰的两侧伴随有一系列的侧峰，一般称这些侧峰为光栅的边模。如将光栅应用于一些对边模的抑制比要求较高的器件如密集波分复用器，这些侧峰的存在是一个不良的因素，它严重影响器件的信道隔离度。为减小光栅边模，人们提出了一种行之有效的办法一切趾所谓切趾，就是用一些特定的函数对光纤光栅