7数字滤波器的设计教程详解.pptVIP

220 1.直接设计； 累试； 极点，峰值；零点，谷值 时域、频域（FIR滤波器的设计）； 最优化； 2. 间接设计； IIR滤波器的设计 246 7.3 模拟滤波器数字化 模拟滤波器数字化 1. 冲击响应不变法（时域） 2. 双线性变换（频域） 注意：性能指标的转换方法必须与数字化方法相对应 （一）冲击响应不变法（时域） 从时域确定S平面和Z平面之间的映射关系 数字性能指标 ? 模拟滤波器的性能指标； （? = ?T） 0 ha(t) h(n) 取样 取样 x(t) ha(t) y(t) x(t) y(t) x(nT) y(nT) x(n)=x(nT) h(n)=Th(nT) y(n)=y(nT) 0 T T T 0 t t t 模拟滤波器数字化（1) 在任一个离散时间点 nT，两个系统是否具有相同的输出？ 模拟系统? 数字系统? 时域瞬态特征相近 ? 频率特征相近？ 247 模拟滤波器数字化（2) 频域分析 xa(t) 取样 xa(nT) 而 ha(t) Ha(Ω) 取样 Tha(nT) 248 模拟滤波器数字化（3) 步骤： 1）对已知的Ha(s) 进行拉氏反变换，求得ha(t)； 2）对ha(t) 进行取样，得ha(nT)； 3）令h(n)=T ha(nT)，以求得h(n)； 4）对h(n) 进行Z 变换，得H（Z）。 即： 模拟滤波器数字化（4) 单阶极点： 模拟滤波器的系统函数可表示为： 而且，一般MN，因此上式可以分解为部分分式形式。 对Ha(s) 两边进行拉氏反变换，可得： 模拟滤波器数字化（5) 单阶极点： 对ha(t) 以周期T进行取样，有： 按冲激响应不变准则，可得： 对上式两边进行Z变换，便得数字滤波器的系统函数： 模拟滤波器数字化（6) 单阶极点： S平面与Z平面之间的映射关系 极点间的映射关系： 由此得S平面与Z平面间的映射关系： 令 则 所以： 模拟滤波器数字化（7) S平面与Z平面之间的映射关系： 1）当σ=0 时，r=1，即S平面虚轴映射为 Z平面的单位圆 ; 2) 当σ0时，r1，即S平面右半平面映射到Z平面的单位圆之外； 3）当σ0时，r1，即S平面左半平面映射到Z平面的单位圆之内； S平面与Z平面之间的映射关系： S平面上 映射为 Z 平面上 -π≤ω≤π 在S平面上，任何一个宽度为（2π/T）的水平带，映射为一Z平面； → H(ejω)为周期函数； 模拟滤波器数字化（8) 模拟滤波器数字化（9) 小结： 1）根据 h(n) =T ha(nT) 从时域完成数字化设计，DF和AF之间具有相近似时域瞬态特征； 2）数字频率 ω 和模拟滤波器频率Ω之间的关系为 ω=ΩT，即两者间呈线性关系； 5）当Ha(Ω)不严格限带，或ha(t) 变化不太平稳（在阻带内存在振荡，例如Chebyshev II 滤波器），而设计性能要求又高时，则不宜采用此法。另外，此方法不能直接设计高通、带阻滤波器。 3）延拓相加与混叠； 4)用该方法设计出的数字滤波器在 附近的频率响应会不同程度的偏离模拟滤波器在 附近的频率响应。因此要求模拟滤波器是一个带限滤波器，而实际的频响不可能真正带限。只要在折叠频率 以上频响衰减加快，用该方法设计出的DF可以较好地满足要求。 模拟滤波器数字化（10) (二）、双线性变换法 原理: 冲激响应不变法存在混迭问题; 双线性变换主要目的: 从根本上解决混叠问题. 具体方法: S平面 Z平面 S1 平面 π/T -π/T Im(z) Re(z) jΩ1 jΩ 0 σ σ1 1) 从S平面到S1平面的映射; 2) 从S1平面到Z平面的映射; AF的Ha(s) 表示为： 令 M = N，则： 时域分析: 再令j = N-i，则： 因此，AF的基本单元是积分器 s-1，而利用某种数字网络代替此单元，就可把AF转变为相应的DF. 模拟滤波器数字化（11) 249 模拟滤波器数字化（12) 设： 则： 设输入信号为Xa(t) 则输出为： 设有0t1t2，则： 时域分析: ∵ ∴ 250 模拟滤波器数字化（13) 数字化, 即令y(n)=ya(nT)，x(n)=xa(nT)，则得差分方程: Z变换,得： 因此, 数字积分器的系统函数： 时域分析: 即： 则： 251 模拟滤波器数字化（14) 设 s=σ+jΩ， 则: S平面同Z平面之间的映射关系: S平面 Z平面 σ0，右半平面 r1，单位圆外 σ=0，虚轴 r=1，单位圆 σ0，左半平面 r1，单位圆内 令σ=0，得: 当 Ω 由 ω 由 -π 到π 变化 252 模拟滤波器数字化（15) 利用双线性变换，把 转换成