数字滤波器原理

1、IIR数字滤波器设计原理利用双线性变换设计IIR滤波器（只介绍巴特沃斯数字低通滤波器的设计），首先要设计出满足指标要求的模拟滤波器的传递函数，然后由通过双线性变换可得所要设计的IIR滤波器的系统函数。如果给定的指标为数字滤波器的指标，则首先要转换成模拟滤波器的技术指标，这里主要是边界频率的转换，对指标不作变化。边界频率的转换关系为。接着，按照模拟低通滤波器的技术指标根据相应设计公式求出滤波器的阶数和截止频率；根据阶数查巴特沃斯归一化低通滤波器参数表，得到归一化传输函数；最后，将代入去归一，得到实际的模拟滤波器传输函数。之后，通过双线性变换法转换公式，得到所要设计的IIR滤波器的系统函数。步骤及

2、内容1 用双线性变换法设计一个巴特沃斯IIR低通数字滤波器。设计指标参数为：在通带内频率低于时，最大衰减小于；在阻带内频率区间上，最小衰减大于。2 以为采样间隔，绘制出数字滤波器在频率区间上的幅频响应特性曲线。3 程序及图形程序及实验结果如下：%iir_1.m%lskyp%rp=1;rs=15;wp=.2*pi;ws=.3*pi;wap=tan(wp/2;was=tan(ws/2;n,wn=buttord(wap,was,rp,rs,'s'z,p,k=buttap(n;bp,ap=zp2tf(z,p,k;bs,as=lp2lp(bp,ap,wap;bz,az=bilinear(

3、bs,as,.5;h,f=freqz(bz,az,256,1;plot(f,abs(h;title('双线性z变换法获得数字低通滤波器，归一化频率轴'xlabel('omega/2pi'ylabel('低通滤波器的幅频相应'grid;figure;h,f=freqz(bz,az,256,100;ff=2*pi*f/100;absh=abs(h;plot(ff(1:128,absh(1:128;title('双线性z变换法获得数字低通滤波器,频率轴取0,pi/2'xlabel('omega'ylabel('低

4、通滤波器的幅频相应'grid on;运行结果：窗函数法设计FIR数字滤波器原理如果所希望的滤波器的理想频率响应函数为，则其对应的单位脉冲响应为窗函数设计法的基本原理是用有限长单位脉冲响应逼近。由于往往是无限长序列，且是非因果的，所以用窗函数将截断，并进行加权处理，得到：就作为实际设计的FIR数字滤波器的单位脉冲响应序列，其频率响应函数为式中，为所选窗函数的长度。用窗函数法设计的滤波器性能取决于窗函数的类型及窗口长度的取值。设计过程中，要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度。设待求滤波器的过渡带用表示，它近似等于窗函数主瓣宽度。因过渡带近似与窗口长度成反比，

5、决定于窗口形式。例如，矩形窗A=4，海明窗A=8等。按照过渡带及阻带衰减情况，选择窗函数形式。原则是在保证阻带衰减满足要求的情况下， 尽量选择主瓣窄的窗函数。这样选定窗函数类型和窗口长度后，求出单位脉冲响应，再求出。是否满足要求，要进行验算。一般在的尾部加零使长度满足的整数次幂，以便用FFT计算。如果要观察细节，补零点数增多即可。如果不满足要求，则要重新选择窗函数类型和长度，再次验算，直至满足要求。如果要求线性相位特性，则还必须满足：根据上式中的正、负号和长度的奇偶性又将线性相位FIR滤波器分成四类。要根据所设计的滤波特性正确选择其中一类。例如，要设计线性相位低通特性，可选择一类， 而不能选一

6、类。1 根据下列技术指标，设计一个线性相位的FIR数字低通滤波器。通带截止频率，通带允许波动；阻带截止频率，阻带衰减。程序如下：%fir_1.m%lskyp%wp=0.2*pi;N=10;b=fir1(N,wp/pi,boxcar(N+1;t=0:10;subplot(211stem(t,b;title('滤波器的单位脉冲响应'h,w=freqz(b,1,256,1;h=abs(h;subplot(212plot(w,htitle('幅频响应'运行结果：2 写出）中所设计的低通滤波器的单位脉冲响应；并调用fir1(函数得到所设计的低通滤波器的单位脉冲响应，调用f

7、ft(函数进行频响验证。打印输出各部分结果。答：单位脉冲响应： 验证程序：%fir_2.m%lskyp%wp=0.2*pi;N=10;b=fir1(N,wp/pi,boxcar(N+1;t=0:10;subplot(211stem(t,b;title('滤波器的单位脉冲响应'h=fft(b,16;h=abs(h;subplot(212f=linspace(0,0.5,8plot(f,h(1:8title('用fft验证'运行结果：3 编程验证窗长和窗形状对实际滤波器性能的影响。如要求用窗函数法设计一个线性相位FIR数字低通滤波器，用理想低通滤波器作为逼近滤波器，

8、截止频率，用四种窗函数（矩形窗，汉宁窗（升余弦窗），哈明窗（改进的升余弦窗），布莱克曼窗）设计该滤波器，选择窗函数的长度两种情况。程序如下：%fir_3.m%lskyp%wc=pi/4;N=15;window=boxcar(N+1;b=fir1(N,wc/pi,window;h,w=freqz(b,1,256,1h=abs(h;subplot(221plot(w,hgrid on;title('矩形窗，N=15'N=33;window=boxcar(N+1;b=fir1(N,wc/pi,window;h,w=freqz(b,1,256,1h=abs(h;subplot(223p

9、lot(w,hgrid on;title('矩形窗,N=33'N=15;window=hanning(N+1;b=fir1(N,wc/pi,window;h,w=freqz(b,1,256,1h=abs(h;subplot(222plot(w,hgrid on;title('汉宁窗,N=15'N=33;window=hanning(N+1;b=fir1(N,wc/pi,window;h,w=freqz(b,1,256,1h=abs(h;subplot(224plot(w,hgrid on;title('汉宁窗,N=33'figure(2wc=pi

10、/4;N=15;window=hamming(N+1;b=fir1(N,wc/pi,window;h,w=freqz(b,1,256,1h=abs(h;subplot(221plot(w,hgrid on;title('汉明窗，N=15'N=33;window=hamming(N+1;b=fir1(N,wc/pi,window;h,w=freqz(b,1,256,1h=abs(h;subplot(223plot(w,hgrid on;title('汉明窗,N=33'N=15;window=blackman(N+1;b=fir1(N,wc/pi,window;h,w=freqz(b,1,256,1h=abs(h;subplot(222plot(w,hgrid on;title(