数字信号处理DSP 课程设计报告new.docVIP

语音信号滤波去噪 ——使用脉冲响应不变法设计的巴特沃斯滤波 摘要 本课程设计主要运用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用一个满足指标的，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析，根据结果和学过的理论得出合理的结论。 MATLAB 1.课程设计目的和要求 1.1 课程设计目的 《数字信号处理》课程设计是在学生完成数字信号处理和MATLAB的结合后的基本实验以后开设的。本课程设计的目的是为了让学生综合数字信号处理和MATLAB并实现一个较为完整的小型滤波系统。这一点与验证性的基本实验有本质性的区别。开设课程设计环节的主要目的是通过系统设计、软件仿真、程序安排与调试、写实习报告等步骤，使学生初步掌握工程设计的具体步骤和方法，提高分析问题和解决问题的能力，提高实际应用水平。 1.2 课程设计的要求 （1）滤波器指标必须符合工程设计。 （2）设计完后应检查其频率响应曲线是否满足指标。 （3）处理结果和分析结论应该一致，而且应符合理论。 （4）独立完成课程设计并按要求编写课程设计报告。 2 .设计原理 用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用一个满足指标的，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析）∞）； 结构上存在着输出到输入的反馈，也就是结构上是递归型的。 IIR滤波器采用递归型结构，即结构上带有反馈环路。可以直接型、型、级联型、并联型四种，都具有反馈回路。IIR数字滤波器在设计上可以借助成熟的模拟滤波器的成果，巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Cauer)滤波器、贝塞尔(Bessel)滤波器等，这些典型的滤波器各有特点。有现成的设计数据或图表可查在设计一个IIR数字滤波器时，我们根据指标先写出模拟滤波器的公式，然后通过一定的变换，将模拟滤波器的公式转换成数字滤波器的公式。 式中，N为正整数，称之为滤波器的阶数，N越大，通带和阻带的近似性越好，过渡带也越陡峭。称为截止频率。当时，有 即 ， 所以又称为巴特沃思低通滤波器的3分贝带宽。 巴特沃思低通滤波器的特点如下： 当时，，即在处衰减。 当时，，，或 ，为通带最大衰减。即不管N为多少，所有的特性曲线都通过－3dB,或说衰减3Db,这就是3dB不变性。 在的通带内有最大平坦的幅度特性，即N阶巴特沃思低 通滤波器在处， 前的阶导数为0，因而巴特沃思滤波器又称为最平幅度特性滤波器。随着由0变到，单调减小，N越大，减小的越慢，也是通内特性越平坦。 当，即在过渡带及阻带中，也随着增加而单调减小， 但是，故比通带内衰减的速度要快的多，N越大，衰减速度越大。当，即频率为阻带截止频率时，衰减为，为阻带最小衰减。 图2.2.1 巴特沃思滤波器幅度特性及其与N的关系 图2.2.2 贝塞尔，契比雪夫，巴特沃思滤波器特性 2.3 脉冲响应不变法 冲激响应不变法是使数字滤波器的单位冲激响应序列模仿模拟滤波器的单位冲激响应。将模拟滤波器的单位冲激响应加以等间隔抽样，使正好等于的抽样值，即满足 其中T是抽样周期。 如果令是的拉普拉斯变换，为的Z变换，根据抽样序列的z变换与模拟信号的拉普拉斯变换关系，得 则可看出，冲激响应不变法将模拟滤波器的s平面变换成数字滤波器的z平面，这个从s到z的变换正是从s平面变换到z平面的变换关系。如图2.3.1所示，s平面上每一条宽度为的横条都将重叠地映射到整个z平面上，而每一条横条的左半边映射到z平面单位圆以内，右半边映射到z平面单位圆以外，而s平面上虚轴（轴）映射到z平面单位圆上，虚轴上每一段长为的线段都映射到z平面单位圆上一周。由于s平面每一条横轴都要重叠地映射到z平面上，这正好反映了是和的周期延拓之间有变换关系,故冲激响应不变法并不相当于从s平面到z平面的简单代数映射关系。 图2.3.1 冲激响应不变法映射关系 设计步骤 3.1设计流程图 语音信号滤波去噪——使用脉冲不变响应法设计的巴特沃斯滤波器的设计流程如图3.1.1所示： 图3.1.1 脉冲响应不变法巴特沃斯滤波器对语音信号去噪流程图 3.2语音信号的采集 点击windows系统桌面的“开始”按钮，点击开始菜单栏里的“附件”，选择“录音机”选项，点击录音机“文件”选项，进入“声音选定”设置，把属性一栏设置成“8000Hz，8位，单声道，7KB/秒”(见图3.2.1)。点击确定，然后开始语言信号的采集，采集时间为1秒左右为最佳。采集的声音文件以“.wav”格式存储（见图3.2.2）。 图3.2.1