使用matlab计算长周期光纤光栅的有效折射率_第1页
使用matlab计算长周期光纤光栅的有效折射率_第2页
使用matlab计算长周期光纤光栅的有效折射率_第3页
使用matlab计算长周期光纤光栅的有效折射率_第4页
全文预览已结束

付费下载

下载本文档

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

文档简介

1、使用matlab计算长周期光纤光栅的有效折射率近几年来,一种在光纤中写入光栅的技术引起了人们的极大兴趣,光纤光栅的出现,给光通信领域带来了极大的变革,尤其是在滤波和传感方面影响尤为深远。按照光纤光栅周期的长短, 通常周期把小于一微米的光纤光栅称为短周期光纤光栅,又称为布拉格光纤光栅或反射光栅,而周期为几十甚至几百微米的光纤光栅称为长周期光纤光 栅,又称为透射光栅。短周期光纤光栅的特点是传输方向相反的模式之间发生耦合,属于反射型带通滤波器,长周期光纤光栅的特点是同向传输的纤芯基模和包层模之间的耦合,没有后向反射,属于投射型带阻滤波器。长周期光纤光栅透射谱的理论模拟在我们研究长周期光纤光栅的各方面

2、属性的过程中,耦合模理论提供了精确有效的解, 从而被我们广泛采纳。长周期光纤光栅的模式耦合主要是指纤芯基模和同向传播的包层模之间的耦合,由光纤的模式理论可以知道在包层中存在很多种包层模,为了方便讨论,我们在这一部分只讨论光纤光栅的纤芯基模和一个包层模之间的耦合。首先,我们定义一下谐振波长,就是在光纤有效折射率调制无穷小的情况下的谐振波长,可以用下面的公式来表示 赤d=AneffA,式子中Aneff表示光纤基模和包层模之间的有效折射 率之差,即 Aneff =neffi neff2。第二步,我们来定义一下,直流耦合系数和交叉耦合系数,非别为直流耦合系数, 'meff ''

3、 neffFeff ;:-,-:(-)'D交流耦合系数二-neffK =第三步,有了上面的定义,我们接着引入直流耦合率和交叉耦合率的概念,分别为交叉耦合率t =IR |2 直流耦合率t =色4 -I R(0)I2Is I2I R(0)I2我们可以验证,直流耦合率和交叉耦合率之和等于一,这说明用直流耦合率和交叉耦合率来表示透射谱和反射谱是恰当的。有了上面一系列的叙述,我们就可以用matlab对长周期光纤光栅进行模仿,在这一过程中,我们对一些列的参数做如下设定,光纤有效折射率之差« =爵1 - Jff2 = 0.042 ,e e e光栅周期A=350Nm,周期数为 N=100,光

4、栅长度L=Ax N =35mm (这相当于公式中的z),由谐振波长计算公式我们可以得到谐振波长九D =1470nm 。我们分别取三个不同的折射率调制分别为 虱=1.2x10 , 虱=2.1妇0 , 虱f = 4.2x10 。以下是用matlab模拟出来的直流耦合率和交叉耦合率的图像14401460 1430 AW 怡 2Q 15 扣 1 涸 1530 1EOO 该长如n口耳口当居&冬器匍- q a nlilt £dit Vlh Xrisetr t ImIsnd w草Cl吕旧匙规方落戚 E3fall Esh 亳 Ya nv Ihiar tQ迁 Lap TlxIqr 2j

5、1;lp140J "20斗归旱!:;京长周期光纤光栅有效折射率的计算由上文的讨论我们可以看出,长周期光纤光栅谱线模拟中一个很重要的因素就是光纤光栅的有效折射率,即纤芯的有效折射率和包层的有效折射率。分别对着两个折射率进行计算。表征弱导光纤纤芯基模的色散方程为v JTEJl(V/)=vTb Kl(V晅 。其中Jo(V J-b) K°(V .b)Jo和Ji是0阶和1阶Bessel函数,Ko和Ki是0阶和1阶修正的Bessel函数。22co 22、22V = (2兀/人)aiWi -山 表示归一化频率;b=(neff 一山)/(山一表示模式的归一化有coco效折射率。neff表示

6、纤芯基模的有效折射率。下面我将计算neff的程序写再下面:首先是写出色散方程公式function f=myfwlc (bj V);f=V*sqrt (1-b) *besselj (lj Vt:sqrt (l-b)/bessel j (Oj V*sqrt (1-b) )-V*sqrt (b) *besselk (1 V*sqrt Cb)/besselk (0, V*sqrt (b);接着再代入参数求得纤芯有效折射率nCffclcclear all;fornat long;al=4150;53345;n2=L 52793;length=1200;for n=l:1:400length1(n)-length;V(n) = (2*pi/length)*al*sqft (”*2-n2“2);b(n)=f zero t tbo) inyfunc Cbo, V(n)j 0. 3) ;%why it must use length%br=f sero (b) inyfuneffnO (b length, V) 5 0. 3);ncoeff (rj =sqrt (b (n). *nl * 2+(l-b (n) ); *n2 * 2);length=lenethl (n) +1;enddlwrit e (' nc

温馨提示

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

评论

0/150

提交评论