数字信号处理(六)IIR数字滤波器的设计
文章目录
- 数字滤波器
- 数字滤波器技术指标
- 数字低通滤波器的幅频响应曲线
- IIR滤波器设计方法
- IIR滤波器的函数模型设计法(间接法)
- 模拟低通滤波器的技术指标
- 模拟滤波器原型介绍
- 1、巴特沃斯模拟低通滤波器
- 介绍
- 巴特沃斯滤波器的设计步骤
- 2、切比雪夫低通滤波器
- 介绍
- 切比雪夫滤波器的设计步骤
- 3、椭圆模拟低通滤波器
- 介绍
- 对比
- 模拟到数字滤波器的转换
- 脉冲响应不变法
- 双线性变换法
- 数字滤波器的频带变换
数字滤波器
什么是数字滤波器?
指输入输出均为数字信号,通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。是离散系统的通用名称。
数字滤波器的分类:
经典滤波器,特点是其输入信号中有用的频率成分和期望滤除的频率成分各占不同的频带,通过一个合适的选频滤波器滤除某个频带或频率成分的干扰,得到纯净信号,达到滤波的目的
现代滤波器,用于信号与干扰的频带互相重叠的情况,现代滤波器是根据随机信号的一些统计特性,在某种最佳准则下,最大限度的抑制干扰,恢复信号,从而达到滤波的目的,如维纳滤波器
、卡尔曼滤波器
等
经典数字滤波器从滤波特性上可分为低通
、高通
、带通
、带阻
和全通
数字滤波器从单位脉冲响应长度分类,可分为无限长单位脉冲响应(IIR)滤波器
和有限长单位脉冲响应(FIR)滤波器
离散线性时不变系统的差分方程为:y(n)=∑i=0Mbix(n−i)−∑k=1Naky(n−k)y(n)=\sum_{i=0}^Mb_ix(n-i)-\sum_{k=1}^Na_ky(n-k)y(n)=∑i=0Mbix(n−i)−∑k=1Naky(n−k),那么其系统函数可表示为:H(z)=∑i=0Mbiz−i1+∑k=1Nakz−kH(z)=\frac{\sum_{i=0}^M b_iz^{-i}}{1+\sum_{k=1}^Na_kz^{-k}}H(z)=1+∑k=1Nakz−k∑i=0Mbiz−i
- 当N≥1,aka_kak中至少有一个非零系数时,该系统中存在反馈回路,其所对应的滤波器称作无限长脉冲响应滤波器(infinite impulse response filter),简称IIR滤波器。N是IIR滤波器的阶数,表示系统中反馈环的个数。一般假定IIR数字滤波器满足M≤N,这时将系统称为N阶的IIR数字滤波器。
- 当aka_kak均为零系数时,其所对应的滤波器称作有限长脉冲响应滤波器(finite impulse response filter),简称FIR滤波器
数字滤波器技术指标
常用的数字滤波器一般属于选频滤波器,假设数字滤波器的频率响应函数H(ejw)H(e^{jw})H(ejw)用下式表示:
H(ejw)=∣H(ejw)∣ejθ(w)H(e^{jw})=|H(e^{jw})|e^{j\theta(w)} H(ejw)=∣H(ejw)∣ejθ(w)
- ∣H(ejw)∣|H(e^{jw})|∣H(ejw)∣称为幅频特性函数,表示信号通过该滤波器后各频率成分振幅衰减情况
- θ(ω)\theta(\omega)θ(ω)为相频特性函数,反映各频率成分通过滤波器后在时间上的延时情况
- 一般选频滤波器的技术指标由幅频特性给出
实际滤波器的通带和阻带中都允许有一定的误差容限,即通带不是完全平的,阻带不是绝对衰减到零,通带和阻带之间还有一定的过渡带
数字低通滤波器的幅频响应曲线
通带内和阻带内允许的衰减一般用dB数表示,归一化后通带衰减αp\alpha_pαp和阻带衰减αs\alpha_sαs分别定义为:
αp=−20lg∣Ha(jωp)∣dB\alpha_p=-20lg|H_a(j\omega_p)|dB αp=−20lg∣Ha(jωp)∣dB
αs=−20lg∣Ha(jωs)∣dB\alpha_s=-20lg|H_a(j\omega_s)|dB αs=−20lg∣Ha(jωs)∣dB
当幅度下降到2/2\sqrt{2}/22/2时,ω=ωc\omega=\omega_cω=ωc,此时αp=3dB\alpha_p=3dBαp=3dB,称为3dB通带截止频率
ωp,ωs,ωc\omega_p,\omega_s,\omega_cωp,ωs,ωc统称为边界频率,它们是数字滤波器的重要参数
IIR滤波器设计方法
直接法: 直接在频域或者时域中设计数字滤波器
间接法: 借助于模拟滤波器的设计原型进行设计的,又称滤波器的函数模型设计法,较常用。
IIR滤波器的函数模型设计法(间接法)
在IIR滤波器设计过程中,通常利用模拟滤波器来设计数字滤波器
借助模拟滤波器的设计方法的基本步骤:
将数字滤波器的技术指标转换成模拟滤波器的技术指标
按转换后技术指标、设计模拟低通滤波器的Ha(s)H_a(s)Ha(s)
然后将Ha(s)H_a(s)Ha(s)按某种方法转换成数字滤波器的系统函数H(z)H(z)H(z)
如果不是低通,则必须先将其转换成低通模拟滤波器的技术指标
设计过程中用到的两种变换:
①频带变换: 将低通转换为高通、带通、带阻滤波器
②变换域变换: 将模拟滤波器转换成数字滤波器,从s->z或者Ha(s)−>H(z)H_a(s)->H(z)Ha(s)−>H(z)的变换
IIR滤波器的间接设计法是以模拟滤波器为基础的,有若干典型的模拟滤波器原型可供选择:
- 巴特沃思滤波器:具有
单调下降
的幅频特性 - 切比雪夫滤波器:幅频特性在通带或者阻带有等波纹特性,可以提高选择性
- 椭圆滤波器:选择性相对前两种是最好的,单通带和阻带内均呈现等波纹幅频特性,相位特性的非线性也稍严重
- 贝塞尔滤波器
模拟低通滤波器的技术指标
模拟低通滤波器的技术指标有αp,Ωp,αs,Ωs\alpha_p,\Omega_p,\alpha_s,\Omega_sαp,Ωp,αs,Ωs,其中Ωp和Ωs\Omega_p和\Omega_sΩp和Ωs分别称为通带截止频率和阻带截止频率,αp\alpha_pαp是通带内的最大允许衰减,αs\alpha_sαs是阻带内最小允许衰减。归一化的αp,αs\alpha_p,\alpha_sαp,αs分别表示如下:
αp=−10lg∣Ha(jΩp)∣2dB\alpha_p=-10lg|H_a(j\Omega_p)|^2dB αp=−10lg∣Ha(jΩp)∣2dB
αs=−10lg∣Ha(jΩs)∣2dB\alpha_s=-10lg|H_a(j\Omega_s)|^2dB αs=−10lg∣Ha(jΩs)∣2dB
Ha(jΩc)=2/2H_a(j\Omega_c)=\sqrt{2}/2Ha(jΩc)=2/2时,−10lg∣Ha(jΩc)∣2=3dB-10lg|H_a(j\Omega_c)|^2=3dB−10lg∣Ha(jΩc)∣2=3dB,故称Ωc\Omega_cΩc为模拟低通滤波器的3dB截止频率
模拟低通滤波器的上述技术指标给定后,需要设计滤波器的系统函数Ha(s)H_a(s)Ha(s),具体步骤如下:
如果我们能由αp,Ωp,αs,Ωs\alpha_p,\Omega_p,\alpha_s,\Omega_sαp,Ωp,αs,Ωs,求出∣Ha(jΩ)∣2|H_a(j\Omega)|^2∣Ha(jΩ)∣2,根据∣Ha(jΩ)∣2=Ha(jΩ)Ha∗(jΩ)=Ha(s)Ha(−s)∣s=jΩ|H_a(j\Omega)|^2=H_a(j\Omega)H_a^*(j\Omega)=H_a(s)H_a(-s)|_{s=j\Omega}∣Ha(jΩ)∣2=Ha(jΩ)Ha∗(jΩ)=Ha(s)Ha(−s)∣s=jΩ,Ha(s)Ha(−s)H_a(s)H_a(-s)Ha(s)Ha(−s)的极点是对称的,Ha(s)H_a(s)Ha(s)必须是因果稳定的,因此极点必须落在s平面的左半平面,相应的Ha(−s)H_a(-s)Ha(−s)的极点必须落在s平面的右半平面。由此,可求出Ha(s)H_a(s)Ha(s)
模拟滤波器原型介绍
1、巴特沃斯模拟低通滤波器
介绍
特点:通带和阻带的幅度响应都是平的
N阶巴特沃斯(Butterworth)模拟低通滤波器的幅频响应的平方函数为:
∣Ha(jΩ)∣2=11+(Ω/Ωc)2N|H_a(j\Omega)|^2=\frac 1{1+(\Omega/\Omega_c)^{2N}} ∣Ha(jΩ)∣2=1+(Ω/Ωc)2N1
其中ΩC\Omega_CΩC是巴特沃斯低通滤波器的3dB截止频率
∣Ha(j0)∣2=1,∣Ha(j∞)∣2=0|H_a(j0)|^2=1,|H_a(j\infty)|^2=0 ∣Ha(j0)∣2=1,∣Ha(j∞)∣2=0
巴特沃斯模拟低通滤波器的性质:
①幅频响应平方函数∣Ha(jΩ)∣2|H_a(j\Omega)|^2∣Ha(jΩ)∣2是 Ω\OmegaΩ的单调减函数 :
d∣Ha(jΩ)∣2dΩ=−2N(Ω/Ωc)2N−1[1+(Ω/Ωc)2N]2<0\frac{d|H_a(j\Omega)|^2}{d\Omega}=-\frac{2N(\Omega/\Omega_c)^{2N-1}}{[1+(\Omega/\Omega_c)^{2N}]^2}<0 dΩd∣Ha(jΩ)∣2=−[1+(Ω/Ωc)2N]22N(Ω/Ωc)2N−1<0
②对于所有的N,当Ω=0\Omega=0Ω=0时,∣Ha(jΩ)∣2=1|H_a(j\Omega)|^2=1∣Ha(jΩ)∣2=1
③Ω=Ωc\Omega=\Omega_cΩ=Ωc时,∣Ha(jΩ)∣=2/2|H_a(j\Omega)|=\sqrt{2}/2∣Ha(jΩ)∣=2/2,即在∣Ha(jΩ)∣=1|H_a(j\Omega)|=1∣Ha(jΩ)∣=1处有3dB的衰减
④Ω>Ωc\Omega>\Omega_cΩ>Ωc时,随Ω\OmegaΩ增大,∣Ha(jΩ)∣|H_a(j\Omega)|∣Ha(jΩ)∣迅速下降,Ωc\Omega_cΩc衰减的速度与阶数N有关,N愈大,衰减速度愉快,过渡带愈窄,N->∞∞∞时,趋向于理想低通滤波器
⑤在Ω=0\Omega=0Ω=0处各阶导数存在且等于0,因此在该点得到最大值
巴特沃斯滤波器的系统函数:
∵\because∵
Ha(s)Ha(−s)=∣Ha(jΩ)∣2=11+(Ω/Ωc)2N∣Ω=s/j=11+(sjΩc)2N=(jΩc)2Ns2N+(jΩc)2NH_a(s)H_a(-s)=|H_a(j\Omega)|^2=\frac 1{1+(\Omega/\Omega_c)^{2N}}|_{\Omega=s/j}=\frac 1{1+(\frac{s}{j\Omega_c})^{2N}}=\frac{(j\Omega_c)^{2N}}{s^{2N}+(j\Omega_c)^{2N}} Ha(s)Ha(−s)=∣Ha(jΩ)∣2=1+(Ω/Ωc)2N1∣Ω=s/j=1+(jΩcs)2N1=s2N+(jΩc)2N(jΩc)2N
∴\therefore∴
Ha(s)Ha(−s)H_a(s)H_a(-s)Ha(s)Ha(−s)的极点为pk=Ωc[ejπ2N(2k+N−1)],k=1,2,...,2Np_k=\Omega_c[e^{j\frac \pi {2N}(2k+N-1)}],k=1,2,...,2Npk=Ωc[ej2Nπ(2k+N−1)],k=1,2,...,2N
由于Ha(s)H_a(s)Ha(s)的因果稳定性,要求极点必须落在s平面左半平面,即极点满足π2<π(2k+N−12N)<3π2\frac \pi2<\pi(\frac{2k+N-1}{2N})<\frac{3\pi}22π<π(2N2k+N−1)<23π
化简得:12<k<N+12\frac 12<k<N+\frac1221<k<N+21,即k=1,2,…,N
因此,N阶巴特沃斯模拟低通滤波器N个极点
pk=Ωc[ejπ2N(2k+N−1)],k=1,2,...,Np_k=\Omega_c[e^{j\frac \pi {2N}(2k+N-1)}],k=1,2,...,N pk=Ωc[ej2Nπ(2k+N−1)],k=1,2,...,N
其系统函数为
Ha(s)=ΩcN∏k=1N(s−pk)H_a(s)=\frac{\Omega_c^N}{\prod_{k=1}^N(s-p_k)} Ha(s)=∏k=1N(s−pk)ΩcN
巴特沃斯滤波器的设计步骤
1、根据模拟滤波器的设计指标αp,Ωp,αs,Ωs\alpha_p,\Omega_p,\alpha_s,\Omega_sαp,Ωp,αs,Ωs确定滤波器的阶数N
N=lg[(100.1αp−1)/(100.1αs−1)]2lg(Ωp/Ωs)N=\frac{lg[(10^{0.1\alpha_p}-1)/(10^{0.1\alpha_s}-1)]}{2lg(\Omega_p/\Omega_s)} N=2lg(Ωp/Ωs)lg[(100.1αp−1)/(100.1αs−1)]
取N为比计算结果大的最小整数,就是巴特沃斯模拟低通滤波器的阶数
2、确定滤波器的3dB截止频率
Ωcp=Ωp100.1αp−12N,Ωcs=Ωs100.1αs−12N\Omega_{cp}=\frac{\Omega_p}{\sqrt[2N]{10^{0.1\alpha_p}-1}},\Omega_{cs}=\frac{\Omega_s}{\sqrt[2N]{10^{0.1\alpha_s}-1}} Ωcp=2N100.1αp−1Ωp,Ωcs=2N100.1αs−1Ωs
实际设计时,Ωc\Omega_cΩc可在Ωcp≤Ωc≤Ωcs\Omega_{cp}\leq \Omega_c \leq \Omega_{cs}Ωcp≤Ωc≤Ωcs范围内选择
3、求出N个极点,pk=Ωc[ejπ2N(2k+N−1)],k=1,2,...,Np_k=\Omega_c[e^{j\frac \pi {2N}(2k+N-1)}],k=1,2,...,Npk=Ωc[ej2Nπ(2k+N−1)],k=1,2,...,N,得到滤波器的系统函数Ha(s)H_a(s)Ha(s)。
2、切比雪夫低通滤波器
切比雪夫模拟低通滤波器的幅频响应在一个频带中具有等波纹特性
- Chebyshev−IfiltersChebyshev-I filtersChebyshev−Ifilters 通带等波纹响应
- Chebyshev−IIfiltersChebyshev-II filtersChebyshev−IIfilters 阻带等波纹响应
相同指标下,切比雪夫滤波器比巴特沃斯滤波器阶数低
介绍
切比雪夫I型滤波器的平方幅度响应为:
∣Ha(jΩ)∣2=11+ϵ2CN2(Ω/Ωp)|H_a(j\Omega)|^2=\frac1{1+\epsilon^2C_N^2(\Omega/\Omega_p)} ∣Ha(jΩ)∣2=1+ϵ2CN2(Ω/Ωp)1
ϵ\epsilonϵ是小于1的正数,表示滤波器通带内幅度波动的程度,ϵ\epsilonϵ愈大,幅度波动也愈大
Ωp\Omega_pΩp是通带截止频率,λ=Ω/Ωp\lambda=\Omega/\Omega_pλ=Ω/Ωp称为对Ωp\Omega_pΩp的归一化频率
CN(x)C_N(x)CN(x)是N阶切比雪夫多项式:
CN(x)={cos[Narccos(x)]∣x∣≤1cosh[Narccosh(x)]∣x∣>1C_N(x)=\begin{cases} cos[Narc cos(x)]&|x|\leq1\\ cosh[Narccosh(x)]&|x|>1 \end{cases} CN(x)={cos[Narccos(x)]cosh[Narccosh(x)]∣x∣≤1∣x∣>1
- 当∣x∣≤1|x|\leq1∣x∣≤1时,∣CN(x)∣≤1|C_N(x)|\leq1∣CN(x)∣≤1。CN(x)C_N(x)CN(x)在-1和+1之间振荡,振荡的次数与N成正比,具有等波纹性
- 当∣x∣≤1|x|\leq1∣x∣≤1时,ϵ2CN2(x)\epsilon^2C_N^2(x)ϵ2CN2(x)在0到ϵ2\epsilon^2ϵ2之间波动,函数1+ϵ2CN2(Ω/Ωp)1+\epsilon^2C_N^2(\Omega/\Omega_p)1+ϵ2CN2(Ω/Ωp)的倒数即是平方幅度响应∣Ha(jΩ)∣2|H_a(j\Omega)|^2∣Ha(jΩ)∣2,其在[0,Ωp][0,\Omega_p][0,Ωp]上x=(Ω/Ωp)<1x=(\Omega/\Omega_p)<1x=(Ω/Ωp)<1,因此具有等波纹波动,最大值为1,最小值为11+ϵ2\frac1{1+\epsilon^2}1+ϵ21
- 当∣x∣>1|x|>1∣x∣>1时,∣CN(x)∣>1|C_N(x)|>1∣CN(x)∣>1,且∣CN(x)∣|C_N(x)|∣CN(x)∣随着∣x∣|x|∣x∣的增加单调上升。在[Ωp,∞)[\Omega_p,\infty)[Ωp,∞)上x=(Ω/Ωp)>1x=(\Omega/\Omega_p)>1x=(Ω/Ωp)>1,1+ϵ2CN2(Ω/Ωp)1+\epsilon^2C_N^2(\Omega/\Omega_p)1+ϵ2CN2(Ω/Ωp)单调上升,因此∣Ha(jΩ)∣2|H_a(j\Omega)|^2∣Ha(jΩ)∣2单调减小
- N相同时,切比雪夫滤波器比巴特沃斯滤波器具有较窄的过渡带
切比雪夫滤波器的设计步骤
1、确定模拟滤波器的设计指标Ωp,Ωs,αp,αs\Omega_p,\Omega_s,\alpha_p,\alpha_sΩp,Ωs,αp,αs
2、确定参数ϵ\epsilonϵ
ϵ2=100.1αp−1\epsilon^2=10^{0.1\alpha_p}-1 ϵ2=100.1αp−1
3、求出滤波器的阶数N,取N为比计算结果大的最小整数,就是切比雪夫模拟低通滤波器的阶数
N=arccosh(100.1αs−1100.1αp−1)1/2/arccosh(Ωs/Ωp)N=arccos h(\frac{10^{0.1\alpha_s}-1}{10^{0.1\alpha_p}-1})^{1/2}/arccosh(\Omega_s/\Omega_p) N=arccosh(100.1αp−1100.1αs−1)1/2/arccosh(Ωs/Ωp)
4、由下式求归一化切比雪夫I型低通滤波器的极点
pnk=σk+jΩk,k=1,2,...,N,其中p_{nk}=\sigma_k+j\Omega_k,k=1,2,...,N,其中 pnk=σk+jΩk,k=1,2,...,N,其中
σk=−sinh(β)sin(2k−1)π2N,Ωk=−cosh(β)cos(2k−1)π2N,β=arsinh(1/ϵ)N\sigma_k=-sinh(\beta)sin\frac{(2k-1)\pi}{2N},\Omega_k=-cosh(\beta)cos\frac{(2k-1)\pi}{2N},\beta=\frac{arsinh(1/\epsilon)}N σk=−sinh(β)sin2N(2k−1)π,Ωk=−cosh(β)cos2N(2k−1)π,β=Narsinh(1/ϵ)
5、求归一化切比雪夫I型低通滤波器的系统函数Hn(s)H_n(s)Hn(s):
N为奇数时,Hn(s)=sinhβ1+sinhβ∏k=1N/2σk2+Ωk2s2−2σks+(σk2+Ωk2)H_n(s)=\frac{sinh\beta}{1+sinh\beta}\prod_{k=1}^{N/2}\frac{\sigma_k^2+\Omega_k^2}{s^2-2\sigma_ks+(\sigma_k^2+\Omega_k^2)}Hn(s)=1+sinhβsinhβ∏k=1N/2s2−2σks+(σk2+Ωk2)σk2+Ωk2
N为偶数时,Hn(s)=11+ϵ2∏k=1N/2σk2+Ωk2s2−2σks+(σk2+Ωk2)H_n(s)=\frac1{\sqrt{1+\epsilon^2}}\prod_{k=1}^{N/2}\frac{\sigma_k^2+\Omega_k^2}{s^2-2\sigma_ks+(\sigma_k^2+\Omega_k^2)}Hn(s)=1+ϵ21∏k=1N/2s2−2σks+(σk2+Ωk2)σk2+Ωk2
3、椭圆模拟低通滤波器
介绍
通带和阻带均具有等波纹响应,幅频响应特性与FIR等波纹滤波器类似
幅频响应平方函数∣Ha(jΩ)∣2=11+ϵ2UN2(Ω/Ωc)|H_a(j\Omega)|^2=\frac1{1+\epsilon^2U^2_N(\Omega/\Omega_c)}∣Ha(jΩ)∣2=1+ϵ2UN2(Ω/Ωc)1
其中UN(.)U_N(.)UN(.)是N阶雅克比椭圆函数
对比
一、当阶数相同时,对相同的通带最大衰减apa_pap,和阻带最小衰减asa_sas:
1、巴特沃斯滤波器具有单调下降的幅频特性,过渡带最宽。两种类型的切比雪夫滤波器的过渡带宽度相等,比巴特沃斯滤波器的过渡带窄,但比椭圆滤波器的过渡带宽
2、切比雪夫I型滤波器在通带具有等波纹幅频特性,过渡带和阻带是单调下降的幅频特性
3、切比雪夫II型滤波器的通带幅频响应几乎与巴特沃斯滤波器相同,阻带是等波纹幅频特性
4、椭圆滤波器的过渡带最窄,通带和阻带均是等波纹幅频特性
二、复杂性:
在满足相同的滤波器幅频响应指标条件下,巴特沃斯滤波器阶数最高,椭圆滤波器的阶数最低,而且阶数差别较大。所以,就满足滤波器幅频响应指标而言,椭圆滤波器的性价比最高,应用较广泛。
模拟到数字滤波器的转换
模拟滤波器Ha(s)H_a(s)Ha(s)到数字滤波器H(z)H(z)H(z)的转换通过s到z的复值映射来实现,需要满足:
- 因果稳定的模拟滤波器变成数字滤波器,仍应是因果稳定的,保证s平面的左半平面要映射到z平面的单位圆内部
- 数字滤波器的频率响应 应模仿模拟滤波器的频率特性,即s平面的虚轴映射为z平面的单位圆,相应的频率之间呈线性关系
Ha(s)H_a(s)Ha(s)从s平面转换到z平面的典型方法有:脉冲响应不变法
、双线性变换法
脉冲响应不变法
基本原理: 使数字滤波器的单位脉冲响应h(n)模仿模拟滤波器的单位冲激响应ha(t)h_a(t)ha(t)
具体步骤:
以理想采样信号h^a(t)\hat h_a(t)h^a(t)为桥梁,推导从模拟滤波器转换到数字滤波器时,s平面和z平面之间的映射关系:
设h^a(t)=∑n=−∞∞ha(t)δ(t−nT)\hat h_a(t)=\sum_{n=-\infty}^{\infty}h_a(t)\delta(t-nT)h^a(t)=∑n=−∞∞ha(t)δ(t−nT),拉氏变换得到:
H^a(s)=∫−∞∞h^a(t)e−stdt=∫−∞∞[∑n=−∞∞ha(t)δ(t−nT)]e−stdt=∑nha(nT)e−snT\hat H_a(s)=\int_{-\infty}^{\infty}\hat h_a(t)e^{-st}dt=\int_{-\infty}^{\infty}[\sum_{n=-\infty}^{\infty}h_a(t)\delta(t-nT)]e^{-st}dt=\sum_nh_a(nT)e^{-snT} H^a(s)=∫−∞∞h^a(t)e−stdt=∫−∞∞[n=−∞∑∞ha(t)δ(t−nT)]e−stdt=n∑ha(nT)e−snT
与H^a(s)=∑nh(n)e−snT=∑nh(n)z−n∣z=esT=H(z)∣z=esT\hat H_a(s)=\sum_nh(n)e^{-snT}=\sum_nh(n)z^{-n}|_{z=e^{sT}}=H(z)|_{z=e^{sT}}H^a(s)=∑nh(n)e−snT=∑nh(n)z−n∣z=esT=H(z)∣z=esT比较得:
理想采样信号的拉氏变换与相应的采样序列的Z变换之间的映射关系为:z=esTz=e^{sT}z=esT,即脉冲响应不变法对应的s和z平面之间的映射关系。
利用Ha(s)H_a(s)Ha(s)直接变换成H(z)H(z)H(z):
H(z)=Ha(s)∣1s−sk=11−eskTz−1,k=1,2,...,N,其中sk是模拟低通滤波器Ha(s)的极点H(z)=H_a(s)|_{\frac1{s-s_k}=\frac1{1-e^{s_kT}z^{-1}}},k=1,2,...,N,其中s_k是模拟低通滤波器H_a(s)的极点 H(z)=Ha(s)∣s−sk1=1−eskTz−11,k=1,2,...,N,其中sk是模拟低通滤波器Ha(s)的极点
脉冲响应不变法设计IIT数字滤波器的步骤:
1、确定采样间隔T
2、根据采样间隔T,将给定的数字滤波器频率指标转换为模拟滤波器的频率指标:Ωp=ωp/T,Ωs=ωs/T\Omega_p=\omega_p/T,\Omega_s=\omega_s/TΩp=ωp/T,Ωs=ωs/T
3、根据指标αp,Ωp,αs,Ωs\alpha_p,\Omega_p,\alpha_s,\Omega_sαp,Ωp,αs,Ωs,设计模拟滤波器Ha(s)H_a(s)Ha(s)
4、用脉冲响应不变法,将模拟滤波器Ha(s)H_a(s)Ha(s)变换为数字滤波器H(z)H(z)H(z):
H(z)=Ha(s)∣1s−sk=11−eskTz−1,k=1,2,...,N,其中sk是模拟低通滤波器Ha(s)的极点H(z)=H_a(s)|_{\frac1{s-s_k}=\frac1{1-e^{s_kT}z^{-1}}},k=1,2,...,N,其中s_k是模拟低通滤波器H_a(s)的极点 H(z)=Ha(s)∣s−sk1=1−eskTz−11,k=1,2,...,N,其中sk是模拟低通滤波器Ha(s)的极点
双线性变换法
脉冲响应不变法是从s平面到z平面的多值映射,会产生频谱混叠。
双线性变换法采用非线性频率压缩的方法,将整个频率轴压缩到±π/T\pm\pi/T±π/T之间,再用z=esTz=e^{sT}z=esT转换到z平面上,由于s平面与z平面是一一对应的关系,从而有效地避免了脉冲响应不变法中多值映射所引起的频谱混叠现象;双线性变换法将s的左半平面映射到z平面的单位圆内,因此也满足了对系统因果稳定性的要求。
双线性变换法介绍:
设模拟滤波器的系统函数是Ha(s),s=jΩH_a(s),s=j\OmegaHa(s),s=jΩ,经过非线性频率压缩后用H^a(s1),s1=jΩ1\hat H_a(s_1),s_1=j\Omega_1H^a(s1),s1=jΩ1表示
用正切变换实现频率压缩:Ω=2Ttan(12Ω1T)\Omega=\frac 2Ttan(\frac12\Omega_1T)Ω=T2tan(21Ω1T),T是采样间隔
当Ω1从−π/T\Omega_1从-\pi/TΩ1从−π/T经过0变化到π/T\pi/Tπ/T时,Ω\OmegaΩ则由−∞-\infty−∞经过0变化到+∞∞∞,实现了s平面上整个虚轴完全压缩到s1s_1s1平面上虚轴的±π/T\pm\pi/T±π/T之间的转换
用正切变换实现频率压缩:Ω=2Ttan(12Ω1T)\Omega=\frac 2Ttan(\frac12\Omega_1T)Ω=T2tan(21Ω1T):
根据欧拉公式的sinθ=ejθ−e−jθ2j,cosθ=ejθ+e−jθ2sin\theta=\frac{e^{j\theta}-e^{-j\theta}}{2j},cos\theta=\frac{e^{j\theta}+e^{-j\theta}}{2}sinθ=2jejθ−e−jθ,cosθ=2ejθ+e−jθ
所以jtanθ=ejθ−e−jθejθ+e−jθjtan\theta=\frac{e^{j\theta}-e^{-j\theta}}{e^{j\theta}+e^{-j\theta}}jtanθ=ejθ+e−jθejθ−e−jθ
因此
jΩ=2Tjtan(12Ω1T)=2TejΩ1T/2−e−jΩ1T/2ejΩ1T/2+e−jΩ1T/2=2T1−e−jΩ1T1+e−jΩ1Tj\Omega=\frac2Tjtan(\frac12\Omega_1T)=\frac2T\frac{e^{j\Omega_1T/2}-e^{-j\Omega_1T/2}}{e^{j\Omega_1T/2}+e^{-j\Omega_1T/2}}=\frac2T\frac{1-e^{-j\Omega_1T}}{1+e^{-j\Omega_1T}} jΩ=T2jtan(21Ω1T)=T2ejΩ1T/2+e−jΩ1T/2ejΩ1T/2−e−jΩ1T/2=T21+e−jΩ1T1−e−jΩ1T
s=jΩ=2T1−e−s1T1+e−s1T=2T1−z−11+z−1,其中z=es1Ts=j\Omega=\frac2T\frac{1-e^{-s_1T}}{1+e^{-s_1T}}=\frac2T\frac{1-z^{-1}}{1+z^{-1}},其中z=e^{s_1T} s=jΩ=T21+e−s1T1−e−s1T=T21+z−11−z−1,其中z=es1T
展开得T2sz+T2s−z+1=0\frac T2sz+\frac T2s-z+1=02Tsz+2Ts−z+1=0
如上图所示,由s平面到s1s_1s1平面的非线性频率压缩,使H^a(s1)\hat H_a (s_1)H^a(s1)带限于π/T\pi/Tπ/T,因此,再用脉冲响应不变法从s1s_1s1平面转换到z平面不可能产生频谱混叠现象
因此,模拟系统Ha(s)H_a(s)Ha(s)和对应的数字系统H(z)的关系为:H(z)=Ha(s)∣s=2T1−z−11+z−1H(z)=H_a(s)|_{s=\frac2T\frac{1-z^{-1}}{1+z^{-1}}}H(z)=Ha(s)∣s=T21+z−11−z−1
模拟频率Ω\OmegaΩ与数字频率ω\omegaω之间的关系:Ω=2Ttanω2\Omega=\frac 2Ttan\frac \omega 2Ω=T2tan2ω, 即s平面上的模拟频率Ω\OmegaΩ与z平面上的数字频率ω\omegaω成非线性正切关系
双线性变换法设计IIR数字滤波器的步骤:
1、确定参数T。为简单起见,可取T=2。通常取ΩcT=1\Omega_cT=1ΩcT=1
2、确定数字低通滤波器的技术指标:通带截止频率ωp\omega_pωp、通带衰减αp\alpha_pαp、阻带截止频率ωp\omega_pωp、阻带衰减αp\alpha_pαp。采用双线性变换法模数频率间的转换关系,将数字低通滤波器的指标转换成模拟低通滤波器的技术指标:Ω=2Ttan(12ω)\Omega=\frac 2Ttan(\frac 1 2\omega)Ω=T2tan(21ω)
3、按照模拟低通滤波器的技术指标设计模拟低通滤波器
4、将模拟滤波器Ha(s)H_a(s)Ha(s),通过下式从s平面转换到z平面,得到数字低通滤波器系统函数H(z):
H(z)=Ha(s)∣s=2T.1−z−11+z−1H(z)=H_a(s)|_{s=\frac2T. \frac{1-z^{-1}}{1+z^{-1}}} H(z)=Ha(s)∣s=T2.1+z−11−z−1
数字滤波器的频带变换
IIR滤波器间接设计法的两种思路:
从原型低通数字滤波器到各种类型的数字滤波器的频带变换也称作Z平面变换法
设原型低通数字滤波器的系统函数为HL(z)H_L(z)HL(z),所需类型数字滤波器的系统函数为H(z)H(z)H(z),这种变换是将HL(z)H_L(z)HL(z)的zLz_LzL平面映射到H(z)H (z)H(z)的z平面
用zL和zz_L和zzL和z来表示变换前后的z平面,假设从zLz_LzL平面到z平面的映射关系为:
zL−1=G(z−1),G(Z−1)必须是z−1的有理函数,这样H(z)才是可实现的z_L^{-1}=G(z^{-1}),G(Z^{-1})必须是z^{-1}的有理函数,这样H(z)才是可实现的 zL−1=G(z−1),G(Z−1)必须是z−1的有理函数,这样H(z)才是可实现的
H(z)=HL(z)∣zL−1=G(z−1),其中,zL平面的单位圆必须映射到z平面的单位圆上zL平面的单位圆内部必须映射到z平面的单位圆内部H(z)=H_L(z)|_{z_L^{-1}=G(z^{-1})}, 其中,z_L平面的单位圆必须映射到z平面的单位圆上\\z_L平面的单位圆内部必须映射到z平面的单位圆内部 H(z)=HL(z)∣zL−1=G(z−1),其中,zL平面的单位圆必须映射到z平面的单位圆上zL平面的单位圆内部必须映射到z平面的单位圆内部
如果能找到适当的从单位圆周到单位圆周的数字变换把原型低通数字的通带频带变换成其他类型数字滤波器所要求的通带频带,就等于完成了这种滤波器的设计
设wL和ww_L和wwL和w分别是zL和zz_L和zzL和z的频率变量,则在各自单位圆zL=ejwL和z=ejwz_L=e^{jw_L}和z=e^{jw}zL=ejwL和z=ejw上有:
∣zL−1∣=∣G(z−1)∣=∣G(e−jw)∣=1,zL是单位圆|z_L^{-1}|=|G(z^{-1})|=|G(e^{-jw})|=1,z_L是单位圆 ∣zL−1∣=∣G(z−1)∣=∣G(e−jw)∣=1,zL是单位圆
因此,函数G的幅频特性对所有的频率都为1,是一个全通函数
全通滤波器一般形式为:
zL−1=G(z−1)=±∏k=1Nz−1−pk∗1−pkz−1,其中N全通函数的阶数,pk是G(z−1)的极点且必须在单位圆内,z_L^{-1}=G(z^{-1})=\pm\prod_{k=1}^N\frac{z^{-1}-p_k^*}{1-p_kz^{-1}},\\其中N全通函数的阶数,p_k是G(z^{-1})的极点且必须在单位圆内, zL−1=G(z−1)=±k=1∏N1−pkz−1z−1−pk∗,其中N全通函数的阶数,pk是G(z−1)的极点且必须在单位圆内,
选择合适的N和pkp_kpk就可以得到各种类型的变换
数字信号处理(六)IIR数字滤波器的设计相关推荐
- 基于MATLAB的数字信号处理(4) IIR数字滤波器设计及软件实现
文章目录 一.实验目的 二.实验原理 三.设计IIR数字滤波器 一.实验目的 熟悉用双线性变换法设计 IIR 数字滤波器的原理与方法 学会调用 MATLAB 信号处理工具箱中滤波器设计函数(或滤波器设 ...
- 数字信号处理(七)FIR数字滤波器的设计
文章目录 FIR滤波器 线性相位FIR滤波器的条件及特点 线性相位FIR滤波器 线性相位条件 线性相位FIR滤波器幅度特性 类型1:h(n)=h(N-n-1),N=奇数 类型2:h(n)=h(N-n- ...
- IIR数字滤波器的设计及应用——MATLAB
一.实验目的 (1)熟悉双线性变换法和双重映射法设计IIR数字滤波器的原理与方法. (2)掌握IIR数字滤波器的MATLAB实现方法设计各种滤波器. (3)观察分析滤波器输入输出数据波形,理解数字滤波 ...
- MATLAB——IIR数字滤波器的设计
1.基础知识 1.1.数字滤波器设计的基本步骤 我们知道模拟滤波器的设计是数字滤波器的设计的基础.在学习数字信号处理的过程中,IIR数字滤波器的设计的步骤是 (1)确定采样间隔Ts或者采样频率fs. ...
- 一阶IIR数字滤波器的设计
一阶IIR数字滤波器的设计 最简单的低通滤波器传递函数入手 对原始滤波器的改造 低通变高通 低通变带通 高通变带阻 从模拟到数字,采用双线性变换,简单方便 频域分析 结论 最简单的低通滤波器传递函数入 ...
- 一阶shelf IIR数字滤波器的设计和实现
一阶shelf IIR数字滤波器的设计和实现 还是从最简单的低通滤波器传递函数入手 第一步:给低通增加一个增益控制旋钮(系数) 第二步:并联一个直通信号 y(n)= x(n) low变high 从模拟 ...
- 数字信号处理5:FIR滤波器设计
文章目录 1. 滤波器初识 2. 最直观的滤波方式:频域滤波 3. 傅里叶变换中的加窗 4. FIR滤波器设计 5. 总结 之前的一系列博客中,详细分解了从卷积到FFT的相关知识,不过那些属于理论,是 ...
- m基于Matlab的fir和iir数字滤波器的设计与仿真
目录 1.算法概述 2.仿真效果预览 3.MATLAB部分代码预览 4.完整MATLAB程序 1.算法概述 MATLAB系统供了许多工具箱(Toolbox),借助于信号处理工具箱(signal pro ...
- Matlab | 数字信号处理:用窗函数法设计FIR数字滤波器
========================================== 博主github:https://github.com/MichaelBeechan 博主CSDN:https:/ ...
最新文章
- 系统学习机器学习之监督学习
- linux l显示详细信息,fdisk -l显示信息详解
- POJ 2421 Constructing Roads MST kruskal
- mysqld_multi stop 不能停掉mysql
- HttpHandler HttpModule入门篇
- 在CentOS 7中安装Jetty服务器
- 数据分析方法(营销模型篇)--最全的9种分析方法
- 微信小程序地图和百度地图定位位置不一样
- 推荐一个220V控制12V的电路板继电器-220v降压控制继电器
- 解决宝塔面板开启HTTS后,Edge不能访问的问题
- ecshop系统前后台出现的几个错误修正升级到PHP5.6后
- 潜伏在前端巅峰 中同出来的前端秘笈
- 建立枚举类型Weekday,假设今日为星期二,计算10天后为星期几。
- 字节跳动小程序对接环信IM遇到的问题
- php google地图,google地图 标记地图代码
- flutter圆形进度和条形进度绘制
- PLC模拟输入滤波10种方法
- 通过网页来打开app指定页面
- 关于弱监督学习的详细介绍——A Brief Introduction to Weakly Supervised Learning
- Mybatis多表查询