数字信号处理95668new.docVIP

数字信号处理 实验一 用FFT进行谱分析 一、实验目的 1.进一步加深对DFT算法原理和基本性质的理解（因为FFT只是DFT的一种快速算法，所以FFT的运算结果必然满足DFT的基本性质）。 2.熟悉FFT算法原理和FFT子程序的应用。 3.学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用FFT。 二、实验仪器 微型计算机 三、实验步骤 1.复习DFT的定义、性质和用DFT作谱分析的有关内容。 2.复习FFT算法原理与编程思想，并对照DIT-FFT运算流图和程序框图，读懂本实验提供的FFT子程序。 3.编制信号产生子程序，产生以下典型信号供谱分析用： （1-1） （1-2） （1-3） （1-4） （1-5） （1-6） 应当注意，如果给出的是连续信号，则首先要根据其最高频率确定采样速率fs 以及由频率分辨率选择采样点数N，然后对其进行软件采样（即计算），产生对应序列。对信号，频率分辨率的选择要以能分辨开其中的三个频率对应的谱线为准则。对周期序列，最好截取周期的整数倍进行谱分析，否则有可能产生较大的分析误差。请实验者根据DFT的隐含周期性思考这个问题。 4.编写主程序。 图1给出了主程序流程图及框图供参考。本实验提供FFT子程序和通用绘图子程序。 5.按实验内容要求，上机实验，并写出实验报告。 四、实验内容 1.对下列所给出的信号逐个进行谱分析。 下面给出针对各信号的FFT变换区间N以及对连续信号的采样速度fs，供实验时参考。 、、、、：N=64； ：fs=64（Hz），N=64 。 2.令=+，用FFT计算64点离散付里叶变换 =DFT[] （1-7） 并根据DFT的对称性,由求出=DFT[] 和=DFT[]，并与1中所得结果比较。 [提示：取N=64时，=，=－。] 3.令=+j，重复2。　 五、思考题 1.在N=8和N=16两种情况下，、的幅频特性会相同吗？为什么？ 2.如果周期信号的周期预先不知道，如何用FFT进行分析？ 六、实验报告要求 1.简述实验原理及目的。 2.结合实验中所得给定典型序列幅频特性曲线，与理论结合比较，并分析说明误差产生的原因以及用FFT作谱分析时有关参数的选择方法。 3.总结实验所得主要结论。 4.简要回答思考题。 附一：程序讲解 1、信号产生子程序： ①定义变量类型： 由于FFT中定义了四个变量（AR、AI、N、M）。因此，信号产生程序中必然有这四个变量。但变量名称不能与子程序相同，所以我们该利用(ar,ai,B,A),其中变量ar,ai为实型数组，B，A为整型变量。 ②B为以2为底A的对数。 ③选择信号种类： 应用键盘输入命令要求输入信号的编号（1~6），例如为1时进入信号1子程序。 ④产生信号：先使用循环语句，n从0循环到N-1。判断n的取值，并按信号要求分别给实部ar(n)、虚部ai(n)赋值。 ⑤调用FFT子程序。 计算|x(n)| 时，用一循环语句n=0, N-1,计算 |x(n)|= （1-8） 就完成了对信号1的FFT变换。信号2~5与信号1相同。特别要注意信号6，由于给出了fs=64Hz，根据T=1/fs ,可得到 Xa(k)=T·DFT[x(n)] （1-9） 所以 1/fs DFT[x(n)]=Xa(k) （1-10） 将中t换成n， 先计算出DFT[x(n)]， 再乘以1/fs ，就可以得到DFT[]。 主程序的框图见图1－1，FFT流程图见图1－2。 2、FFT子程序： FFT子程序中所需变量为AR、AI、M、N，其中AR为信号实部；AI为信号虚部；N= ，N表示变换序列长度。先确定变量类型，N、M为整型；AR、AI为实型。通过学习课本知道，要进行FFT变换必须先进行倒序，方法及框图如1－3。 把N/2→LH，LH→J，N-2→N 。I从1循环到N。 ①如果IJ，把LH赋给K； ②否则的话X(I)→X(J)→A(I)，T→(J)后执行①。再判断J是否小于K，如果是则J+K→J，继续循环；否则J-K→J，K