控制理论拾遗:截止频率与闭环带宽

截止频率

在频域分析中,开环系统的截止频率和闭环系统的带宽都是常用的性能指标。教材通常会指出:闭环带宽越大,系统响应越快;开环截止频率越高,系统响应通常也越快。两个指标分别由开环和闭环频率特性定义,为什么都能反映系统的快速性?二者之间又有什么关系?

这里先说明一下术语。对于一般负反馈系统,前向通道和反馈通道的传递函数分别为G(s)G(s)和H(s)H(s),开环传递函数为G(s)H(s)G(s)H(s)。满足

|G(jωc)H(jωc)|=1 |G(j\omega_c)H(j\omega_c)|=1

的频率ωc\omega_c,国内教材中常见「截止频率」「剪切频率」「幅值穿越频率」等称呼,英文教材通常称为增益穿越频率(gain crossover frequency)。为了与国内传统自动控制原理教材保持一致,本文统一称ωc\omega_c为截止频率。

对于典型的低通型闭环系统,幅频特性从低频平坦段首次下降至低频值的1/21/\sqrt2,即下降3dB3\ \mathrm{dB}时对应的频率ωB\omega_B,通常称为闭环带宽。它和开环截止频率是两个不同的概念。

闭环带宽与响应速度之间的关系很好理解。带宽越大,系统能够在更宽的频率范围内保持较小的幅值衰减,输入中较快变化的频率成分也能够传递到输出,因此时域响应通常更快。

但截止频率ωc\omega_c是根据开环频率特性定义的。为什么一个开环指标越大,闭环响应通常也越快?这个结论背后的原因并不像闭环带宽那样显而易见。

反馈作用的频率范围

为了把这个问题的基本机理讲清楚,下面先考虑最常见的单位负反馈系统,即H(s)=1H(s)=1。此时开环传递函数就是G(s)G(s),输入、输出和误差信号分别为R(s)R(s)、C(s)C(s)和E(s)E(s),有

Φ(s)=C(s)R(s)=G(s)1+G(s),E(s)R(s)=11+G(s) \Phi(s)=\frac{C(s)}{R(s)}=\frac{G(s)}{1+G(s)},\qquad \frac{E(s)}{R(s)}=\frac{1}{1+G(s)}

截止频率的定义也相应简化为|G(jωc)|=1|G(j\omega_c)|=1。理解截止频率为什么与快速性有关,关键在于弄清楚|G(jω)||G(j\omega)|究竟意味着什么。

设输入中含有一个频率为ω\omega的正弦分量。如果希望输出很好地跟踪它,那么输出中也必须产生一个大小相近的同频分量。在单位负反馈系统中,这个输出不是凭空产生的,而是由误差信号E=R−CE=R-C经过前向通道G(s)G(s)产生,即

C(jω)=G(jω)E(jω) C(j\omega)=G(j\omega)E(j\omega)

这里存在一个很基本的矛盾:误差必须保留到足以通过前向通道产生所需输出的程度,但从跟踪性能来看,误差又应该尽可能小。

因此,一个好的反馈系统真正希望做到的是:在需要跟踪的频率上,只需很小的误差,就能驱动出足够大的输出。这里说的误差和输出,都是指该频率上的同频分量。而这恰恰要求该频率上的|G(jω)||G(j\omega)|足够大。从这个角度看,|G(jω)||G(j\omega)|的物理意义就很直观了:它表示在频率ω\omega处,单位幅值的误差同频分量经过前向通道后,能够产生多大的输出同频分量。需要说明的是,这里强调的是动态驱动关系,若开环含积分环节,在直流稳态下所需误差可以趋于零,此时输出由积分器状态维持,而非由当前误差瞬时产生。

如果|G(jω)|=100|G(j\omega)|=100,那么幅值仅为0.010.01的误差,经过前向通道后就能够产生幅值约为1的输出。对于幅值为1的输入来说,只需要留下很小的误差,就足以维持所需的输出。此时

|ER|=|11+G|≈1|G|=0.01 \left|\frac{E}{R}\right| = \left|\frac{1}{1+G}\right| \approx\frac{1}{|G|} =0.01

同时|Φ(jω)|≈1|\Phi(j\omega)|\approx1。也就是说,在这个频率上,输入幅值为1时,误差只占大约1%1\%,而输出几乎能够完整跟随输入。

反过来,如果|G(jω)|=0.01|G(j\omega)|=0.01,即使误差的幅值已经达到1,也只能通过前向通道产生幅值约为0.010.01的输出。此时有|E/R|≈1|E/R|\approx1、|C/R|≈|G|=0.01|C/R|\approx|G|=0.01。换句话说,输出几乎没有跟上输入,绝大部分输入都表现为误差。

对于常见的低通型控制系统,开环幅值总体上会随着频率升高而下降。在非零频率上,低频段通常|G|≫1|G|\gg1,反馈能把同频误差压得很小;随着频率升高,|G||G|逐渐减小,系统为了产生同样幅值的输出,需要容许越来越大的误差;到了|G|≪1|G|\ll1的高频段,即使误差已经接近输入本身,产生的输出仍然很小,反馈实际上已经难以控制这么快的动态变化。

而截止频率正好满足|G(jωc)|=1|G(j\omega_c)|=1,位于上述两种状态的过渡区域。因此,ωc\omega_c可以理解为一个很有意义的频率尺度:

它大致标志着反馈还能有效控制到多高的频率。

截止频率越高,意味着较大的开环幅值能够维持到更高频率,反馈能够有效作用的频率范围通常也越宽。

二者的数学联系

上面的分析说明了截止频率为什么能够反映反馈作用的频率范围。进一步看,截止频率与闭环带宽之间还有非常直接的数学联系。

对于单位负反馈系统,有Φ=G/(1+G)\Phi=G/(1+G)。低频处|G|≫1|G|\gg1时Φ≈1\Phi\approx1,高频处|G|≪1|G|\ll1时Φ≈G\Phi\approx G,因此闭环幅频特性从通带进入衰减区的过渡自然发生在|G|≈1|G|\approx1附近,而截止频率ωc\omega_c正是|G(jωc)|=1|G(j\omega_c)|=1的频率。

在ωc\omega_c处,若开环相位为ϕc\phi_c,由于开环幅值恰好为1,可以写成G(jωc)=ejϕcG(j\omega_c)=e^{j\phi_c}。于是

|Φ(jωc)|=1|1+ejϕc|=12|cos⁡(ϕc/2)| |\Phi(j\omega_c)| = \frac{1}{|1+e^{j\phi_c}|} = \frac{1}{2|\cos(\phi_c/2)|}

相位裕度γ=180∘+ϕc\gamma=180^\circ+\phi_c,对于通常的正相位裕度,有

|Φ(jωc)|=12sin⁡(γ/2) |\Phi(j\omega_c)| = \frac{1}{2\sin(\gamma/2)}

这个关系说明,在开环截止频率处,闭环幅值还取决于相位裕度。

最典型的情况是γ=90∘\gamma=90^\circ。此时ϕc=−90∘\phi_c=-90^\circ,有|Φ(jωc)|=1/2|\Phi(j\omega_c)|=1/\sqrt2,正好对应−3dB-3\ \mathrm{dB}。如果闭环低频增益近似为1,那么此时恰有ωB=ωc\omega_B=\omega_c。如果相位裕度不是90∘90^\circ,二者一般不会严格相等。例如γ=60∘\gamma=60^\circ时,有|Φ(jωc)|=1|\Phi(j\omega_c)|=1,说明在ωc\omega_c处闭环幅值尚未下降到−3dB-3\ \mathrm{dB}。对于典型的低通型闭环系统,其−3dB-3\ \mathrm{dB}带宽通常还会高于ωc\omega_c。

对于相位裕度通常小于90∘90^\circ且频率响应比较规整的低通型系统,工程上常用

ωc≤ωB≲2ωc \omega_c\leq\omega_B\lesssim2\omega_c

来估计闭环带宽,即闭环带宽通常略高于截止频率,但一般不会相差太多。有些材料将上式写成ωc≤ωB≤2ωc\omega_c\leq\omega_B\leq2\omega_c,但它应理解为典型系统中的经验范围,而不是对任意反馈系统都成立的严格数学界限。

所以,截止频率并不等于闭环带宽,但对于频率响应比较规整且稳定裕度合理的系统,两者通常处于相近的频率范围。这就是为什么开环截止频率能够作为闭环响应速度的重要近似指标。

最简单的例子

考虑开环传递函数

G(s)=ωcs G(s)=\frac{\omega_c}{s}

它在ω=ωc\omega=\omega_c时恰好满足|G(jωc)|=1|G(j\omega_c)|=1,所以ωc\omega_c就是截止频率。

单位负反馈后,

Φ(s)=ωcs+ωc \Phi(s)=\frac{\omega_c}{s+\omega_c}

这是一个标准一阶系统,其时间常数为τ=1/ωc\tau=1/\omega_c,闭环带宽恰好也是ωc\omega_c,而10%∼90%10\%\sim90\%上升时间约为tr≈2.2/ωct_r\approx2.2/\omega_c。因此,在这个最简单的系统中,提高截止频率就是直接减小闭环时间常数,响应速度随之提高。

开环和闭环Bode图

这个例子还解释了一个经典的频域设计经验。实际系统通常希望开环幅频曲线在截止频率附近以约−20dB/dec-20\ \mathrm{dB/dec}的斜率穿过0dB0\ \mathrm{dB},而积分环节1/s1/s恰好具有−20dB/dec-20\ \mathrm{dB/dec}的幅值斜率和−90∘-90^\circ的相位。如果系统在截止频率附近具有一段较宽而规整的−20dB/dec-20\ \mathrm{dB/dec}中频段,其局部特性就很接近上面的简单模型。

非单位反馈系统

上面的讨论以单位负反馈为例。对于一般负反馈系统,前向通道和反馈通道分别为G(s)G(s)和H(s)H(s),开环传递函数为G(s)H(s)G(s)H(s),截止频率由|G(jωc)H(jωc)|=1|G(j\omega_c)H(j\omega_c)|=1定义。

因此,前面对反馈作用强弱的讨论可以相应推广。这里误差信号定义为E=R−HCE=R-HC,且有E/R=1/(1+GH)E/R=1/(1+GH)。在相应频率,当|GH|≫1|GH|\gg1时,有|E/R|≪1|E/R|\ll1,反馈能够有效压低误差,当|GH|≪1|GH|\ll1时,有|E/R|≈1|E/R|\approx1,反馈作用很弱。需要注意,在非单位反馈系统中,误差很小意味着HC≈RHC\approx R,即C≈R/HC\approx R/H,而不一定是输出CC直接跟踪输入RR。单位反馈H=1H=1正是其中最简单的特殊情况。

不过此时输入到输出的闭环传递函数为

Φ(s)=G(s)1+G(s)H(s)=1H(s)G(s)H(s)1+G(s)H(s) \Phi(s) = \frac{G(s)}{1+G(s)H(s)} = \frac{1}{H(s)} \frac{G(s)H(s)}{1+G(s)H(s)}

其中G(s)H(s)/(1+G(s)H(s))G(s)H(s)/(1+G(s)H(s))具有与单位反馈系统完全相同的结构,而前面的1/H(s)1/H(s)还会进一步影响输出端的闭环频率特性。

若H(s)H(s)在关注频段内近似为常数,则

Φ(s)≈1H⋅G(s)H(s)1+G(s)H(s). \Phi(s)\approx \frac{1}{H}\cdot\frac{G(s)H(s)}{1+G(s)H(s)}.

其中1/H1/H只是恒定缩放,不改变−3dB-3\ \mathrm{dB}带宽,因此前面关于截止频率与闭环带宽的解释仍基本成立。如果H(s)H(s)本身具有明显的频率特性,则1/H(s)1/H(s)还会改变闭环幅频特性的形状,此时不能简单地根据截止频率判断输出端的闭环带宽。

截止频率的权衡

截止频率越高,意味着希望反馈一直作用到更快的动态,但这并不是没有代价的。

随着频率升高,被控对象通常会出现越来越大的相位滞后,同时模型也越来越不准确。执行器和传感器的未建模动态以及测量噪声等问题在高频往往更加突出;时间延迟和右半平面零点等非最小相位因素,也会进一步限制能够安全实现的截止频率。

如果一味提高ωc\omega_c,把0dB0\ \mathrm{dB}交越点推入相位滞后严重的频段,相位裕度就可能明显下降,闭环出现较大的谐振峰和超调,甚至失稳。即使系统保持稳定,更高的带宽通常也会使系统对高频测量噪声和未建模动态更加敏感,并可能需要更大的控制作用。

所以更准确的说法不是「截止频率越大越好」,而是:

在保持足够稳定裕度和健壮性的前提下,较高的截止频率通常对应较宽的闭环带宽和更快的响应。

这也是为什么频域设计不能只关心开环幅频曲线在哪里穿过0dB0\ \mathrm{dB},还必须同时关心它以怎样的斜率和相位穿过去。

总结

闭环带宽回答的是「输出还能跟到多快」,开环截止频率回答的是「反馈环路在什么频率附近开始力不从心」。二者分别从结果和机制两个角度描述了系统从通带到衰减的过渡段,并不是同一个量。

对单位负反馈系统,远低于截止频率ωc\omega_c时开环增益|G|≫1|G|\gg1,只需较小的同频误差就能驱动出所需的同频输出,反馈作用强。高于ωc\omega_c后开环增益回落,反馈对快速变化的控制能力逐步下降,ωc\omega_c就是反馈控制能力开始退化的分界点。

对于幅频特性规整、相位裕度合理的系统,这个分界点和闭环−3dB-3\ \mathrm{dB}带宽频率相近,具体差值主要受ωc\omega_c处的相位裕度影响。一般负反馈系统要以开环传递函数G(s)H(s)G(s)H(s)为基准来判断。若反馈通道H(s)H(s)自身频率特性显著,闭环带宽需单独分析。

参考资料

  1. S. Skogestad, I. Postlethwaite, Multivariable Feedback Control: Analysis and Design, Chapter 2, Section 2.4.5 “Bandwidth and Crossover Frequency”. PDF
  2. K. J. Åström, R. M. Murray, Feedback Systems: An Introduction for Scientists and Engineers, Chapter 12: Frequency Domain Design. PDF
  3. R. M. Murray, Frequency Domain Design, CDS 101/110a Lecture 8-1, California Institute of Technology. PDF