数字信号处理CH7.pptVIP

7.0 设计过程 由应用定指标 由特性定类型：IIR简单，FIR有线性相位 求满足指标的H(z) 设计实现所用结构、字长、算法 仿真确认 本章重点讨论求H(z)，有现成的软件(如MATLAB)辅助求解滤波器系数，关键了解各种算法特性 理想选频滤波器 对不同频率分量幅度加以改变都统称滤波器 例：二维图像的高低通滤波效果 例: 二维图像的带通带阻滤波效果 原图 频谱 带阻 输出 带通输出 常用滤波器指标 通带增益 通带纹波δ1 通带截止频率ωp 阻带衰减 阻带纹波δ2 阻带截止频率ωs 过渡带宽 Δω=ωs - ωp 7.1 IIR滤波器设计：连续时间滤波器转换法 典型连续时间模拟滤波器 Butterworth：通、阻带内单调 Chebyshev Type I：通带等纹波、阻带单调 Type II：通带单调、阻带等纹波 椭圆函数：通、阻带内皆等纹波 相同指标时，允许纹波将减少阶数，降低系统运算复杂度 算法见参考书附录，可用MATLAB实现 典型模拟滤波器幅频特性 例：复习巴特沃斯低通设计 巴特沃斯平方函数：|Hc(jΩ)|2=1/[1+(Ω/Ωc)2N] 特性：单调，Ωc为半功率点，N加大过渡带减小 算法：利用单调性，用通带、阻带截止频率处条件可解出二待定因子Ωc、阶数N N必须为整数 IIR转换法I：脉冲响应不变法 由应用确定H(ejω)指标，一般为对|H(ejω)|的要求 由Ω=ω/Td得Hc(jΩ)指标 用模拟滤波器设计法求出hc(t) 用Td采样最终得：h[n]=Tdhc(nTd) 尽管Td并非真实采样周期T，在来回转换过程中作用抵消，可取任意值或1，但hc(t)必须带限：Hc(jΩ)=0，|Ω|≥π/Td。故该方法不适宜作高通！ 对实际非理想带限，提高带外衰减指标后忽略 对应系统函数的转换公式 拉氏变换系统函数：Hc(s)= Ak/(s-sk) 脉冲响应：hc(t) = Akesktu(t) 则：h[n]= TdAkesknTd u[n] = TdAk(eskTd)n u[n] H(z)= TdAk/(1 - eskTd z-1) 对比Hc(s)与H(z)，部分分式项系数乘Td S平面极点sk对应Z平面极点zk= eskTd 例：数字巴特沃斯LPF 给定指标：0.89≤|H(ejω)|≤1， 0≤|ω|≤0.2π |H(ejω)|≤0.18，0.3≤|ω|≤π 令Td=1，由Ω=ω/Td=ω，得连续时间系统指标： 0.89≤|Hc(jΩ)|≤1， 0≤|Ω|≤0.2π |Hc(jΩ)|≤0.18，0.3π≤|Ω|≤π 通带条件：|Hc(j0.2π)|=0.89 阻带条件：|Hc(j0.3π)|=0.18 代入巴特沃斯平方函数，得：N=5.88，Ωc=0.705 例：数字巴特沃斯LPF(续) 放大N至整数6，此时通、阻带指标都略有提高 注意该LPF非带限，为减小混叠，提高放给阻带, 用N=6及通带条件重算，最终得 Ωc=0.703 为因果稳定，S平面|Ω|≤π/T 左半对应Z平面单位圆内|ω|≤π 用S平面左半三对极点，可得 Hc(s)，进而用前述转换关系求出 H(z) 例：数字LPF(续) 6阶数字巴特沃斯LPF 对数幅频，幅频与群延迟 由于混叠，在|ω|=π处幅频增加一倍，对数幅频增加6dB 对|ω|=0.2π、0.3π截止频率处混叠影响很小，系统满足给定指标 IIR转换法II：双线性变换法 S平面至Z平面非线性变换： jΩ轴，-∞≤Ω≤∞ ? Z平面单位圆，-π≤ω≤π 频率变换关系：设z=ejω在单位圆上求值，代入得：s=σ+jΩ=(2/Td)(1-ejω)/(1+ejω)=j(2/Td)tg(ω/2) 虚实部对应，得σ=0，证明jΩ轴与单位圆对应，且：Ω = (2/Td)tg(ω/2) 双线性变换中频率非线性变换 频率轴非线性压缩，无混叠问题，适应分段恒幅，包括非带限的高通 线性相位、幅度(如微分、残余边带)经双线性变换将失去线性，可用脉冲响应不变法 S平面与Z平面的转换关系 稳定因果，S平面左半空间对应Z平面单位圆内 S平面jΩ轴对应Z平面单位圆上 截止频率的预畸变 Ω = (2/Td)tg(ω/2) Ωp = (2/Td)tg(ωp/2) Ωs = (2/Td)tg(ωs/2) Td同样任意可取1 例：用双线性变换法设计数字巴特沃斯LPF 给定指标：0.89≤|H(ejω