应用FFT实现信号频谱分析

电信类课程试验报告

2.2

clear all

N=100;

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=100') subplot(1,2,2)

k=0:length(magXK)-1; title('x(n)N=100') subplot(1,2,2)

k=0:length(magXK)-1; k=k*(2/100)

stem(k,magXK,'.'); xlabel('k');ylabel('|X(k)|');

2.3复合函数

clear all

N=100;

n=0:N-1;

xn=0.9*sin(2*pi*n/N)+0.6*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=100')

subplot(1,2,2)

k=0:length(magXK)-1;

title('x(n)N=100')

subplot(1,2,2)

k=0:length(magXK)-1;

k=k*(2/100)

stem(k,magXK,'.');

xlabel('k');ylabel('|X(k)|');

title('X(k)N=100');

3.3 N=64