第七章无限脉冲响应数字滤波器的.ppt

数字滤波器设计方法 １） IIR滤波器设计方法 借助模拟滤波器设计方法： 计算机辅助设计：直接在频域或时域设计 ２）FIR滤波器设计方法 窗函数法 频率采样法 等波纹逼近法——需计算机辅助设计 ３）线性相位滤波器设计方法 FIR滤波器：常用（相位特性严格线性，这是AF无法达到的） IIR滤波器：必须使用全通网络对其非线性相位特性进行相位校正 其设计步骤是：   所以，已知H（S）通过变换 可以得到DF的系统函数H（Z）。? ???映射关系： 当σ0时，r1 说明S平面的左半平面映射到Z平面的单位圆周内当σ=0时，r=1 说明S平面的jΩ轴映射到Z平面的单位圆周当σ0时，r1 说明S平面的右半平面映射到Z平面的单位圆周外。 ???所以多个S值映射为单个Z值。S平面的jΩ轴映射为Z平面的单位圆周，所以ΩT每增加或减少2?，即Ω增加或减少2?/T,对应于单位圆周上逆时针旋转一周。S平面上每一条宽为2?/T的横带部分都将重叠地映射为整个Z平面，每一条横带的左、右半平面分别映射为单位园内、外。如图7.2.1所示：? ????????? 以上分析说明：映射的是H（Z）与H（S）周期延拓的关系，而不是H（Z）与H（S）本身的关系。即： ? ?????  图7.2.2 脉冲响应不变法的频率混叠现象可看出：如果Ha(jΩ) 中有些分量的频率大于Ωs/2(=π/T),则DF的H(ejw)中必然产生混叠现象,设计出的DF的响应就会产生失真，所以要求：AF的频响Ha(jΩ)的最高频率Ωh≤π/T。即：? 总结：用脉冲响应不变法设计IIR DF的步骤为：1、通过w=ΩT将给定的DF的技术要求转化为AF的 技术要求2 、设计此AF，得到Ha(S) 3 、对H(s)进行Laplace反变换求得ha(t) 4 、令5、求H（Z）。? 实际上问题可简化，一般Ha(s)的极点si是一个复数，且以共轭成对的形式出现，将一对复数共轭极点放在一起，形成一个二阶基本节。如果模拟滤波器的二阶基本节的形式为:  ? 如果Ha(s)的极点si是一个实数，则AF是一阶基本节的形式: 则DF为：?? 例7.2.1 已知模拟滤波器的传输函数Ha(s)为?用脉冲响应不变法将Ha(s)转换成数字滤波器的系统函数H(z)。解：?? 以T=1s代入，用H1(z)表示，T=0.1s时用H2(z)表示，则 分别画出它们的幅频特性如图7.2.3所示。其中图7.2.3(a)表示模拟滤波器的幅频特性。????? 图7.2.3 例7.2.1的幅度特性  由  可求得 ??? ??∴ 当σ0时，|Z|1 说明S左半平面映射到Z平面的单位圆 周内当σ=0时,|Z|=1 说明S平面的jΩ轴映射到Z平面的单位圆周当σ0时,|Z|1 说明S右半平面映射到Z平面的单位圆周外。? 令s=jΩ, z=e jω，有： ? 即在给定DF的ωs,ωp后，在设计AF滤波器时，并不直接按这个数据设计，而是先通过进行预变形，根据Ωp,Ωs进行AF的设计，求出H（S），再用 代入，即可求出DF的H（Z）。?  (3)按照模拟滤波器的技术指标设计模拟滤波器Ha(s)。  (4)将模拟滤波器Ha(s)，从s平面转换到z平面，得到数字滤波器系统函数H(z)。 ? 为简化设计，将模拟滤波器各系数和经双线性变换法得到的数字滤波器的各系数之间关系，列成表格供设计时使用。? ? 例7．3．1：设有一数字处理系统，它的抽样频率fsa为2000Hz，希望在此系统中设计一个一阶低通数字滤波器，使其通带中允许的最大衰减为3dB，通带上限频率fp为400 Hz ? ???数字低通系统函数： ?? 解: (1) 数字高通的技术指标为 ωp=0.8πrad, αp=3dB; ωs=0.44πrad, αs=15dB (3)模拟原型低通滤波器的技术指标计算如下： 查表6.2.1，得到归一化模拟低通传输函数G(p)为 ?(6)用双线性变换法将模拟高通H (s)转换成数字高通H(z)： 解 (1)数字带通滤波器技术指标为 通带上截止频率