专业课程设计-光纤光栅的特性分析.doc_第1页
专业课程设计-光纤光栅的特性分析.doc_第2页
专业课程设计-光纤光栅的特性分析.doc_第3页
专业课程设计-光纤光栅的特性分析.doc_第4页
专业课程设计-光纤光栅的特性分析.doc_第5页
已阅读5页,还剩13页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

专业课程设计报告设计题目:光纤光栅的特性分析专业班级:姓 名:班内序号:指导教师:地 点:时 间: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

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论