exam_5_1.m

来自「程序利用快速傅立叶变换及逆变换模拟了脉冲在光纤中传输时三阶色散对脉冲形状的影响」· M 代码 · 共 28 行

M=15;
tend=10;          %单位为皮秒ps
T=80;             %窗口，单位为皮秒ps
N=2^M;            %采样点总数
dt=T/N;
n=-N/2:(N/2-1);
t=n*dt;
A=zeros(size(t,2),1);                %初始化A
Tow=find(t<tend&t>-tend);
%1ps高斯脉冲
T0=1;                              %T0与脉冲半极大宽度成正比 1ps
z=2;                                   %传输距离km
A(Tow,1)=exp(-1/2.*((Tow-N/2-1)*dt./T0)'.^2);  %guass波
[AW,f]=cftbyfft(A,t,0);              %AW为A傅立叶变换后的  为频域中
w=2*pi.*f';                          %频率到圆频率的转换
[BW,t]=icftbyifft(N.*AW.*exp(i./6.*0.1*z.*w.^3),f,0);
plot(t,BW,'b')
hold on

T0=1;                              %T0与脉冲半极大宽度成正比 1ps
z=4;                                   %传输距离km
A(Tow,1)=exp(-1/2.*((Tow-N/2-1)*dt./T0)'.^2);  %guass波
[AW,f]=cftbyfft(A,t,0);              %AW为A傅立叶变换后的  为频域中
w=2*pi.*f';                          %频率到圆频率的转换
[BW,t]=icftbyifft(N.*AW.*exp(i./6.*0.1*z.*w.^3),f,0);
plot(t,BW,'r')
hold on