实验一用FFT进行信号谱分析报告
- 格式:doc
- 大小:680.50 KB
- 文档页数:8
实验一用FFT进行信号谱分析
一.实验目的:
A.进一步加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快速算法,
所以FFT的运算结果必然使DFT得基本性质)。
B.熟悉FFT算法原理和FFT子程序的应用。
C.学习用FFT对连续信号和时域离散信号进行谱分析的算法,了解可能出现的分析误
差及其原因,以便在实际中正确应用FFT。
二.实验原理:
DFT是离散傅里叶级数(DFS)在时域和频域的两个主值序列构成的变换对:
用DFT来分析有限长序列的频谱特性。
FFT是DFT的快速算法,常用的两种快速算法有按时间抽取的FFT算法和按频率抽取的FFT算法,其中按时间抽取的FFT算法如下式:
按频率抽取的FFT算法如下:
三.实验数据
1、及其8点和16点DFT
2、及其8点和16点DFT
3、及其8点和16点DFT
4、的8点和16点波形及其DFT
5、的8点和16点波形及其DFT
6、的16点、32点和64点采样序列波形及其DFT
选7时,①计算并图示和
及其DFT。
②程序自动计算并绘图验证DFT的共轭对称性。
当N=16时,,。
即为的共轭对称分量,而是的共轭反对称分量。
根据DFT的共轭对称性,应有以下结果:
的8点和16点波形及其DFT
绘出和的模。
它们正是图中16点的和。
选8时,①计算并图示和
及其DFT。
②程序自动计算并绘图验证DFT的共轭对称性的第二种形式:如果,,则,。
其中,。
的8点和16点DFT
程序计算结果如下:
及
,正好与图中
的
16点
及
相同。
及
()()()j
k X IDFT n x n x o i /][885==,正好与图中16点
的及
相同。
四. 思考题
(1)当N=8时,
和
的幅频特性会相同吗?为什么?N=16 呢?
(2)对于周期序列,如果周期不知道,如何用FFT 进行谱分析? 答:(1)当N=8时,
和的幅频特性会相同. 当N=16时,和
的幅频特性会不相同。
(2)、如果
的周期预先不知道,可截取M 点进行DFT ,即
再将截取长度扩大1倍,截取
比较和,如果两者的主谱差别满足分析误差要求,则以或近似表示的频谱,否则,继续将截取长度加倍,直至前后两次分析所得主谱频率差别满足误差要求。
设最后截取长度为,则表示点的谱线强度。
通过实验,我知道了用FFT对信号作频谱分析是学习数字信号处理的重要容。
经常需要进行谱分析的信号是模拟信号和时域离散信号。
对信号进行谱分析的重要问题是频谱分辨率D和分析误差。
频谱分辨率直接和FFT的变换区间N有关,因为FFT能够实现的频率分辨率是2л/N≤D。
可以根据此式选择FFT的变换区间N。
误差主要来自于用FFT作频谱分析时,得到的是离散谱,而信号(周期信号除外)是连续谱,只有当N较大时,离散谱的包络才能逼近于连续谱,因此N要适当选择大一些。
周期信号的频谱是离散谱,只有用整数倍周期的长度作FFT,得到的离散谱才能代表周期信号的频谱。
如果不知道信号周期,可以尽量选择信号的观察时间长一些。
对模拟信号进行频谱分析时,首先要按照采样定理将其变成时域离散信号。
如果是模拟周期信号,也应该选取整数倍周期的长度,经过采样后形成周期序列,按照周期序列的谱分析进行。