FIR音频滤波.doc

音频滤波：去噪 FIR系统的系统函数 系数b0,b1,…,bM即为系统的单位抽烟响应h(0)，h(1)，…，h(M)，且当nM时，h(n)=0。 FIR滤波器的设计主要是建立在对理想滤波器频率特性作某种近似的基础上，这些近似方法有窗函数法、频率抽样法以及等波纹切比雪夫逼近法。 窗函数法设计的基本思想是把给定的频率响应通过，求得脉冲响应，然后利用加窗函数对它进行截断和平滑，以实现一个物理可实现且具有线性相位的 FIR 数字滤波器的设计目的。其核心是从给定的频率特性，通过加窗确定有限长单位取样响应h(n)。频率采样法设计的基本思想是，把给出的理想频率响应进行取样，通过 IDFT 从频谱样点直接求得有限脉冲响应。等波纹切比雪夫逼近法，利用 MATLAB 提供的 remez 函数实现Remez算法，设计滤波器逼近理想频率响应。 一、窗函数法 1.1 设计原理 设计FIR数字滤波器的最简单的方法是窗函数法，通常也称之为傅立叶级数法。 窗函数法设计的基本思想是，加窗截取理想滤波器的单位脉冲响应的一段作为所要设计的滤波器的脉冲响应。其核心是从给定的频率特性，通过加窗确定有限长单位取样响应h(n)。设计过程如下： 加窗的作用是通过把理想滤波器的无限长脉冲响应hd(n)乘以窗函数w(n)来产生一个被截断的脉冲响应，即h(n)= hd(n)w(n)，并且对频率响应进行平滑。 用窗函数法设计FIRDF的具体设计步骤如下： （1）选择窗函数类型和长度，写出才窗函数的表达式。 根据阻带最小衰减选择窗函数的的类型，再根据过渡带宽度确定所选窗函数的长度。用窗函数法设计的FIRDF通带波纹幅度近似等于阻带波纹幅度。一般阻带最小衰减达到40dB以上，则通带最大衰减就小于0.1dB。所以用窗函数法设计FIRDF时，通常只考虑阻带最小衰减就可以了。 （2）构造希望逼近的频率响应函数 。 根据设计需要，一般选择线性相位理想滤波器（理想低通、理想高通、理想带通、理想带阻）。理想滤波器的截止频率近似为最终设计的FIRDF的过渡带中心频率，幅度函数衰减一半（约-6dB）。所以一般取 ，和分别为通带边界频率和阻带边界频率。 （3）计算。 （4）加窗得到设计结果：。 常见的窗函数有矩形窗、bartlett窗、三角窗、Hamming窗、Hanning窗、Blackman窗、Chebyshey窗、Kaiser窗等等。矩形窗具有最窄的主瓣，但也有最大的边瓣峰值和最慢的衰减速度。Hamming窗和Hanning窗主瓣稍宽，但有较小的边瓣和较大的衰减速度，是较为常用的窗函数。 1.2 实验流程 1.2.1 获取wav文件并加噪 wav音乐文件的获取，可以自行采集，也可以从网上下载。 然后在Matlab软件平台下，利用函数wavread对语音信号进行采样，记住采样频率和采样点数。通过wavread函数的使用，让我们很快理解了采样频率、采样位数等概念。采集完成后在信号中加入一个单频噪声。(对应程序见附录5.1) 接下来，输出显示原始信号和加噪信号的时域图和频谱图。(程序见附录5.2) 图1-1 原始信号和加噪信号的时域图和频谱图 1.2.2 带阻滤波器概念 根据前面所加的噪声知道，这里需要使用带阻滤波器来滤除加上的噪声。 数字带阻滤波器也具有频率响应的周期性，频率变量以数字频率ω来表示(ω=ΩT=Ω/fs，Ω为模拟角频率，T为抽样时间间隔，fs为抽样频率)。所以，数字滤波器设计中，必须给出抽样频率。下图显示数字带阻滤波器的理想幅频响应，这样的理想频率响应是不太可能实现的，原因是频带之间幅度响应是突变的，因而其单位抽样响应是非因果的。 一般来说，滤波器的性能要求往往以频率响应的幅度特性的允许误差来表征。以低通滤波器为例，如下图所示，频率响应有通带、过渡带和带阻三个范围(非理想)。 在通带内幅度响应以误差α1逼近于1，即 在带阻中，幅度响应以误差α2逼近于0，即 其中ωc，ωst分别为通带截止频率和阻带截止频率。为了逼近理想低通滤波器的特性，还必须有一个非零宽度(ωc-ωst)的滤波带，在这个过渡带内的频率响应平滑地从通带下降到阻带。 1.2.3 滤波器设计 在该滤波器的设计中，我们给出该滤波器的性能指标如下： fpd=1800; fsd=1850; fsu=1950; fpu=2000; Rp=1; As=40; 截止频率也可以任意自选，在单频噪声干扰附近即可。在这里，很重要的是通带截止db值的设置。这个值一定要根据我们使用的设计滤波器的方法来设定。因为我使用的是NUTTALLWIN窗法，NUTTALLWIN窗函数中，滤波器的过渡带宽为15.4，最小阻带衰减为108db。所以，一定要将通带截止db值设置的小于108，所以，我将其设置为100db。 在这里我是使用窗函数法设计上面要求