fs=200 #采样频率
t=LinRange(-0.5,0.5,fs)#时间轴长度及分度
τ=0.05  #高电平持续时间  
        #=矩形脉冲函数rectpuls(t,τ),其中T为高电平持续时间，分别给出三种不同宽度的脉冲，
         宽度分别为τ、0.5τ、1.5τ=#
x = rectpuls(t,τ) 
x1=rectpuls(t,0.5*τ)
x2=rectpuls(t,1.5*τ)
        #绘制时域波形
ax1=subplot(2, 1, 1)
ax2=subplot(2, 1, 2)
plot(ax1,t,x,t,x1,t,x2)
xlabel("Time (s)")
ylabel("Magnitude")
title("矩形脉冲的傅里叶变换演示")
        #对时域信号做快速傅里变换获得频域数据
y=fft(x)
y1=fft(x1)
y2=fft(x2)
         #由于在理论上我们获得的是双边谱，即正负频率都有。
         #因此需要对fft后的信号做移动将0频信号移到中心，可以采用fftshift()函数
yshift=fftshift(y)
y1shift=fftshift(y1)
y2shift=fftshift(y2)
         #确定绘制频谱图的横轴
fshift=LinRange(-1,1,length(yshift))
f1shift=LinRange(-1,1,length(y1shift))
f2shift=LinRange(-1,1,length(y2shift))
         #绘制频谱图
plot(ax2,fshift,abs.(yshift),f1shift,abs.(y1shift),f2shift,abs.(y2shift))
xlabel("Frequency (Hz)")
ylabel("Magnitude")