ch4_3脉冲响应不变法讲述.ppt

脉冲响应不变法 脉冲响应不变法 近代数字信号处理(Advanced Digital Signal Processing) 电子信息工程学院 信号与图像处理研究室 IIR数字滤波器设计的基本思想 模拟低通滤波器设计 模拟域频率变换 脉冲响应不变法 双线性变换法 IIR数字滤波器的基本结构 利用MATLAB设计IIR DF 问题的提出 脉冲响应不变法的基本原理 脉冲响应不变法设计DF的步骤 脉冲响应不变法 (Impulse Invariance) 问题的提出 如何将模拟滤波器转变为数字滤波器? 1. 脉冲响应不变法 2. 双线性变换法 Wp,Ws wp,ws H(s) H(z) 频率 变换 设计模拟 滤波器 AF到DF 的转换 脉冲响应不变法的基本原理 对模拟滤波器的单位冲激响应h(t)等间隔抽样来获得数字滤波器的单位脉冲响应h[k] 脉冲响应不变法由H(s)获得H(z)步骤: 1．对H(s)进行Laplace反变换获得h(t)。 2．对h(t)等间隔抽样得到h[k]。 3．计算h[k]的z变换得到H(z)。 脉冲响应不变法的基本原理 脉冲响应不变法由H(s)获得H(z)步骤: H(s) h(t) h[k] H(z) 拉氏反变换 抽样t=kT z变换 设H(s)只有一阶极点，即 对h(t)等间隔抽样得 脉冲响应不变法的基本原理 脉冲响应不变法由H(s)获得H(z) —— 单极点情况 H(s) h(t) h[k] H(z) 拉氏反变换 抽样t=kT Z变换 脉冲响应不变法的基本原理 H(ejW)和H(jw)的关系 无混叠时： 数字滤波器在W点的频率响应和模拟滤波器在w 点的频率响应只差一个常数因子1/T 数字频率W与模拟频率w的关系为 W =wT 脉冲响应不变法的基本原理 脉冲响应不变法的优缺点 缺点：存在频谱混叠，故不能用脉冲响应不变法设 计高通、带阻等滤波器。 优点：数字滤波器和模拟滤波器的频率关系为线性 W =wT 脉冲响应不变法设计DF的步骤 1. 将数字滤波器的频率指标{Wk}转换为 模拟滤波器的频率指标{wk} 2. 由模拟滤波器的指标设计模拟滤波器的H(s)。 3. 利用脉冲响应不变法，将H(s)转换H(z)。 Wp,Ws wp,ws H(s) H(z) w=W/T 设计模拟 滤波器 脉冲响应不变 脉冲响应不变法设计DF的步骤 [numd,dend] = impinvar(num,den,Fs) num,den：AF分子、分母多项式的系数向量 Fs=1/T：抽样频率 numd,dend：DF分子、分母多项式的系数向量 脉冲响应不变法的MATLAB实现 例: 利用BW型模拟低通滤波器和脉冲响应不变法设计满足指标Wp=p/3，Ap=3dB，N=1的数字低通滤波器。 解： 在数字滤波器的设计过程中，参数T可以被抵消。故常取T=1 (1) 将数字低通指标转换成模拟低通指标w=W / T. (2) 设计wp=Wp/T,Ap=3dB的一阶BW型模拟低通滤波器 (3) 将模拟低通滤波器转换成数字低通滤波器 为了消除模、数滤波器频率响应幅度中的1/T，常将TH(s)转化成H(z) 例：利用AF-BW filter及脉冲响应不变法设计一DF，满足 Wp=0.2p, Ws=0.6p, Ap?2dB, As?15dB。 解： (1) 将数字低通指标转换成模拟低通指标，取 T=1. wp=0.2p, ws=0.6p, Ap?2dB, As?15dB (2) 设计模拟低通滤波器 （BW型） =2 =0.801 3 例：利用AF-BW filter及脉冲响应不变法设计一DF，满足 Wp=0.2p, Ws=0.6p, Ap?2dB, As?15dB 。 解： (3) 将模拟低通滤波器转换成数字低通滤波器 极点为s1= -0.567 8 + 0.565 4j, s2= -0.567 8 - 0.565 4j 利用 可得DF的系统函数为 例：利用AF-BW filter及脉冲响应不变法设计一DF，满足 Wp=0.2p, Ws=0.6p, Ap?2dB, As?15dB 。 %Design DF BW low-pass filter using impulse invariance %DF BW LP specfication Wp=0.2*pi; Ws=0.6*pi; Ap=2; As=15; Fs=1;