语音信号滤波去噪报告书课程：数字信号处理指导老师:完成组员：完成日期：2013.01.05摘要本课程设计主要是下载一段语音信号，绘制其波形并观察其频谱。

然后在该语言信号中加一个噪音，利用布莱克曼和矩形窗窗设计一个FIR滤波器，对该语音信号进行虑噪处理，然后比较滤波前后的波形与频谱。

在本课程设计中，是用MATLAB的集成环境完成一系列的设计。

首先对加噪的语音信号进行虑波去噪处理，再比较滤波前后的频率响应曲线，若一样则满足所设计指标，否则不满足。

也可以调用函数sound听滤波前后其语音信号是否带有噪声。

若无噪声也说明该滤波器的设置也是成功的。

关键词语音信号；MATLAB；FIR滤波器；滤波去噪；1 引言人们在语音通信的过程中将不可避免的会受到来自周围环境的干扰，例如传输媒介引入的噪声，通信设备内部的电噪声，乃至其他讲话者的话音等。

正因为有这些干扰噪声的存在，接受者接受到的语音已不是原始的纯净语音信号，而是受噪声干扰污染的带噪声语音信号。

而本课程设计就是利用MATLAB集成环境用布莱克曼窗的方法设计一个FIR滤波器，对语音信号进行滤波去噪处理，并将虑噪前后的频谱图进行对比。

1.1 课程设计目的数字信号处理课程设计是数字信号处理课程的重要实践性环节，是学生在校期间一次较全面的工程师能力训练，在实现学生总体培养目标中占有重要地位。

综合运用本课程的理论知识进行频谱分析以及滤波器设计，通过理论推导得出相应结论，并利用MATLAB 作为编程工具进行计算机实现，从而复习巩固了课堂所学的理论知识，提高了对所学知识的综合应用能力，并从实践上初步实现了对数字信号的处理。

本课程设计能使学生对通信工程领域各种技术的DSP实现的设计有较熟练的掌握。

且通过自身的实践，对DSP的设计程序、内容和方法有更深入的掌握，提高实际运用的能力。

并可综合运用这些知识解决一定的实际问题，使学生在所学知识的综合运用能力上以及分析问题、解决问题能力上得到一定的提高。

1.2课程设计的要求（1）、录制一段个人自己的语音信号，并对录制的信号进行采样，画出采样后语音信号的时域波形和频谱图。

（2）、给定滤波器的性能指标，采用窗函数法和双线性变换设计滤波器，并画出滤波器的频率响应。

（3）、用自己设计的滤波器对采集的信号进行滤波，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化并回放语音信号；（4）、通过利用各种不同的开发工具实现语音信号的滤波去噪，掌握数字信号的分析方法和处理方法。

而且通过课程设计能够培养学生严谨的科学态度，认真的工作作风和团队协作精神。

（5）、在老师的指导下，要求独立完成课程设计的全部内容，并按要求编写课程设计学年论文，能正确阐述和分析设计和实验结果。

1.3设计平台MATLAB名称是有两个英文单词Matrix和Laboratory的前三个字母组成。

MATLAB7.0是美国MathWorks公司开发的优秀计算软件MATLAB的最新版本。

MATLAB自20世纪80年代面世以来，以其强大的数值计算能力、优秀的绘图功能以及与其他软件良好的交互功能在众多的数学计算软件中独领风骚，特别是它源代码的开放性使用户可以二次开发，受到了广大使用者的格外赞赏。

MATLAB是一个为科学和工程计算机专门设计的交互式大型软件，是一个可以完成各种精确计算和数据处理的、可视化的、强大的计算工具。

它集图和精确计算与一身，在应用数学、物理、化工、机电工程、医药、金融和其他需要进行复杂数值计算的领域得到了广泛应用。

它不仅是一个在各类工程设计中便于使用的计算工具，在世界各地的高等院校中十分流行，在各类工业应用中更有不俗的表现。

MATLAB可以几乎所有的PC机和大型计算机上运行，适用于Window、UNIX等多种系统平台。

本课程设计我们就可以直接诶使用MATLAB提供的模块，实现模拟通信系统的仿真。

MATLAB软件有很强的开放性和适应性。

在保持内核不变的情况下，MATLAB可以针对不同的应用学科推出相应的工具箱，目前已经推出了图像处理工具箱、信号处理工具箱、小波工具箱、神经网络工具箱以及通信工具箱等多个学科的专用工具箱，极大的方便了不同学科的研究工作。

国内已有越来越多的科研和技术人员认识到MABLAB的强大作用，并在不同领域内使用MATLAB来快速实现科研构想和提高工作效率。

2 设计原理2.1 FIR滤波器FIR(Finite Impulse Response)滤波器：有限长单位冲激响应滤波器，是数字信号处理系统中最基本的元件，它可以在保证任意幅频特性的同时具有严格的线性相频特性，同时其单位抽样响应是有限长的，因而滤波器是稳定的系统。

因此，FIR 滤波器在通信、图像处理、模式识别等领域都有着广泛的应用。

FIR 数字滤波器以其良好的线性特性被广泛应用于现代电子通信系统中，是数字信号处理的重要内容之一。

在实际信号处理中，往往要求系统兼具实时性和灵活性，而已有的一些软件或硬件实现方案(如DSP)则难以同时达到这两方面的要求。

使用具有并行处理特性的FPGA 来实现FIR 滤波器，既有很强的实时性，又兼顾了灵活性，为数字信号处理提供了一种很好的解决方案。

FIR 滤波器系数计算较为繁琐，在设计时借助Matlab 工具箱，选择合适的窗函数，可以方便地计算滤波器系数，并分析其幅频、相频特性。

有限长单位冲激响应（FIR ）滤波器有以下特点：(1)、系统的单位冲激响应h (n)在有限个n 值处不为零；(2)、系统函数H (z)在|z|>0处收敛，极点全部在z = 0处（因果系统）；(3)、 结构上主要是非递归结构，没有输出到输入的反馈，但有些结构中（例如频率抽样结构）也包含有反馈的递归部分。

设FIR 滤波器的单位冲激响应h (n)为一个N 点序列，0 ≤ n ≤ N —1，则滤波器的系统函数为：H(z) = ∑-=10)(N n z n h -n (2-1)就是说，它有（N —1）阶极点在z = 0处，有（N —1）个零点位于有限z 平面的任何位置因此H(z)是永远稳定的。

稳定和香味特性是FIR 滤波器突出的优点。

FIR 滤波器有以下几种基本结构：横截型（卷积型、直接型）、级联型、频率抽样型、快速卷积结构。

FIRDF 的设计方法主要分为两类：第一类是基于逼近理想滤波器特性的方法，包括窗函数法、频率采样法和等波纹最佳逼近法；第二类是最优设计法。

2.2窗口设计法数字信号处理的主要数学工具是博里叶变换．而傅里叶变换是研究整个时间域和频率域的关系。

不过，当运用计算机实现工程测试信号处理时，不可能对无限长的信号进行测量和运算，而是取其有限的时间片段进行分析。

做法是从信号中截取一个时间片段，然后用观察的信号时间片段进行周期延拓处理，得到虚拟的无限长的信号，然后就可以对信号进行傅里叶变换、相关分析等数学处理。

无线长的信号被截断以后，其频谱发生了畸变，原来集中在f(0)处的能量被分散到两个较宽的频带中去了(这种现象称之为频谱能量泄漏)。

为了减少频谱能量泄漏，可采用不同的截取函数对信号进行截断，截断函数称为窗函数，简称为窗。

窗函数设计法的基本思路是用FIRDF逼近希望的滤波特性。

设希望逼近的滤波器的频率响应函数为h d(e jω)，其单位脉冲响应用h d(n)表示。

为了设计简单方便，通常选择h d(e j ω)为具有片段常数特性的理想滤波器。

因此h d(n)是无限长非因果序列，不能直接作为FIRDF的单位脉冲响应。

窗函数设计法就是截取h d(n)为有限长的一段因果序列，并用合适的窗函数进行加权做为FIRDF的单位脉冲响应h(n)。

下面介绍窗函数设计法的基本设计过程。

窗口设计法的主要工作是计算h d(n)和w(n)，但当H d(e jω)较为复杂时，h d(n)就不容易由反付里叶变换求得。

这时一般可用离散付里叶变换代替连续付里叶变换，求得近似值。

窗口法的设计步骤如下：（1）、通过傅里叶变换忽的理想滤波器的单位脉冲响应h d(n)。

（2）、根据指标选择窗口形状、大小和位置。

确定窗口类型的主要依据是过渡带宽和阻带最小衰耗的指标。

（3）、给定理想频响由h d (e j ω)和h d (n)，加窗得h(n)=w(n)h d (n)。

（4）、检验滤波器的性能。

由h （n ）求H (e j ω) 是否在误差容限之内。

如果不满足，则返回第（2）步。

以上步骤中h d (n)、H(e j ω)的计算可采用傅氏变换的现成公式和程序，窗函数w(n)也是现成的。

但整个设计过程不能一次完成，因为窗口类型和大小的选择没有解析公式可一次算出。

整个设计可用计算机编程来做。

窗口法的优点是简单，有闭合的公式可用，性能及参数都有表格资料可查，计算程序简单，较为实用。

缺点是当H d (e j ω)较为复杂时，h d (n)就不容易由反付里叶变换求得。

边界频率因为加窗的影响而不易控制。

窗口函数对理想特性的影响：改变了理想频响的边沿特性，形成过渡带，宽为N 4π , 等于WR(ω)的主瓣宽度；过渡带两旁产生肩峰和余振（带内、带外起伏），取决于 WR(ω)的旁瓣，旁瓣多，余振多；旁瓣相对值大，肩峰强，与 N 无关；N 增加,过渡带宽减小,肩峰值不变。

因主瓣附近(2-2)其中x=N ω/2,所以N 的改变不能改变主瓣与旁瓣的比例关系，只能改变W R （ω）的绝对值大小和起伏的密度，当N 增加时，幅值变大，频率轴变密，而最大肩峰永远为8.95%，这种现象称为吉布斯（Gibbs ）效应。

肩峰值的大小决定了滤波器通带内的平稳程度和阻带内的衰减，所以对滤波器的性能有很大的影响。

改变窗函数的形状，可改善滤波器的特性，窗函数有许多种，但要满足以下两点要求:窗谱主瓣宽度要窄，以获得较陡的过渡带；相对于主瓣幅度，旁瓣要尽可能小，使能量尽量集中在主瓣中，这样就可以减小肩峰和余振，以提高阻带衰减和通带平稳性。

但实际上对同样长度的窗这两点不能兼得，一般总是通过增加主瓣宽度来换取对旁瓣的抑x x N N N N N W R sin 2/)2/sin()2/sin()2/sin()(=≈=ωωωωω制。

2.3 布莱克曼窗布莱克曼窗的时域形式可表示为：R N n N n n )]14cos(08.0)12(cos 5.042.0[)(w -+--=ππN (n) (2-3) 它的频域特性为：W 0.42W )(=ωR 0.25[W )(+ωR W N +--)12(πωR ]12(-+N πω 0.04[W +R )14(--N πωW +R )]14(-+N πω (2-4) 其中)(ωR W 为矩形窗函数的幅度频率特性。

增加一个二次谐波余弦分量，可进一步降低旁瓣，但主瓣宽度进一步增加，为N π12。