基于Matlab的FFT算法

clf;
fs=100;N=128;
n=0:N-1;t=n/fs;
x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t);
y=fft(x,N);  
mag=abs(y);
f=n*fs/N;


subplot(2,2,1),plot(f,mag);
ylabel('振幅');title('N=128');grid on;


subplot(2,2,2),plot(f(1:N/2),mag(1:N/2));
ylabel('振幅');title('N=128');grid on;
N=64;
subplot(2,2,3),plot(f,mag);
ylabel('振幅');title('N=1024');grid on;


subplot(2,2,4),plot(f(1:N/2),mag(1:N/2));

ylabel('振幅');title('N=1024');grid on;