传输矩阵法计算啁啾光栅反射谱和.m

传输矩阵法计算啁啾光栅反射谱 本节中给出传输矩阵法计算啁啾光栅反射谱和 时延的完整M at lab 程序实例, 本段代码在 M at lab513 和W indow s2000 环境中运行通过。

PI=3.142;
%定义光栅长度、有效折射率、啁啾系数=0.05nm/cm、折射率调制深度、光栅所分段数
L=0.04;Neff=1.45;C=10*1e-9;deltaN=0.0001;SectionN=50;
Lamda=[1549.5:0.001:1550.6]*1e-9; %计算波长范围
for k=1:1101 %1100=(1549.5-1550.6)/0.001为总计算点数
F=[1,0;0,1]; %初始化传输矩阵
for i=1:SectionN %计算相应波长下的传输矩阵
deltaN=0.00005*exp((-56*(-L/2+i*L/SectionN)^4)/(L)^4);%切趾
LamdaD=(1550-C*L/2+C*i*L/SectionN)*1e-9; %计算每段光栅本地中心波长
delta=2*PI*Neff*(1/Lamda(k)-1/LamdaD)+2*PI*deltaN/Lamda(k)+(4*PI*Neff)*C*(-L/2+i*L/SectionN)/(LamdaD^2); 
Kac=PI*deltaN/Lamda(k);
RB=(Kac^2-delta^2)^0.5;
F=F*[cosh(RB*L/SectionN)-j*(delta/RB)*sinh(RB*L/SectionN),-j*(Kac/RB)*sinh(RB*L/SectionN);j*(Kac/RB)*sinh(RB*L/SectionN),cosh(RB*L/SectionN)+j*(delta/RB)*sinh(RB*L/SectionN)];
end
Y1(k)=(abs(-F(3)/F(1)))^2;
Q(k)=phase(-F(3)/F(1));
end
Y(1)=Q(1);Y(2)=Q(2);Y(3)=Q(3);
for i=4:1101
    if(abs(Q(i-1)-Q(i)<=1)) %对相位不连续处求导的简单处理
    Y(i)=(1549.5+i*0.001)^2*1e-18/(2*PI*3e-4)*((Q(i-1)-Q(i))/(0.001e-9));
    else
        Y(i)=(1549.5+i*0.001)^2*1e-18/(2*PI*3e-4)*((Q(i-3)-Q(i-2))/(0.001e-9));
    end
end
hold on;
plot(Lamda*1e9,Y,'b'); %绘制时延曲线
plot(Lamda*1e9,Y1*500,'b'); %绘制放大500倍的反射谱图型