《06IIR(2)_脉冲响应不变法》-课件.ppt

数字信号处理 习题 5.1 习题 5.3 习题 5.4 回顾 回顾 MATLAB函数 MatLab提供的函数： [N,Wc]=buttord(Wp,Ws,Rp,Rs,’s’) 模拟滤波器 [N,Wc]=buttord(Wp,Ws,Rp,Rs) 数字滤波器 [B,A]=butter(N,Wc,’ftype’,’s’) 模拟滤波器 [B,A]=butter(N,Wc,’ftype’) 数字滤波器 阅读 P161 【例6.2.2】 自编函数 函数 1：由Ωc和N求分子、分母的系数。 function [b,a] = u_buttap(N,Omegac); [z,p,k] = buttap(N); % 归一化巴特沃斯函数【见教材P160】 p = p*Omegac; k = k*Omegac^N; B = real(poly(z)); % 多项式 b0 = k; b = k*B; a = real(poly(p)); 例 1 （续） 用函数计算： 代码如下： N=3; Omegac=2; [b a]=u_buttap(N,Omegac) 结果为： b = 8.0000 a = 1.0000 4.0000 8.0000 8.0000 例 2（续） 由Wp,Ws,Rp,As求分子、分母的系数。 function [b,a] = afd_butt(Wp,Ws,Rp,As); if Wp = 0 error(通带边缘必须大于 0) end if Ws = Wp error(阻带边缘必须大于通带边缘) end if (Rp = 0) | (As 0) error(通带波动或阻带衰减必须大于0) end N = ceil((log10((10^(Rp/10)-1)/(10^(As/10)-1)))/(2*log10(Wp/Ws))); OmegaC = Wp/((10^(Rp/10)-1)^(1/(2*N))); [b,a]=u_buttap(N,OmegaC); 例 2（续） 【例2】的代码 Wp = 0.2*pi; Ws = 0.3*pi; Rp = 7; As = 16; Ripple = 10 ^ (-Rp/20); Attn = 10 ^ (-As/20); [b,a] = afd_butt(Wp,Ws,Rp,As); [C,B,A] = sdir2cas(b,a) 结果为： C = 0.1238 B = 0 0 1 A = 1.0000 0.4985 0.2485 0 1.0000 0.4985 脉冲响应不变法 P177 第6.3节 脉冲响应不变法 脉冲响应不变法 脉冲响应不变法 脉冲响应不变法 脉冲响应不变法 例 4 —— P182 【例6.3.2】 【例 6.3.2】设计一个数字低通滤波器， 指标如下： 代码如下： T=1; %采样周期 Wp=0.2*pi/T; Rp=1; Ws=0.35*pi/T; Rs=10; [N,Wc] = buttord(Wp, Ws, Rp, Rs, s); [B,A]=butter(N, Wc, ‘s’); % 求Ha(s) [Bz,Az] = impinvar(B,A,1/T); % 用脉冲响应不变法求H(z) [c,b,a]=dir2par(Bz,Az) % 由直接型转化为并联型 结果为： c = 0 b = 0.7719 -0.0680 -0.7719 0.2398 a = 1.0000 -0.8774 0.2136 1.0000 -1.0410 0.5276 练习 1 P194 第9题 题意： 拟设计一低通数字滤波器，其技术指标如下： 采样间隔： 要求：（1）采用巴特沃斯模拟低通滤波器模型； （2）采用脉冲响应不变法进行转换。 脉冲响应不变法 练习 2 P193 第1题 小结 可见，N为奇数时，实轴上有极点；N为偶数时，实轴上没有极点。但极点决不会落在虚轴上，这样滤波器才有可能是稳定的。 为形成稳定的滤波器，Ha(s)Ha(-s)的2N个极点中只取S左半平面的N个极点为Ha(s)的极点，而右半平面的N个极点构成Ha(-s)的极点。 Ha(s)的表示式为 回顾 分子系数ΩcN由Ha(s)的低频特性决定，代入Ha(0)=1，可求得； 而sk为 k=1, 2, …, N 回顾 模拟低通滤波器指标： 由参数Ap