专业课程设计-光纤光栅的特性分析.doc

专业课程设计报告设计题目：光纤光栅的特性分析专业班级：姓 名：班内序号：指导教师：地 点：时 间：2012.6.042012.6.15 电子科学与技术教研室光纤光栅的特性分析摘要 光纤光栅是利用石英光纤的紫外线光敏特性将光波导结构直接做在光纤上形成的光纤波导器件。根据特定的光栅结构，光纤光栅可以做成滤波器、反射器、色散补偿器等。由于光纤光栅器件易于与光纤连接，因此，在光纤通信中有着重要的应用。我们通过传输矩阵的方法利用入射光的反射谱和透射谱来研究不同光纤光栅的特性。我们用Matlab软件对sinc切趾光纤光栅、均匀光纤光栅、取样光纤光栅的不同特性进行仿真，然后我们研究光栅长度对均匀光纤光栅特性的影响，最后我们分析光栅周期对均匀光纤光栅特性的影响。关键字：传输矩阵 均匀光纤光栅 sinc切趾光纤光栅 取样光纤光栅 反射谱 透射谱 光栅长度 光栅周一、引言由于入射光的反射谱和透射谱是反映光纤光栅特性的重要参数，我们通过传输矩阵的方法研究入射光的反射谱和透射谱来反映不同光纤光栅的特性。所谓传输矩阵法，是利用分层的方法把光纤光栅的折射率量化，得出一个表示光纤光栅本身特性的传输矩阵。光纤光栅的模型如下图：.LAYERINTERFACE 图中的分层思想是把光纤光栅的折射率量化为i层，则光纤光栅的传输矩阵计算方方法如下：两层不同折射率相邻表面传输矩 当折射率为的层内传输矩阵，其中，是真空中的波矢量，为折射率为的层距。则整段光纤光栅的传输矩阵可表示为：根据这个计算式，任何光纤光栅的传输矩阵M可以通过折射率n(z)求得，并且取样长度越长，计算精度越大。 我们用表示光纤光栅的反射率，表示光纤光栅的透射率，则： 2、 设计方案1.1 研究不同光纤光栅的反射谱和折射谱1.1.1 均匀光纤光栅折射率随光纤方向的表达式为：，其中，为纤芯折射率值，为调制深度，为光纤光栅的周期。下面我们用Matlab进行仿真，研究其特性。Matlab仿真程序中，我们用表示入射光的波长源程序为： clear; nn=15000; a=10*10(-3)/nn;di=a;i=1;for z=0:a:10*10(-3)n(i)=1.468+5*10(-5)*(cos(2*pi*z/(527.9*10(-9); i=i+1;endi=1;t=1;for k0=1548*10(-9):0.02*10(-9):1553*10(-9) M=1 0;0 1; for i=1:1:nn M1=n(i)+n(i+1),n(i)-n(i+1);n(i)-n(i+1),n(i)+n(i+1)/(2*n(i); M2=exp(j*(2*pi)/k0*n(i)*di),0;0,exp(j*(-1)*(2*pi)/k0*n(i)*di); M=M*M2*M1; endr=M(2,1)/M(1,1);R(t)=(abs(r)2; e=1/M(1,1);T(t)=(abs(e)2; t=t+1; endk0=1548*10(-9):0.02*10(-9):1553*10(-9);plot(k0,R); figure;k0=1548*10(-9):0.02*10(-9):1553*10(-9);plot(k0,T);仿真结果： 反射谱：透射谱：1.1.2 sic 切趾函数光纤光栅折射率随光纤方向的表达式为：，其中，为纤芯折射率值，为纤芯折射率的平均增加值， 为纤芯折射率的最大变化量，为均匀光纤光栅的周期,是切趾函数的参数。下面我们用Matlab进行仿真，研究其特性。Matlab仿真程序中，我们用表示入射光的波长。源程序为： clear;nn=15000;a=0.014/nn;di=a;i=1;for z=0:a:1.4*10(-2); n(i)=1.452+7.5*10(-4)+7.5*10(-4)*sinc(z-0.014/2)/0.001)*(sin(2*pi*z/(535*10(-9);i=i+1;endi=1;t=1;for k0=1550*10(-9):0.02*10(-9):1557*10(-9) M=1 0;0 1; for i=1:1:nn M1=n(i)+n(i+1),n(i)-n(i+1);n(i)-n(i+1),n(i)+n(i+1)/(2*n(i); M2=exp(j*(2*pi)/k0*n(i)*di),0;0,exp(j*(-1)*(2*pi)/k0*n(i)*di); M=M*M2*M1; endr=M(2,1)/M(1,1);R(t)=(abs(r)2;e=1/M(1,1);T(t)=(abs(e)2;t=t+1;endk0=1550*10(-9):0.02*10(-9):1557*10(-9);plot(k0,R)figure;k0=1550*10(-9):0.02*10(-9):1557*10(-9);plot(k0,T);仿真结果为：透射率：反射谱：1.1.3 sic 取样光纤光栅折射率随光纤方向的表达式为：，其中，光栅周期，为切趾函数的参数。下面我们用Matlab进行仿真，研究其特性。Matlab仿真程序中，我们用表示入射光的波长。源程序：clear; nn=15000;a=0.014/nn;di=a;i=1; for z=0:a:0.014n(i)=1.452+4*10(-4)+4*10(-4)*sinc(z-0.014/2)/1.29*10(-4)*(sin(2*pi*z/(535*10(-9);i=i+1;endt=1;for k0=1550*10(-9):0.01*10(-9):1557*10(-9) M=1 0;0 1; for i=1:1:nn M1=n(i)+n(i+1),n(i)-n(i+1);n(i)-n(i+1),n(i)+n(i+1)/(2*n(i); M2=exp(j*(2*pi)/k0*n(i)*di),0;0,exp(j*(-1)*(2*pi)/k0*n(i)*di); M=M*M2*M1; end r=M(2,1)/M(1,1); R(t)=(abs(r)2; e=1/M(1,1); T(t)=(abs(e)2; t=t+1; end k0=1550*10(-9):0.01*10(-9):1557*10(-9) plot(k0,R)figure; k0=1550*10(-9):0.01*10(-9):1557*10(-9)plot(k0,T)仿真结果：透射谱：反射谱：1.2 切趾函数不同，均匀光纤光栅的特性分析折射率随光纤方向的表达式为：，其中，为纤芯折射率值，为纤芯折射率的平均增加值， 为纤芯折射率的最大变化量，为均匀光纤光栅的周期，则为切趾函数。1.2.1 sinc切趾函数光纤光栅当时，为sic切趾光纤光栅，其中，为光纤光栅的长度，是切趾函数的参数。下面我们用Matlab进行仿真，研究其特性。Matlab仿真程序中，我们用表示入射光的波长源程序为：1.2.2 升余弦切趾光纤光栅当时，为升余弦切趾函数光纤光栅。下面我们用Matlab进行仿真，研究其特性。Matlab仿真程序中，我们用表示入射光的波长源程序为： clear;nn=15000;a=0.014/nn;di=a;i=1;for z=0:a:1.4*10(-2); n(i)=1.452+7.5*10(-4)+7.5*10(-4)*(1/2)*(1-(cos(2*pi*z/0.014)*(sin(2*pi*z/(535*10(-9); i=i+1;endi=1;t=1;for k0=1550*10(-9):0.02*10(-9):1557*10(-9) M=1 0;0 1; for i=1:1:nn M1=n(i)+n(i+1),n(i)-n(i+1);n(i)-n(i+1),n(i)+n(i+1)/(2*n(i); M2=exp(j*(2*pi)/k0*n(i)*di),0;0,exp(j*(-1)*(2*pi)/k0*n(i)*di); M=M*M2*M1; endr=M(2,1)/M(1,1);R(t)=(abs(r)2;e=1/M(1,1);T(t)=(abs(e)2;t=t+1; endk0=1550*10(-9):0.02*10(-9):1557*10(-9);plot(k0,R)1.2.3 正弦切趾函数光纤光栅当时，为正弦切趾函数光纤光栅下面我们用Matlab进行仿真，研究其特性。Matlab仿真程序中，我们用表示入射光的波长。源程序为： clear;nn=15000;a=0.014/nn;di=a;i=1;for z=0:a:1.4*10(-2); n(i)=1.452+7.5*10(-4)+7.5*10(-4)*sin(pi*z/0.014)*(sin(2*pi*z/(535*10(-9);i=i+1;endi=1;t=1;for k0=1550*10(-9):0.02*10(-9):1557*10(-9) M=1 0;0 1; for i=1:1:nn M1=n(i)+n(i+1),n(i)-n(i+1);n(i)-n(i+1),n(i)+n(i+1)/(2*n(i); M2=exp(j*(2*pi)/k0*n(i)*di),0;0,exp(j*(-1)*(2*pi)/k0*n(i)*di); M=M*M2*M1; endr=M(2,1)/M(1,1);R(t)=(abs(r)2;e=1/M(1,1);T(t)=(abs(e)2;t=t+1; endk0=1550*10(-9):0.02*10(-9):1557*10(-9);plot(k0,R,-.g)仿真结果：由上图分析可知，不同切趾函数对光纤光栅反射率谱的旁瓣滤除作用不同，比较可得，升余弦切趾函数和正弦切趾函数的反射率谱的带宽相同，但升余弦切趾函数比正弦切趾函数的反射率峰值低。sinc切趾函数的带宽要比升余弦切趾函数和正弦切趾函数大，但是反射率有多个波峰且峰值显著降低。1.3 均匀光纤光栅的长度对其反射谱的影响折射率随光纤方向的表达式为：，其中，为纤芯折射率值，为调制深度，为光纤光栅的周期。我们分别取光纤光栅的长度为：、下面我们用Matlab进行仿真，研究其特性。Matlab仿真程序中，我们用表示入射光的波长。源程序： clear; nn=15000; a=11*10(-3)/nn;di=a;i=1;for z=0:a:11*10(-3) n(i)=1.468+5*10(-5)*(cos(2*pi*z/(527.9*10(-9); i=i+1;endi=1;t=1;for k0=1548*10(-9):0.02*10(-9):1553*10(-9) M=1 0;0 1; for i=1:1:nn M1=n(i)+n(i+1),n(i)-n(i+1);n(i)-n(i+1),n(i)+n(i+1)/(2*n(i); M2=exp(j*(2*pi)/k0*n(i)*di),0;0,exp(j*(-1)*(2*pi)/k0*n(i)*di); M=M*M2*M1; endr=M(2,1)/M(1,1);R(t)=(abs(r)2;t=t+1; endk0=1548*10(-9):0.02*10(-9):1553*10(-9);plot(k0,R,r)仿真结果：理论分析：在折射率调制深度不变的情况下，从图中可看出随长度的增加反射谱越来越尖锐，峰值越来越大。1.4 均匀光纤光栅折射率的调制深度对其反射谱的影响折射率随光纤方向的表达式为：，其中，为纤芯折射率值，为光纤光栅的周期，为光纤光栅的长度。我们分别取调制深度、下面我们用Matlab进行仿真，研究其特性。Matlab仿真程序中，我们用表示入射光的波长。源程序： clear; nn=15000; a=5*10(-3)/nn;di=a;i=1;for z=0:a:5*10(-3) n(i)=1.468+4*10(-4)*(cos(2*pi*z/(527.9*10(-9); i=i+1;endi=1;t=1;for k0=1548*10(-9):0.02*10(-9):1553*10(-9) M=1 0;0 1; for i=1:1:nn M1=n(i)+n(i+1),n(i)-n(i+1);n(i)-n(i+1),n(i)+n(i+1)/(2*n(i); M2=exp(j*(2*pi)/k0*n(i)*di),0;0,exp(j*(-1)*(2*pi)/k0*n(i)*di); M=M*M2*M1; endr=M(2,1)/M(1,1);R