长沙理工大学-Read.DOC

语音信号滤波去噪——使用脉冲响应不变法设计的巴特沃斯滤波器 学生姓名：熊隽迪 指导老师：胡双红 摘 要 本课程设计主要运用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用一个满足指标的，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析，根据结果和学过的理论得出合理的结论。MATLAB 1 引 言 本课程设计主要利用麦克风采集一段语音信号，1.2 课程设计的要求 （1）滤波器指标必须符合工程设计。 （2）设计完后应检查其频率响应曲线是否满足指标。 （3）处理结果和分析结论应该一致，而且应符合理论。 （4）独立完成课程设计并按要求编写课程设计报告。 1.3 设计平台 MATLAB7.1 MATLAB是美国MathWorks公司生产的一个为科学和工程计算专门设计的交互式大型软件，是一个可以完成各种精确计算和数据处理的、可视化的、强大的计算工具。它集图示和精确计算于一身，在应用数学、物理、化工、机电工程、医学、金融和其他需要进行复杂数值计算的领域得到了广泛应用。它不仅是一个在各类工程设计中便于使用的计算工具，而且也是一个在数学、数值分析和工程计算等课程教学中的优秀的数学工具，在世界各地的高和大型计算机上运行，适用于Windows、UNIXMATLAB作为一种科学计算的高级语言之所以受欢迎，就是因为它有丰富的函数资源和工具箱资源，编程人员可以根据自己的需要选择函数，而无需再去编写大量繁琐的程序代码，从而减轻了编程人员的工作负担，被称为第四代编程语言的MATLAB最大的特点就是简洁开放的程序代码和直观实用的开发环境[1]。 2 设计原理 用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用一个满足指标的，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析IIR数字滤波器具有无限宽的冲激响应，与模拟滤波器相匹配，所以IIR滤波器的设计可以采取在模拟滤波器设计的基础上进一步变换的方法。其设计方法主要有经典设计法、直接设计法和最大平滑滤波器设计法。FIR数字滤波器的单位脉冲响应是有限长序列。它的设计问题实质上是确定能满足所要求的转移序列或脉冲响应的常数问题，设计方法主要有窗函数法、频率采样法和等波纹最佳逼近法等。IIR滤波器具有无限长脉冲响应，因此能够与模拟滤波器相匹敌；因此，IIR滤波器设计的基本方法是利用复值映射将大家熟知的模拟滤波器变换为数字滤波器。这一方法的优势在于各种模拟滤波器设计（AFD）表格和映射在文献中普遍都能获得。这个基本方法称为A/D（模拟-数字）滤波器变换。然而，AFD表格仅对低通滤波器适用，而同时要得到设计其他频率选择性滤波器（高通、带通、带阻等等）。为此，需要对低通滤波器实行频带变换。这种IIR滤波器设计的基本方法存在两种途径： 途径1： 期望的IIR滤波器 途径2： 期望的IIR滤波器 我们将研究途径二的设计基本方法，在这种IIR滤波器设计方法将按下列步骤进行： （1）设计模拟低通滤波器。 （2）研究并实行滤波器变换以得到数字低通滤波器。 （3）研究并实行频带变换以便从数字低通滤波器得到其他数字滤波器[2]。 2.2 巴特沃斯滤波器 巴特沃斯滤波器是电子滤波器的一种。巴特沃斯滤波器的特点是通频带的频率响应曲线最平滑。这种滤波器最先由英国工程师斯替芬·巴特沃斯(Stephen Butterworth)在1903发表英国《无线电工程》期刊的一篇论文中提出的。 巴特沃斯滤波器的特性是通频带内的频率响应曲线最大限度平坦，没有起伏，而在组频带则逐渐下降为零。在振幅的对数对角频率的波得图上，从某一边界角频率开始，振幅随着角频率的增加而逐步减少，趋向负无穷大。 其振幅平方函数具有如2-1式： （2-）2.1所示： 图2.1 巴特沃兹filter 振幅平方函数 过渡带：通带→阻带间过渡的频率范围，Ωc：截止频率。 理想滤波器的过渡带为Ω，阻带|H(jΩ)|=0，通带内幅度|H(jΩ)|=常数，H(jΩ)线性相?位。通带内，分母Ω/Ωc1,相应(Ω/Ωc)2N随N的增加而趋于0，A(2)→1，/Ωc1，N的增加，Ωe/Ωc1，所以A(Ω2)快速下降。 Ω=Ωc时， ，幅度衰减 ，相当于3bd衰减点。 振幅平方函数的极点可写成如式2-2： Ha(-s).Ha(s)= （2-2） 可分解为2N个一次因式令分母为零，→ 可见，Butterworth 滤波器的振幅平方函数有2N个极点，它们均匀对称地分布在|s|=Ωc的圆周上。 2.3 脉冲响应不变法 如果从模拟到数字滤波器我们想要保留脉冲响应的形状，那么就得到一种方法称为脉冲不