实验三 IIR滤波器的设计实验.doc

实验三、 无限长单位脉冲响应(IIR)滤波器的设计方法实验目的无限冲激响应（ IIR ）数字滤波器的常用指标和设计过程。截止)。其中，Wp 和Ws分别是滤波器的通带截止频率和阻止截止频率，单位为rad/s；Rp和Rs分别是通带最大衰减系数和阻带最小衰减系数，单位为dB。 （2）用于计算N阶巴特沃思归一化（=1）模拟低通滤波器系统函数的零、极点和增益因子，返回长度为N的向量z和p分别给出N个零点和极点，G是滤波器增益。得到的滤波器系统函数形式如下： 其中，和分别是向量z和p的第k个元素。如果要从零、极点得到系统函数的分子和分母多项式系数向量B和A，可以调用结构转换函数 （3），结构转换后系统函数的形式为 其中，M 是向量B的长度，N是向量A的长度，分别是向量B和A的第k个元素。 .（3） 用于计算巴特沃思模拟滤波器系统函数中分子和分母多项式系数向量B和A，其中N和分别是滤波器的阶和3dB截止频率，返回向量B和A中的元素和分别是上面的表示式中的分母和分子系数。ftype缺省时，设计低通滤波器;ftype=high时，设计高通滤波器；ftype=stop时，设计带阻滤波器，此时为向量，且ftype缺省时，设计带通滤波器，带通的频率区间为。S缺省时，设计数字滤波器。 例如：设计一个满足下列指标要求的巴特沃思模拟低通滤波器。指标 Matlab程序： 运行结果如下： N=6 Wc=2.2496 B.= 0 0 0 0 0 0 129.5917 A=1.0000 8.6916; 37.7720; 104.0667;191.1447 ;222.5973 ; 129.5917 即： 幅频响应曲线如图6.3.4所示。 2、Matlab中设计数字滤波器的函数都是采用双线性变换法，将模拟滤波器转换为数字滤波器。这些函数及其调用格式如下（巴特沃思数字滤波器）： （1） 该格式用于计算巴特沃思数字滤波器的阶N和3dB截止频率的归一化值（关于归一化）。调用参数Wp 和Ws分别是数字滤波器的通带截止频率和阻带截止频率的归一化值（关于归一化），要求和，其中1表示数字频率（对应模拟频率，为采样频率）。和分别是通带最大衰减和阻带最小衰减，单位dB。当时，设计高通滤波器；当和是二元向量时，设计带通（）或带阻（）滤波器，这时返回参数也是二元向量。 （2） 该格式用于计算巴特沃思模拟滤波器系统函数中分子和分母多项式系数向量B和A。调用参数N 和WC分别是巴特沃思数字滤波器的阶和3dB截止频率归一化值（关于归一化）。当ftype缺省时，设计低通滤波器；当ftype=high时，设计高通滤波器；ftype=stop时，设计带阻滤波器，此时为二元向量，和分别是带阻滤波器的通带3dB下，上截止频率的归一化值（关于归一化）；为向量，且ftype缺省时，设计带通滤波器，带通的频率区间为。注意设计出的带通和带阻数字滤波器是2N阶的，这是因为带通滤波器可表示为一个N阶低通滤波器与一个N阶高通滤波器的级联。 由函数的返回向量B和A可写出数字滤波器的系统函数为： 其中，分别是向量B和A的第k个元素。 例如：用双线性变换法设计巴特沃思数字低通滤波器，指标要求 通带截止频率，通带最大衰减； 阻带截止频率，阻带最小衰减 解 设计步骤如下： （1）给出数字滤波器的指标。 （2）计算巴特沃思数字滤波器的阶N和3dB截止频率。 ； （3）用双线性变换法设计巴特沃思数字低通滤波器。 运行结果如下： =0.0335 0.1006 0.1006 0.0335 =1.000 -1.4245 0.8827 -1.1900 设计出的数字滤波器系统函数为 程序如下： 例：利用脉冲响应不变法，把系统函数为的模拟滤波器变换成等价的数字滤波器，采样间隔 手算： Matlab: B=[1,1]; A=[1,5,6]; T=0.1; Fs=1/T; [Bz,Az]=impinvar(B,A,Fs);%用脉冲响应不变法将模拟滤波器变换成数字滤波器 运行结果： Bz=0.1000 -0.0897 Az=1.000 -1.5595 0.6065 1.设采样周期T=250μs（采样频率fs =4kHz），用脉冲响应不变法和双线性变换法设计一个三阶巴特沃兹滤波器，其3dB边界频率为fc =1kHz。 [B,A]=butter(3,2*pi*1000,s); [Bz,Az]=impinvar(B,A,4000); [h1,w]=freqz(Bz,Az); [B,A]=butter(3,2/0.00025,s); [Bz1,Az1]=bilinear