用FFT对信号作频谱分析.doc

实验三：用FFT对信号作频谱分析
1．实验目的
学****用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析
误差及其原因，以便正确应用FFT。
2. 实验原理
用FFT对信号作频谱分析是学****数字信号处理的重要内容。经DFT的幅频特性图
title('(1b)16点DFT[x_1(n)]');xlabel('ω/π');ylabel('幅度');
axis([0,2,0,*max(abs(X1k16))])
figure(2)
n=(0:length(X2k8)-1)/((length(X2k8))/2)
subplot(3,2,3);stem(n,abs(X2k8),'.'); %绘制8点DFT的幅频特性图
title('(2a) 8点DFT[x_2(n)]');xlabel('ω/π');ylabel('幅度');
axis([0,2,0,*max(abs(X2k8))])
n=(0:length(X2k16)-1)/((length(X2k16))/2)
subplot(3,2,4);stem(n,abs(X2k16),'.'); %绘制16点DFT的幅频特性图
title('(2b)16点DFT[x_2(n)]');xlabel('ω/π');ylabel('幅度');
axis([0,2,0,*max(abs(X2k16))])
n=(0:length(X3k8)-1)/((length(X3k8))/2)
subplot(3,2,5);stem(n,abs(X3k8),'.'); %绘制8点DFT的幅频特性图
title('(3a) 8点DFT[x_3(n)]');xlabel('ω/π');ylabel('幅度');
axis([0,2,0,*max(abs(X3k8))])
n=(0:length(X3k16)-1)/((length(X3k16))/2)
subplot(3,2,6);stem(n,abs(X3k16),'.'); %绘制16点DFT的幅频特性图
title('(3b)16点DFT[x_3(n)]');xlabel('ω/π');ylabel('幅度');
axis([0,2,0,*max(abs(X3k16))])
%实验内容(2) 周期序列谱分析==================================
N=8;n=0:N-1; %FFT的变换区间N=8
x4n=cos(pi*n/4);
x5n=cos(pi*n/4)+cos(pi*n/8);
X4k8=fft(x4n); %计算x4n的8点DFT
X5k8=fft(x5n); %计算x5n的8点DFT
N=16;n=0:N-1; %FFT的变换区间N=16
x4n=cos(pi*n/4);
x5n=cos(pi*n/4)+cos(pi*n/8);
X4k16=fft(x4n); %计算x4n的16点DFT
X5k16=fft(x5n); %计算x5n的16点DFT
figure(3)
n=(0:length(X4k8)-1)/((length(X4k8))/2)