均匀布拉格光栅的原理及MATLAB反射谱仿真

1、均匀布拉格光栅的原理及 MATLAB射谱仿真张睿一、 前言光纤光栅是纤芯折射率受到周期性微扰而形成的一种全光纤无源器件，自问世以来，由于其与光纤通信系统兼容、体积小、插入损耗低、结构简单、成本低等等，广泛应用于光纤通信、光纤传感、光信息处理等领域，对于光纤光栅的分析，通常主要用耦合模理论、付立叶变换理论与传输矩阵理论，本文主要得用传输矩阵理论对均匀布拉格光栅的反射谱进行理论分析和仿真。均匀布拉格光栅的原理假设光纤光栅的模型如下：八_ _B -z iZi+1图光纤光栅的输入与输出如图可知输入为：A z、B zi 1 ;输出为：B zi、A zi 1 ，但是为了表示方便，输入为：A z、B zi

2、,输出为：A zi 1、B zi 1 。利用麦克斯韦方程组可以得到光波在光波导中的耦合模方程:dAdzdBdz* Ae j(2z)其中：一由边界条件：可以得到相移光栅的传输矩阵:zi(4)ziziziFzizi 1A ziB zi其中:FZiZi 1sl1Sl2(6)S11cosh(s(ziS12S21s21s22Zi)j -sinh(s(4 1 sj sinh(s(4 1 s4)j -sinh(s(Zi 1zi)scosh(s(zi 1为光纤的耦合系数。整个相移光栅的传输矩阵可以表示为:其反射率可以表示成:MATLA皈射谱仿真Zi)Zi) jsinh(s(41s砧1.Fz2Z1(9)仿真所用

3、的参数为：布拉格光栅的中心波长1、在kl相同情况下的反射谱4)(8)1550nm,光栅有效折射率。在kl=5的情况下，当l=2mn 5mm 10mm 20mmW,其反射谱如下: 101ngj曰07ta=5 ；l=2rmil=lOnnrnl=2ftnm由以上图片可知2kSrnm皑怅ir i.j!100%昨5k2QI巧EE但是边缘波长的反射率也越来越高，因而选择较低的耦合系数，可以得到较窄的带宽的波长。kl不同情况下的反射谱l越长的情况下，光栅的反射谱越来越窄，在中心波长的反射率最高由图可知，k越大，即耦合系数越高时，反射谱的带宽越宽，中心波长的反射率都能达到特性对于光纤通信器件，光纤传感器件具有

4、很高的实用价值。02Hie Reilejtwfe! urFBGOb23匚如=一_=加三国四、MATLA魄程序% 传输矩阵的功能函数 function f =tansmit_fiber(L,kappa,num,n_eff,lambda,lambda_Brag)delta(num)=2*pi* n_eff *(1./lambda(num)-1./lambda_Brag);s(num)=sqrt(kappa.A2-delta(num).A2);s11(num尸cosh(s(num)*L)-i*(delta(num)/s(num)*sinh(s(num)*L);s12(num尸-i*(kappa/s(

5、num)*sinh(s(num)*L);s21(num尸i*(kappa/s(num)*sinh(s(num)*L);s22(num尸cosh(s(num)*L)+i*(delta(num)/s(num)*sinh(s(num)*L);f=s11(num) s12(num);s21(num) s22(num);endfunction FBG TransferMatix(n eff,L,lambda Brag)变量定义% n_eff% l % lambda_Brag % kappa % delta光栅有效折射率光栅长度光栅中心波长光波模式互耦合系数模间失谐量%disp('输入的参数如下)

6、；disp('光栅有效折射率为：');disp(n_eff);disp('布拉格光栅的中心波长为');disp(lambda_Brag);disp('布拉格光栅的长度为：');disp(L)lambda=1e-9*linspace(1308,1312,500);kappa_L=5;kappa=kappa_L/L;%交流耦合系数F=1 0;0 1;for num=1:500f=tansmit_fiber(L,kappa,num,n_eff,lambda,lambda_Brag);f=f*F;r3(num尸f(2,1)/f(1,1);R3(num)=(abs(-r3(num)A2;%反射率endsubplot(2,1,1)p