《信号分析与处理》_第 9 章 滤波器设计.pptVIP

9.4FIR数字滤波器设计3.窗函数法设计步骤用窗函数设计FIR滤波器可分为如下几个步骤。(1)首先由给定滤波器的理想频率响应Hd(ejω),利用傅里叶逆变换,求理想滤波器的单位样值响应,即(2)根据过渡带以及阻带衰减的要求,参考表9--3的基本参数,选择窗函数的形式w[n],并估计窗口的大小N;(3)计算滤波器的单位样值响应h[n]=w[n]hd[n];(4)检验设计出的数字滤波器是否满足技术指标,若不满足所需性能指标,重复步骤(2)至(4),直到满足指标为止。上一页下一页返回9.4FIR数字滤波器设计9.4.4频率采样法设计FIR滤波器频率采样法是从频域出发,把给定的理想频率响应Hd(ejω)等间隔采样N个点得到再对Hd(k)作IDFT,得到为了保证所设计的FIR滤波器具有线性相位,采样频率采样法设计FIR滤波器时需满足如表9-4所列的约束条件。上一页下一页返回9.4FIR数字滤波器设计9.4.5FIR滤波器与IIR滤波器的比较FIR滤波器和IIR滤波器两种滤波器在其自身特点和设计方法上都是各有所长,在此对他们做简单比较,可作为实际设计滤波器时的参考。从滤波器设计来看,IIR滤波器的最大优点是可以利用模拟滤波器的设计技术,而模拟滤波器的设计有大量图表可查,设计技术成熟,所以IIR滤波器的设计比较简单,主要用于设计具有片断常数特性的数字滤波器,但是它在设计过程中往往脱离不了模拟滤波器的模式影响,此外IIR滤波器的不足之处是不容易实现线性相位。FIR滤波器在满足单位样值响应对称性的情况下可以实现严格的线性相位。上一页下一页返回9.4FIR数字滤波器设计从实现成本来看,IIR滤波器可以用较低的阶数实现较高的技术指标,设计成本较低。对于同样的滤波器设计指标,FIR滤波器要求的阶数通常要比IIR滤波器的阶数高5至10倍,因此FIR滤波器的设计成本较高。从滤波器结构上看,IIR滤波器一般采用递归结构,为使系统稳定,极点必须全部位于单位圆内,这种结构在有限字长处理是容易引起较大误差。FIR滤波器一般不采用递归结构,总是稳定的系统,因此误差很小。从运算角度看,FIR滤波器可以利用FFT算法实现,运算速度很快,运算量很小。IIR滤波器则不可能这样计算。总之,FIR滤波器和IIR滤波器各有所长,在实际运用时需多方面考虑。在对相位特性要求不高的场合可以偏向于选择IIR滤波器,对于要求线性相位的场合则最后选择FIR滤波器。上一页返回9.5本章相关的MATLAB命令及应用9.5.1IIR数字滤波器设计【例96】设计一个butterworth低通滤波器,通带截止频率为100Hz,阻带截止频率为200Hz,通带衰减rp小于2dB,阻带衰减rs大于15dB,取样频率为fs=1500Hz。(1)请用冲激响应不变法设计,并作出其频响特性图;(2)请用双线性法设计,并作出其频响特性图;1)冲激响应不变法程序clc;clear;closefp=100;fss=200;fs=1500;rp=2;rs=15;wp1=2*pi*fp/fs;%求出通带数字角频率wp1wss1=2*pi*fss/fs;%求出阻带数字角频率wss1下一页返回9.5本章相关的MATLAB命令及应用[N,wc]=buttord(2*pi*fp,2*pi*fss,rp,rs,s)%求出最小阶数N和截止频率wc[z,p,k]=buttap(N);[ba,aa]=zp2tf(z,p,k)%求出归一化模拟滤波器的分子分母系数ba和aa[b,a]=lp2lp(ba,aa,wc)%求出截止频率为wc的模拟滤波器的分子分母系数b和a[bz,az]=impinvar(b,a,fs)[H,f]=freqz(bz,az);ma=20*log10(abs(H));mp=angle(H)*180/pi;subplot(2,1,1);plot(f/pi,ma)subplot(2,1,2);plot(f/pi,mp)上一页下一页返回9.5本章相关的MATLAB命令及应用2)双线性法程序clc;clear;closefp=100;fss=200;fs=1500;rp=2;rs=15;wp1=2*pi*fp/fs;%求出通带数字角频率wp1w