创建一个幅值为10，脉冲宽度为0.5，时长为0-6秒的矩形脉冲

t=0:0.001:6;
Y=(t>0&t<2.75)*0 + (t>=2.75&t<3.25)*10 + (t>=3.25&t<6)*0;
MyFourierT(Y);
title('矩形脉冲的频谱');

function MyFourierT(x)
X=abs(fft(x));
f=(0:length(X)-1)*1000/length(X);
plot(f,X);
axis([0 200,-1,1100]);
xlabel('频率/Hz');
ylabel('振幅');
grid on
end