数字信号处理第6章无限脉冲响应数字滤波器的设计-(精品课件).pptVIP

解 考虑到通带和过渡带的重要，在0~0.2π区间，每隔0.01π取一点ωi值，在0.2π~π区间每隔0.1π取一点ωi值，并增加一点过渡带，在ω=0.1π处 |Hd(e jω)|=0.5。 1.0,ω=0,0.01π,0.02π,:,0.09π 0.5,ω=0.1π 0.0,ω=0.11π,0.12π,:,0.19π 0.0,ω=0.2π,0.3π,:,π N=29，取k=1，系统函数为 待求的参数是A，a1,b1,c1,d1。设初始值θ=(0000.25)T经过90 次迭代，求得E=1.2611，系统函数零、极点位置为 零点0j0 极点0j1.3213916 为使滤波器因果稳定，将极点按其倒数搬入单位圆内，再进行62次优化迭代，求得结果为 ?零点0j0 极点0j0 Ag=0E=0.56731 误差函数用下式表示： (6.6.12) 3. 在时域直接设计IIR数字滤波器 设我们希望设计的IIR数字滤波器的单位脉冲响应为 hd(n)，要求设计一个单位脉冲响应h(n)充分逼近hd(n)。下 面我们介绍这种设计方法。 设滤波器是因果性的，系统函数为 (6.6.13) 谢谢！ 图6.4.7 例6.4.2图——用脉冲响应不变法设计的数字低通滤波器的幅度特性 (2) 用双线性变换法设计数字低通滤波器。 ① 数字低通技术指标仍为 ωp=0.2πrad,αp=1dB; ωs=0.3πrad,αs=15dB ② 模拟低通的技术指标为 ③ 设计巴特沃斯低通滤波器。阶数N计算如下： 取N=6。为求Ωc,将Ωs和αs代入(6.2.18)式中，得 到Ωc=0.7662rad/s。这样阻带技术指标满足要求，通 带指标已经超过。 根据N=6，查表6.2.1得到的归一化传输函数Ha(p)与脉冲响应不变法得到的相同。为去归一化，将p=s/Ωc代入Ha(p)，得实际的Ha(s)， ④ 用双线性变换法将Ha(s)转换成数字滤波器H(z)： 图6.4.8 例6.4.2图——用双线性变换法设计的数字低通滤波器的幅度特性 6.5 数字高通、带通和带阻滤波器的设计 例如高通数字滤波器等。具体设计步骤如下： (1) 确定所需类型数字滤波器的技术指标。 (2) 将所需类型数字滤波器的技术指标转换成所需类型模拟滤波器的技术指标，转换公式为 (3)将所需类型模拟滤波器技术指标转换成模拟低通滤波器技术指标(具体转换公式参考本章6.2节)。 (4)设计模拟低通滤波器。 (5)将模拟低通通过频率变换，转换成所需类型的模拟滤波器。 (6)采用双线性变换法，将所需类型的模拟滤波器转换成所需类型的数字滤波器。 例6.5.1 设计一个数字高通滤波器，要求通带截止频率ωp=0.8πrad，通带衰减不大于3dB，阻带截止频率ωs=0.44πrad,阻带衰减不小于15dB。希望采用巴特沃斯型滤波器。 解 (1)数字高通的技术指标为 ωp=0.8πrad,αp=3dB; ωs=0.44πrad,αs=15dB (2) 模拟高通的技术指标计算如下： 令T=1，则有 (3)模拟低通滤波器的技术指标计算如下： 将Ωp和Ωs对3dB截止频率Ωc归