数字信号处理实验指导书技巧.docVIP

电 子 信 息 学 院 《数 字 信 号 处 理》 实验指导书 适用专业：通信工程、网络工程 贵 州 大 学 二OO 七年八月 前言 通过本课程的学习要求学生学习和掌握数字信号处理中的一些基本理论和处理方法：离散时间系统、离散傅立叶变换及其性质、IIR数字滤波器的原理与设计、FIR数字滤波器的原理与设计。为了使学生更好的理解和深刻掌握以上知识，培养学生对数字信号的分析和处理能力，设置了以下几个实验项目： 实验一 离散傅立叶变换的性质及应用 实验二 因果性数字系统的时域实现 实验三 数字巴特沃思滤波器的设计 实验四 用凯塞窗设计线性相位FIR滤波器 其中实验四为综合设计性实验。 学生应认真阅读《数字信号处理》教材中的与实验相关的章节内容，提前做好实验预习，做到每个实验前明确实验目的、掌握实验的基本内容及操作方法；在实验中正确使用实验设备，认真观察实验结果；实验后根据要求做好总结，上交实验报告。 目 录 实验一：离散傅立叶变换的性质及应用 1 实验二：因果性数字系统的时域实现 4 实验三：数字巴特沃思滤波器的设计 6 实验四：用凯塞窗设计线性相位FIR滤波器 8 实验报告的基本内容及要求 11 贵州大学实验报告 12 实验一：离散傅立叶变换的性质及应用 实验学时：2 实验类型：验证 实验要求：必修 一、实验目的 1、了解DFT的性质及其应用。 2、熟悉MATLAB编程特点。 二、实验内容 1、用三种不同的DFT程序实现一维数字信号的傅立叶变换。 2、利用DFT实现两序列的卷积运算，并研究DFT点数与混叠的关系。 3、研究高密度频谱与高分辨率频谱。 4、序列的内插和抽取时所对应的傅立叶变换。 三、实验原理 1、DFT变换 正变换： 反变换： 2、序列卷积 设序列的长度为N，序列的长度为M。则分别对两个序列作点的DFT得到和，则两序列的线性卷积等于。即时域卷积频域为相乘关系。 四、实验组织运行要求 集中授课的教学形式 五、实验条件 计算机一台； Matlab软件 六、实验步骤 1、用三种不同的DFT程序计算的傅立叶变换，并比较三种程序计算机的运行时间。 （1）编制用for循环语句的M函数文件dft1.m，用循环变量逐点计算； （2）编写用MATLAB矩阵运算的M函数文件dft2.m，完成下列矩阵运算： （3）调用FFT库函数，直接计算; （4）分别利用上述三种不同方式编写的DFT程序计算序列的傅立叶变换，并画出相应的幅频和相频特性，再比较各个程序的计算机运行时间。 2、利用DFT实现两序列的卷积运算，并研究DFT点数与混叠的关系。 给定，。利用圆周卷积计算线性卷积（快速卷积，频域）；计算不同DFT点数下的圆周卷积输出；并用函数stem(n,y)画出相应图形。选择不同的DFT点数进行对比，观察其混叠效应。 3、研究高密度频谱与高分辨率频谱 设有连续信号 以采样频率对该信号采样，分析下列三种情况的幅频特性。 （1）采集数据长度点，做点的DFT，并画出幅频特性。 （2）采集数据长度点，补零到256点的DFT，并画出幅频特性。 （3）采集数据长度点，做点的DFT，并画出幅频特性。 观察三幅不同频率特性图，分析和比较它们的特点以及形成的原因。 4、实现序列的内插和抽取所对应的傅立叶变换。 给定序列，做128点的傅立叶变换，并求 和 为整数 对应的傅立叶变换（128点）。比较这三个计算结果得到的幅频特性图，分析其差别产生的原因。选择不同的插值倍数和抽样倍数对比其幅频的变化。 七、思考题 1、直接计算点DFT与用FFT计算256点，理论上速度差别应有多少？ 2、什么是高密度频谱？什么是高分辨率频谱？实验步骤3中为区分这三个频率分量至少应采集几个样本？ 3、实验步骤4中，分别求和与原的关系，并说明为什么抽取后能分辨两个频率分量？有无任何条件的限制？ 八、实验报告 1、实验预习：前期实验和教材相关内容 2、列出本实验编写的所有文件及各项实验结果曲线，加注必要的说明。 3、对实验结果作理论计算，解释实验结果。 4、总结实验体会及实验中存在的问题。 实验二：因果性数字系统的时域实现 实验学时：2 实验类型：验证 实验要求：必修 一、实验目的 通过本实验的学习使学生熟悉和掌握数字信号处理系统中因果性数字系统的时域实现。 二、实验内容 1、编程实现FIR滤波 2、编程实现纯递归IIR滤波； 差分方程：