数字信号处理实验报告11121.doc

《数字信号处理》实验报告 专 业 学 号 姓 名 实验一 利用FFT实现快速卷积 一、实验目的 1．加深理解FFT在实现数字滤波（或快速卷积）中的重要作用，更好的利用FFT进行数字信号处理。 2．掌握循环卷积和线性卷积两者之间的关系。 二、实验原理 用FFT来快速计算有限长度序列的线性卷积。 这种方法就是先将输入信号x(n)通过FFT变换为它的频谱采样值，然后再和FIR滤波器的频响采样值H(k)相乘，H(k)可事先存放在存储器中，最后再将乘积H(k)X(k)通过快速傅里叶变换（简称IFFT）还原为时域序列，即得到输出y(n)。 现以FFT求有限长序列的卷积及求有限长度序列与较长序列间的卷积为例来讨论FFT的快速卷积方法。 序列x(n)和h(n)的长差不多。设x(n)的长为N1，h(n)的长为N2，要求 用FFT完成这一卷积的具体步骤如下： ①为使两有限长序列的线性卷积可用其循环卷积代替而不发生混叠，必须选择循环卷积长度，若采用基2-FFT完成卷积运算，要求（为整数）。 ②用补零方法使x(n)和h(n)变成列长为N的序列。 ③用FFT计算x(n)和h(n)的N点离散傅里叶变换。 ④完成X(k)和H(k)乘积， ⑤用FFT计算 的离散傅里叶反变换得 三、主要实验仪器及材料 微型计算机、Matlab6.5教学版。 四、实验内容 1．数字滤波器的脉冲响应为。 输入序列x(n)可选下列几种情况 ①。 ②。 ③。 2．实验前，预先编制一个应用FFT实现数字滤波器的通用程序。 3．上机独立调试，并打印或记录实验结果。 4．将实验结果与预先笔算的结果比较，验证其正确性。 五、思考题 1． 六、实验报告要求 1.简述实验原理及目的。 2.记录调试运行情况及所遇问题的解决方法。 3. 给出实验结果，并对结果作出分析。 4．简要回答思考题。 实验二 数字滤波器的结构 实验目的 加深对数字滤波器分类与结构的了解； 明确数字滤波器的基本结构及其相互间的转换方法； 掌握用MATLAB进行数字滤波器各种结构相互间转换的子函数及程序编写方法。 实验原理 一个离散LSI系统可用系统函数来表示； 也可用差分方程来表示： 当至少有一个不为0时，则在有限z平面上存在极点，表示一个IIR数字滤波器；当全都为0时，系统不存在极点，表示一个FIR系统。 IIR数字滤波器的基本结构分为直接Ⅰ型、直接Ⅱ型、级联型和并联型。 FIR数字滤波器的基本结构分为横截型、级联型、并联型、、线性相位型和频率抽样型。 实验仪器 微型计算机、MATLAB 实验内容 已知一个IIR系统的系统函数为 将其从直接型转换为级联型和并联型结构，并画出各种结构的流程图。 已知一个FIR系统的系统函数为 将其从横截型转换为级联型结构，并画出各种结构的流程图。 实验程序 试验结果分析 实验三 用双线性变换法设计IIR数字滤波器 实验目的 加深对双线性变换法设计FIR数字滤波器的基本原理的理解； 掌握用双线性变换法设计数字低通、高通、带通滤波器的设计； 了解MATLAB有关双线性变换法设计的常用子函数。 实验原理 1. 设计思想 s平面与z平面之间满足以下映射关系： s平面的虚轴单值地映射于z平面的单位圆上，s平面的左半平面完全映射到z平面的单位圆内。双线性变换不存在混叠问题。 双线性变换时一种非线性变换，这种非线性引起的幅频特性畸变可通过预畸而得到校正。 2. 设计步骤 （(1)确定数字滤波器的性能指标：通带临界频率、阻带临界频率；通带内的最大衰减；阻带内的最小衰减； ?(2)确定相应的数字角频率； ? (3)计算经过预畸的相应模拟低通原型的频率； ? (4)根据计算模拟低通原型滤波器的阶数N，并求得低通原型的传递函数； ? (5)用上面的双线性变换公式代入，求出所设计的传递函数； ? (6)分析滤波器特性，检查其指标是否满足要求。 实验仪器 微型计算机、MATLAB 实验内容 (1) 采用双线性变换法设计一个巴特沃斯数字低通滤波器，要求：通带截止频率，阻带截止频率，滤波器的采样频率，并画出滤波器的幅频响应曲线和相频响应曲线以及零极点图。 实验程序 试验结果分析 思考题 使用双线性变换法时，数字频率与模拟频率有何关系？会带来什么影响？如何解决？ 实验四 用窗函数法设计FIR数字滤波器 实验目的 加深对窗函数法设计FIR数字滤波器的基本原理的理解； 学习用MATLAB语言的窗函数