实验三--用FFT对信号作频谱分析.doc

文档介绍：实验三--用FFT对信号作频谱分析
实验三：用FFT对信号作频谱分析
1．实验目的
学****用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析
误差及其原因，以便正确应用FFT。
2. 实验原理('(1b)16µãDFT[x_1(n)]');
xlabel('¦Ø/¦Ð ');ylabel('·ù¶È');
axis([0,15,0,*max(abs(X1k16))]);figure(2)
subplot(2,2,1);mstem(X2k8);title('(2a)8µãDFT[x_2(n)]');xlabel('¦Ø/¦Ð');ylabel('·ù¶È');
axis([0,7,0,*max(abs(X2k8))])
subplot(2,2,2);mstem(X2k16);
title('(2b)16µãDFT[x_2(n)]');xlabel('¦Ø/¦Ð ');ylabel('·ù¶È');axis([0,15,0,*max(abs(X2k16))])
subplot(2,2,3);mstem(X3k8);title('(3a)8µãDFT[x_3(n)]');xlabel('¦Ø/¦Ð ');ylabel('·ù¶È');
axis([0,7,0,*max(abs(X3k8))])
subplot(2,2,4);mstem(X3k16);title('(3b)16µãDFT[x_3(n)]');
xlabel('¦Ø/¦Ð ');ylabel('·ù¶È');axis([0,15,0,*max(abs(X3k16))])
N=8;n=0:N-1;x4n=cos(pi*n/4);x5n=cos(pi*n/4)+cos(pi*n/8);X4k8=fft(x4n);X5k8=fft(x5n);
N=16;n=0:N-1;
x4n=cos(pi*n/4);x5n=cos(pi*n/4)+cos(pi*n/8);X4k16=fft(x4n);X5k16=fft(x5n);figure(3)
subplot(2,2,1);mstem(X4k8);title('(4a)8µãDFT[x_4(n)]');xlabel('¦Ø/¦Ð ');ylabel('·ù¶È');
axis([0,7,0,*max(abs(X4k8))])
subplot(2,2,3);mstem(X4k16);title('(4b)16µãDFT[x_4(n)]');
xlabel('¦Ø/¦Ð');ylabel('·ù¶È');axis([0,15,0,*max(abs(X4k16))])
subplot(2,2,2);mstem(X5k8);
title('(5a)8µãDFT[x_5(n)]');xlabel('¦Ø/¦Ð ');ylabel('·ù¶È');axis([0,7,0,*max(abs(X5k8))])
subplot(2,2,4);mstem(X5k16);title('(5b)16µãDFT[x_5(n)]');xlabel('¦Ø/¦Ð ');ylabel('·ù¶È');
axis([0,15,0,*max(abs(X5k16))])
figure(4)
Fs=64;T=1/Fs;N=16;n=0:N-1;
x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T);
X6k16=fft(x6nT);
Tp=N*T;
F=1/Tp;
k=-N/2:N/2-1;
fk=k*F;
subplot(3,1,1);
stem(fk,abs(X6k16),'.');
box on
title('(6a)16|点DFT[x_6(nT)]');
xlabel('f(Hz)');
ylabel('幅度');
axis([-N*F/2-1,N*F/2-1,0,*max(abs(X6k16))])
N=32;n=0:N-1;
x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T);
X6k32=fft(x6nT);
X6k32=fftshift(X6k32);
Tp=N*T;F=1/Tp;
k=-N/2:N/2-1;fk=k*F;
subplot(3,1,2);
stem(fk,abs(X6k32),'.');
box on
title('(6b)32点|DFT[x_6(nT)]');
xlabel('f(Hz)');
ylabel('幅度');
axi