滤波器设计之巴特沃斯滤波器

巴特沃斯滤波器简介

巴特沃斯低通滤波器的幅度平方函数定义为:
∣ H ( j λ ) ∣ = 1 1 + C 2 λ 2 N \left | H(j\lambda ) \right |= \frac{1}{1+C^{2}\lambda ^{2N}} H(jλ)=1+C2λ2N1
其中C为一常数参数,N为滤波器阶数, λ \lambda λ为归一化低通截止频率。
λ = Ω Ω p \lambda = \frac{\Omega }{\Omega{p}} λ=ΩpΩ

巴特沃斯滤波器设计步骤

  1. 设计巴特沃斯低通滤波器指标
    Ω p \Omega_{p} Ωp:通带截止频率;
    α p \alpha_{p} αp : 通带最小衰减,单位dB;
    Ω s \Omega_{s} Ωs :阻带开始频率;
    α s \alpha_{s} αs :阻带最大衰减,单位dB;
  2. 计算归一化频率
    λ p = Ω p Ω p = 1 , λ s = Ω s Ω p \lambda_{p} = \frac{\Omega_{p}}{\Omega _{p}} = 1 ,\lambda _{s}= \frac{\Omega _{s}}{\Omega _{p}} λp=ΩpΩp=1,λs=ΩpΩs
    α = 3 d B \alpha = 3dB α=3dB时, Ω p = Ω C \Omega _{p}=\Omega _{C} Ωp=ΩC为通常意义上的截止频率。
  3. 根据设计要求求出阶次N和参数C
    C 2 = 1 0 α p 10 − 1 C^{2}=10^{\frac{\alpha _{p}}{10}}-1 C2=1010αp1
    N = l g ( a ) l g ( λ s ) N = \frac{lg (a)}{lg(\lambda _{s})} N=lg(λs)lg(a)其中
    a = 1 0 α s 10 − 1 1 0 α p 10 − 1 a=\sqrt{\frac{10^{\frac{\alpha _{s}}{10}}-1}{10^{\frac{\alpha _{p}}{10}}-1}} a=1010αp11010αs1
    注意当 α p = 3 d B \alpha _{p} = 3dB αp=3dB 时 C=1
  4. 利用N查表
    通过N值查表获得归一化巴特沃斯低通滤波器的系统函数;
    在这里插入图片描述
  5. 变换成需要的滤波器
    低通滤波器: H ( s ) = H ( p = s Ω p ) H(s)=H(p=\frac{s}{\Omega _{p}}) H(s)=H(p=Ωps)

巴特沃斯低通滤波器设计实例

设计一个巴特沃斯低通滤波器,要求截止频率 f p = 5000 H z f_{p} = 5000Hz fp=5000Hz,通带最大衰减 α p = 3 d B \alpha _{p}=3dB αp=3dB,阻带起始频率 f s = 10000 H z f_{s} = 10000Hz fs=10000Hz,阻带最小衰减 α s = 30 d B \alpha _{s}=30dB αs=30dB

解:已知 Ω p = 2 π f p = 2 π × 5000 , α p = 3 d B , Ω s = 2 π f s = 2 π × 10000 , α s = 30 d B \Omega _{p}=2\pi f_{p} = 2\pi \times 5000,\alpha _{p} = 3dB,\Omega _{s}=2\pi f_{s}=2\pi \times 10000,\alpha _{s} = 30dB Ωp=2πfp=2π×5000,αp=3dB,Ωs=2πfs=2π×10000,αs=30dB

  1. 计算归一化频率
    λ p = Ω p Ω p = 1 , λ s = Ω s Ω p = 2 \lambda _{p}=\frac{\Omega _{p}}{\Omega _{p}}=1,\lambda _{s}=\frac{\Omega _{s}}{\Omega _{p}}=2 λp=ΩpΩp=1,λs=ΩpΩs=2
  2. 计算巴特沃斯滤波器阶次N和参数C
    C 2 = 1 0 α P 10 − 1 = 1 0 0.3 − 1 = 1 C^{2}=10^{\frac{\alpha _{P}}{10}}-1=10^{0.3}-1=1 C2=1010αP1=100.31=1
    a = 1 0 α s 10 − 1 1 0 α p 10 − 1 = 1 0 3 − 1 1 0 0.3 − 1 = 31.637 a=\sqrt{\frac{10^{\frac{\alpha _{s}}{10}}-1}{10^{\frac{\alpha _{p}}{10}}-1}}=\sqrt{\frac{10^{3}-1}{10^{0.3}-1}}=31.637 a=1010αp11010αs1 =100.311031 =31.637
    N = l g ( a ) l g ( λ s ) = l g ( 31.637 ) l g ( 2 ) = 4.982 N=\frac{lg(a)}{lg(\lambda _{s})}=\frac{lg(31.637)}{lg(2)}=4.982 N=lg(λs)lg(a)=lg(2)lg(31.637)=4.982
    所以N选择5;
  3. 利用N查表获得归一化巴特沃斯低通滤波器的系统函数:
    H ( p ) = 1 p 5 + 3.2361 p 4 + 5.2361 p 3 + 5.2361 p 2 + 3.2361 p + 1 H(p)=\frac{1}{p^{5}+3.2361p^{4}+5.2361p^{3}+5.2361p^{2}+3.2361p+1} H(p)=p5+3.2361p4+5.2361p3+5.2361p2+3.2361p+11
    4.去掉归一化
    H ( s ) = H ( p ) → p = s Ω p H(s)=H(p)\overset{p=\frac{s}{\Omega _{p}}}{\rightarrow} H(s)=H(p)p=Ωps
    = 1 0 20 π 5 s 5 + 3.236 × 1 0 4 π s 4 + 5.236 × 1 0 8 π 2 s 3 + 5.236 × 1 0 12 π 3 s 2 + 3.236 × 1 0 16 π 4 s + 1 0 20 π 5 =\frac{10^{20}\pi ^{5}}{s^{5}+3.236\times 10^{4}\pi s^{4}+5.236\times 10^{8}\pi ^{2}s^{3}+5.236\times 10^{12}\pi ^{3}s^{2}+3.236\times 10^{16}\pi ^{4}s+10^{20}\pi ^{5}} =s5+3.236×104πs4+5.236×108π2s3+5.236×1012π3s2+3.236×1016π4s+1020π51020π5

注意事项

Ω p \Omega_{p} Ωp:通带截止频率;
α p \alpha_{p} αp : 通带最小衰减,单位dB;
Ω s \Omega_{s} Ωs :阻带开始频率;
α s \alpha_{s} αs :阻带最大衰减,单位dB;
这四个变量中,衰减其实真实值是负数,例如-3dB,-20dB,但是在计算之中已经将符号考虑在内,所以不用代符号计算。

高通滤波器低通滤波器之间存在一个简单的转换函数关系,后续将会继续跟进。

已标记关键词 清除标记
相关推荐
©️2020 CSDN 皮肤主题: 大白 设计师:CSDN官方博客 返回首页