使用matlab计算长周期光纤光栅的有效折射率

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

2、的过程中，耦合模理论提供了精确有效的解，从而被我们广泛采纳。长周期光纤光栅的模式耦合主要是指纤芯基模和同向传播的包层模之间的耦合，由光纤的模式理论可以知道在包层中存在很多种包层模，为了方便讨论，我们在这一部分只讨论光纤光栅的纤芯基模和一个包层模之间的耦合。首先，我们定义一下谐振波长，就是在光纤有效折射率调制无穷小的情况下的谐振波长，可以用下面的公式来表示=An,式子中An表示光纤基模和包层模之间的有效折射Deffeff率之差，即An=n一n。effeff1eff2第二步，我们来定义一下，直流耦合系数和交叉耦合系数，非别为.An+6nAn直流耦合系数交流耦合系数G=兀(eff一ef)I丿D兀nK

3、=eff第三步，有了上面的定义，我们接着引入直流耦合率和交叉耦合率的概念，分别为|R|2直流耦合率t二S=IR|2(0)|s|2交叉耦合率t二严IR|2(0)我们可以验证，直流耦合率和交叉耦合率之和等于一，这说明用直流耦合率和交叉耦合率来表示透射谱和反射谱是恰当的。有了上面一系列的叙述，我们就可以用matlab对长周期光纤光栅进行模仿，在这一过程中，我们对一些列的参数做如下设定，光纤有效折射率之差An二n-n二0.042，effeff1eff2光栅周期,=350pm,周期数为N=100，光栅长度L=,N二35mm（这相当于公式中的z),由谐振波长计算公式我们可以得到谐振波长二1470nm。我们

4、分别取三个不同的D折射率调制分别为，n=1.2x105，,n=2.1x105，,n=4.2x105。effeffeff以下是用matlab模拟出来的直流耦合率和交叉耦合率的图像ZiltEditVjwInsertIodIeZLCLdM0口长周期光纤光栅有效折射率的计算由上文的讨论我们可以看出，长周期光纤光栅谱线模拟中一个很重要的因素就是光纤光栅的有效折射率，即纤芯的有效折射率和包层的有效折射率。分别对着两个折射率进行计算。7J(V1b)“7K(Vb)表征弱导光纤纤芯基模的色散方程为V1b1丿=Vb1。其中J(V1b)K(Vb)00J和J是0阶和1阶Bessel函数，K和K是0阶和1阶修正的Bes

5、sel函数。0101V=(2/)an2一n表示归一化频率；=(nco2-n2)/n2一n2表示模式的归一化有112eff212效折射率。nco表示纤芯基模的有效折射率。下面我将计算nco的程序写再下面：effeff首先是写出色散方程公式接着再代入参数求得纤芯有效折射率需clcclearall;formatlong:al=4150;nl=l.53345;nHl.52793;length=1200;forn=l:1:4001engthl(n.)=lengizh;V(nO=(2*pi/lengthjVa1*sqrt(n.12-ri22);b(n)=zer0-(bo)uiyfunc(bojV(nJ)0.3);SSvhyitmu.stuseLength%br=fzero(tOznyfuneffnO(bj.lengthyV).j0.jj;ncoeff(n)sqxt(b(n).*nl2+(l-b.(n)；.*n2岔;length=lengthl(ri)+1;enddliRwrite(znu