语音信号滤波去噪——使用双线性变换法设计切比雪夫II型滤波器精选.docVIP

语音信号滤波去噪——使用双线性变换法设计切比雪夫II型滤波器 学生姓名：李谷越 指导老师：胡双红 摘 要 本课程设计的主要目的石设计一个切比雪夫II型滤波器，并对语音信号进行滤波去噪。开发平台为MATLAB7.1,设计方法为双线性变换法。用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用一个满足指标的II型IIR滤波器，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析 关键词 课程设计；滤波去噪；IIR滤波器；双线性变换法；MATLAB7.1 1 引 言 用麦克风采集一段8000Hz8k的单声道语音信号，绘制波形并观察其频谱，给定Hz，阻带截止频率为2100Hz，通带波纹为1dB，阻带波纹为60dB，用一个满足指标的II型IIR滤波器，对该语音信号进行滤波去噪处理 课程设计目的 综合运用数字信号处理的理论知识进行频谱分析和滤波器设计，通过理论推导得出相应结论，再利用 MATLAB 作为编程工具进行计算机实现，从而加深对所学知识的理解，建立概念。本次使学生加深对理论知识理解的同时增强其逻辑斯维的能力，另一方面对课堂教学中的理论知识做一个总结和补充。熟悉和巩固模拟滤波器的设计方法和原理 掌握双线性变换方法设计方法 实现滤波器设计的有关经典算法 熟练掌握使用MATLAB语言设计各种要求的数字滤波器 1.2课程设计的要求 （1）学会 MATLAB 的使用，掌握 MATLAB 的程序设计方法 （2）滤波器指标必须符合工程实际根据模拟滤波器的性能指标，确定数字滤波器指标； （3）采用双线性变换法，设计满足上述性能指标要求的ChebyshevI型数字低通滤波器； 4）设计完后应检查其频率响应曲线是否满足指标 （5）处理结果和分析结论应该一致，而且应符合理论 （6）独立完成课程设计并按要求编写课程设计报告书 1.3设计平台 MATLAB 2 设计原理 IIR滤波器具有无限长脉冲响应，因此能够与模拟滤波器相匹敌；因此，IIR滤波器设计的基本方法是利用复值映射将大家熟知的模拟滤波器变换为数字滤波器。这一方法的优势在于各种模拟滤波器设计（AFD）A/D（模拟-数字）滤波器变换。然而，AFD表格仅对低通滤波器适用，而同时要得到设计其他频率选择性滤波器（高通、带通、带阻等等）。为此，需要对低通滤波器实行频带变换。这种IIR滤波器设计的基本方法存在两种途径： 途径1： 期望的IIR滤波器 途径2： 期望的IIR滤波器 我们将研究途径二的设计基本方法，在这种IIR滤波器设计方法将按下列步骤进行： （1）设计模拟低通滤波器。 （2）研究并实行滤波器变换以得到数字低通滤波器。 （3）研究并实行频带变换以便从数字低通滤波器得到其他数字滤波器。 设IIR滤波器的输入序列为，则IIR滤波器的输入序列和输出序列之间的关系可用下列方程式表示： (2-1) 其中和是滤波器的系数，且中至少有一个非零。与上述差分方程对应的传递函数为： (2-2) 由传递函数可以发现无限长单位冲激响应滤波器有如下特点： 单位冲激响应是无限长的。 （2）系统传递函数在有限z平面上有极点存在。 （3）结构上存在着输出到输入的反馈，也就是结构上是递归型的。 设是某个模拟滤波器的频率响应，那么低通滤波器在幅度平方响应上的技术指标给出为： （2-3） （2-4） 式中是通带波纹参数，是通带截止频率以rad/s（弧度每秒）计，A是阻带衰减参数，以及是阻带截止频率以rad/s计。 图2.1 模拟低通滤波器技术指标要求 由图可见，必须满足 （2-5） （2-6） 参数和A是分别与以dB计的参数和有关的，这些关系是 （2-7） （2-8） 波纹和的绝对标尺是通过下式与和A有关的： 利用幅度平方响应给出的模拟滤波器要求（2-1）和（2-2）式不包含任何相位信息。现在，为了求s域的系统函数，考虑 那么有 （2-12） 或者有 （2-9） 轴是以镜像对称方式分布的【2】。 2.2 切比雪夫II型滤波器 切比雪夫滤波器是在通带或阻带上频率响应幅度等波纹波动的滤波器。在通带波动的为“I型切比雪夫滤波器”，在阻带波动的为“II型切比雪夫滤波器”。切比雪夫滤波器在过渡带比巴特沃斯滤波器的衰减快，但频率响应的幅频特性不如后者平坦。切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小，但是在通频带内存在幅度波动。 双线性变换由下式定义： （2-11） 这里T是一个参数。也叫线性分式变换，因为当乘开之后得到 这在每个变量上都是线性的，如果另