%创建一个幅值为5，20HZ的正弦函数

t=0:0.001:1;
Y=5*sin(2*pi*20*t);
MyFourierT(Y);
title('正弦的频谱');

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