基于MATLAB的IIR滤波器的语音信号去噪要点.docVIP

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