2016年北航数字信号处理上机实验二实验报告.docxVIP

离散时间信号处理实验报告实验二FIR数字滤波器的设计与实现班 级 学 号 姓 名 日 期实验二FIR数字滤波器的设计与实现滤波器设计与实现是数字信号处理中最基本的内容之一。我们知道滤波器分IIR滤波器（无限冲激响应滤波器）和FIR滤波器（有限冲激响应滤波器）两种。FIR滤波器的设计方法以直接逼近所需离散设计系统的频率响应为基础。FIR滤波器可以很容易地获得线性相移特性，不存在不稳定的问题，是实际系统中广为采用的一种数字滤波器。FIR滤波器的设计，通常有窗函数设计法、频率抽样设计法和最佳逼近设计法。窗函数设计法比较简单，它的频率特性是理想滤波器频谱与窗的频谱的卷积，因而，其频率特性取决于窗的类型和长度。频率抽样设计法比较直观，但由于频域的采样会造成时域的混叠，从而滤波器叠性能不可能很高，为提高滤波器的性能，可以在过渡带加上0～1之间的过渡点。本实验主要采用窗函数设计法。FIR滤波器在Matlab中的实现包括时域卷积和filter滤波等方法，通过本实验可以对两种方法的实现过程和实现效果进行分析比较。相关Matlab函数说明1、有限冲激响应数字滤波器设计函数fir1和fir2Matlab函数fir1和fir2可以用来设计加窗的有限冲激响应数字滤波器。两个函数均产生一个线性相位设计。函数fir1可用于设计常规的低通、高通、带通和带阻线性相位有限冲激响应滤波器。对于抽样频率为2Hz的情况，命令b = fir1(N,Wn)在向量b中返回以的升幂排列的N阶低通或带通滤波器的冲激响应系数。对于低通设计，归一化截止频率由标量Wn给定，它是在0和1之间的一个数。对于带通设计，Wn是包含指定通带边界的一个双元素向量[Wn1,Wn2]，其中0Wn1Wn21。命令b = fir1(N,Wn,’high’)其中N为一个偶数，用于设计高通滤波器。命令b = fir1(N,Wn,’stop’)其中Wn是一个双元素向量，用于设计带阻有限冲激响应滤波器。若没有指定，用海明窗作为默认值。命令b = fir(N,Wn,taper)使用在向量taper中给定的长度为N+1的窗系数。然而，窗系数必须事先用一个适当的MATLAB函数来产生，如blackman，hamming，hanning，chebwin或kaiser。用到的命令为taper = blackman(N) taper = hamming(N) taper = hanning(N) taper = chebwin(N) taper = kaiser(N,beta)函数fir2可用于设计带有任意形状幅度响应的线性相位有限冲激响应滤波器。在其基本形状中，命令是b = fir2(N,fpts,mval)它在长度为N+1的向量b中返回以的升幂排列的冲激响应系数。fpts是给定的频率点向量，在0到1范围内以递增的顺序排列，其中第一个频率点是0，而最后的频率点是1。和前面一样，抽样频率被假定为2Hz。mval是给定频率点的给定幅度值的向量，因此也必须和fpts的长度一样，默认使用海明窗。为使用其他窗口，所用到的命令是b = fir2(N,fpts,mval,taper)其中向量taper包含给定的窗系数。2、FIR滤波器的实现(1) FIR滤波器的时域卷积实现滤波器是一种能够让某些频率分量通过而拒绝其他频率成分或对某些频率成分进行修正的系统。对于一个线性时不变系统，其时域输入和输出关系可用卷积的方法描述，。对时域而言，，，均为数字信号，可在数字系统中实现。滤波器按照单位脉冲响应在某一有限区间外是否为零分为有限冲击响应（FIR）滤波器和无限冲击响应（IIR）滤波器。若单位脉冲响应在某一有限区间以外都是零，称为有限冲击响应FIR系统；若单位脉冲响应在某一有限区间外不是零，则称为无限冲击响应IIR系统。在线性时不变系统的卷积描述中，对于FIR滤波器，为有限项，范围从0到M，故。即通过卷积运算，可实现对输入信号的滤波处理。(2) FIR滤波器的filter滤波实现MATLAB中有滤波函数filter，用法为y=filter(b,a,x)。其中b和a为差分方程系数，x为输入序列，y为滤波输出序列。对于FIR滤波器，a设置为1即可。MATLAB实现主要包括两个部分，即滤波器设计和滤波。代码如下：%%fir滤波器系数设计fl_kaiser = [750 1150];fl_mag = [1 0];fl_dev = [0.01 0.05];[fl_n_kaiser,fl_wn,fl_beta,fl_ftype] = kaiserord(fl_kaiser,fl_mag,fl_dev,fs);h = fir1(fl_n_kaiser,fl_wn,fl_ftype,kaiser(fl_n_kaise