无限长单位脉冲响应滤波器的设计方法.PPTVIP

4、N阶特性阶数N等于通带内最大和最小值个数的总和。可由幅频特性中看出N阶数。且当：N=奇数，则Ω=0处有一最大值，N=偶数，则Ω=0处有一最小值。N阶公式N值是根据阻带的边界条件来确定的切比雪夫Ⅰ型与巴特沃斯低通的A2(Ω)曲线对比3.3从模拟滤波器低通原型到各种数字滤波器的频率变换实际应用中的数字滤波器有低通、高通、带通和带阻等类型。设计各种类型的数字滤波器通常可以把一个归一化的原型模拟低通滤波器经模拟频带变换成所需类型的模拟滤波器，再通过脉冲响应不变法或双线性变换法转换为所需类型的数字滤波器。如图所示：归一化模拟低通模拟低通、高通、带通、带阻数字低通、高通、带通、带阻模拟-模拟频带变换双线性变换或脉冲响应不变法模拟-模拟-数字归一化模拟低通数字低通、高通、带通、带阻模拟-数字频带变换模拟-数字设计IIR滤波器的两种频率变换法3.3.1低通变换（模拟低通-数字低通）首先将数字低通滤波器的性能要求转换为与之相对应的模拟低通滤波器的性能要求，根据此性能要求设计模拟低通滤波器。然后通过脉冲响应不变法或双线性变换法，将此模拟低通滤波器Ha(s)数字化为所需的数字低通滤波器。可见，H(z)与采样周期T有关，T越小，H(z)的相对增益越大，这是不希望的。为此，实际应用脉冲响应不变法时稍作一点修改，即求出H(z)后，再乘以因子T,使H(z)只与fc和fs的相对值有关，而与采样频率fs无直接关系。例如，与的数字滤波器具有相同的传递函数，这一结论适合于所有的数字滤波器设计。最后得：Matlab程序%巴特沃思滤波器[B,A]=butter(3,2*pi*1000,s);%[b,a]=butter(n,Wn,‘s’)，n为滤波器的阶数Wn为边界频率，按s的降幂排列%脉冲响应不变法[num1,den1]=impinvar(B,A,4000);%4000为采样频率[h1,w]=freqz(num1,den1);[B,A]=butter(3,2/0.00025,‘s’);％2/0.00025预畸变模拟滤波器边界频率[num2,den2]=bilinear(B,A,4000);%双线性法[h2,w]=freq(num2,den2);f=w/pi*2000;plot(f,abs(h1),-.,f,abs(h2),-);grid;xlabel(频率/Hz);ylabel(幅值/dB);频率/Hz三阶巴特沃兹滤波器的频率响应幅度/dB上图为两种设计方法所得到的频响，对于双线性变换法，由于频率的非线性变换，使截止区的衰减越来越快，最后在折叠频率处形成一个三阶传输零点,这个三阶零点正是模拟滤波器在处的三阶传输零点通过映射形成的。因此,双线性变换法使过渡带变窄,对频率的选择性改善,而脉冲响应不变法存在混淆,且没有传输零点。3.3.2高通变换（模拟低通-数字高通）fswsfpwp1ApAsfw|H(ejw)|通带截止频率：fp(wp)又称为通带下限频率。通带衰减：Ap阻带截止频率：fs(ws)又称阻带上限截止频率。阻带衰减：As1、高通滤波器的性能指标由于倒数关系不改变模拟滤波器的稳定性，因此，也不会影响双线变换后的稳定条件映射到即映射到即图1高通变换频率关系这一曲线的形状与双线性变换时的频率非线性关系曲线相对应，只是将坐标倒置，因而通过这一变换后可直接将模拟低通变为数字高通,如图2。1.01.003.2常用的模拟低通滤波器特性1、由幅度平方函数确定模拟滤波器的系统函数2阶2阶例设已知求对应的零极点分布取左半平面的零极点构成Ha(s)3.2.1巴特沃思滤波器不管