哈工大数字信号处理实验介绍.docVIP

Harbin Institute of Technology 数字信号处理 学生姓名： 江世凯 学号： 1122110307 班级： 1221103 专业： 电子科学与技术 任课教师： 李杨 所 在 单 位： 电子工程系 2014年11月 实验一、用FFT作谱分析 一、实验目的 (1) 进一步加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快速算法， 所以FFT的运算结果必然满足DFT的基本性质)。 (2) 熟悉FFT算法原理和FFT子程序的应用。 (3) 学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用FFT。 二、实验内容 (1) 编制信号产生子程序， 产生以下典型信号供谱分析用： (2) 画出1 中所给出的信号，并逐个进行谱分析。 下面给出针对各信号的FFT变换区间N以及对连续信号x6(t)的采样频率fs， 供实验时参考。 x1(n), x2(n), x3(n), x4(n), x5(n): N=8, 16 x6(t): fs=64(Hz), N=16, 32, 64（n=0:1:69） (3) 令x(n)=x4(n)+x5(n)， 用FFT计算 8 点和 16 点离散傅里叶变换， X(k)=DFT［x(n)］ (4) 令x(n)=x4(n)+jx5(n)， 重复(2)。 三、程序框图 图1.实验程序框图 四、实验过程 (1) 复习DFT的定义、 性质和用DFT作谱分析的有关内容。 (2) 复习FFT算法原理与编程思想， 并对照DIT-FFT运算流图和程序框图， 读懂本实验提供的FFT子程序。 (3) 编制信号产生子程序， 产生以下典型信号供谱分析用： (4) 编写主程序。 (5) 按实验内容要求， 上机实验， 并写出实验报告。 五、实验结果及分析 1.实验结果 （1）对x1(n)进行FFT变换（N=8和N=16） 图2 x1(n)进行FFT变换结果 （2）对x2(n)进行FFT变换（N=8和N=16） 图3 x2(n)FFT变换结果 （3）对x3(n)进行FFT变换（N=8和N=16） 图4 x3(n) FFT变换结果 （4）对x4(n)进行FFT变换（N=8和N=16） 图5 x4(n) FFT变换结果 （5）对x5(n)进行FFT变换（N=8和N=16） 图6 x5(n) FFT变换结果 （6）对x6(n)进行FFT变换（N=16、N=32和N=64） 图7 x6(n) FFT变换结果 （7）对x(n)=x4(n)+x5(n)进行FFT变换（N=8和N=16） 图8 x4(n)+x5(n) FFT变换结果 （8）对x(n)=x4(n)+jx5(n)进行FFT变换（N=8和N=16） 图9 x4(n)+jx5(n)FFT变换结果 （9）GUI设计结果 2.分析： (1) 在N=8时，x2(n)和x3(n)的幅频特性会相同吗? 为什么? N=16呢? 答：在N=8的时候，x3(n)只是x2(n)以它的长度8为周期，将其延拓成周期 序列然后加以移位，最后截取主值区间（n=0到8）的序列值，因此x3(n)是x2(n) 进行循环位移后的结果。由于序列的循环位移性质，因此，循环位移只影响到DFT后的相频特性，而不影响幅频特性，因此x2(n)和x3(n)的幅频特性会相同。若N=16，此时，x2(n)和x3(n)做DFT即为它们分别进行末尾补零后再进行的DFT，则此时两个经过补零以后的序列就不满足循环位移的性质，因 此x2(n)和x3(n)的幅频特性就会发生变化。 (2) 如果周期信号的周期预先不知道， 如何用FFT进行谱分析 答：在用FFT对信号进行频谱分析时，首先要对信号进行采样，使之变成离散信号，然后就可按照前面的方法用FFT对离散信号进行频谱分析。按照采样定理，采样频率fs应当大于信号最高频率的2倍。若预先不知道周期信号的周期，应先适当截取M点进行FFT，再将截取的长度扩大1倍重新截取，比较