信号处理课程设计基于MATLAB的FIR滤波器语音信号去噪Word下载.docx

信号处理课程设计基于MATLAB的FIR滤波器语音信号去噪Word下载.docx

《信号处理课程设计基于MATLAB的FIR滤波器语音信号去噪Word下载.docx》由会员分享，可在线阅读，更多相关《信号处理课程设计基于MATLAB的FIR滤波器语音信号去噪Word下载.docx（28页珍藏版）》请在冰豆网上搜索。

MATLAB7.0

前言

语音信号在数字信号处理中占有极其重要的地位，因此选择通过对语音信号的研究来巩固和掌握数字信号处理的基本能力十分具有代表性。

对数字信号处理离不开滤波器，因此滤波器的设计在信号处理中占有极其重要的地位。

而MATLAB软件工具箱提供了对各种数字滤波器的设计。

本设计“在MATLAB平台上实现对语音信号的去噪研究与仿真”综合运用了数字信号处理的各种基本知识，进而对不带噪语音信号进行谱分析以及带噪语音信号进行谱分析和滤波处理。

通过理论推导得出相应的结论，再通过利用MATLAB作为编程工具来进行计算机实现比价已验证推导出来的结论。

在设计过程中，通过设计FIR数字滤波器来完成滤波处理。

在设计过程中，运用了MATLAB对整个设计中的图形的绘制和一些数据的计算以及仿真。

此次课程中主要的要求是用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用kaiser窗设计一个满足指标的FIR滤波器，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析，根据结果和学过的理论得出合理的结论。

与不同信源相同滤波方法的同学比较各种信源的特点，与相同信源不同滤波方法的同学比较各种滤波方法性能的优劣。

第一章设计要求及环境

1.1课程设计的要求

1、学会MATLAB的使用，掌握MATLAB的程序设计方法；

2、掌握windows环境下语音信号采集方法；

3、掌握数字信号处理的基本概念、基本理论和方法；

4、掌握MATLAB设计FIR数字滤波器的方法；

5、学会用MATLAB对信号进行分析和处理。

1.2工作平台简介

课程设计的主要设计平台式MATLAB7.0。

如下图1-1所示：

MATLAB的名称源自MatrixLaboratory，它是美国MathWorks公司生产的一个为科学和工程计算专门设计的交互式大型软件，是一个可以完成各种精确计算和数据处理的、可视化的、强大的计算工具。

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

它不仅是一个在各类工程设计中便于使用的计算工具，而且也是一个在数学、数值分析和工程计算等课程教学中的优秀的教学工具，在世界各地的高等院校中十分流行，在各类工业应用中更有不俗的表现。

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

总的来说，该软件有三大特点。

一是功能强大。

具有数值计算和符号计算、计算结果和编程可视化、数学和文字统一处理、离线和在线计算等功能；

二是界面友善、语言自然。

MATLAB以复数处理作为计算单元，指令表达与标准教科书的数学表达式相近；

三是开放性强。

当学好MATLAB的同时，会更好的帮助自己去就解决一些难题，而且MATLAB拥有非常好的发展前途，对我们未来的帮助也是不可限量的。

图1-1MATLAB7.0设计平台

第二章、设计原理

2.1数字信号处理简介

数字信号处理（DigitalSignalProcessing，简称DSP）是一门涉及许多学科而又广泛应用于许多领域的新兴学科。

20世纪60年代以来，随着计算机和信息技术的飞速发展，数字信号处理技术应运而生并得到迅速的发展。

在过去的二十多年时间里，数字信号处理已经在通信等领域得到极为广泛的应用。

数字信号处理是利用计算机或专用处理设备，以数字形式对信号进行采集、变换、滤波、估值、增强、压缩、识别等处理，以得到符合人们需要的信号形式。

随着信息技术的迅猛发展，数字信号处理已成为一个极其重要的学科和技术领域。

在通信、语音、图像、自动控制和家用电器等众多领域得到了广泛的应用。

数字滤波是数字信号处理的重要环节，它在数字信号处理中占有着重要的地位,它具有可靠性好、精度高、灵活性大、体积小、重量轻等优点。

随着数字技术的发展，数字滤波器越来越受到人们的重视，广泛地应用于各个领域。

数字滤波器的输入输出信号都是数字信号，它是通过一定的运算过程改变输入信号所含频率成分的相对比例或者滤除某些频率成分来实现滤波的，这种运算过程是由乘法器、加法器和单位延迟器组成的。

数字滤波器是数字信号处理技术的重要内容，其对数字信号进行的最常见处理是保留数字信号中的有用频率成分和去除信号中的无用频率成分。

按照时间域的特性，数字滤波器可以分为无限冲激脉冲响应数字滤波器（IIR滤波器）和有限冲激脉冲响应数字滤波器（FIR滤波器）。

2.2FIR滤波器

有限长单位脉冲响应数字滤波器（FiniteImpulseResponseDigitalFilter,缩写FIRDF）：

有限长单位冲激响应滤波器，是数字信号处理系统中最基本的元件，最大优点是可以实现线性相性滤波，它可以在保证任意幅频特性的同时具有严格的线性相频特性，同时其单位抽样响应是有限长的，因而滤波器是稳定的系统。

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

FIR滤波器的设计方法主要分为两类：

第一类是基于逼近理想滤波器器特性的方法包括窗函数法、频率采样法、和等波纹最佳逼近法；

第二类是最优设计法。

设FIRDF的单位脉冲响应

的长度为

，则其频率响应函数为

（2-1）

一般将

表示成如下形式：

（2-2）

式中，

是

的实函数（可以去负值）。

与前面的表示形式，即

相比，

与

不同。

与

为了区别于幅频响应函数

和相频响应函数

，称

为幅频特性函数，称

为相频特性函数。

第一类线性相位FIRDF的相位特性函数是

的严格线性函数：

（2-3）

第二类线性相位FIRDF的相位特性函数如下：

（2-4）

是常数，

是起始相位。

在信号处理中很有实用价值（如希伯尔特变换器），这是FIRDF除了线性相位滤波外，还具有真正交变换作用。

2.3窗口设计法

FIR滤波器的设计方法有许多种，如窗函数设计法、频率采样设计法和最优化设计法等。

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

设希望逼近的滤波器的频率响应为

，其单位脉冲响应用

表示。

为了设计简单方便，通常选择

为具有片段常数特性的理想滤波器。

因此

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

窗口设计法就是截取

为有限长的一段因果序列，并用合适的窗口函数进行加权作为FIRDF的单位脉冲响应

。

窗口设计法基本步骤如下：

（1）构造希望逼近的频率响应函数

以低通线性相位FIRDF设计为例，一般选择

为线性理想低通滤波器，即

（2-5）

（2）求出

对

进行IFT得到

（2-6）

（3）加窗得到FIRDF的单位脉冲响应

，

（2-7）

称为窗口函数，其长度为

如果要求第一类线性相位FIRDF，则要求

关于

点偶对称。

而

点偶对称，所

，同时要求

常见的窗函数，可以分为以下主要类型：

（1）幂窗--采用时间变量某种幂次的函数，如矩形、三角形、梯形或其它时间（t）的高次幂；

（2）三角函数窗--应用三角函数，即正弦或余弦函数等组合成复合函数，例如汉宁窗、海明窗、kaiser窗等；

（3）指数窗--采用指数时间函数，例如高斯窗等。

其性能如表2-1所示：

表2-1常见窗函数性能表

名称

滤波器

过渡带宽

最小阻带衰减

矩形

1.8π/M

21dB

PARZENWIN

6.6π/M

56dB

巴特利特

6.1π/M

25dB

FLATTOPWIN

19.6π/M

108dB

汉宁

6.2π/M

44dB

GAUSSWIN

5.8π/M

60dB

汉明

51dB

BARTHANNWIN

3.6π/M

40dB

布莱克曼

11π/M

74dB

BLACKMANHARRIS

16.1π/M

109dB

BOHMANWIN

51.5dB

CHEBWIN

15.2π/M

113dB

NUTTALLWIN

15.4π/M

TUKEYWIN

2.4π/M

22dB

2.4kaiser窗

人们研究窗函数主要是它可以有效地抑制频谱泄漏，本论文中采用的加窗插值FFT方法，窗函数的选择非常重要。

在间谐波检测的频谱分析时，最好的是窗函数主瓣窄、旁瓣低且跌落快，但对于同一窗函数，以上两个要求又是互相矛盾的。

因为要是增加主瓣的宽度，旁瓣就会降低，反之，若想主瓣能变得又高又窄，旁瓣就会增高。

在选择窗的时候，应根据间谐波特征和研究目的来选择。

目前，已经有20多种窗函数，在电力系统谐波与间谐波检测中常用的窗函数有矩形窗、海宁窗（Hanning）、汉明窗（Hamming）和布来克曼窗（Blackman）。

数字信号处理领域较早就提出了凯泽窗（Kaiser），它通过改变参数可以达到不同的性能，正由于其优良的窗函数特性，广泛应用于高通、低通、带通、带阻等各种滤波器的设计。

凯泽窗（Kaiser）是利用贝塞尔函数来逼近需要的理想窗，其时域函数形式如式（2-8）。

（2-8）

为零阶第一类修正的贝赛尔函数，可用式（2-9）的级数表示。

（2-9）

凯泽窗（Kaiser）的幅度公式为：

（2-10）

由公式可知，凯泽窗的β值与贝赛尔函数级数项数n是两个独立的参数，但都会对凯泽窗造成影响，为进一步的了解和研究凯泽窗，有必要分别对其分别介绍。

随着β增大，主瓣加宽，旁瓣幅度减少。

图2-1给出了不同

值所对应的凯泽窗（Kaiser）时域和频率图形（采样点取1024）。

图2-1凯泽窗（Kaiser）时域和频域图形

时域图中，从上往下,

分别为0，1.231，2.341，3.440，4.394，5.440，6.451，7.321，8.500，9.412，10.399，为了方便比较，分别摘取

为0、4.394、8.500、10.399来进行幅度特性对比，明显看出

为一个可自由选择的凯泽窗函数形状参数，能调节主瓣和旁瓣的宽度[12]。

利用MATLAB的求得不同的

值的凯泽窗（Kaiser）详细数据（采样点数定为1024），总结出表2-2。

表2-2不同β值凯泽窗（Kaiser）的详细数据（采样点数为1024）

Kaiser窗函数

泄露系数（%）

旁瓣相对衰减（dB）

主瓣-3dB带宽（rad）

β=0

9.15

-13.3

0.001709

β=1.231

5.48

-15.4

β=2.341

1.65

-20.1

0.0019531

β=3.44