第章无限脉冲响应数字滤波器的设计IIR.ppt

* Ch.6 Design of IIR * 函数butter的调用格式 函数butter的调用格式为 [b, a]= butter(n, Wc,) % 设计数字Butterworth滤波器 [b, a]= butter(n, Wc, ftype) % 设计模拟Butterworth滤波器 其中，n为滤波器阶数，Wc为截止频率。 * Ch.6 Design of IIR * 函数cheby1的调用格式 函数cheby1的调用格式为 [b, a]= cheby1(n, Rp, Wc) % 设计数字Chebyshev滤波器 [b, a]= cheby1(n, Rp, Wc, ftype) % 设计模拟Chebyshev滤波器 其中，n为滤波器阶数， Rp为通带内的纹波系数，Wc为截止频率。 * Ch.6 Design of IIR * 例：设计ButterWorth低通滤波器 例5.8 设计一模拟butterworth低通滤波器，通带截止频率300Hz，通带最大衰减2dB，阻带截止频率800Hz，阻带最小衰减30dB。 解 滤波器的阶数和截止频率可由式和确定，程序段为 Wp= 2*pi*300;Ws= 2*pi*800;Rp= 2;Rs= 30; N= ceil((log10((10^(0.1*Rs)-1)/(10^(0.1*Rp)-1)))/(2*log10(Ws/Wp))); Wc= Wp/((10^(Rp/10)-1)^(1/(2*N))); [b,a]= butter(N, Wc, s); freqs(b,a) * Ch.6 Design of IIR * 例5.8 程序运行结果 运行程序，得到N= 4，Wc= 2.0157e+003。幅频特性和相频特性如图所示。 * Ch.6 Design of IIR * 5.6.2 模拟滤波器到数字滤波器的转换 设模拟滤波器系统函数为 数字滤波器的系统函数为 从模拟滤波器到数字滤波器的转换有两种方法，即脉冲响应不变法和双线性变换法。 * Ch.6 Design of IIR * 脉冲响应不变法 脉冲响应不变法: 用代换Ha(s)中的(s-sk)即可得到H(z)，从而将模拟滤波器转换为数字滤波器格式。 可用函数impinvar实现，调用格式为 [bz, az]= impinvar(b, a, fs) 其中，fs为取样频率。 * Ch.6 Design of IIR * 双线性变换法 双线性变换法: 用代换Ha(s)中的s即可得到H(z)，从而将模拟滤波器转换为数字滤波器格式。 可用函数bilinear实现，调用格式为 [zd, pd, kd]= bilinear(z, p, k, fs) 其中，z，p，k和zd，pd，kd分别为s域和z域系统函数的零点、极点和增益。 * Ch.6 Design of IIR * 例5.9 利用impinvar将一模拟低通滤波器变换成数字滤波器(取样频率为10Hz)， 程序段为 [b,a]= butter(4,.3,’s’); [bz,az]= impinvar(b,a,10); 程序运行结果为 bz=1.0e-006 * -0.0000 0.1324 0.5192 0.1273 0 az=1.0000 -3.9216 5.7679 -3.7709 0.9246 例：模拟滤波器转换数字滤波器 * * * * * 5.预畸变 例5-11：设计一个一阶巴特沃思数字低通滤波器，其通 带截止频率为200Hz，通带最大衰减为3dB，采样频率为 1000Hz。 * 例5-11：巴特沃思数字低通，通带截止频率为200Hz，通 带最大衰减为3dB，采样频率为1000Hz。 * 预畸变就是将临界模拟频率事先加以畸变，然 后经变换后正好映射到所需要的数字频率上。 * 例5-12 已知fp=0.3kHz，Rp=1dB，fs=0.2kHz， Rs=20dB，T=1ms, 利用双线性变换法设计一个Chebyshev I型数字高通滤波器。 解 MATLAB实现程序如下： Rp=1;Rs=20;T=0.001;fp=300;fs=200; wp=2*pi*fp*T;ws=2*pi*fs*T; %预畸变 wp1=(2/T)*tan(wp/2);ws1=(2/T)*tan(ws/2); * %设计模拟滤波器 [n,wn]=cheb1ord(wp1,ws1,Rp,Rs,s); [b,a]=cheby1(n,Rp,wn,high,s