pulse1.m

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

m=1000;
t=0:0.5e-11:50e-9
Pin1=5e-3;    %mw
Pin2=0;       %mw
for k=1:1:m
    Pin1=Pin1+10e-3*sin(pi/2*k)/k/pi*exp(i*k*pi*1e9/5*t);
    Pin2=Pin2+10e-3*sin(pi/2*k)/k/pi*exp(-i*k*pi*1e9/5*t);
end
Pin=Pin1+Pin2;
plot(t/1e-9,abs(Pin)*1000)