基于MATLAB的FIR滤波器的程序滤波,课程设计-基于Matlab的FIR滤波器设计与仿真

基于Matlab的FIR滤波器设计与仿真

实验地点：曹光彪信息楼219 实验日期：7月13日至7月18日

一、实验目的

1、掌握基本的MATLAB编程方法；

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

3、学会用MATLAB来编程实现FIR滤波器；

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

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

二、实验内容

1、制作数据源：

用电脑采集或用软件截取5至10秒的语音(如“宁波大学”...)(注意用wavread函数将语音文件读入到matlab时，有“右声道”和“左声道”两个信号，只要将其中一个声道作为信号就可以。)

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

事先取一个参考值为0，再取噪声的方差，方差的取值是根据原始语音的幅度取的噪声，噪声强度不宜过大或过小，适中即可。然后通过 randn(size(a)) .* sqrt(noise_var) + noise_mu;语句产生与原始语音长度相同的随机噪声，最后把这个噪声添加到原始语音中，得到加噪语音信号。

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

有限长单位冲激响应( FIR) 数字滤波器具有严格的线性相位,又具有任意的幅频特性。同时FIR 系统只有零点,系统是稳定的,因而容易实现线性相位和允许实现多通道滤波器。只要经过一定的延时,任何非因果有限长序列都能变成因果的有限长序列, 因而总能用因果系统来实现。FIR 滤波器由于单位冲激响应是有限长的,可以用快速傅立叶变换( FFT) 算法来实现过滤信号,从而大大提高运算效率。由于FIR 滤波器具有以上优点,在信号处理和数据传输中得到了广泛的应用。运用窗函数设计滤波器是FIR 数字滤波器设计的主要方法之一,由于运算简单,又有闭合形式的公式可循,因而

很实用。

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

是窗函数设计法。

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

窗,汉宁窗, 海明窗, 布拉克曼窗, 凯泽

窗(β =7. 865) ,各种窗函数基本参数比

较如下表1 。由表1得矩形窗虽然简便，

但存在明显的吉布斯效应，主瓣和第一旁

瓣之比只有13dB,为了克服这些缺陷，设

计时可根据给定的滤波器技术指标，选择

滤波器长度M和窗函数w(n)，使其具有最

窄宽度的主瓣和最小旁瓣