数字信号处理实验报告

实

验 报 告

姓名： 0907090117 0907090118 0907090122 孙飞姜万宝张战东

实验二：用FFT作谱分析 一、实验目的：

1、进一步加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快速算法，所以FFT的运算结果必然满足DFT的基本性质)。 2、熟悉FFT算法原理和FFT子程序的应用。

3、学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用FFT。 二、实验环境（软件、硬件及条件）：

Matlab6.5

三、实验内容及其步骤：

1、复习DFT的定义、性质和用DFT作谱分析的有关内容。 2 、复习FFT算法原理与编程思想，并对照DIT-FFT运算流图和程序框图，读懂本实验提供的FFT子程序。

3、编制信号产生子程序，产生以下典型信号供谱分析用： （1）、x1(n)?R4(n) clear

m=input('FFT点数='); n=m-1;

subplot(1,2,1); i=0:0.5:3;

x1=sign(sign(4-i)+1); stem(i,x1,'.'); xlabel('n'); ylabel('x1(n)'); title('x1(n)'); if(m==8) subplot(1,2,2); xb=fft(x1,m); i=0:n;

stem(i,abs(xb),'.'); xlabel('k'); ylabel('X1(k)');

title('x1(n)的8点FFT'); else

subplot(1,2,2); xa=fft(x1,m); i=0:n;

stem(i,abs(xa),'.');

xlabel('k'); ylabel('X1(k)');

title('x1(n)的16点FFT'); end

x1(n)的8点FFT：

x1(n)10.90.80.70.6)(n10.5x0.40.30.20.100123nx1(n)的16点FFT：

x1(n)的8点FFT7654321002468kX1(k)