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

* 双线性变换法适合分段常数特性滤波器的设 计。分段常数型模拟滤波器经变换后仍为分段常数 型数字滤波器，但临界频率点产生畸变。 这种频率的畸变，可通过频率的预畸加以校正。 * * 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); %双线性变换 [bz,az]=bilinear(b,a,1/T); [db,mag,pha,grd,w]=freqz_m(bz,az); plot(w/pi,db);axis([0,1,-30,2]) * * 5.利用AF设计IIR数字低通滤波器的步骤 1 .确定数字低通滤波器（LPDF）的性能指标; 2 .将LPDF指标转换为LPAF的性能指标（主要指 边界频率）; 脉冲响应不变法： 双线性变换法： 3.设计LPAF; 4.将 * 5.5 数字高通、带通及带阻滤波器的设计 2.方法1的一般实现步骤 可利用MATLAB中的buttord、cheb1ord、cheb2ord、ellipord等函数来估计满足性能指标的模拟低通原型阶数和截止频率。 可利用MATLAB中的buttap、cheb1ap、cheb2ap、ellipap等函数来实现。 可利用MATLAB中的lp2lp（LP→LP）、lp2hp、lp2bp、lp2bs等函数来实现。 * 例5-13 用双线性变换法设计一个ChebyshevⅠ型 数字带通滤波器，设计指标为： Rp=1dB，wp1=0.4πrad，wp2=0.6πrad，Rs=40dB，ws1=0.2πrad，ws2=0.8πrad，T=1ms。 解：MATLAB程序如下： %确定所需类型数字滤波器的技术指标 Rp=1;Rs=40;T=0.001; wp1=0.4*pi;wp2=0.6*pi;ws1=0.2*pi;ws2=0.8*pi; * * 一般以 为 参考频率进行 归一化 4.巴特沃思滤波器的图表法设计 （1）将 、 进行归一化处理，得到 、 （2）由图5-5，得阶数N 通 带 内 衰 减 阻 带 内 衰 减 * （3）查表，得H(p)的分母多项式 * （4）去归一化 例5-3 利用图表法设计例5-2所述的巴特沃斯滤波器。 * 5.巴特沃思滤波器的MATLAB实现 在MATLAB信号处理工具箱中，与Butterworth滤波器 相关的有三个函数： （1）[n,wn]=buttord(wp,ws,Rp,Rs,‘s’)，求出给定通带 截止频率wp、阻带截止频率ws、通带最大衰减Rp、阻带最 小衰减Rs条件下，所需要的最小Butterworth模拟滤波器阶 数n和3dB截止频率wn。 （2）模拟低通原型滤波器设计函数buttap，调用格式 为[z,p,k]=buttap(n)，其中参数z、p、k分别为滤波器的零 点、极点和增益。 * （3）Butterworth滤波器完全设计函数，调用格式为 [b,a]=butter(n,wn,‘ftype’,‘s’)，其中’ftype’为滤波器类型，缺省时表示低通或带通滤波器。b、a分别为滤波器传输函数 分子、分母多项式系数向量。 例5-4 利用MATLAB实现例5-2所述巴特沃斯滤波器。 解 MATLAB实现程序如下： wp=2*pi*5000; ws=2*pi*10000; Rp=3; Rs=20; %计算满足性能指标的滤波器阶数n和3dB截止频率wn [n,wn]=buttord(wp,ws,Rp,Rs,s) * %设