Differences between four kinds of filters

摘要:本文研究了低通、高通、带通、带阻模拟滤波器和数字滤波器的频率响应特性,单位脉冲响应以及零极点分布图,并定性解释了零极点分布图与频率响应特性的关系。

1.Analog Filters

a.Low Pass

$w_p=0.4\pi; w_s=0.6\pi; A_p=3dB ; A_s=10dB $
H ( s ) = 2.232 s 3 + 2.614 s 2 + 3.416 s + 2.232 H(s)=\frac {2.232}{s^3+2.614s^2+3.416s+2.232} H(s)=s3+2.614s2+3.416s+2.2322.232

b.High Pass

w p = 0.6 π w s = 0.4 π ; A p = 3 d B ; A s = 10 d B ; w_p=0.6\pi w_s=0.4\pi;A_p=3dB;A_s=10dB; wp=0.6πws=0.4π;Ap=3dB;As=10dB;
H ( s ) = s 3 s 3 + 3.625 s 2 + 6.569 s + 5.953 H(s)=\frac {s^3}{s^3+3.625s^2+6.569s+5.953} H(s)=s3+3.625s2+6.569s+5.953s3

c.Band Pass

w p 1 = 0.4 π ; w p 2 = 0.6 π ; w s 1 = 0.3 π ; w s 2 = 0.7 π ; A p = 3 d B ; A s = 10 d B ; w_{p_1}=0.4\pi ;w_{p_2}=0.6\pi; w_{s_1}=0.3\pi;w_{s_2}=0.7\pi;A_p=3dB;A_s=10dB; wp1=0.4π;wp2=0.6π;ws1=0.3π;ws2=0.7π;Ap=3dB;As=10dB;
H ( s ) = 0.4196 s 2 s 4 + 0.9161 s 3 + 5.157 s 2 + 2.17 s + 5.611 H(s)=\frac {0.4196s^2}{s^4+0.9161s^3+5.157s^2+2.17s+5.611} H(s)=s4+0.9161s3+5.157s2+2.17s+5.6110.4196s2

d.Band Stop

w p 1 = 0.3 π ; w p 2 = 0.7 π ; w s 1 = 0.4 π ; w s 2 = 0.6 π ; A p = 3 d B ; A s = 10 d B ; w_{p_1}=0.3\pi;w_{p_2}=0.7\pi; w_{s_1}=0.4\pi;w_{s_2}=0.6\pi;A_p=3dB;A_s=10dB; wp1=0.3π;wp2=0.7π;ws1=0.4π;ws2=0.6π;Ap=3dB;As=10dB;
H ( s ) = s 4 + 4.737 s 2 + 5.611 s 4 + 1.539 s 3 + 5.922 s 2 + 3.646 s + 5.611 H(s)=\frac {s^4+4.737s^2+5.611}{s^4+1.539s^3+5.922s^2+3.646s+5.611} H(s)=s4+1.539s3+5.922s2+3.646s+5.611s4+4.737s2+5.611

Discussion

(1)比较四个滤波器的频率响应 ∣ H ( j ω ) ∣ |H(j\omega)| H(jω) ω ⊆ [ − 2 π , 2 π ] \omega\subseteq[-2\pi,2\pi] ω[2π,2π],横坐标坐标归一化 π \pi π
4types_Hw.jpg

(2)比较四个滤波器的单位脉冲响应 h ( t ) h(t) h(t)
4types_ht.jpg
(3)比较四个滤波器的零极点分布图
4types_pzmap.jpg
(4)定性分析零极点与频率响应 ∣ H ( j ω ) ∣ |H(j\omega)| H(jω)的关系。
考虑 s = j ω s=j\omega s=jω

  • 当低通时, ω \omega ω很小,所以分子不能有 s s s项,且必须有常数项,所以低通滤波器没有零点。
  • 当高通时, ω \omega ω很大,所以分母要和分子同阶,故零极点的数目相同才能保证高频时 ∣ H ( j ω ) ∣ |H(j\omega)| H(jω)接近1。
  • 分子分母同阶还有一个有趣的现象是,化简 H ( s ) H(s) H(s)时,会出现常数1,对应在 h ( t ) h(t) h(t) δ ( t ) \delta(t) δ(t),表明在高频时, h ( t ) h(t) h(t)的其他分量都与信号卷积为0,只有 δ ( t ) \delta(t) δ(t)的贡献,所以可以通过。这也是为什么 h ( t ) h(t) h(t)第2、4图中的 t = 0 t=0 t=0的位置为负的原因, δ ( t ) \delta(t) δ(t)画不出来,只能画出其他部分。
  • 带阻滤波器具有低通和高通的性质,所以有常数项,也有最高项;而带通滤波器没有低通和高通的性质,故没有常数项和最高项,但有平方项,故在某 ω \omega ω处存在极值,即为通过的频率带。

2.Digital Filters

a.Low Pass

w p = 0.4 π ; w s = 0.6 π ; A p = 3 d B ; A s = 10 d B ; w_p=0.4\pi ; w_s=0.6\pi;A_p=3dB;A_s=10dB; wp=0.4π;ws=0.6π;Ap=3dB;As=10dB;
H ( z ) = 0.2292 + 0.4584 z − 1 + 0.2292 z − 2 1 − 0.2675 z − 1 + 0.1843 z − 2 H(z)=\frac {0.2292+0.4584z^{-1}+0.2292z^{-2}}{1-0.2675z^{-1}+0.1843z^{-2}} H(z)=10.2675z1+0.1843z20.2292+0.4584z1+0.2292z2

b.High Pass

w p = 0.6 π ; w s = 0.4 π ; A p = 3 d B ; A s = 10 d B ; w_p=0.6\pi ; w_s=0.4\pi;A_p=3dB;A_s=10dB; wp=0.6π;ws=0.4π;Ap=3dB;As=10dB;
H ( z ) = 0.2292 − 0.4584 z − 1 + 0.2292 z − 2 1 + 0.2675 z − 1 + 0.1843 z − 2 H(z)=\frac {0.2292-0.4584z^{-1}+0.2292z^{-2}}{1+0.2675z^{-1}+0.1843z^{-2}} H(z)=1+0.2675z1+0.1843z20.22920.4584z1+0.2292z2

c.Band Pass

w p 1 = 0.4 π ; w p 2 = 0.6 π ; w s 1 = 0.3 π ; w s 2 = 0.7 π ; A p = 3 d B ; A s = 10 d B ; w_{p_1}=0.4\pi ;w_{p_2}=0.6\pi ; w_{s_1}=0.3\pi;w_{s_2}=0.7\pi;A_p=3dB;A_s=10dB; wp1=0.4π;wp2=0.6π;ws1=0.3π;ws2=0.7π;Ap=3dB;As=10dB;
H ( z ) = 0.0995 − 0.1989 z − 2 + 0.0995 z − 4 1 + 0.9316 z − 2 + 0.3294 z − 4 H(z)=\frac {0.0995-0.1989z^{-2}+0.0995z^{-4}}{1+0.9316z^{-2}+0.3294z^{-4}} H(z)=1+0.9316z2+0.3294z40.09950.1989z2+0.0995z4

d.Band Stop

w p 1 = 0.3 π ; w p 2 = 0.7 π ; w s 1 = 0.4 π ; w s 2 = 0.6 π ; A p = 3 d B ; A s = 10 d B ; w_{p_1}=0.3\pi ;w_{p_2}=0.7\pi ; w_{s_1}=0.4\pi;w_{s_2}=0.6\pi;A_p=3dB;A_s=10dB; wp1=0.3π;wp2=0.7π;ws1=0.4π;ws2=0.6π;Ap=3dB;As=10dB;
H ( z ) = 0.4733 + 0.9467 z − 2 + 0.4733 z − 4 1 + 0.6469 z − 2 + 0.2465 z − 4 H(z)=\frac {0.4733+0.9467z^{-2}+0.4733z^{-4}}{1+0.6469z^{-2}+0.2465z^{-4}} H(z)=1+0.6469z2+0.2465z40.4733+0.9467z2+0.4733z4

Discussion

(1)比较四个滤波器的频率响应 ∣ H ( j Ω ) ∣ |H(j\Omega)| H(jΩ) Ω ⊆ [ − 2 π , 2 π ] \Omega\subseteq[-2\pi,2\pi] Ω[2π,2π],横坐标坐标归一化 π \pi π
4types_Hz.jpg
(2)比较四个滤波器的单位脉冲响应 h [ k ] h[k] h[k]
4types_hk.jpg
(3)比较四个滤波器的零极点分布图
4types_zplane.jpg
(4)定性分析零极点与频率响应 ∣ H ( j Ω ) ∣ |H(j\Omega)| H(jΩ)的关系。
由零极点可以得到
H ( j Ω ) = ( z − p 1 ) ( z − p 2 ) ( z − z 1 ) ( z − z 2 ) H(j\Omega)=\frac{(z-p_1)(z-p_2)}{(z-z_1)(z-z_2)} H(jΩ)=(zz1)(zz2)(zp1)(zp2)
的形式。考虑 z = e j Ω z=e^{j\Omega} z=ejΩ,不同的 Ω \Omega Ω实际对应零极点图中单位圆上的不同位置:
4types_zplane2.jpg

∣ H ( j Ω ) ∣ |H(j\Omega)| H(jΩ)的大小实际上等价于
绿 色 线 段 的 长 度 乘 积 红 色 线 段 的 长 度 乘 积 \frac{绿色线段的长度乘积}{红色线段的长度乘积} 线绿线
两个零点表示乘两次。因此容易看出:

  • Ω \Omega Ω为低频0附近,图1图4的结果较大,表示低通和带阻滤波器可以低频通过。而图2图3的结果为零,表示高通和带通滤波器不能低频通过。
  • Ω \Omega Ω为高频 π \pi π附近,图1图3的结果为0,表示低通和带通滤波器不能高频通过。而图2图4的结果较大,表示高通和带阻滤波器可以高频通过。
  • 我们可以很容易从零点分布看出哪些频率的波不能通过。极点与x轴的位置会影响 A p , A s A_p,A_s Ap,As的衰减情况。
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值