第5章IIR数字滤波器的设计分析.ppt

二、切比雪夫滤波器的设计方法 右图示出了阶N=0,4,5时的切比雪夫多项式特性。 由图可见： (1)切比雪夫多项式的过零点在|x|≤1的范围内； (2)当|x|1时，|CN(x)|≤1,在|x|1范围内具有等波纹性； (3)当|x|1时，CN(x)是双曲线函数，随x单调上升。 2、切比雪夫滤波器（chebyshev）特点： 三、椭圆滤波器（考尔滤波器） 四、小结： 一、由模拟滤波器设计IIRDF步骤： 二、由Ha(s)转换为Hn(z) 遵循原则 5.3.1 脉冲响应不变法（标准z变换法） 2、步骤： 频率非线性变换影响： 5、反畸变（预畸） 总结：双线性变换法设计数字滤波器步骤 方法二：数字域映射 1、频率变换设计原则 满足条件： 2、变换关系导出： 全通函数定义： 3、低通到低通转换 4、低通到高通转换 5、低通到带通转换 6、低通到带阻转换 例题低通到各种滤波器转换 例题：设计一个数字带通滤波器 5.5 数字高通、带通及带阻滤波器的设计 1.数字高通、带通及带阻滤波器的设计方法 2.方法1的一般实现步骤 重点 可利用MATLAB中的buttord、cheb1ord、cheb2ord、ellipord等函数来估计满足性能指标的模拟低通原型阶数和截止频率。 可利用MATLAB中的buttap、cheb1ap、cheb2ap、ellipap等函数来实现。 可利用MATLAB中的lp2lp（LP→LP）、lp2hp、lp2bp、lp2bs等函数来实现。 例5-10 用双线性变换法设计一个ChebyshevⅠ 型数字带通滤波器，设计指标为： Rp=1dB，wp1=0.35π，wp2=0.65π，Rs=40dB， ws1=0.2π，ws2=0.8π，T=1ms。 解：MATLAB程序如下： %确定所需类型数字滤波器的技术指标 Rp=1;Rs=40;T=0.011; wp1=0.4*pi;wp2=0.6*pi;ws1=0.2*pi;ws2=0.8*pi; %将所需类型数字滤波器的技术指标转换成模拟滤 波器的技术指标 wp3=(2/T)*tan(wp1/2);wp4=(2/T)*tan(wp2/2); ws3=(2/T)*tan(ws1/2);ws4=(2/T)*tan(ws2/2); %将所需类型模拟滤波器技术指标转换成模拟低通 滤波器技术指标 wp=[wp3,wp4];ws=[ws3,ws4]; [n,wn]=cheb1ord(wp,ws,Rp,Rs,s); %设计模拟低通滤波器 [z,p,k]=cheb1ap(n,Rp);[b,a]=zp2tf(z,p,k); %频率变换 w0=sqrt(wp3*wp4);Bw=wp4-wp3; [b1,a1]=lp2bp(b,a,w0,Bw); %双线性变换法 [bz,az]=bilinear(b1,a1,1/T); [db,mag,pha,grd,w]=freqz_m(bz,az); plot(w/pi,db); axis([0,1,-50,2]); 仿真波形 以数字高通滤波器为例，给出设计步骤： 1）将数字高通滤波器的技术指标转换为模拟高通 滤波器指标。 2）将模拟高通滤波器的技术指标转换为低通原型 滤波器指标。 3）设计模拟低通滤波器。 4）将模拟低通滤波器映射为数字低通滤波器。 5）将数字低通原型滤波器映射为数字高通滤波器。 将数字低通滤波器转换成所需的各种截止频率的低通、高通、带通和带阻的数字滤波器实质是： 实现从一个z平面到另一个z平面的映射。 N为全通的阶次。 N阶全通 一阶全通 一阶全通 二阶全通 4．频率混叠效应 根据时域采样理论，采样序列h(n)的z变换与模 拟信号ha(t)的拉普拉斯变换之间满足如下关系 将 和 代入得数字滤波器的频率响应 数字滤波器的频响是模拟滤波器频响的周期延 拓。 与时域采样定理吻合！ 图5-11 脉冲响应不变法中的频率混叠现象 只有当模拟滤波器的频响是带限的，且带限于 折叠频率以内，即 有 此时，数字滤波器的频率响应重现模拟滤波器 的频率响应而不产生混叠失真。 但是，任何一个实际的模拟滤波器，其频响都 不可能是真正带限的，这就不可避免的会产生混叠 失真。 答：当模拟滤波器的频响在折叠频率以上处衰 减越大时，混叠失真就越小。 问题：要求混叠 失真小，对模拟滤波 器的频响有何要求？ 5、实际应用中修正 6．优缺点 优点： 1）数字滤波器的单位脉冲响应完全模仿模拟 滤波器的冲激响应，所以