chirp analysis.m

来自「这是几个仿真利用半导体激光器的交叉相位调制进行全光波长转换的程序」· M 代码 · 共 42 行

te=2.2e-9;                    % carrier life     ns 
tp1=1.6e-12;                  % model1 photon life      ps
tp2=1.2e-12;                  % model2 photon life     ps
a=2.5e-20;                    % gain coefficient        m^2 
Nt=1.0e27;                    % transparency carrier density      m^(-3)
e=6.8e-23;                    % saturation coefficient            m^3
vg=0.75e8;                    % group velocity                    m/s
T=0.4;                        % confinement factor
q=1.6e-19;                    % C
L=250e-6;                     % m
d=0.2e-6;                     % m
w=2.0e-6;                     % m
bsp=0.004;                    % auto-spontaneous emission factor    
bc=5;                         % line-width enhancement factor
%%%%%%%%%%%%%%%%%%%%%%%%%%%输入脉冲函数%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
m=1000;
t=0:0.5e-9:50e-9;
pin1=5e-3;                    %mw
pin2=0;                       %mw
for k=1:1:m
   pin1=pin1+10*sin(pi/2*k)/k/pi*exp(i*k*pi*1e9/5*t);
   pin2=pin2+10*sin(pi/2*k)/k/pi*exp(-i*k*pi*1e9/5*t);
end
Pin=pin1+pin2;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
dt=1e-12;                     % ps
P(1)=0;
N(1)=0;
I=20e-3;                      % A
for J=2:50000
    N(J)=N(J-1)+(I/(q*L*w*d)-N(J-1)/te-vg*a*(N(J-1)-Nt)*(1-e*P1(J-1)-e*(P2(J-1)+Pin(J-1)))*P1(J-1)-vg*a*(N(J-1)-Nt)*(1-e*P1(J-1)-e*(P2(J-1)+Pin(J-1)))*(P2(J-1)+Pin(J-1)))*dt;
    P1(J)=P1(J-1)+(T*vg*a*(N(J-1)-Nt)*(1-e*P1(J-1)-e*(P2(J-1)+Pin(J-1)))*P1(J-1)-P1(J-1)/tp1+T*bsp*N(J-1)/te)*dt;
    P2(J)=P2(J-1)+(T*vg*a*(N(J-1)-Nt)*(1-e*P1(J-1)-e*(P2(J-1)+Pin(J-1)))*(P2(J-1)+Pin)-P2(J-1)/tp2+T*bsp*N(J-1)/te)*dt;
    Fi(J)=Fi(J-1)+(1/2*bc(T*vg*a*(N(J-1)-Nt)-1/tp1))*dt;
end
plot(t,N);
plot(t,P1);
plot(t,P2);
plot(t,Fi);