广西师范大学电子工程学院ouqibiao中国矿业大学subaishun1第1 章归于三层模型的长周期光纤光栅透射谱的仿真1.1 芯层导模和包层模式有效折射率的确定一、芯层导模有效折射率的确定芯层导模LP01（HE11）的色散方程如下[1]：)()()1()1(10101bVKbVKbVbVJbVJbV=−−−(1-1)式中，22211/2nnaV−⋅=λπ是光纤在波长λ下的V 参数，2221222nnnnbcoeff−−=为光纤导模的归一化有效折射率，coeffn表示纤芯导模的有效折射率。给定光纤参数和光波长，应用二分法即可实现对芯层导模的有效折射率的求解。求解芯层导模的有效折射率程序如下（程序共有三个文件，这三个文件放在同一文件夹中，运行 CoreNeffSolving.m）：文件 1：求解芯层模式有效折射率的程序CoreNeffSolving.m%---------------------------------------------------------------------------------------%%文件名：CoreNeffSolving.m%%描述：模拟一段范围内的长周期光纤光栅的透射谱图%输入：n1,n2,n3分别为光栅芯层、包层和环境层的折射率%...a1,a2分别为光纤芯层半径和包层半径%%date:2011.07.24%%ouqibiao%%结果：%%-----------------------------------------------------------------------------------------%clearallformatlongn1=1.4681;n2=1.4628;n3=1.0;a1=4.15e-6;a2=62.5e-6;epsilon=1e-10;%计算误差wl=1550e-9;neffcore=CoreDichotamy(n2,n1,wl,n1,n2,a1,epsilon);文件 2：二分法的程序 CoreDichotamy.mfunctionc=CoreDichotamy(a,b,wl,n1,n2,a1,err)fun=@CoreFormula;f1=feval(fun,wl,a,n1,n2,a1);f2=feval(fun,wl,b,n1,n2,a1);iff1.*f2>0disp('Note:f1*f2>0')elsemax=real(1+round((log(a-b)-log(err))/log(2)));fork=1:max广西师范大学电子工程学院ouqibiao中国矿业大学subaishun2c=(a+b)./2;f3=feval(fun,wl,c,n1,n2,a1);iff3==0a=c;b=c;elseiff2.*f3>0b=c;f2=f3;elsea=c;f1=f3;endifabs(b-a)
