「课程设计-基于Matlab的FIR滤波器设计与仿真」.pdfVIP

基于Ｍatlａb的FIR滤波器设计与仿真

实验地点：实验日期：7月13日至7月18日

一、实验目的

1、掌握基本的ＭATＬAB编程方法；

2、理解FIR滤波器的设计原理；

3、学会用MＡTLAB来编程实现FIR滤波器;

4、掌握基本的simulink交互式仿真，并对FIR滤波器模型并进行仿真；

5、学会对所得的结果进行分析。

二、实验内容

１、制作数据源：

用电脑采集或用软件截取5至10秒的语音(如“宁波大学”.．.）(注意用wavrｅａｄ

函数

将语音文件读入到ｍａtlab时，有“右声道”和“左声道”两个信号,只要将其中一个声道

作为信号就可以。)

２、信号中混入随机噪声(注意信噪比，噪声强度不要太大):

事先取一个参考值为0,再取噪声的方差，方差的取值是根据原始语音的幅度取的噪声,

噪声强度不宜过大或过小,适中即可。然后通过rａｎdn（size(a）).*sqｒt(noise_

ｖａr)＋nｏise_mu；语句产生与原始语音长度相同的随机噪声，最后把这个噪声添加到

原始语音中，得到加噪语音信号。

3、设计一个FＩＲ滤波器对前面的混合信号进行处理、降低噪声，并对性能进行分析:

有限长单位冲激响应(FＩR)数字滤波器具有严格的线性相位，又具有任意的幅频特性。

同时FIR系统只有零点,系统是稳定的，因而容易实现线性相位和允许实现多通道滤波器。

只要经过一定的延时,任何非因果有限长序列都能变成因果的有限长序列,因而总能用因

果系统来实现。FIR滤波器由于单位冲激响应是有限长的,可以用快速傅立叶变换(FFT)

算法来实现过滤信号,从而大大提高运算效率。由于ＦIR滤波器具有以上优点,在信号处理

和数据传输中得到了广泛的应用。运用窗函数设计滤波器是FIR数字滤波器设计的主要方法

之一,由于运算简单,又有闭合形式的公式

可循,因而很实用。

有限长FIR数字滤波器的设计方法主要

是窗函数设计法。

常用的窗函数有以下几种,矩形窗，三

角窗，汉宁窗,海明窗,布拉克曼窗,凯

泽窗(β=7.865),各种窗函数基本参

数比较如下表１。由表１得矩形窗虽然

简便，但存在明显的吉布斯效应，主瓣和

第一旁瓣之比只有1３dB,为了克服这些

缺陷,设计时可根据给定的滤波器技术指

标，选择滤波器长度M和窗函数w(n),使

其具有最窄宽度的主瓣和最小旁瓣

4、用Simuｌink设计FIR滤波器模型并进行仿真，所得结果进行分析：

Sｉｍulink创建模型并寻找解决实际问题的方法，可以用它轻松地搭建一个系统模型,

并设置模型参数和仿真参数。由于Simulink是交互式的应用程序，因此在仿真过程中，可

以在线修改仿真参数,并立即观察到改变后的仿真结果。

本文将介绍在Ｓimulink环境下应用数字滤波器设计模块进行信号处理（滤波）的可视

化设计技术。由于调用了Ｍaｔｌaｂ高度集成化和可视化的计算机仿真环境，用户可以不用

编写任何程序语句，仅仅通过鼠标拖放和点击操作,就可完成复杂的数字滤波器的设计工作。

下面通过一个应用示例介绍这种基于可视化的信号滤波技术,用数字滤波器设计模块（Diｇ

itaｌFiltｅrＤesignBlock)设计、分析和实现滤波器。

三、实验结果与分析

1、制作数据源：

我用电脑录制了同学的一段语音,内容是：现代语音处理技术。

该语音大概有4，5秒之长。然后用wａvread函数将语音文件读入到ｍatlab中，

其实现的语句如下：

[ａ,fｓ］＝ｗａｖreａd(现代语音处理技术.wａv);％读入语音

a=a(:,1）;

从中我们知道fｓ语音的采样频率。

2、信号中混入随机噪声：

ｎoｉｓe＿mｕ=0;

noｉｓe_vａr=０.005；

ａ０=ｒandn(s