数字信号处理-IIR系统设计.pptVIP

44阶remez滤波器 FIR系统与IIR系统的对比 6阶elliptic滤波器 MATLAB中的设计演示工具：filtdemo * * 第八章 IIR系统的设计 IIR系统的设计概述 原型滤波器与模拟系统设计 模拟系统与数字系统的变换 数字信号处理 模拟系统与数字系统的变换方式 从数字滤波器频谱到模拟滤波器频谱； 从模拟滤波器系统到数字滤波器系统。 冲激响应不变法 双线性变换法 冲激响应不变法 冲激响应不变法 冲激响应不变法的变换示例 设模拟系统函数为： 设采样周期T=1，对应数字系统函数为 冲激响应不变法的特点 从时域角度进行对应，可以保障系统暂态响应时间不变； 可以将s平面左边的极点对应到z平面单位圆内，保障系统的稳定性； 数字滤波器频率响应与模拟滤波器频率响应保持良好线性关系； 冲激响应不变法的局限 由于采样定理限制，模拟系统的频率响应必须具有带限特点，否则会导致频率混迭，因此冲激响应不变法只适用于阻带没有纹波的低通或带通滤波器。 冲激响应不变法的系统转换函数 [bz,az]=impinvar(b,a) （b,a）为模拟滤波器系数 （bz,az）为数字滤波器系数 冲激响应不变法的设计仿真 [z,p,k]=buttap(10);[b,a]=zp2tf(z,p,k); %[z,p,k]=cheb2ap(10,25);[b,a]=zp2tf(z,p,k); t=0:0.1:100;[H,w]=freqs(b,a);h=impulse(b,a,t); subplot(2,2,1),plot(t,h);axis([0,50,-0.3,0.4]); subplot(2,2,2),plot(w,abs(H));axis([0,5,0,1.2]); [bz,az]=impinvar(b,a);N=0:100 hz=dimpulse (bz,az,N);[Hz,wz]=freqz(bz,az); subplot(2,2,3),stem(N,hz);axis([0,50,-0.3,0.4]); subplot(2,2,4),plot(wz,abs(Hz));axis([0,5,0,1.2]); Butterworth滤波器的设计仿真 Chebyshev2型滤波器的设计仿真 双线性变换法 利用非线性函数将数字频率区间对应到模拟频率区间： 双线性变换法 完成模拟系统设计后，再进行反变换，从模拟系统函数得到数字系统函数： 双线性变换法的变换示例 设模拟系统函数为： 数字系统函数为 双线性变换法的特点 没有采样过程，不存在频率混叠问题，适应于各类滤波器的变换； 变换将s平面的虚轴对应到z平面的单位圆，可以保持系统稳定性不变； 变换在高频区域体现出强烈的非线性。 双线性变换法的特点 为了减少非线性关系的影响，实现正确的变换，可以利用参数k对非线性区的范围进行调节；在初步设计时，通常可以先将k值选择在最高模拟截止频率附近，再根据仿真结果进行调整。 双线性变换法的系统转换函数 [bz,az]=bilinear(b,a,k) （b,a）为模拟滤波器系数， （bz,az）为数字滤波器系数， k为非线性调节系数。 双线性变换法的设计仿真 [z,p,k]=cheb2ap(10,25);[b,a]=zp2tf(z,p,k); t=0:0.1:100;[H,w]=freqs(b,a);h=impulse(b,a,t); subplot(2,2,1),plot(t,h);axis([0,50,-0.3,0.3]); subplot(2,2,2),plot(w,abs(H));axis([0,5,0,1.2]); [bz,az]=bilinear(b,a,0.9);N=0:100; hz=dimpulse (bz,az,N);[Hz,wz]=freqz(bz,az); subplot(2,2,3),stem(N,hz);axis([0,50,-0.3,0.3]); subplot(2,2,4),plot(wz,abs(Hz));axis([0,5,0,1.2]); 双线性变换法的设计仿真 从原型滤波器到数字滤波器的变换 模拟变换法的MATLAB设计仿真函数 函数使用说明： 用于频率选择滤波器的设计； 采用双线性变换法进行运算； 数字频率区间需要归一到[0，1]； [b,a]=butter(n,wp,‘ftype) [b,a]=cheby1(n,Rp,wp,‘ftype’) [b,a]=cheby2(n,Rs,ws,‘ftype’) [b,a]=ellip(n,Rp,Rs,wp,‘ftype’) 模拟变换法的MATLAB设计仿真函数 1 选择原型滤波器； 2 将设计