语音信号滤波噪.docVIP

语音信号滤波去噪——使用双线性变换法设计的椭圆滤波器 学生姓名：王哲 指导老师：胡双红 摘 要 本课程设计主要利用MATLAB。设计一个语音信号滤波去噪——使用双线性变换法设计的椭圆滤波器，用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用双线性变换法设计的一个满足指标的椭圆IIR滤波器，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析，根据结果和学过的理论得出的结论合理。 关键词 课程设计；滤波去噪；IIR滤波器；双线性变换法；椭圆滤波器；MATLAB 1 引言 本课程设计主要利用MATLAB，设计一个语音信号滤波去噪——使用双线性变换法设计的椭圆滤波器。DSP的应用主要是在DSP处理器上或者软件中执行的一种算法，这就需要一定的编程量，就可以利用MATLAB实现。而对滤波器和频谱分析仪的结构的分析和设计是DSP的最重要的内容之一，所以本课程设计要对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析。 1.1课程设计目的 用麦克风采集一段语音信号，绘制波形并观察其频谱，给定技术指标，用双线性变换法设计一个椭圆低通滤波器，对语音信号进行滤波处理，比较滤波前后的波形和频谱并进行分析【1】。 1.2课程设计的要求 （1）滤波器指标必须符合工程实际。 （2）设计完后应检查其频率响应曲线是否满足指标。 （3）处理结果和分析结论应该一致，而且应符合理论。 （4）独立完成课程设计并按要求编写课程设计报告书。 1 .3设计平台 美国Mathwork公司于1967年推出了“Matrix Laboratory”（缩写为Matlab）软件包，并不断更新和扩充。目前最新的5.x版本（windows环境）是一种功能强、效率高便于进行科学和工程计算的交互式软件包。其中包括：一般数值分析、矩阵运算、数字信号处理、建模和系统控制和优化等应用程序，并集应用程序和图形于一便于使用的集成环境中。在此环境下所解问题的Matlab语言表述形式和其数学表达形式相同，不需要按传统的方法编程。 Matlab作为一种新的计算机语言，要想运用自如，充分发挥它的威力，也需先系统地学习它。但由于使用Matlab编程运算与人进行科学计算的思路和表达方式完全一致，所以不象学习其它高级语言--如Basic、Fortran和C等那样难于掌握。实践证明，你可在几十分钟的时间内学会Matlab的基础知识，在短短几个小时的使用中就能初步掌握它.从而使你能够进行高效率和富有创造性的计算。 Matlab大大降低了对使用者的数学基础和计算机语言知识的要求，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一高效的科研助手。自推出后即风行美国，流传世界。 2 设计原理 2.1 IIR滤波器 根据IIR数字滤波器的设计原理，提出了IIR数字滤波器的快速设计方法，并在Matlab环境下做出了实现快速设计IIR 数字滤波器的设计系统.在该系统中，只需要将数字滤波器的技术性能指标根据指定的或是需要的设计方法（如双线性变换法或脉冲响应不变法等），转换为模拟滤波器的技术性能指标，根据指定的模拟滤波器设计出相应的数字滤波器，该系统中指定的模拟滤波器为巴特沃斯滤波器和切比雪夫滤波器，该系统为快速、高效地设计IIR数字滤波器提供了一个可靠而有效的工作平台。 2.2 双线性变换法 脉冲响应不变法的主要缺点是频谱交叠产生的混淆，这是从S平面到Z平面的标准变换z＝esT的多值对应关系导致的,为了克服这一缺点，设想变换分为两步： 第一步：将整个S平面压缩到S1平面的一条横带里； 第二步：通过标准变换关系将此横带变换到整个Z平面上去。 由此建立S平面与Z平面一一对应的单值关系，消除多值性，也就消除了混淆现象。 图 2.1 双线性变换法中的复平面映射 为了将s平面的jΩ轴压缩到s1平面jΩ轴上的- 一段上，可通过以下的正切变换实现： 这里C是待定常数，下面会讲到用不同的方法确定C，可使模拟滤波器的频率特性与数字滤波器的频率特性在不同频率点有对应关系。 经过这样的频率变换，当Ω1由 时 Ω由 即映射了整个jΩ轴。 ??? 将这一关系解析延拓至整个s平面，则得到s平面 平面的映射关系： 再将s1平面通过标准变换关系映射到z平面，即令 通常取C=2/T 最后得S平面与Z平面的单值映射关系：  双线性换法的主要优点是S平面与Z平面一一单值对应，S平面的虚轴(整个jΩ)对应于Z平面单位圆的一周，S平面的Ω=0处对应于Z平面的ω=0处， 对应 即数字滤波器的频率响应终 止于折迭频率处，所以双线性变换不存在混迭效应。 2.3 椭圆低通滤波器 椭圆模拟低通原型滤波器的平方幅值响应函数为： 式中，为小于1的正数，表示波