《数字信号处理实验实验二.docVIP

数字信号处理实验 实验2 用FFT作谱分析 专业： 学号： 姓名： 日期： 一、实验目的： 1、 进一步加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快 速算法，所以FPT的运算结果必然满足DFT的基本性质)。 2、 熟悉FFT算法原理和FFT子程序的应用。 3、 学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现 的分析误差及其原因，以便在实际中正确应用FFT。 二，实验原理 由于直接计算DFT需要很大的运算量，在一些要求对信号进行实时处理的场合或需要处理的信息量比较大的地方实不可取，为此设计了快速傅里叶变换，即FFT。它分为按时间抽选与按频率抽选两种方式，由于它仅仅是计算傅里叶的一种方法，所以其运算过程必然满足傅里叶变换基本性质。对于一个N点序列，一般设N为2的整数倍，不足补零。即N=2^L,则按频率抽选的FFT需要进行L运算，每级有N/2个蝶形运算，每个蝶形运算需要乘法1次，负数加法2次，总共需要N/2*L次乘法，N*L次加法，比直接计算，大大节省了运算量。如果输入为倒位序，则结果输出为自然序。 下面是8点FFT算法流图 从图中我们可以看到计算FFT一般步骤为先把序列进行倒位运算，然后在计算相应的蝶形结。此方法具有原址运算，所需存储空间小等优点。 三，实验内容及步骤 1），对2中所给出的信号逐个进行谱分析。下面给出各个信号的FFT变换区间N以及连续信号的采样频率，供实验时参考。 ，：N=8，16 上图是x5(n)的波形 上图是x5(n)的FFT后频谱图。我们可以看出时域是离散的周期的，频域也是离散周期的。 上图是x4(n)时域波形图 上图是x4(n)频域的FFT频谱图 2）令，用FFT计算8点和16点 从时域到频域的波形及频谱我们可以，时域相加，对应频域的相加，满足DFT线性性质。 3）令，重复（2）。 x4(n)与x5(n)是两个序列，通过组合变成一个序列，通过一次FFT运算即可求出这两个序列的傅里叶变换。在这里x4(n)的FFT变换对应着x(n)的FFT的共轭对称序列，jx5(n)的FFT对应着其共轭反对称序列，得到x(n)的傅里叶变换后，通过公式就可以求出x4(n)和x5(n)的FFT值来。 上述频域图形既是恢复后的x4(n)与x5(n)的FFT变换，可以看到结果是正确的。 三、思考题： 1、在N=8时，和的幅频特性会相同吗？为什么？N=16呢？ 它们的幅频特性相同，我们可以对x2(n)以8为周期进行周期延拓，然后再取其主值区间序列，则是x2(n),x3(n)既是x2(n)经周期延拓后得到的序列再进行圆周移位，向左移或向右移4个单位，再取其主值区间而得到x3(n),我们知道序列在时域的圆周移位会在频域产生相位的移位，但是幅度不会变化。而当N=16时，x3(n)则无法通过x2(n)的移位得到。所以其幅频特性就不会相同。 四，试验体会 通过本次试验，让我更加熟悉了Matlab的功能及用途，编程也更加得 心应手，同时我对FFT算法有了更深入的理解。