数字信号处理实验三 用双线性变换法设计IIR数字滤波器.pptVIP

实验三 用双线性变换法设计 IIR 数字滤波器 测控技术与仪器研究室 指导教师：林晓钢 实验目的 ? 1. 掌握用双线性变换法设计 IIR DF 的原理 及具体设计方法，熟悉用双线性变换法设 计 IIR DF 的计算机编程。 ? 2. 观察用双线性变换法设计的 DF 的频响特 性，了解双线性变换法的特点。 ? 3. 熟悉用双线性变换法设计 BW 和 CB 型 DF 的 全过程。 实验原理及方法 ? 为了克服冲激响应不变法产生的频率混叠现象，这是从 S 平面到 Z 平面的标准变换 z ＝ e sT 的多值对应关系导致的，为了 克服这一缺点，产生了双线性变换法。 0 ? -1 1 ] Im[ z j ? j ] Re[ z 0 s 平面 z 平面 双线性变换法的映射函数： 1 1 1 2 1 2 1 1 2 T s z s z T T z s ? ? ? ? ? ? ? ? ? 分 式 展 开 1 0 2 2 T T sz s z ? ? ? ? 固定其中一个变量，则另一个是线性的。或者说，这种变 换对于 s 和 z 是双线性的。 实验原理及方法 实验原理及方法 ? （ 1 ）把 s = σ + j Ω 带入上式得： 因此： （ 2 ）整个左半平面映射成单 位圆的内部，因此这是一个 稳定的变换。 （ 3 ）虚轴以一一对应的方式 映射成单位圆，因此在频域 中不会出现混叠。 1 1 2 2 2 2 T T T T z j j ? ? ? ? ? ? ? ? ? 1 2 2 0 1 1 2 2 1 2 0 1 1 2 1 2 2 0 1 1 2 2 T T j z T T j T j z T j T T j z T T j ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 0 ? -1 1 ] Im[ z j ? j ] Re[ z 0 s 平面 z 平面 实验原理及方法 1 2 0 1 2 j T j e T j ? ? ? ? ? ? ? ? ? ? z 由于幅度为 1 ，把 ω 看做 Ω 的函数求解，得到： 1 2 2 2 2 tan tan T T ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 这说明 Ω 和 ω 的关系是非线性的，但是没出现混叠。在把 Ω 变换为 ω 时产生了非线性畸变。为了克服它带来的问题，通常要使 Ω 按上式预 修正，以抵消畸变的影响。 （ 4 ? ? ? ? 1 1 2 1 1 a z s T z H z H s ? ? ? ? ? ? IIR 数字滤波器设计流程图 求最小阶数 N ， ω c Buttord,cheblord Cheb2ord,ellipord 合为一步的设计函数 Butter ， cheb1 ， cheb2 ， ellip ， besself 模拟低通滤波器原型设计 buttap,cheb1ap,cheb2ap, besselap,ellipap 函数 频率变换（变为高通、 带通、带阻等） lp2lp,lp2hp,lp2b p,lp2bs 模拟数字变换 bilinear ， impinvar 滤波器 系数 B ， A 滤波器系数 B ， A 典型设计 直接设计 设计指标 IIR 数字滤波器设计流程 ? 设计 IIR 数字滤波器的一般步骤： ? 把给出的数字滤波器的性能指标转换为模拟滤波器的 性能指标 ? 根据转换后的性能指标，通过滤波器阶数选择函数， 来确定滤波器的最小阶数 N 和固有频率 Wn ? 由最小阶数 N 得到低通滤波器原型 ? 由固有频率 Wn 把模拟低通滤波器原型转换为低通、高 通、带通、带阻滤波器 ? 运用脉冲响应不变法或双线性变换法把模拟滤波器转 换成数字滤波器 实验程序流程（ Matlab ） ? 低通滤波器，采样频率为 1000Hz ，通带临界频率 f l =200Hz ，通带内衰减小于 1dB （ α p =1 ）；阻带临界频率 f h =300Hz ，阻带内衰减大于 25dB （ α s =25 ）。设计一个数字滤 波器满足以上参数。 ? FS=1000; ? Fl=200;Fh=300; % 通带、阻带截止频率 ? Rp=1;Rs=25; ? wp=Fl*2*pi; % 临界频率采用角频率表示 ? ws=Fh*2*pi; % 临界频率采用角频率表示 ?