基于matlab的iir滤波器的语音信号去噪要点.docxVIP

摘 要 滤波器设计在数字信号处理中占有极其重要的地位，本次课程设计主要是录制一段语音 信号对其进行加噪处理，然后利用 IIR 低通滤波器对加有随机噪声的语音信号进行滤波处理 及时频谱分析，画出滤波之后的频谱图与时域波形，并对信号滤波处理前后进行分析比较， 分析信号的变化。通过对对所设计滤波器的仿真和频率特性分析，由仿真结果可以看出，所 设计的滤波器能够实现对语音信号的语音有效去噪，并对滤波前后的语音信号进行对比。 关键词： 去噪；滤波器； MATLAB 1 目 录 一 语音信号去噪的设计任务 错误！未定义书签。 二 语音信号去噪的基本原理 3 2.1 数字滤波器的基本设计方法 3 2.2 双线性变换法 4 2.3 数字滤波器设计基本思想 5 2.4 数字滤波器的设计步骤 6 2.5 采样定理 7 三 基于 MATLAB的仿真结果及结果分析 10 3.1 IIR 高通滤波器的仿真 10 3.2 原始语音信号的录制 10 3.3 语音信号的时频域分析 11 3.4 加随机噪声后的时频域分析 12 3.5 滤波前后的时频域比较 15 总结 17 参考文献 18 致谢 19 附录 20 2 一 基本原理 1.1 数字滤波器的基本设计方法 IIR 数字滤波器的设计一般有两种方法： 一个是借助模拟滤波器的设计方法进行。 其设计 步骤是，先设计模拟滤波器，再按照某种方法转换成数字滤波器。这种方法比较容易一些， 因为模拟滤波器的设计方法已经非常成熟， 不仅有完整的设计公式， 还有完善的图表供查阅；另外一种直接在频率或者时域内进行， 由于需要解联立方程， 设计时需要计算机做辅助设计。其设计步骤是： 先设计过渡模拟滤波器得到系统函数 H a (s) ，然后将 H a (s) 按某种方法转换成数字滤波器的系统函数 H (z) [1] 。为了保证转换后的 H ( z) 稳定且满足技术指标要求， 对转换关 系提出两点要求 : （1）因果稳定的模拟滤波器转换成数字滤波器，仍是因果稳定的。 （2）数字滤波器的频率相应模仿模拟滤波器的频响特性， s平面的虚轴映射为 z平面的单 位圆，相应的频率之间呈线性关系。 利用模拟滤波器成熟的理论设计 IIR 数字滤波器的过程是： (1) 确定数字低通滤波器的技术指标：通带边界频率 p 、通带最大衰减 p 、阻带截止频 率 s 、阻带最小衰减 s 。 将数字低通滤波器的技术指标转换成相应的模拟低通滤波器的技术指标。 按照模拟低通滤波器的技术指标设计过渡模拟低通滤波器。 用所选的转换方法，将模拟滤波器 H a (s) 转换成数字低通滤波器系统函数 H (z) 。 IIR 数字滤波器的设计流程图 2-1 如下： 数字滤波器 变换 模拟滤波器 技术指标 Ω=ｇ（ ω ） 技术指标 Ω 模 拟 滤 波 器 变换 设计方法 数字滤波器 模拟滤波器 H (Z )S=f （Z） Ha (S) [1] 图 2-1 IIR 数字滤波器的设计步骤流程图 成熟的模拟滤波器设计方法主要有脉冲响应不变法和双线性变换法。 3 2.2 双线性变换法 脉冲响应不变法的主要缺点是产生频率响应的混叠失真。这是因为从 S 平面到Ｚ平面是 多值的映射关系所造成的。为了克服这一缺点，可以采用非线性频率压缩方法，将整个频率 轴上的频率范围压缩到 - π/ T～π/ T 之间，再用 z=esT 转换到 Z 平面上。也就是说，第一步先 将整个 S 平面压缩映射到 S1 平面的 - π/ T～π/ T 一条横带里；第二步再通过标准变换关系 z=es1T 将此横带变换到整个 Z 平面上去。这样就使 S平面与 Z 平面建立了一一对应的单值关系， 消除了多值变换性，也就消除了频谱混叠现象，映射关系如图 2-2 所示 [2] 。 j j 1 jIm[z] / T －1 1 o o 1 oRe[z] / T S平面 S 平面 Z平面 1 [3] 图 2-2 双线性变换的映射关系 为了将 S 平面的整个虚轴 j Ω 压缩到 S1 平面 j Ω 1 轴上的 - π / T 到 π/ T 段上，可以通过 以下的正切变换实现 2 tan 1T （2-1 ） T 2 式中 , T 仍是采样间隔。 当Ω 1 由- π/ T 经过 0 变化到 π/ T 时，Ω由- ∞经过 0 变化到 +∞，也即映射了整个 j Ω轴。将式（ 2-1 ）写成 j 2 . e j 1T / 2 e T e j 1T / 2 e  j T / 2 j 1T / 2 （2-2 ） 将此关系解析延拓到整个 S 平面和 S1 平面，令 j Ω=s，j Ω1=s1，则得 s 2 . esi T / 2 e Tesi T / 2 e  si T / 2 si T / 2  2 tanh s1T