Et=1*(t>0.1&t<0.5)表示方波，然后用fft() fftshift()函数画出傅里叶变换后的频域振幅

clc;
clear;

%时域光谱
delta=500;
t = 0:1/delta:1;
f0=20;
Et=1*(t>0.1&t<0.5);
figure(1);
plot(t,Et);

%频域光谱
Ew=fft(Et);
Ew=fftshift(Ew);%若需要作fftshift,则必先做fft，不能直接做fftshift或只做fft而不做fftshift,都会得到错误的结果
n=length(t);
x=(-n/2:n/2-1)*(delta/n);
figure(2);%振幅
plot(x,abs(Ew));%实函数傅里叶变换完之后一般是复函数，所以用abs()取模
figure(3);%实部
plot(x,abs(real(Ew)));
figure(4);%虚部
plot(x,abs(imag(Ew)));

%把变换后的画成极坐标图
figure(5);
theta=angle(Ew);
rho=abs(Ew);
polar(theta,rho)

结论：方波FFT的意义是，使用一系列不同频率的正弦函数表示了方波，各个正弦函数有不同的振幅，就是频域纵坐标的值，纵坐标越大，表示原信号中这个正弦函数频率分量越强