光纤光栅的特性分析
- 格式:doc
- 大小:610.50 KB
- 文档页数:8
一、 实验要求
对光纤光栅进行特性分析;
分析,光栅长度、分层数、谱宽等参数对反射光谱的影响;
利用MATLAB进行程序设计,通过软件仿真的形式实现
二、 实验原理
光纤光栅是利用光纤中的光敏性而制成的。
光敏性是指当外界入射的紫外光照射到纤芯中掺锗的光纤时,光纤的折射率将随光强而发生永久性改变。人们利用这种效层内折射率看成不变利用00exp()0=0exp()iiLAYERiijkndMjknd。最后利用矩阵的叠成得到光纤光栅总的传输矩阵应可在几厘米之内写入折射率分布光栅,称为光纤光栅。
光纤光栅最显著的优点是插入损耗低,结构简单,便于与光纤耦合,而且它具有高波长选择性。
光纤光栅有很多种分析法,但目前技术都不太理想。由于反射率是反映光纤光栅特性的一个重要参数。这里利用分层的思想将光纤光栅分层处理,每一层看做折射率n恒定不变,层与层之间折射率不同利用,nnnnijijINTERFACEijijnnMnn11122122MMMM。光纤光栅的反射系数121111rEzMEzM,反射率R=2r。根据不同的入射光波长有不同的反射率,最后绘出反射率与入射光波长的图谱。
以此实现对光纤光栅的特性分析。
三、 实验方案
我们取得是48.645*10^(-4)的光纤长度,15000的分层数,350个点。
1、程序:
clear;
nn=15000;
a=48.636*10^(-4)/nn;
di=a;
i=1;
for z=0:a:48.636*10^(-4)
n(i)=1.452+0.75*10^(-3)*((sin(pi*z/(535*10^(-9))))^2);
i=i+1;
end
wl=1.5541*10^(-6);
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;
end
r=M(2,1)/M(1,1);
R(t)=(abs(r))^2;
t=t+1;
end
plot(R)
2、结果截图:
图一、按步进画图的结果
图二、按波长画图的结果
四、 数据分析
通过对参数的修改我们可得到以下结论:
1.反射率与光栅长度的关系
反射率是光纤光栅的一个重要参数2.14和2.15直接描述了反射率R和光栅长度L的关系。下面图3.1,3.2,.3.3分别描述了不同耦合系数(即不同n)时候,R和L的关系。光栅中心波长nm5.827,V=2.405,)11(*2VnnC折射率扰动n分别为444410*4,10*3,10*2,10*1。
反射率与光栅长度的关系
可见对折射率扰动大的光栅,长度较短也可以达到高的反射率。
描述n分别为334410*2,10*1,10*8,10*6时,反射率与光栅长度的关系。
反射率与光栅长度的关系
描述n分别为555510*4,10*3,10*2,10*1时,反射率与光栅长度的关系。
反射率与光栅长度的关系
2.有效长度cL与折射率扰动的关系
取反射率R=0.9时,光栅长度为有效长度cL,可得有效长度cL与n的关系。
n从0变化到410*5,其他参数仍照上面选取,可以得到如下曲线:
光栅有效长度和折射率扰动的关系
可见在反射率一定的情况下,折射率扰动越大,光栅的长度可以做的越短。图3.5,3.6描述了n从0变化到310*5,0变化到510*5 时候cL与n的关系。
光栅有效长度和折射率扰动的关系
光栅有效长度和折射率扰动的关系
3.谱线宽度
光栅的另一个重要特性是谱线宽度,我们取半峰谱线宽度为光栅线宽。图.3.7描述了n变化对的影响。折射率扰动大会加宽谱线带宽,光栅的谱线宽度还与光栅长度L有关。图3.8描述了410*1n时,线宽和光栅长度L的关系。
根据公式2122])()2[(LnnBFWHM,我们取中心波长mb610*5497.1, n=1.462,710*3.5,mL410*6,510*5~0n
线宽与折射率的关系
线宽与光栅长度的关系
4.光纤光栅反射光谱特性
根据公式:
)(cosh)(sinh)(sinh)0()0(),(222222211SLSSLSLCaaLR…………………………… (2.14)
当0,即n2时,满足相位匹配条件,2.9可以化为: )(tanh2maxCLR
当22C时,入射光的反射率
QLkQLCaaLR22222211cos)(sin)0()0(),(…………………………………………(2.15)
其中222CQ
我们假设光纤各项参数为:mb610*5497.1,n=1.462,710*3.5,mL410*6,310*4n,V=2.405
得到3.9光栅反射光谱特性曲线
光栅反射光谱特性曲线
从上图我门可以得出2个结论:
(1):存在峰值反射率。当δβ=0 时,有峰值反射率;当δβ≠0 时,反射谱有边带存在,边带的反射率大大降低。δβ= 0 时有λ= 2nΛ= B,这称为光纤光栅的Bragg条件,其中B为Bragg 波长。即在一阶Bragg 波长2 nΛ=B 处,有最大反射率)(tanh2maxCLR。
(2): λ=B 时,由上式可以看出:耦合系数愈高,峰值反射率愈高,愈接近于1 ,反射谱边带的峰值反射率也相应增大。
五、 实验结论
反射率在布拉格波长处有个最大值,边带衰减很快,值很小几乎为0;呈正态分布变化。
由图我们可以知道存在一个布拉格波长λ=1.54*10^(-6)m,此时当光栅长度、入射光波长一定时有一个最大的反射率我们可以发现,分层厚度一定的条件下,反射率受光栅长度的影响极大,最大可以接近1。之后反射率就不再随着光栅长度的变化而变化了。
而且我们可以得到光栅长度一定,分的层数越多,实验结果越精确,得到的反射率一般也相对的大一些。实现矩阵连乘时,要注意矩阵之间是左乘还是右乘,注意是从哪边入射。我们得到的实验数据是在均匀分层的情况下得到的,假如是不均匀分层的话,折射率是一个与Z有关的函数,每层不在相当于一个不变的常数。但仍可以利用矩阵的方法分析,得到反射率与入射光波长的图谱,只不过折射率的处理方法有点改变。
六、 心得体会
通过实验,我又掌握了一门新的设计语言MATLAB。它是矩阵实验室之意。除具备卓越的数值计算能力外,它还提供了专业水平的符号计算,文字处理,可视化建模仿真和实时控制等功能。MATLAB的基本数据单位是矩阵,它的指令表达式与数学,工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完相同的事情简捷得多。这对以后的学习和工作都有很大帮助。
七、 参考书目
《光纤通信》 方强、梁孟等 北邮出版社
《光纤光栅原理与应用》 饶云江 科学出版社
《光纤光栅理论基础与传感技术》 张自嘉 科学出版社