数字信号处理实验报告案例.docVIP

实 验 报 告 课程名称： 数字信号处理实验 实验项目： 1．用FFT对信号作频谱分析 2．FIR数字滤波器设计与软件实现 3．IIR数字滤波器设计与软件实现 实验人员： *** 实验班级： 13电子信息工程2班 总 份 数： 1份 指导教师： 朱小松 电子与信息 学院 实验楼409 实验室 二〇一五 年 十一 月 实验三 用FFT对信号作频谱分析 一 实验目的 加深对DFT算法原理和基本性质的理解； 熟悉FFT算法原理和FFT子程序的应用； 学习用FFT对连续信号和时域信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用FFT。 二 实验设备 计算机，CCS 5.3 版软件，实验箱，DSP仿真器，短接块，导线。 三 基本原理 （1）离散傅立叶变换DFT的定义：将时域的采样变换成频域的周期性离散函数，频域的采样也可以变换成时域的周期性离散函数，这样的变换称为离散傅立叶变换，简称DFT。 （2）FFT是DFT的一种快速算法，将DFT的N2步运算减少为（N/2）log2N步，极大的提高了运算的速度。快速傅氏变换（FFT）是离散傅氏变换的快速算法，它是根据离散傅氏变换的奇、偶、虚、实等特性，对离散傅立叶变换的算法进行改进获得的。 ( 3)旋转因子的变化规律：N点DIT―FFT运算流图中，每级都有N/2个蝶形。每个蝶形都要乘以因子WpN，称其为旋转因子， （4）蝶形运算规律：DFT运算称为蝶形运算,而整个FFT就是由若干级迭代的蝶形运算组成,而且这种算法采用原位运算,故只需N个存储单元第一列只有一种类型的蝶形运算，系数是 ,以后每列的蝶形类型，比前一列增加一倍，到第是N/2个蝶形类型，系数是，共N/2个。由后向前每推进一列，则用上述系数中偶数序号的那一半，例如第列的系数则为参加蝶形运算的两个数据点的间距，则是最末一级最大，其值为N/2,向前每推进一列，间距减少一半。 四 实验步骤 1. 复习DFT的定义、性质和用DFT作谱分析的有关内容； 2. 复习FFT算法原理与编程思想，并对照DIT-FFT运算流程图和程序框图，了解本实验提供的FFT子程序； 3. 运行CCS软件，对样例程序进行跟踪，分析结果；记录必要的参数。 （1） 启动CCS 5.3，在D盘新建一个文件夹并命名为lz，单击Project/Import CCS5，选择select search –directory,选择D盘，在选择DSP实验及材料，最后点击DSPprogram 里面的dhlab39 FFT ,单击确定。注意，一定要选择copy projects into才算完成导入程序。打开“dhlab.39 FFT”工程文件；选择SRC文件夹里面的lab-FFT.； （2）编译程序，出现错误，则更改查找路径，再次编译程序则不会出现错误； （3）调试程序，按下debug ,在按下Resume再次运行； （4）查看图像，点击Tools里面的graph/single Time,并设置参数。结果如下； 五 实验收获 用FFT作谱分析时有关参数的选择方法。 答：如果用FFT对模拟信号进行谱分析，首先要把模拟信号转换成数字信号，转换时要求知道模拟信号的最高截止频率，以便选择满足采样定理的采样频率。用FFT作谱分析时，要求做FFT的点数服从２的整数幂。如果要进行谱分析的模拟信号是周期信号，最好选择观测时间是信号周期的整数倍。 实验四 IIR数字滤波器设计与软件实现 一、实验目的 （1） 熟悉用双线性变换法设计IIR数字滤波器的原理与方法； （2） 学会调用MATLAB信号处理工具箱中滤波器设计函数（或滤波器设计分析工具FDATool）设计各种IIR数字滤波器，学会根据滤波需求确定滤波器指标参数。 （3） 掌握IIR数字滤波器的MATLAB实现方法。 （4） 通过观察滤波器输入、输出信号的时域波形及其频谱，建立数字滤波的概念。 二、实验设备 计算机，CCS 2.0 版软件，DSP仿真器，实验箱，短接块，导线。 三、实验原理 设计IIR数字滤波器一般采用间接法（脉冲响应不变法和双线性变换法），应用最广泛的是双线性变换法。基本设计过程是： ① 将给定的数字滤波器的指标转换成过渡模拟滤波器的指标； ② 设计过渡模拟滤波器；③ 将过渡模拟滤波器系统函数转换成数字滤波器的系统函数。MATLAB信号处理工具箱中的各种IIR数字滤波器设计函数都是采用双线性变换法。第6章介绍的滤波器设计函数butter、cheby1 、cheby2 和ellip可以分别被调用来直接设计巴特沃斯、切比雪夫1、