数字信号处理第4章IIR数字滤波器的3.pptVIP

IIR数字滤波器设计的基本思想 模拟低通滤波器设计 模拟域频率变换 脉冲响应不变法 双线性变换法 利用MATLAB设计IIR DF 问题的提出 脉冲响应不变法的基本原理 脉冲响应不变法的基本原理 脉冲响应不变法的基本原理 脉冲响应不变法的基本原理 脉冲响应不变法的基本原理 例：利用AF-BW filter及脉冲响应不变法设计一DF，满足 Wp=0.2p, Ws=0.6p, Ap?2dB, As?15dB 例：利用AF-BW filter及脉冲响应不变法设计一DF，满足 Wp=0.2p, Ws=0.6p, Ap?2dB, As?15dB 例：利用AF-BW filter及脉冲响应不变法设计一DF，满足 Wp=0.2p, Ws=0.6p, Ap?2dB, As?15dB 例：利用AF-BW filter及脉冲响应不变法设计一DF，满足 Wp=0.2p, Ws=0.6p, Ap?2dB, As?15dB 例：利用AF-BW filter及脉冲响应不变法设计一DF，满足 Wp=0.2p, Ws=0.6p, Ap?2dB, As?15dB * * 问题的提出 脉冲响应不变法的基本原理 脉冲响应不变法设计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 例: 设一阶模拟低通滤波器的系统函数为 利用脉冲响应不变法求H(z)，并分别画出AF与DF的幅度响应。 解： 利用单极点H(s) 与H(z)的映射关系，可得 AF与DF的频率响应分别为 其中W =wT，抽样频率为50，200Hz的幅度响应如下 例: 设一阶模拟低通滤波器的系统函数为 利用脉冲响应不变法求H(z)，并分别画出AF与DF的幅度响应。 解： 0 5 10 15 20 25 0 0.2 0.4 0.6 0.8 1 Hz DF AF fsam=50 Hz 例: 设一阶模拟低通滤波器的系统函数为 利用脉冲响应不变法求H(z)，并分别画出AF与DF的幅度响应。 解： 0 20 40 60 80 100 0 0.2 0.4 0.6 0.8 1 Hz DF AF fsam=200 Hz 脉冲响应不变法的优缺点 缺点：存在频谱混叠，故不能用脉冲响应不变法设 计高通、带阻等滤波器。 优点：数字滤波器和模拟滤波器的频率关系为线性 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实现 解： (1) 将数字低通指标转换成模拟低通指标w=W / T. wp=0.2p/T, ws=0.6p/T, Ap?2dB, As?15dB (2) 设计模拟低通滤波器 （BW型） =2 =0.801