我是靠谱客的博主 会撒娇火车,最近开发中收集的这篇文章主要介绍滤波器设计之巴特沃斯滤波器巴特沃斯滤波器简介,觉得挺不错的,现在分享给大家,希望可以做个参考。

概述

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

  • 巴特沃斯滤波器简介
    • 巴特沃斯滤波器设计步骤
    • 巴特沃斯低通滤波器设计实例
    • 注意事项

巴特沃斯滤波器简介

巴特沃斯低通滤波器的幅度平方函数定义为:
∣ H ( j λ ) ∣ = 1 1 + C 2 λ 2 N left | H(jlambda ) 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}=2pi f_{p} = 2pi times 5000,alpha _{p} = 3dB,Omega _{s}=2pi f_{s}=2pi 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.236times 10^{4}pi s^{4}+5.236times 10^{8}pi ^{2}s^{3}+5.236times 10^{12}pi ^{3}s^{2}+3.236times 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,但是在计算之中已经将符号考虑在内,所以不用代符号计算。

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

最后

以上就是会撒娇火车为你收集整理的滤波器设计之巴特沃斯滤波器巴特沃斯滤波器简介的全部内容,希望文章能够帮你解决滤波器设计之巴特沃斯滤波器巴特沃斯滤波器简介所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。

本图文内容来源于网友提供,作为学习参考使用,或来自网络收集整理,版权属于原作者所有。
点赞(37)

评论列表共有 0 条评论

立即
投稿
返回
顶部