《实验指导书 数字信号处理 2011电子、电信11275341》由会员分享，可在线阅读，更多相关《实验指导书 数字信号处理 2011电子、电信11275341（17页珍藏版）》请在金锄头文库上搜索。
1、数字信号处理实实 验验 指指 导导 书书王晓宁 编沈阳大学信息工程学院1目录目录实验一实验一 信号的频谱分析2实验二：LPF、HPF、BPF的幅频特性6实验三：DF设计8【附录A】 MATLAB简介12【附录B】 m文件简介142课程编号：课程编号：1127534111275341 课程类别：基础必修课课程类别：基础必修课 适用层次：本科适用层次：本科 适用专业：电子信息工程、电子科学与技术适用专业：电子信息工程、电子科学与技术 课程总学时：课程总学时：3232 适用学期适用学期: :第 5 学期 实验学时：实验学时：1212 开设实验项目数：开设实验项目数：3 3 撰写人：王晓宁撰写人：王晓
2、宁 审核人：周昕审核人：周昕 教学院长：王晖教学院长：王晖实验一实验一 信号的频谱分析一、实验目的与要求通过实验确实了解信号的时域、频域特征，对这两个信号空间有较真实的了解熟悉 MATLAB 的操作界面；熟悉常用 MATLAB 常用指令的用法。了解语音信号的采集、回放处理。加深对 DFT 算法原理和基本性质的理解；熟悉 FFT 算法原理和 FFT 子程序的应用；学习用 FFT 对连续信号和时域信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用 FFT。二、实验类型验证（4 学时）三、实验原理及说明 （1）利用计算机OS中的播放软件播放歌曲，所播放出来的音乐通过其声卡的输
3、出端（主机箱背面的绿色音频输出端口）输出模拟音频信号。（2）此模拟音频信号经黑色音频线被送达计算机的音频输入端（主机箱背面的粉色音频输出端口），并由计算机内部的声卡进行声音信号的采集。（3）在MATLAB软件环境中导入数字音乐信号，仔细观察声音信号时间域（频率域）波形（曲线）体现出的周期特征。离散傅立叶变换 DFT 的定义：将时域的采样变换成频域的周期性离散函数，频域的采样也可以变换成时域的周期性离散函数，这样的变换称为离散傅立叶变换，简称DFT。利用 DFT 将所采集的音频信号进行时、频域转换，并观察二者之间的联系。FFT是 DFT 的一种快速算法，将 DFT 的 N2步运算减少为步，极大的
4、提高了运算N2log2N的速度。旋转因子的变化规律。蝶形运算规律。3基 2FFT 算法。四、实验仪器序 号名 称主要用途 1计算机DSP 硬件支持2MATLAB6.5 版以上软 件DSP 软件支持五、实验内容和步骤【准备工作】1 用音频对录线，连接计算机主机箱的的音频“麦克输入”接口（粉色）和主机的音频输出口（绿色） ；2 将音乐文件“东邪西毒.mp3”或同学们自己准备的音乐文件复制到计算机中，并用系统中的音频播放器（如“千千静听”或“media player”等）播放（注意，为了实验方便，最好将播放方式设置成“循环”方式） ；（一）时频分析1 将实验程序 au_sample2.m 复制到装有
5、 Matlab 软件的磁盘的相应 Matlab 目录中的“work”文件夹中；2 双击桌面上 Matlab 软件的快捷图标，打开 Matlab 软件；3 在菜单栏中选择 “FileopenM-File”或是点击快捷按钮，在对话窗中选定本次试验程序 au_sample2.m，在 m 文件编辑器中打开；4 仔细阅读程序，了解程序执行过程；5 在程序代码第 14 行的“-”位置双击鼠标左键，建立断点（如图所示红点所在位置） ；6 在菜单栏中选择 “DebugRun”或点击快捷按钮、或按 F5 运行程序；（结果可以用虚拟示波器观察或用耳机听）实现语音信号的回放；47 单击采用步进方式运行程序，观察语音
6、采集的（双声道）各帧信号波形和频谱； 8 记录波形，写实验报告。9 修改程序（或调用程序 au_sample3.m） ，对语音信号做重低音加重（断点设置如下图） ，重新观察（听）音频信号的变化，并加以记录，将修改过程、所观察到的结果记录在实验报告上。（二）FFT1、 复习 DFT 的定义、性质和用 DFT 作谱分析的有关内容；2、 将实验程序 au_fft.m 复制到装有 Matlab 软件的磁盘的相应 Matlab 目录中的“work”文件夹中；3、 双击桌面上 Matlab 软件的快捷图标，打开 Matlab 软件；在菜单栏中选择 “FileopenM-File”或是点击快捷按钮，在对话窗
7、中选定本次试验程序 au_fft.m，在 m 文件编辑器中打开；54、 仔细阅读程序，了解程序执行过程；5、 在程序代码第 14 行的“-”位置双击鼠标左键，建立断点（如上图所示红点所在位置） ；6、 在菜单栏中选择 “DebugRun”或点击快捷按钮、或按 F5 运行程序；（结果可以用虚拟示波器观察或用耳机听）实现语音信号的回放；7、 单击采用步进方式运行程序，观察语音采集的（双声道）各帧信号波形和FFT 幅频特性； 8、 记录波形，写实验报告。六、实验数据处理与分析在 MAtlab 环境，同步观察输入语音信号波形及其频谱结果，记录，定性地分析时、频域之间的对应关系七、注意事项1. 实验说明
8、：实验说明：在设置断点时，要注意调试过程，即：（1）先在源程序中相应的位置断点，再运行；（2）若要修改程序，要点击，保存程序后再运行，否则新设置的内容不会直接体现在运行结果中；2. 音量的影响：音量的影响：实验时，输入及输出音量应视情况做适当调整，以避免采集信号因强度过大造成溢出或是强度太小造成信号幅度波动太小，不利于观察和记录。3. 主要函数说明：主要函数说明：close all 关闭所有（绘图）窗口6clear all 清除工作台所有变量clc 清屏input(，s); %在命令窗口中显示内容“y”，并将用户输入的内容以字符串（s）形式存入变量“x”for i=1:fr %固定次数循环，循
9、环变量i从1按步长1增至fr时循环结束endy=wavrecord(fr_l ,fs ,ch ,dtype);%将声卡采集的音频信号数据存入变量“y”，其中fr_l为音频信号采样点的个数（即帧长），fs为采样频率，ch为声道数，dtype为保存采样值时的数据类型（“double”指双精度类型）subplot(2,2,1) %在绘图窗口按两行两列划分子图，并指定其中的第一幅（即第一行、第一列）plot(y(:,1) %在制定图区将数组“y”的所有行（用“：”表示）、第一列以曲线形式绘出num2str(i) %将数字i转换为相应数字的字符串数据类型title(zh1) %在图区上方显示标题zh1f
10、loor(fr_l/2); %计算fr_l/2后向下取整Y=fft(y); %计算数组“y”的FFT并将结果赋给变量“Y”YY=abs(Y); %计算数组“Y”的绝对值或模并将结果赋给变量“YY”wavplay(y) %将数组“y”按音频播放h=bartlett(fr_l); %生成长度为fr_l的三角窗并赋给变量hout=zeros(fr*fr_l,ch); %生成fr*fr_l行ch列的零矩阵outY1(:,1)=h1.*Y(:,1); %将数h1和Y（：，1）的对应元素相乘（即求两信号乘积），结果写入Y1(:,1);y1=real(ifft(Y1); %求Y1的IFFT并将其实部赋给数组
11、y1wavwrite(out,fs,name1) %将二维数组out按照采样频率fs写入音频文件name1，保存格式为“*.wav”soundview(name1) %将音频文件name1播放出来并在播放窗口上显示该音频文件的波形八、预习与思考题1.音频信号在时域和频域的表现特征有什么区别？2.音频信号在时域和频域的表现特征有什么区别？3.对于不同的 N，幅频特性会相同吗？为什么？实验二：LPF、HPF、BPF 的幅频特性一、实验目的与要求了解各种不同类别数字滤波器在时频域的异同 熟悉 Matlab 软件的 Simulink 开发环境二、实验类型设计（4 学时）三、实验原理及说明利用 Matl
12、ab 软件提供的 Simulink 仿真功能对混合信源进行滤波并观察结果。其混合信源由正弦基波信号和二次谐波、四次谐波信号相叠加构成，对此混合信号源分别利用 LPF、HPF、BPF 进行滤波，借以观察 LPF、HPF、BPF 滤波器的滤波性能（主要是幅频特性） 。7四、实验仪器序 号名 称主要用途 1计算机DSP 硬件支持2MATLAB6.5 版以上软 件DSP 软件支持五、实验内容和步骤实验内容：利用 Matlab 的仿真功能观察 LPF 的时频域特性利用 Matlab 的仿真功能观察 HPF 的时频域特性利用 Matlab 的仿真功能观察 BPF 的时频域特性观察比较不同 DF 的滤波效果
13、实验步骤：本实验通过调用 Simulink 中的功能模块构成数字滤波器的仿真框图，在仿真过程中，可以双击各功能模块，随时改变参数，获得不同状态下的仿真结果。1. 构造以基波为主的原始信号，(其中线性加)()()(210ncSinnbSinnaSin权值 a、b、c 可以自己设置)，并计算出信源信号中各分量的数字角频率；2. 通过 Simulink 环境下的 Digital Filter Design（数字滤波器设计）模块进行参数设置，分别根据信源信号中各分量的数字角频率设计出 LPF、BPF、HPF，将基波分量、二次谐波分量和四次谐波分量分别滤出，并观察结果；Simulink 仿真图及滤波效果
14、图（例如：可以从上图看到经过数字滤波后分离出了四次谐波分量 scope1。）8六、实验数据处理与分析记录实验结果，定性地分析几种不同滤波器之间的异同。七、注意事项对滤波器的仿真要设定好仿真参数，尤其要注意原始信号数字角频率的计算和LPF、BPF、HPF 三种滤波器的截止角频率的设置。八、预习与思考题数字角频率与模拟角频率之间的关系是什么？以上几种不同类别的滤波器的作用分别是什么？实验三：DF 设计一、实验目的与要求熟悉设计 IIR 数字滤波器的原理与方法（双线性变换） ；掌握数字滤波器的计算机仿真方法；通过观察对实际信号的滤波作用，获得对数字滤波的感性认识。掌握用窗函数法设计 FIR 数字滤波
15、器的原理和方法；熟悉线性相位 FIR 数字滤波器特性；了解各种窗函数对滤波特性的影响。二、实验类型综合（4 学时）三、实验原理及说明无限冲击响应数字滤波器的基础理论；巴特沃斯低通滤波器（AF）设计方法；双线性变换法 DF 设计步骤。有限冲击响应数字滤波器的基础理论有限冲击响应数字滤波器具有严格的线性相位特性，其单位脉冲响应仅含有有限个非零值，是因果的有限长序列。目前关于有限冲击响应数字滤波器的设计方法主要有三种：窗函数法，频率取样法和切比雪夫等波纹逼近的最优化设计方法。各种方法都有其优缺点，使用时应根据需要进行选择。一般应用较多的是第一和第三种方法。因为窗函数法比较简单，可应用现成的窗函数公式，在技术指标要求不高时比较方便灵活的。最优化设计法必须借助计算机计算。切比雪夫等波纹线性相位滤波器的机助设计程序已经比较完善，采样瑞梅兹迭代算法，设计效率也很高，在应用中越来越占优势。窗函数设计法原理9使设计的滤波器逼近理想单位取样响应序列，可采取如下方法；一是用有限项和来逼近无限项和；二是将无限长的响应序列进行有限延时，从而由非因果的系统
