yy.m

fft做谱分析各种连续信号 用MATLAB实现

%x6(n)=cos(8пt)+cos(16пt)+cos(20пt) fs=64kHz,N=16,32,64
m=input('FFT点数=');
n=0:(m-1);
subplot(2,2,1);
x6=cos(pi*n/8)+cos(pi*n/4)+cos(pi*n*5/16);
stem(n,x6,'.');
xlabel('n');
ylabel('x6(n)');
title('x6(n)=cos(8пt)+cos(16пt)+cos(20пt)的函数');
subplot(2,2,2);
xa=fft(x6,16);
i=0:15;
stem(i,abs(xa),'.');
xlabel('k');
ylabel('x6(n)');
title('x6(n)的16点FFT');
subplot(2,2,3);
xb=fft(x6,32);
i=0:31;
stem(i,abs(xb),'.');
xlabel('k');
ylabel('x6(n)');
title('x6(n)的32点FFT');
subplot(2,2,4);
xc=fft(x6,64);
i=0:63;
stem(i,abs(xc),'.');
xlabel('k');
ylabel('x6(n)');
title('x6(n)的64点FFT');