滤波器基础系列博客,传送门:

滤波器基础01——滤波器的种类与特性

滤波器基础02——滤波器的传递函数与性能参数

滤波器基础03——Sallen-Key滤波器、多反馈滤波器与Bainter陷波器

滤波器基础04——全通滤波器

滤波器基础05——巴特沃斯、切比雪夫与贝塞尔滤波器

滤波器基础06——滤波器设计软件


一. 电路元件的S域模型

通过引入拉式变化,将电路复杂的微分方程关系转化为简单的阻抗分压关系,大大简化了对电路频域性能的分析。

1.1 电阻R

电阻时域的伏安特性关系如下:
U R ( t ) = R i R ( t ) U_R\left( t \right) =Ri_R\left( t \right) UR(t)=RiR(t)
取拉式变换有:
U R ( s ) = R i R ( s ) U_R\left( s \right) =Ri_R\left( s \right) UR(s)=RiR(s)

这意味着电感在复频域或者说s域中保持不变,仍为R

1.2 电感L

流过电感的电流不能突变,对于含有初始值 i L ( 0 − ) i_L\left( 0_- \right) iL(0)的电感L,有如下关系:
u L ( t ) = L d i L ( t ) d t ↔ U L ( s ) = s L I L ( s ) − L i L ( 0 − ) ↔ I L ( s ) = U L ( s ) s L + i L ( 0 − ) s u_L\left( t \right) =L\frac{di_L\left( t \right)}{dt}\leftrightarrow U_L\left( s \right) =sLI_L\left( s \right) -Li_L\left( 0_- \right) \leftrightarrow I_L\left( s \right) =\frac{U_L\left( s \right)}{sL}+\frac{i_L\left( 0_- \right)}{s} uL(t)=LdtdiL(t)UL(s)=sLIL(s)LiL(0)IL(s)=sLUL(s)+siL(0)

对于电流初始值为0的电感,其在s域中变为sL

1.3 电容C

电容两端的电压不能突变,对应含有初始值 u C ( 0 − ) u_C\left( 0_- \right) uC(0)的电容C,有如下关系:
i C ( t ) = C d u C ( t ) d t ↔ I C ( s ) = s C U C ( s ) − C u C ( 0 − ) ↔ U C ( s ) = I C ( s ) s C + u C ( 0 − ) s i_C\left( t \right) =C\frac{du_C\left( t \right)}{dt}\leftrightarrow I_C\left( s \right) =sCU_C\left( s \right) -Cu_C\left( 0_- \right) \leftrightarrow U_C\left( s \right) =\frac{I_C\left( s \right)}{sC}+\frac{u_C\left( 0_- \right)}{s} iC(t)=CdtduC(t)IC(s)=sCUC(s)CuC(0)UC(s)=sCIC(s)+suC(0)
对于电压初始值为0的电容,其在s域中变为1/sC


二. 滤波器的传递函数、幅频特性、相频特性与Bode图

2.1 传递函数

滤波器的传递函数定义为输出电压的拉式变换与输入电压拉式变换的比值,为:
H ( s ) = U o ( s ) U i ( s ) H\left( s \right) =\frac{U_o\left( s \right)}{U_i\left( s \right)} H(s)=Ui(s)Uo(s)
以一阶低通有源滤波器为例,其电路如下图所示:

其传递函数为:
H ( s ) = − R 2 R 1 R 2 C 1 s + 1 H\left( s \right) =\frac{-\frac{R_2}{R_1}}{R_2C_1s+1} H(s)=R2C1s+1R1R2
A m = − R 2 / R 1 , n 1 = R 2 C 1 A_m=-{{R_2}\Big/{R_1}}\text{,}n_1=R_2C_1 Am=R2/R1n1=R2C1,传递函数可写为:
H ( s ) = A m 1 + n 1 s H\left( s \right) =\frac{A_m}{1+n_1s} H(s)=1+n1sAm
这就是一阶低通滤波器的通用传递函数。对于更加复杂的滤波器,都可以写成如下通用形式

H ( s ) = A m × 1 + m 1 s + m 2 s 2 + . . . + m m s m 1 + n 1 s + n 2 s 2 + . . . + n n s n H\left( s \right) =A_m\times \frac{1+m_1s+m_2s^2+...+m_ms^m}{1+n_1s+n_2s^2+...+n_ns^n} H(s)=Am×1+n1s+n2s2+...+nnsn1+m1s+m2s2+...+mmsm
其中, n ⩾ m n\geqslant m nm传递函数分母上s的最高次幂称为滤波器的阶数

根据上式我们可以写出如下多种滤波器的通用形式,如:

s = 0 s=0 s=0 s → ∞ s\rightarrow \infty s代入传递函数,就能理解低通、高通、带通、带阻的含义了。

2.2 幅频特性与相频特性

传递函数可写成实部加虚部的形式:
H ( j ω ) = X ( ω ) + j Y ( ω ) H\left( j\omega \right) =X\left( \omega \right) +jY\left( \omega \right) H()=X(ω)+jY(ω)
输出信号与输入信号幅值比 A ( ω ) A\left( \omega \right) A(ω)称为幅频特性,有:
A ( ω ) = A o / A i = ∣ H ( j ω ) ∣ = X 2 ( ω ) + Y 2 ( ω ) A\left( \omega \right) =A_o/A_i=\left| H\left( j\omega \right) \right|=\sqrt{X^2\left( \omega \right) +Y^2\left( \omega \right)} A(ω)=Ao/Ai=H()=X2(ω)+Y2(ω)

输出信号与输入信号相位之差 φ ( ω ) \varphi \left( \omega \right) φ(ω)为相频特性,有:
φ ( ω ) = φ o − φ i = ∠ H ( j ω ) = { a r c tan ⁡ ( Y ( ω ) / X ( ω ) ) , X ( ω ) > 0 180 ° +    a r c tan ⁡ ( Y ( ω ) / X ( ω ) ) , X ( ω ) < 0 90 ° , X ( ω ) = 0 , Y ( ω ) > 0 − 90 ° , X ( ω ) = 0 , Y ( ω ) < 0 \varphi \left( \omega \right) =\varphi _o-\varphi _i=\angle H\left( j\omega \right) =\begin{cases} \mathrm{arc}\tan \left( Y\left( \omega \right) /X\left( \omega \right) \right) ,& X\left( \omega \right) >0\\ 180\degree +\,\,\mathrm{arc}\tan \left( Y\left( \omega \right) /X\left( \omega \right) \right),& X\left( \omega \right) <0\\ 90\degree,& X\left( \omega \right) =0\text{,}Y\left( \omega \right) >0\\ -90\degree,& X\left( \omega \right) =0\text{,}Y\left( \omega \right) <0\\ \end{cases} φ(ω)=φoφi=H()= arctan(Y(ω)/X(ω)),180°+arctan(Y(ω)/X(ω)),90°,90°,X(ω)>0X(ω)<0X(ω)=0Y(ω)>0X(ω)=0Y(ω)<0

2.3 Bode图

横坐标按 lg ⁡ ω \lg \omega lgω分度,单位为rad/s。幅频曲线纵坐标按:
d b A ( w ) = 20 lg ⁡ ( A ( ω ) ) dbA\left( w \right) =20\lg \left( A\left( \omega \right) \right) dbA(w)=20lg(A(ω))
线性分度,单位是分贝(dB)。对数相频曲线纵坐标按 φ ( ω ) \varphi \left( \omega \right) φ(ω)线性分度,单位为°。由此构成的坐标系称为半对数坐标系。构成的图称为Bode图。

根据传递函数就能画出滤波器的伯德图,可利用similink快速绘制,具体步骤可参考:


三. 滤波器的特征频率 F 0 F_0 F0、品质因数Q与阻尼系数 α \alpha α

3.1 特征频率

在对滤波器长久的研究中,人们发现传递函数上存在一特殊的频率点对滤波器特性影响很大,称其为特征频率,符号为 ω 0 \omega _0 ω0,单位为rad/s,对应以Hz为单位的特征频率为 F 0 F_0 F0。这两者都称为特征频率,表示的是同一个物理量,只是使用的单位不同。

对于一阶、三阶等奇数阶滤波器,特征频率是使得传递函数分母的实部等于虚部的点

对于二阶、四阶等偶数阶滤波器,特征频率是使得传递函数分母的实部为0的点

奇数阶以一阶低通滤波器为例,其传递函数为:
H ( s ) = A m × 1 1 + n 1 s H\left( s \right) =A_m\times \frac{1}{1+n_1s} H(s)=Am×1+n1s1
奇数阶滤波器特征频率 ω 0 \omega _0 ω0是满足传递函数的分母实部等于虚部的点,有:
1 = n 1 ω 0 ⇒ ω 0 = 1 / n 1 1=n_1\omega _0\Rightarrow \omega _0=1/n_1 1=n1ω0ω0=1/n1
偶数阶以二阶低通滤波器为例,其传递函数为
H ( s ) = A m × 1 1 + n 1 s + n 2 s 2 H\left( s \right) =A_m\times \frac{1}{1+n_1s+n_2s^2} H(s)=Am×1+n1s+n2s21
代入 s = j ω s=j\omega s=到上式,有:
H ( j ω ) = A m × 1 1 − n 2 ω 2 + j n 1 ω H\left( j\omega \right) =A_m\times \frac{1}{1-n_2\omega ^2+jn_1\omega} H()=Am×1n2ω2+jn1ω1
偶数阶滤波器特征频率 ω 0 \omega _0 ω0是满足传递函数的分母实部为0的点,有:
1 − n 2 ω 0 2 → ω 0 = 1 n 2 1-n_2{\omega _0}^2\rightarrow \omega _0=\frac{1}{\sqrt{n_2}} 1n2ω02ω0=n2 1

3.2 品质因数

另一个影响滤波器特性的参数称为滤波器品质因素,用字母Q表示,它定义为传递函数在特征频率处的幅值与 A m A_m Am的比值。

对于上述二阶低通滤波器,有:
Q = A m n 1 ω 0 / A m = 1 n 1 ω 0 Q=\frac{A_m}{n_1\omega _0}/A_m=\frac{1}{n_1\omega _0} Q=n1ω0Am/Am=n1ω01

3.3 阻尼系数

滤波器阻尼系数 α \alpha α与品质因数Q表示的是同一个物理量,它们是互为倒数的关系,有:
α = 1 Q \alpha =\frac{1}{Q} α=Q1


四. 归一化滤波器

还是以低通二阶滤波器为例,在已知特征频率 ω 0 \omega _0 ω0和品质因素Q之后,可得:
n 1 = 1 Q ω 0 , n 2 = 1 ω 0 2 n_1=\frac{1}{Q\omega _0}\text{,}n_2=\frac{1}{{\omega _0}^2} n1=Qω01n2=ω021
可见,在引入特征频率与品质因数后,就能将传递函数分母的系数用这两个有实际意义的物理量表示出来,有:

H ( j ω ) = A m × 1 1 + 1 Q ω 0 j ω + 1 ω 0 2 ( j ω ) 2 = A m × 1 1 + 1 Q j ω ω 0 + ( j ω ω 0 ) 2 H\left( j\omega \right) =A_m\times \frac{1}{1+\frac{1}{Q\omega _0}j\omega +\frac{1}{{\omega _0}^2}\left( j\omega \right) ^2}=A_m\times \frac{1}{1+\frac{1}{Q}j\frac{\omega}{\omega _0}+\left( j\frac{\omega}{\omega _0} \right) ^2} H()=Am×1+Qω01+ω021()21=Am×1+Q1jω0ω+(jω0ω)21
x = ω / ω 0 x={{\omega}\Big/{\omega _0}} x=ω/ω0,那么上式可以改写为:
H ( j ω ) = A m × 1 1 + 1 Q j x + ( j x ) 2 H\left( j\omega \right) =A_m\times \frac{1}{1+\frac{1}{Q}jx+\left( jx \right) ^2} H()=Am×1+Q1jx+(jx)21
x为频率与特征频率的比值,称为相对频率,是一个无量纲数。x = 1就表示频率等于特征频率点,这就是滤波器的归一化传递函数,通过引入相对频率,使频率的变化转换为相对频率的变化,此时滤波器特性转化为仅和品质因数Q有关。

滤波器的归一化在设计滤波器,确定不同截止频率的元件参数值时十分有用,这里不再展开。


五. 滤波器的截止频率、中心频率与带宽

5.1 截止频率

滤波器的截止频率定义为幅值比通带幅值低3dB时的频率,记为 ω c \omega _c ωc f c f _c fc 。以一阶低通滤波器为例,其传递函数为:
H ( s ) = 1 1 + s H\left( s \right) =\frac{1}{1+s} H(s)=1+s1
通带幅值为:
d B A ( 0 ) = 20 lg ⁡ ( ∣ H ( 0 ) ∣ ) = 20 lg ⁡ ( 1 ) = 0 d B dBA\left( 0 \right) =20\lg \left( \left| H\left( 0 \right) \right| \right) =20\lg \left( 1 \right) =0dB dBA(0)=20lg(H(0))=20lg(1)=0dB
所以,截止频率处幅值的dB值为-3dB,有:
d B A ( ω c ) = − 3 d B = 20 lg ⁡ ( ∣ H ( j ω c ) ∣ ) = 20 lg ⁡ ( 1 1 + ω c 2 ) = − 10 lg ⁡ ( 1 + ω c 2 ) ⇒ 1 + ω c 2 = 1 0 0.3 ⇒ ω c = 1 0 0.3 − 1 = 1 dBA\left( \omega _c \right) =-3dB=20\lg \left( \left| H\left( j\omega _c \right) \right| \right) =20\lg \left( \frac{1}{\sqrt{1+{\omega _c}^2}} \right) =-10\lg \left( 1+{\omega _c}^2 \right) \\ \Rightarrow 1+{\omega _c}^2=10^{0.3}\Rightarrow \omega _c=\sqrt{10^{0.3}-1}=1 dBA(ωc)=3dB=20lg(H(jωc))=20lg(1+ωc2 1)=10lg(1+ωc2)1+ωc2=100.3ωc=100.31 =1
显然 20 ∗ lg ⁡ ( 0.707 ) = − 3 20*\lg \left( 0.707 \right) =-3 20lg(0.707)=3,所以截止频率处的幅值为通带幅值的0.707倍。功率比 = 幅值比的平方,截止频率处信号功率为通带功率的 0.707 X 0.707 = 0.5倍,所以截止频率也可以描述为功率下降为通带功率一半的频率点

一般来说,截止频率与特征频率是不同的频率点,但对于上述一阶低通滤波器,两频率是一个点。

实际设计滤波器时,通常会先预设好截止频率,所以,明确截止频率、特征频率、品质因数这三者的关系就能在确定截止频率后,快速确定后两者的数值,从而确定滤波器参数。

定义截止频率与特征频率的比值为K,有:
K = ω c ω 0 K=\frac{\omega _c}{\omega _0} K=ω0ωc
根据二阶低通滤波器传递函数:
H ( j ω ) = A m × 1 1 + 1 Q j ω ω 0 + ( j ω ω 0 ) 2 H\left( j\omega \right) =A_m\times \frac{1}{1+\frac{1}{Q}j\frac{\omega}{\omega _0}+\left( j\frac{\omega}{\omega _0} \right) ^2} H()=Am×1+Q1jω0ω+(jω0ω)21

∣ H ( j ω c ) ∣ = ∣ A m × 1 1 + 1 Q j ω c ω 0 + ( j ω c ω 0 ) 2 ∣ = ∣ A m × 1 1 − K 2 + j K Q ∣ = 1 2 ∣ A m ∣ ⇒ ( 1 − K 2 ) 2 + ( K Q ) 2 = 2 ⇒ K = 4 Q 2 − 2 + 4 − 16 Q 2 + 32 Q 4 2 Q \left| H\left( j\omega _c \right) \right|=\left| A_m\times \frac{1}{1+\frac{1}{Q}j\frac{\omega _c}{\omega _0}+\left( j\frac{\omega _c}{\omega _0} \right) ^2} \right|=\left| A_m\times \frac{1}{1-K^2+j\frac{K}{Q}} \right|=\frac{1}{\sqrt{2}}\left| A_m \right| \\ \Rightarrow \sqrt{\left( 1-K^2 \right) ^2+\left( \frac{K}{Q} \right) ^2}=\sqrt{2}\Rightarrow K=\frac{\sqrt{4Q^2-2+\sqrt{4-16Q^2+32Q^4}}}{2Q} H(jωc)= Am×1+Q1jω0ωc+(jω0ωc)21 = Am×1K2+jQK1 =2 1Am(1K2)2+(QK)2 =2 K=2Q4Q22+416Q2+32Q4

所以,对于二阶低通滤波器,在Q确定时,K的值也就确定了即截止频率与特征频率的比值就确定了。通过这一点,我们先预设截止频率,然后确定Q值,这时就能根据公式推出特征频率,进而确定传递函数。以下给出一些典型的Q值对应的K值。

Q0.40.50.60.70710.811.21.52.0
K0.42780.64360.827111.11461.27201.35901.42991.5510

5.2 中心频率

中心频率是仅在带通或带阻滤波器才存在的概念,用 ω c e n t e r \omega _{center} ωcenter f c e n t e r f_{center} fcenter表示,对于低通或高通滤波器不存在中心频率的概念。

以带通滤波器为例。带通滤波器有两个截止频率,均是功率降为通带功率0.5倍的频率点,低频率点记为 f L f_L fL,高频率点记为 f H f_H fH。带通滤波器的中心频率指的是通带的中心点,有:
中心频率 = f H + f L 2 \text{中心频率}=\frac{f_H+f_L}{2} 中心频率=2fH+fL
同理,带阻滤波器也有两个截止频率,中心频率定义为阻带的中心点,计算公式同上。

5.3 带宽

滤波器带宽指的是通带的宽度,英文BandWidth,记为BW。

对于低通滤波器,通带宽度等于截止频率。

对于高通滤波器,没有带宽的概念。如果是有源高通滤波器,受限与运放的带宽,频率太高的信号是无法通过的,此时高通实际上变为了带通,这时才有带宽的概念。

对于带通滤波器, 带宽 = f H − f L \text{带宽}=f_H-f_L 带宽=fHfL,即两个截止频率的差值。

对于带阻滤波器,不存在通带带宽的概念,只存在阻带带宽的概念。


六. 滤波器的增益与衰减

滤波器增益指的是通带内输出幅值与输入幅值的比,用A表示。

对于低通滤波器,增益为传递函数在s = 0这一点的幅值。
二阶低通滤波器: H ( s ) = A m × 1 1 + n 1 s + n 2 s 2 \text{二阶低通滤波器:}H\left( s \right) =A_m\times \frac{1}{1+n_1s+n_2s^2} 二阶低通滤波器:H(s)=Am×1+n1s+n2s21

A = ∣ H ( 0 ) ∣ = ∣ A m × 1 1 + n 1 0 + n 2 0 2 ∣ = ∣ A m ∣ A=\left| H\left( 0 \right) \right|=\left| A_m\times \frac{1}{1+n_10+n_20^2} \right|=\left| A_m \right| A=H(0)= Am×1+n10+n2021 =Am

对于高通滤波器,增益为传递函数在 s → ∞ s\rightarrow \infty s时的幅值。
二阶高通滤波器: H ( s ) = A m × s 2 1 + n 1 s + n 2 s 2 \text{二阶高通滤波器:}H\left( s \right) =A_m\times \frac{s^2}{1+n_1s+n_2s^2} 二阶高通滤波器:H(s)=Am×1+n1s+n2s2s2

A = ∣ H ( ∞ ) ∣ = ∣ A m × ∞ 2 1 + n 1 ∞ + n 2 ∞ 2 ∣ = ∣ A m n 2 ∣ A=\left| H\left( \infty \right) \right|=\left| A_m\times \frac{\infty ^2}{1+n_1\infty +n_2\infty ^2} \right|=\left| \frac{A_m}{n_2} \right| A=H()= Am×1+n1+n222 = n2Am

对于带通滤波器,增益为传递函数在频率处于两个截止频率[ f L f_L fL, f H f_H fH]范围内的幅值。可用传递函数在频率为中心频率 f c f_c fc这一点的幅值来表示。
二阶带通滤波器: H ( s ) = A m × m 1 s 1 + n 1 s + n 2 s 2 \text{二阶带通滤波器:}H\left( s \right) =A_m\times \frac{m_1s}{1+n_1s+n_2s^2} 二阶带通滤波器:H(s)=Am×1+n1s+n2s2m1s

ω c e n t e r = ω H + ω L 2 \omega _{center}=\frac{\omega _H+\omega _L}{2} ωcenter=2ωH+ωL

在中心频率点处,传递函数分母实部为0,所以对于二阶带通滤波器,中心频率=特征频率,有:
A = ∣ H ( ω c e n t e r ) ∣ = ∣ A m × m 1 j ω c e n t e r 1 + n 1 j ω c e n t e r + n 2 ( j ω c e n t e r ) 2 ∣ = ∣ A m m 1 ∣ A=\left| H\left( \omega _{center} \right) \right|=\left| A_m\times \frac{m_1j\omega _{center}}{1+n_1j\omega _{center}+n_2\left( j\omega _{center} \right) ^2} \right|=\left| A_mm_1 \right| A=H(ωcenter)= Am×1+n1jωcenter+n2(jωcenter)2m1jωcenter =Amm1

对于带阻滤波器,不存在增益的说法,只会说阻带衰减。其传递函数如下:
二阶带阻滤波器: H ( s ) = A m × 1 + m 2 s 2 1 + n 1 s + n 2 s 2 \text{二阶带阻滤波器:}H\left( s \right) =A_m\times \frac{1+m_2s^2}{1+n_1s+n_2s^2} 二阶带阻滤波器:H(s)=Am×1+n1s+n2s21+m2s2

1 + m 2 s 2 = 0 ⇒ ω = 1 / m 2 1+m_2s^2=0\Rightarrow \omega =1/\sqrt{m_2} 1+m2s2=0ω=1/m2

显然,分子为0时对应的衰减是最大的,直接衰减到0,但实际情况下,因为元件都不是理想的,衰减永远不可能到真的到0,只能近似到0。


七. 滤波器的相移、延时与群延时

输出正弦波和输入正弦波之间,会存在相移,也就相应的存在延时。当输入波形为复合波形,包含多个不同频率的正弦波,且它们的相移没有特殊的规律,就可能出现不同的延时——这会引起输出波形与输入波形形态不同,造成线性失真。

群延时,Group delay,是一个随频率变化的量,用以描述一个滤波器对不同频率的输入信号,产生的延时是否一致。其数学定义是:
G D = − d φ ( ω ) d ω GD=-\frac{d\varphi \left( \omega \right)}{d\omega} GD=dωdφ(ω)
即相移对角频率求导的负数。之所以取负数,是因为一般都将延时量定义为正值,而这个延时由小于 0 的滞后相移引起。

对于每一个输入信号来说,滤波器产生的相移与滤波器产生的延时之间存在以下关系:
t d e l a y ( ω ) = φ ( ω ) 2 π × T = φ ( ω ) 2 π f = φ ( ω ) ω t_{delay}\left( \omega \right) =\frac{\varphi \left( \omega \right)}{2\pi}\times T=\frac{\varphi \left( \omega \right)}{2\pi f}=\frac{\varphi \left( \omega \right)}{\omega} tdelay(ω)=2πφ(ω)×T=2πfφ(ω)=ωφ(ω)
为保证信号不失真,就需要在频率增大时同步增大相移,则延时就会相等。但频率是无限增加了而相移最大就是超前180°或是滞后180°,不能无限增加,所以延时无法一直相等。因此,我们只关心较低频率处的相移与延时特性。


八. 滤波器的灵敏度

滤波电路由许多元件构成,每个元件参数值的变化都会影响滤波器的性能。滤波器某一性能指标y对某一元件参数x变化的灵敏度记作Sxy,定义为: S x y = ( d y / y ) / ( d x / x ) Sxy=(dy/y)/(dx/x) Sxy=(dy/y)/(dx/x)

该灵敏度与测量仪器或电路系统灵敏度不是一个概念,该灵敏度越小,标志着电路容错能力越强,稳定性也越高


九. 参考

《ADI新概念模拟电路》—— 4.运放电路的频率特性和滤波器,作者杨建国

《ADI滤波器设计教程》—— ADI智库


徐晓康的博客持续分享高质量硬件、FPGA与嵌入式知识,软件,工具等内容,欢迎大家关注。

Logo

旨在为数千万中国开发者提供一个无缝且高效的云端环境,以支持学习、使用和贡献开源项目。

更多推荐