实验6数字滤波器设计..docVIP

实验6 数字滤波器设计 设计传输函数G（z），使其频率响应逼近于给定的频率响应指标，这个过程称为数字滤波器的设计。在得到G（z）之后，就以一种适当的滤波器结构的形式实现它。在前面的实验练习中，已经考虑过了有限冲激响应和无限冲激响应传输函数的实现。 在这个实验练习中，我们将学会设计满足某个给定幅度或增益响应的无限冲激响应或有限冲激响应数字滤波器。 无限冲激响应滤波器设计 最常用于设计无限冲激响应滤波器的方法是基于原型模拟传输函数的双线性变换。模拟传输函数通常是下列类型之一：巴特沃斯型，切比雪夫I ，切比雪夫II和椭圆传输函数。这些滤波器类型之间的差别可通过考虑模拟低通滤波器来说明。 （1）无限冲激响应滤波器的阶数的估计 滤波器设计过程中的第一步是，选择接近所使用的滤波器的类型，然后由滤波器指标来估计 传输函数的阶数。用来估计巴特沃斯型滤波器的阶数的MATLAB命令是： [N,Wn]=buttord(Wp,Ws,Rp,Rs) 其中输入参数是归一化通带边界频率wp,归一化阻带频率ws、单位为Db 的通带波纹RP和单位为Db 的J最小阻带衰减纹RS。由于抽样频率被假定为2HZ，WP和WS均必须是0和1之间的一个数。输出数据是满足指标的最低阶数N和归一化截止频率Wn。若RP=3（Db）,则WN=WP。Buttord也可用于估计高通，带通和带阻的巴特沃斯型滤波器的阶数。 （2）无限冲激响应滤波器设计 在选择了滤波器类型并估计了其阶数后，下一步是确定滤波器的传输函数，为此，MATLAB对所有的四种类型都提供了函数。（实验指导书，92） [num,den]=butter(N,Wn) 其中输入参数N和Wn可以使用函数buttord确定，输出是向量num和den，它们分别是以z-1的升幂排列的传输函数的分子和分母多项式的系数。若Wn是一个标量，butter返回一个N阶的低通传输函数，基Wn是一个双元素向量，它返回一个2N阶的带通传输函数。设计N阶巴特沃斯型数字高通滤波器的MATLAB命令是： [num,den]=butter(N,Wn,’high’) 然而，若Wn是一个双元素向量，命令 [num,den]=butter(N,Wn,’stop’) 返回一个2N阶巴特沃斯型带阻滤波器的传输函数。 程序6.1 说明巴特沃斯型带阻滤波器的设计 %巴特沃斯型带阻滤波器的设计 Ws=[0.4 0.6];Wp=[0.3 0.7];Rp=0.4;Rs=50; %估计滤波器阶数 [N1,Wn1]=buttord(Wp,Ws,Rp,Rs); %设计滤波器 [num,den]=butter(N1,Wn1,stop); %显示传输函数 disp(分子系数是);disp(num); disp(分母系数是);disp(den); 分子系数是 Columns 1 through 8 0.0330 0.0000 0.2972 0.0000 1.1889 0.0000 2.7741 0.0000 Columns 9 through 16 4.1611 0.0001 4.1611 0.0000 2.7741 0.0000 1.1889 0.0000 Columns 17 through 19 0.2972 0.0000 0.0330 分母系数是 Columns 1 through 8 1.0000 0.0000 2.6621 0.0000 4.1451 0.0001 4.1273 0.0001 Columns 9 through 16 2.8977 0.0000 1.4381 0.0000 0.5027 0.0000 0.1178 0.0000 Columns 17 through 19 0.0167 0.0000 0.0011 低通巴特沃斯模拟滤波器设计 要求指标 通带截止频率：fp=3400hz,通带最大衰减：Rp=3Db; 阻带截止频率：fs=4000hz,阻带最大衰减：Rs=40db. Matlab程序如下 fp=3400;fs=4000;Rp=3;Rs=40;%输入滤波器的指标 [N,fc]=buttord(fp,fs,Rp,Rs,s);%计算阶数N和3DB截止频率fc [B,A]=butter(N,fc,s);%设计低通巴特沃斯模拟滤波器/传输函数的系数 [hf,f]=freqs(B,A,1024); %计算模拟滤波器的频率响应 plot(