应用FFT实现信号频谱分析.doc

文档介绍：clear all
N=128;
n=0:N-1;
xn=.^n;
XK=fft(xn,N);
magXK=abs(XK);
phaXK=angle(XK);
subplot(1,2,1);
plot(n,xn);
xlabel('n');ylabel('x(n)');
title('x(n) N=128');
subplot(1,2,2);
k=0:length(magXK)-1;
stem(k,magXK,'.');
xlabel('k');ylabel('|X(k)|');
title('X(k) N=128');
clear all
N=128;
n=0:N-1;
xn=cos(2*pi*n/N);
XK=fft(xn,N);
magXK=abs(XK);
phaXK=angle(XK);
subplot(1,2,1);
plot(n,xn);
xlabel('n');ylabel('x(n)');
title('x(n) N=128');
subplot(1,2,2);
k=0:length(magXK)-1;
stem(k,magXK,'.');
xlabel('k');ylabel('|X(k)|');
title('X(k) N=128');
clear all
N=128;
n=0:N-1;
xn=*sin(2*pi*n/N)+*sin(2*pi*n/N/3);
XK=fft(xn,N);
magXK=abs(XK);
phaXK=angle(XK);
subplot(1,2,1);
plot(n,xn);
xlabel('n');ylabel('x(n)');
title('x(n) N=128');
subplot(1,2,2);
k=0:length(magXK)-1;
stem(k,magXK,'.');
xlabel('k');ylabel('|X(k)|');
title('X(k) N=128');
当n=64时
clear all
N=64;
n=0:N-1;
xn=.^n;
XK=fft(xn,N);
magXK=abs(XK);
phaXK=angle(XK);
subplot(1,2,1);
plot(n,xn);
xlabel('n');ylabel('x(n)');
title('x(n) N=64');
subplot(1,2,2);
k=0:length(magXK)-1;
stem(k,magXK,'.');
xlabel('k');ylabel('|X(k)|');
title('X(k) N=64');
当n=128时
clear all
N=128;
n=0:N-1;
xn=.^n;
XK=fft(xn,N);
magXK=abs(XK);
phaXK=angle(XK);
subplot(1,2,1);
plot(n,xn);
xlabel('n');ylabel('x(n)');
title('x(n) N=128');
subplot(1,2,2);
k=0:length(magXK)-1;
stem(k,magXK,'.');
xlabel('k');ylabel('|X(k)|');
title('X(k) N=128');
clear all
N=256;
n=0:N-1;
xn=.^n;
XK=fft(xn,N);
magXK=abs(XK);
phaXK=angle(XK);
subplot(1,2,1);
plot(n,xn);
xlabel('n');ylabel('x(n)');
title('x(n) N=256');
subplot(1,2,2);
k=0:length(magXK)-1;
stem(k,magXK,'.');
xlabel('k');ylabel('|X(k)|');
title('X(k) N=256');
clear all
N=128;
n=0:N-1;
xn=cos(2*pi*n/N);
XK=fft(xn,2*N);
magXK=abs(XK);
phaXK=angle(XK);
subplot(1,2,1);
plot(n,xn);
xlabel('n');ylabel('x(n)');
title('x(n) N=128');
subplot(1,2,2);
k=0:length(magXK)-1;
stem(k,magXK,'.');
xlabel('k')