第24章用窗函数法的设计FIR数字滤波器.pptVIP

一、实验目的　　(1)加深对窗函数法设计FIR数字滤波器的基本原理的理解。　　(2)学习用MATLAB语言的窗函数法编写设计FIR数字滤波器的程序。　　(3)了解MATLAB有关窗函数法设计的常用子函数。 二、实验涉及的MATLAB子函数　　1.boxcar　　功能：矩形窗。　　格式：w＝boxcar(n)　　说明：boxcar(n)函数可产生一长度为n的矩形窗函数。 　　2.triang　　功能：三角窗。　　格式：w＝triang(n)　　说明：triang(n)函数可得到n点的三角窗函数。三角窗系数为：当n为奇数时， 　　当n为偶数时，　 3.bartlett　　功能：bartlett(巴特利特)窗。　　格式：w＝Bartlett(n)　　说明：bartlett(n)可得到n点的bartlett窗函数。bartlett窗函数系数为 　　4.hamming　　功能：hamming(哈明)窗。　　格式：w＝hamming(n)　　说明：hamming(n)可产生n点的hamming窗。hamming窗函数系数为　　 　　5.hanning　　功能：hanning(汉宁)窗。　　格式：w＝hanning(n)　　说明：hanning(n)可产生n点的hanning窗。hanning窗函数系数为　　 　　6.blackman　　功能：blackman(布莱克曼)窗。　　格式：w＝blackman(n)　　说明：blackman(n)可产生n点的blackman窗。blackman窗函数系数为　　与等长度的hamming和hanning窗相比，blackman窗的主瓣稍宽，旁瓣稍低。 　　7.chebwin　　功能：chebyshev(切比雪夫)窗。　　格式：w＝chebwin(n，r)　　说明：chebwin(n，r)可产生n点的chebyshev窗函数，其傅里叶变换后的旁瓣波纹低于主瓣r dB。注意：当n为偶数时，窗函数的长度为n＋1。 　　8.kaiser　　功能：kaiser(凯塞)窗。　　格式：w＝kaiser(n，beta)　　说明：kaiser(n，beta)可产生n点的kaiser窗函数，其中，beta为影响窗函数旁瓣的b参数，其最小的旁瓣抑制a与b之间的关系为　　增加b可使主瓣变宽、旁瓣的幅度降低。 　　9.fir1　　功能：基于窗函数的FIR数字滤波器设计——标准频率响应，以经典方法实现加窗线性相位FIR滤波器设计，可设计出标准的低通、带通、高通和带阻滤波器。　　格式：　　※b＝fir1(n，Wn)；设计截止频率为Wn的hamming(哈明)加窗线性相位滤波器，滤波器系数包含在b中。当0≤Wn≤1(Wn＝1相应于0.5fs)时，可得到n阶低通FIR滤波器。　　当Wn＝［W1W2］时，fir1函数可得到带通滤波器，其通带为w1ww2。 　　※b＝fir1(n，Wn，￠ftype￠)；可设计高通和带阻滤波器，由ftype决定：　　·当ftype＝high时，设计高通FIR滤波器；　　·当ftype＝stop时，设计带阻FIR滤波器。　　在设计高通和带阻滤波器时，fir1函数总是使用偶对称N为奇数(即第一类线性相位FIR滤波器)的结构，因此当输入的阶次为偶数时，fir1函数会自动加1。 　　※b＝fir1(n，Wn，Window)；利用列矢量Window中指定的窗函数进行滤波器设计，Window长度为n＋1。如果不指定Window参数，则fir1函数采用hamming窗。　　※b＝fir1(n，Wn，￠ftype￠，Window)；可利用ftype和Window参数设计各种加窗的滤波器。　　由fir1函数设计的FIR滤波器的群延迟为n/2。 三、实验原理　　1.运用窗函数法设计FIR数字滤波器　　FIR数字滤波器的系统函数为　　这个公式也可以看成是离散LSI系统的系统函数： 分母a0为1，其余ak全都为0时的一个特例。由于极点全部集中在0点，稳定和线性相位特性是FIR滤波器突出的优点，因此在实际中具有更广泛的使用价值。　　FIR滤波器的设计任务是选择有限长度的h(n)，使传输函数H(ejw)满足技术要求。主要设计方法有窗函数法、频率采样法和切比雪夫等波纹逼近法等。本实验主要介绍用窗函数法设计FIR数字滤波器。 　　用窗函数法设计FIR数字滤波器的基本步骤是：　　(1)根据过渡带和阻带衰减设计指标选择窗函数的类型，估算滤波器的阶数N。　　(2)由数字滤波器的理想频率响应H(ejw)求出其单位冲激响应hd(n)。　　对于理想的数字低通滤波器频率响应，有下列的子程序可以实现(该程序名为idea