门爱东老师DSP讲义第4章2.ppt

设θ和ω分别是 z 平面和 Z 平面的数字频率变量，即各自单位园，令 4.6.1 IIR DF 设计：数字频率变换 任何全通函数一般可表示为 4.6.1 IIR DF 设计：低通-低通数字频率变换 低通 HLP(z) — 低通 Hd(Z) 变换 对于这种情况，HLP(z) 和 Hd(Z) 都是低通数字滤波器，只不过截止频率不同或带宽不同。 当原型 HLP(z) 的数字频率θ从 0 变化到π时，所要求的 Hd(Z) 的数字频率ω也是从 0 变化到 π，而此时全通函数 G(Z-1) 的相角φ(ω) 的变化量为 Nπ。因此，全通函数 G(Z-1) 的阶数 N＝1。 4.6.1 IIR DF 设计：低通-低通数字频率变换 设原型和变换后的低通滤波器的截止频率分别为θc 和ωc，代入上式： 4.6.1 IIR DF 设计：低通-高通数字频率变换 4.6.1 IIR DF 设计：低通-带通数字频率变换 N=2，则全通函数为 4.6.1 IIR DF 设计：低通-带通数字频率变换 4.6.1 IIR DF 设计：低通-带通数字频率变换 方法一：数字频率变换 DF Hd(Z) 频率指标变换为数字低通 HLP(z) 频率指标。为数字LPF 任意选定一个通带截止频率θc，利用表4.8 中变换公式得到相应的数字 LPF 的阻带起始频率θS； 低通 DF 指标变换为低通 AF 指标。利用双线性变换的频率预畸变公式 Ω= (2/T) tan(θ/2) ，把数字LPF HLP(z) 的数字频率指标转换为相应的模拟 LPF HLP(s) 的模拟频率指标； 设计模拟低通滤波器原型 HLP(p) ，得到模拟低通滤波器的传输函数 HLP(s)|p=s/Ωc； 利用双线性变换得到相应的数字低通滤波器传输函数 利用数字频率变换将数字低通滤波器传输函数 HLP(z) 转换成所要求的数字滤波器传输函数 Hd(z)。 4.6.2 IIR DF 设计： 利用数字频率变换 设计模拟原型 Chebyshev I 低通滤波器 数字滤波器设计步骤： 把 DF 技术指标转换为相应的 AF 技术指标； 设计原型低通模拟滤波器 HLP(s)，例如 Butterworth、Chebyshev I、Chebyshev II 或者 Ellip ； 利用脉冲响应不变法或双线性变换法，由模拟滤波器原型 HLP(s) 得到所要求的数字滤波器的传输函数 HLP(z)； 对于其它类型滤波器，采用模拟域 S 或数字域 Z 中完成频率转换。 Matlab 实现 上面的步骤，每一步都有相应的 Matlab 函数或语句； Matlab 还提供了综合性更高的一些函数实现上面的这些步骤； 这些函数可用于设计低通、高通、带通和带阻滤波器。 这些函数都采用了双线性变换，并在 s 域中完成频率转换。 求滤波器阶数 N 和截止频率 wn Matlab 函数 [n,Wn] = buttord(Wp,Ws,Rp,Rs) [n,Wn] = cheb1ord(Wp,Ws,Rp,Rs) [n,Wn] = cheb2ord(Wp,Ws,Rp,Rs) [n,Wn] = ellipord(Wp,Ws,Rp,Rs) 参数 Wp 和 Ws 限制： 对于低通滤波器： Wp Ws； 对于高通滤波器： Wp Ws； 对于带通滤波器：Wp 和 Ws 是具有两个元素的矢量，Wp = [Wp1,Wp2] 和 Ws = [Ws1,Ws2]，并且 Ws1 Wp1 Wp2 Ws2； 对于带阻滤波器：Wp 和 Ws 是具有两个元素的矢量，Wp = [Wp1,Wp2] 和 Ws = [Ws1,Ws2]，并且 Wp1 Ws1 Ws2 Wp2。 数字滤波器设计函数 buttrworth 数字滤波器 [b,a] = butter(n,Wn) 低通、带通 [b,a] = butter(n,Wn,ftype) 高通、带阻 其中的 Wn 为截止频率，在butterworth 滤波器中为 3dB 处的频率，等于函数 [n,Wn]=buttord(Wp,Ws,Rp,Rs) 返回的 Wn。 chebyshev I 数字滤波器 [b,a] = cheby1(n,Rp,Wn) 低通、带通 [b,a] = cheby1(n,Rp,Wn,ftype) 高通、带阻 其中 Wn 为截止频率，等于滤波器的幅度响应为 Rp(dB) 处的频率值，即 Wn = Wp，或者是函数[n,Wn] = cheb1ord(Wp,W