chirp.m

啁啾脉冲的源代码

clear all
format long e
% 定义光栅长度、有效折射率、折射率调制深度、光栅所分段数
L=0.015;
Neff=1.458; 
deltaN=0.000075;
SectionN=60;
Lamda=[1547.00:0.01:1553.00]*1e-9;
C=0.13;
% 计算波长范围
for k=1:601 % 1100 =(1549.5 -1550.6)/0.001 为总计算点数
  F=[1,0;0,1]; % 初始化传输矩阵
   for i=0:1:SectionN % 计算相应波长下的传输矩阵-SectionN/2
      LamdaD=1550e-9*(1+C*(-0.0075+i*L/SectionN)); %0.445*(0.025^2-3*(i*L/SectionN+0.008)^2);    %0.291     0.365   0.441  0.515   0.588
    %计算每段光栅本地中心波长KB
     delta=2*pi*Neff*(1/Lamda(k)-1/LamdaD)-4*pi*Neff/LamdaD.^2*1550e-9*C*(-0.0075+i*L/SectionN); %%-2*pi*Neff*0.216*C*(-0.0075+i*L/SectionN)*(1/Lamda(k))
     Kac=pi*deltaN/Lamda(k)*exp(-4.*log(2)*(i*L/SectionN-L/2).^2/0.04^2);
     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
%Y1=10*log(Y1);
hold on;
ylim([0 0.5]);
plot(Lamda*1e9,Y1,'-r');
%axis('tight')
xlabel('Wavelength \lambda (nm)')
ylabel('Reflectivity')