滤波器设计
(空军工程大学2008研|杭州电子科技大学2017研)
某一低通数字滤波器的各种指标和参量要求如下:
-
巴特沃斯频率响应,采用双线性变换法设计;
-
当0<f<2.5 Hz时,衰减小于3dB;
-
当f<50 Hz时,衰减大于或等于40dB;
-
抽样颇率f=200 Hz
试确定系统函数 H ( z ) H(z) H(z)
注:巴特沃斯归一化滤波器分母多项式的系数 s N + a N − 1 s N − 1 + ⋯ + a 1 s + 1 s^N+a_{N-1}s^{N-1}+\cdots+a_1s+1 sN+aN−1sN−1+⋯+a1s+1
N N N | a 1 a_1 a1 | a 2 a_2 a2 | a 3 a_3 a3 |
---|---|---|---|
1 | 1 | ||
2 | 1.414 | ||
3 | 2 | 2 | |
4 | 2.613 | 3.414 | 2.613 |
解:
N
=
l
g
1
0
0.1
δ
p
−
1
1
0
0.1
δ
s
t
−
1
l
g
(
Ω
p
Ω
s
t
)
N=\frac{lg\sqrt{\frac{10^{0.1\delta_{p}}-1}{10^{0.1\delta_{st}}-1}}}{lg(\frac{{\Omega}_p}{\Omega_{st}})}
N=lg(ΩstΩp)lg100.1δst−1100.1δp−1
(1) 先把模拟频率转化为数字频率
ω p = 2 π f p t f c = π 40 \displaystyle \omega_{p}=\frac{2\pi f_{pt}}{f_c}=\frac{\pi}{40} ωp=fc2πfpt=40π
ω s = 2 π f s t f = π 2 \displaystyle \omega_s=\frac{2\pi f_{st}}{f_{}}=\frac{\pi}{2} ωs=f2πfst=2π
(2) 计算 N N N
Ω p Ω s t = t a n π 80 t a n π 4 = 0.03929 \frac{{\Omega}_p}{\Omega_{st}}=\frac{tan\frac{\pi}{80}}{tan\frac{\pi}{4}}=0.03929 ΩstΩp=tan4πtan80π=0.03929
故, N ≈ 1.43 N\approx1.43 N≈1.43 [ N ] = 2 [N]=2 [N]=2
(3) 计算 H ( s ) H(s) H(s)
因为3dB的截止频率为 Ω c \Omega_c Ωc
G ( s ) = 1 s 2 + 1.414 s + 1 G(s)=\frac{1}{s^2+1.414s+1} G(s)=s2+1.414s+11
H ( s ) = 1 ( s Ω c ) 2 + 1.414 ( s Ω c ) + 1 , Ω c = 15.7 H(s)=\frac{1}{(\frac{s}{\Omega_c})^2+1.414(\frac{s}{\Omega_c})+1},\Omega_c=15.7 H(s)=(Ωcs)2+1.414(Ωcs)+11,Ωc=15.7
H ( z ) = H ( s ) ∣ s = 2 T 1 − z − 1 1 + z − 1 H(z)=H(s)|_{s=\frac{2}{T}\frac{1-z^{-1}}{1+z^{-1}}} H(z)=H(s)∣s=T21+z−11−z−1