实验4 基于MATLAB的FIR数字滤波器设计.docVIP

. .. 实验4 基于MATLAB的FIR数字滤波器设计 实验目的：加深对数字滤波器的常用指标和设计过程的理解。 实验原理：低通滤波器的常用指标： 通带边缘频率，阻带边缘频率 ，通带起伏， 通带峰值起伏， 阻带起伏，最小阻带衰减。 数字滤波器有IIR和FIR两种类型，它们的特点和设计方法不同。 在MATLAB中，可以用b=fir1(N,Wn,’ftype’,taper) 等函数辅助设计FIR数字滤波器。N代表滤波器阶数；Wn代表滤波器的截止频率(归一化频率)，当设计带通和带阻滤波器时，Wn为双元素相量；ftype代表滤波器类型，如’high’高通，’stop’带阻等；taper为窗函数类型，默认为海明窗，窗系数需要实现用窗函数blackman, hamming,hanning chebwin, kaiser产生。 例1 用凯塞窗设计一FIR低通滤波器，通带边界频率 ，阻带边界频率 ，阻带衰减 不小于50dB。 解 首先由过渡带宽和阻带衰减 来决定凯塞窗的N和 ， ， ? 上图给出了以上设计的频率特性，(a) 为N=30直接截取的频率特性(b)为凯塞窗设计的频率特性。凯塞窗设计对应的MATLAB程序为： wn=kaiser(30,4.55); nn=[0:1:29]; alfa=(30-1)/2; hd=sin(0.4*pi*(nn-alfa))./(p