学习目标
- 理解传递函数的基本性质,能把输入、系统、输出分开看。
- 掌握极点(pole)、零点(zero)和特征方程(characteristic equation)的定义。
- 能在 s 平面上画出极点位置,并判断系统的稳定性。
- 掌握 Routh 稳定性判据的构造方法,能在不求解特征方程的前提下判断稳定性。
- 知道 Routh 判据中”必要条件”和”充分条件”的区别,避免做题时漏步骤。
先用人话理解本章在讲什么
控制理论中最核心的问题其实就一个:这个系统稳不稳定?
一个系统不管结构多复杂,拿到手后第一件事不是去算输出波形,而是先判断它会不会”炸掉”——即输出会不会随时间无限增长。如果系统不稳定,后面做任何控制策略都是白费功夫。
本章的逻辑链条是这样的:
- 传递函数是描述系统的代数工具,它把微分方程变成了多项式之比。
- 传递函数的分母多项式叫特征方程,它的根——也就是极点——决定了系统的行为。
- 如果所有极点都在 s 平面的左半边,系统就是稳定的;只要有一个极点跑到右半边,系统就不稳定。
- 问题是:高阶多项式的根不好求。Routh 判据提供了一种不需要真的去解根、只靠系数就能判断稳定性的方法。
初学者最容易卡住的地方:
- 把极点和零点搞混。极点看分母,零点看分子,不要记反。
- 以为 Routh 判据”必要条件”(系数全正)过了就稳定——这远远不够,还要构造 Routh 表。
- 以为虚轴上的极点是”稳定”的——严格来说叫”临界稳定”,在工程上通常不接受。
核心概念
3.1 传递函数(Transfer Function)
一句话理解:
传递函数就是”系统对输入做了什么”的代数描述——把微分方程里的求导运算全部变成 s 的乘法,整个系统就变成了一个多项式之比。
正式定义:
对一个线性时不变(LTI)系统,假设初始条件为零,输出的 Laplace 变换与输入的 Laplace 变换之比就是传递函数:
G(s)=R(s)C(s)
其中 R(s) 是输入的 Laplace 变换,C(s) 是输出的 Laplace 变换。
直观例子:
考虑一阶微分方程:
dtdc(t)+2c(t)=r(t)
两边做 Laplace 变换(零初始条件):
sC(s)+2C(s)=R(s)
整理:
G(s)=R(s)C(s)=s+21
整个系统的动态特性,就浓缩在 s+21 这一个分式里。
传递函数的一般形式:
G(s)=ansn+an−1sn−1+⋯+a1s+a0bmsm+bm−1sm−1+⋯+b1s+b0
分子是 m 次多项式,分母是 n 次多项式,通常 n≥m(物理可实现系统的分母阶数不低于分子)。
容易混淆的点:
- 传递函数只描述系统本身,和输入信号无关。换一个输入,输出会变,但 G(s) 不变。
- 传递函数的前提是零初始条件,有非零初始条件时要另外处理。
- 串联的子系统可以直接把传递函数相乘,并联的可以相加。这是传递函数最方便的地方——把微分方程的卷积变成了代数乘法。
3.2 极点(Pole)与零点(Zero)
一句话理解:
极点是让传递函数”爆炸”(绝对值趋于无穷)的 s 值,零点是让传递函数”消失”(等于零)的 s 值。
正式定义:
对于传递函数
G(s)=ansn+⋯+a0bmsm+⋯+b0
- 零点:使分子多项式等于零的 s 值,即 bmsm+⋯+b0=0 的根。在 s 平面上用 ∘ 标记。
- 极点:使分母多项式等于零的 s 值,即 ansn+⋯+a0=0 的根。在 s 平面上用 × 标记。
直观理解:
极点决定了系统”怎么动”——是衰减、振荡还是发散;零点影响”动多少”——它可以放大或抵消某些输入模式。做稳定性分析时,极点是主角。
3.3 特征方程(Characteristic Equation)
一句话理解:
特征方程就是把传递函数的分母令为零。
ansn+an−1sn−1+⋯+a1s+a0=0
这个方程的根就是极点。系统的全部动态特性(稳定性、响应速度、振荡频率)都由这个方程的根决定。
为什么叫”特征”?
因为它反映了系统本身的固有属性——不管你给系统施加什么输入,系统的固有响应模式都由这些根决定。就像一个人的性格不会因为遇到什么人而改变一样。
3.4 s 平面(Complex Plane)
一句话理解:
s 平面就是画极点和零点的坐标纸。横轴是实部 σ,纵轴是虚部 jω。
关键规则:
- 极点在左半平面(σ<0):系统模态是衰减的,比如 e−3t,会越来越小。
- 极点在右半平面(σ>0):系统模态是发散的,比如 e2t,会越来越大。
- 极点在虚轴上(σ=0):系统模态是等幅振荡的,比如 ejωt,不衰减也不增长。
判断口诀:左稳右不稳,虚轴临界。
3.5 LTI 系统的稳定性
一句话理解:
系统稳定 = 输入去掉之后,输出最终会归零。
正式定义:
LTI 系统稳定当且仅当所有极点都在 s 平面的左半平面(即所有极点的实部都为负)。
三种情况:
| 极点位置 | 实部 σ | 系统模态 | 稳定性 |
|---|
| 左半平面 | σ<0 | eσt 衰减 | 稳定 |
| 右半平面 | σ>0 | eσt 发散 | 不稳定 |
| 虚轴上 | σ=0 | 等幅振荡 | 临界稳定 |
对于复共轭极点对 σ±jω:稳定性只看实部 σ,虚部 ω 只影响振荡频率,不影响稳定性。
常见错误:
- 不稳定条件只要一个极点在右半平面就够了,不需要全部都在右边。
- 虚轴上的极点(包括原点处的极点)在工程上通常视为不稳定,因为任何微小扰动都可能让系统漂移到右半平面。
3.6 极点位置练习:质量-弹簧-阻尼系统
给定系统的传递函数:
G(s)=Y(s)X(s)=s2+2ζωns+ωn22ζωns+ωn2
分母 s2+2ζωns+ωn2=0 的根就是极点。
情况 a:过阻尼(ζ=1.25,ωn=4)
s2+2(1.25)(4)s+42=s2+10s+16=0
用求根公式:
s=2−10±100−64=2−10±6
极点:s1=−2,s2=−8
两个极点都是实数且在左半平面——系统稳定,响应无振荡。
情况 b:欠阻尼(ζ=0.4,ωn=4)
s2+2(0.4)(4)s+42=s2+3.2s+16=0
s=2−3.2±10.24−64=2−3.2±−53.76=−1.6±j3.68
极点是复共轭对:s=−1.6±j3.68
实部为负,系统稳定;但有虚部,响应会有衰减振荡。
关键结论: 同一个系统,改变阻尼比 ζ,极点从左半平面的实轴(过阻尼)变成左半平面的复共轭(欠阻尼)。极点位置直接反映系统行为。
核心公式与推导
4.1 传递函数的极零点表达式
将传递函数分解为因式形式:
G(s)=K⋅(s−p1)(s−p2)⋯(s−pn)(s−z1)(s−z2)⋯(s−zm)
其中 z1,…,zm 是零点,p1,…,pn 是极点,K 是增益。
这个形式的好处:一眼看出极点和零点的位置,便于分析。
4.2 Routh 稳定性判据
目的: 不解特征方程,直接从系数判断所有极点是否都在左半平面。
步骤一:写出特征方程并归一化
sn+a1sn−1+a2sn−2+⋯+an−1s+an=0
确保最高次项系数为 1(如果原方程最高次系数不是 1,先两边同除以它)。
步骤二:必要条件检查(快筛)
如果系数 a1,a2,…,an 中任何一个为零或为负,则系统不稳定或临界稳定。
但这只是必要条件,不是充分条件。 系数全正不代表一定稳定。
步骤三:构造 Routh 表
Routh 表的行数 = n+1(从 sn 到 s0),每行元素从上两行计算得出。
前两行直接填:
| 行 | 列1 | 列2 | 列3 | 列4 | … |
|---|
| sn | 1 | a2 | a4 | a6 | … |
| sn−1 | a1 | a3 | a5 | a7 | … |
| sn−2 | b1 | b2 | b3 | … | |
| sn−3 | c1 | c2 | c3 | … | |
| ⋮ | | | | | |
| s0 | | | | | |
注意前两行取系数的规律:sn 行取第 1、3、5、7… 个系数,sn−1 行取第 2、4、6、8… 个系数。这是最容易搞错的地方之一。
从第三行开始,每个元素的计算公式:
b1=a1a1⋅a2−1⋅a3=a1a1a2−a3
b2=a1a1⋅a4−1⋅a5
b3=a1a1⋅a6−1⋅a7
通用公式:设相邻两行为
| 行 | 第1列 | 第2列 | 第3列 | … |
|---|
| 上一行 | α1 | α2 | α3 | … |
| 下一行 | β1 | β2 | β3 | … |
则下下行为:
γi=β1β1⋅αi+1−α1⋅βi+1
记忆口诀: “主元乘右减上行乘右,除以主元”。主元(pivot)就是紧挨着的上一行第 1 列的元素(也就是 β1)。
步骤四:判断稳定性
看 Routh 表第一列所有元素的符号:
- 第一列全部同号(通常全正) ⇒ 系统稳定
- 第一列符号变化的次数 = 右半平面极点的个数
例如第一列为 {1,3,−2,5,1},从 3 到 −2 变了一次,从 −2 到 5 又变了一次,共两次符号变化 → 有两个极点在右半平面 → 系统不稳定。
常见错误:
- Routh 表第一行的元素顺序是”跳着取”的:1,a2,a4,…,不是 1,a1,a2,…。
- 算某个 bi 时,一定要用相邻两行对应列的元素,不要把列搞错。
- 遇到某行全为零的情况,需要构造辅助方程并对其求导后代入(下面例题 2 会演示)。遇到第一列出现零但整行不为零,用 ϵ(趋近于 0 的小量)替代后继续计算。
配套例题
例题 1:从微分方程求传递函数
题目: 已知系统微分方程
dt2d2c(t)+3dtdc(t)+2c(t)=r(t)
求传递函数 G(s)=C(s)/R(s),并求极点。
解:
两边 Laplace 变换(零初始条件):
s2C(s)+3sC(s)+2C(s)=R(s)
C(s)(s2+3s+2)=R(s)
G(s)=R(s)C(s)=s2+3s+21=(s+1)(s+2)1
特征方程:s2+3s+2=0
极点:s1=−1,s2=−2
两个极点都在左半平面 → 系统稳定。
本题关键点: 注意分子是常数 1(输入 r(t) 的 Laplace 变换在方程右边,没有导数项),所以分子多项式没有零点(严格说,m=0,分子是常数)。
例题 2:Routh 表中遇到全零行
题目: 特征方程为
s4+3s3+3s2+3s+2=0
用 Routh 判据判断稳定性。
解:
第一步:确认所有系数为正(1, 3, 3, 3, 2),必要条件满足,继续。
第二步:构造 Routh 表。
取系数:1,3,3,3,2。注意最高次 n=4,所以:
- s4 行:取第 1, 3, 5 个系数 → 1,3,2
- s3 行:取第 2, 4 个系数 → 3,3,0
| 行 | 列1 | 列2 | 列3 |
|---|
| s4 | 1 | 3 | 2 |
| s3 | 3 | 3 | 0 |
s2 行的计算:
b1=33×3−1×3=39−3=2
b2=33×2−1×0=36=2
| 行 | 列1 | 列2 | 列3 |
|---|
| s4 | 1 | 3 | 2 |
| s3 | 3 | 3 | 0 |
| s2 | 2 | 2 | |
s1 行的计算:
c1=22×3−3×2=26−6=0
c2=22×0−3×0=0
这里 s1 行全为零!这说明特征方程有对称于原点的根(可能是纯虚根或实根 ±α)。
遇到全零行的标准处理方法:用上面一行(s2 行)构造辅助方程:
A(s)=2s2+2=0⟹s2=−1⟹s=±j
A(s) 求导:dsdA=4s,系数为 4。
用 4 替代 s1 行第 1 列的 0,继续计算 s0 行:
s0 行=44×2−2×0=2
| 行 | 列1 | 列2 |
|---|
| s1 | 4(替换0) | |
| s0 | 2 | |
第一列:1,3,2,4,2 ——全部为正,没有符号变化。
但是,辅助方程告诉我们 s=±j 是特征方程的根——这是虚轴上的极点。
结论: Routh 表第一列无符号变化,说明右半平面没有极点;但辅助方程解出虚轴上的极点,系统临界稳定。
本题的教训: Routh 表第一列全正只说明”右半平面没有极点”,但出现全零行时要额外检查辅助方程——辅助方程的根是特征方程的根的子集,可能是实根或纯虚根。
例题 3:六阶系统的 Routh 表(Lecture 原题)
题目: 特征方程
a(s)=s6+4s5+3s4+2s3+s2+4s+4=0
构造 Routh 表并判断稳定性。
解:
系数:1,4,3,2,1,4,4(n=6)
填前两行: 第 1 行取奇数位(1, 3, 5, 7),第 2 行取偶数位(2, 4, 6)。
| 行 | 列1 | 列2 | 列3 | 列4 |
|---|
| s6 | 1 | 3 | 1 | 4 |
| s5 | 4 | 2 | 4 | 0 |
计算 s4 行:(主元 = s5 行第 1 列 = 4)
b1=44×3−1×2=412−2=2.5
b2=44×1−1×4=44−4=0
b3=44×4−1×0=416=4
计算 s3 行:(主元 = s4 行第 1 列 = 2.5)
c1=2.52.5×2−4×0=2.55−0=2
c2=2.52.5×4−4×0=2.510−0=4
计算 s2 行:(主元 = s3 行第 1 列 = 2)
d1=22×0−2.5×4=20−10=−5
d2=22×4−2.5×0=28−0=4
计算 s1 行:(主元 = s2 行第 1 列 = −5)
e1=−5−5×4−2×4=−5−20−8=−5−28=5.6
计算 s0 行:(主元 = s1 行第 1 列 = 5.6)
f1=5.65.6×4−(−5)×0=4
完整 Routh 表第一列:
| 行 | 第1列 |
|---|
| s6 | 1 |
| s5 | 4 |
| s4 | 2.5 |
| s3 | 2 |
| s2 | −5 |
| s1 | 5.6 |
| s0 | 4 |
分析符号变化:
- s3 行的 2 → s2 行的 −5:变号(正到负),第 1 次
- s2 行的 −5 → s1 行的 5.6:变号(负到正),第 2 次
共 2 次符号变化 → 有 2 个极点在右半平面 → 系统不稳定。
例题 4:四组极点图的稳定性分析
题目: 给定四组极点在 s 平面上的位置,分别判断系统的稳定性。
| 配置 | 极点位置描述 |
|---|
| (a) | 两个极点都在左半平面 |
| (b) | 两个极点都在虚轴上(关于实轴对称) |
| (c) | 一个极点在原点,一个在左半平面 |
| (d) | 一个极点在左半平面,一个在右半平面 |
分析:
- (a) 所有极点在左半平面 → 稳定。暂态响应随时间衰减至零。
- (b) 极点在虚轴上 → 临界稳定。响应是等幅正弦振荡,既不衰减也不增长。
- (c) 原点处的极点对应常数模态(不衰减),但也不增长 → 临界稳定。工程上通常视为不稳定。
- (d) 右半平面有极点 → 不稳定。哪怕只有一个极点在右半平面,响应就会随时间指数增长。
重点难点总结
极点与零点的区分
| 极点(Pole) | 零点(Zero) |
|---|
| 位置 | 传递函数分母的根 | 传递函数分子的根 |
| 效果 | 使 $ | G(s) |
| 对稳定性的影响 | 直接决定系统稳定性 | 影响响应形状,不影响稳定性 |
| 画法 | × | ∘ |
稳定性判断三步法
- 看极点位置:全部在左半平面 = 稳定。
- 不想解高阶方程时用 Routh 判据。
- Routh 判据先检查系数全正(必要条件),再构造表(充分条件)。
Routh 判据做题清单
- 特征方程最高次系数化为 1。
- 检查所有系数是否存在零或负数,存在则直接判定不稳定。
- 正确填写前两行——sn 行取第 1、3、5… 个系数,sn−1 行取第 2、4、6… 个系数。
- 逐行计算,注意公式的分子是”交叉相乘再相减除以主元”。
- 如果某行全为零,构造辅助方程,求其根并求导后代入。
- 最后统计第一列符号变化次数。
典型坑
- 分母阶数小于分子阶数:物理上不可实现,考试中出现说明需要先检查化简。
- Routh 表某列出现零但整行不为零:用 ϵ(趋近于 0+ 的小量)替代,然后正常计算,最后判断 ϵ→0+ 时第一列的符号。
- 忘记归一化:如果特征方程最高次系数不是 1(如 2s3+4s2+3s+1=0),先两边除以 2 再构造 Routh 表。
- 辅助方程的根要代回原方程验证:辅助方程的根是原特征方程根的子集,但不能直接用辅助方程判断所有极点的稳定性。
自测题与答案
自测题 1:极点计算
求传递函数
G(s)=s3+5s2+8s+4s+3
的极点和零点,并判断稳定性。
答案
零点: 分子 s+3=0⇒s=−3。
极点: 分母 s3+5s2+8s+4=0。尝试 s=−1:
(−1)3+5(−1)2+8(−1)+4=−1+5−8+4=0
所以 s=−1 是一个根。做多项式除法:
s3+5s2+8s+4=(s+1)(s2+4s+4)=(s+1)(s+2)2
极点:s1=−1,s2=s3=−2(二重极点)。
全部极点在左半平面,系统稳定。
自测题 2:Routh 判断(含全零行)
特征方程:
s3+2s2+4s+8=0
用 Routh 判据判断稳定性。
答案
系数:1,2,4,8,全部为正,必要条件满足。
Routh 表:
| 行 | 列1 | 列2 |
|---|
| s3 | 1 | 4 |
| s2 | 2 | 8 |
| s1 | 22×4−1×8=0 | |
| s0 | | |
s1 行为零!用 s2 行构造辅助方程:A(s)=2s2+8=0
s2=−4⟹s=±j2
求导:dsdA=4s,用 4 替代 s1 行第 1 列。
继续:s0 行 = 44×8−2×0=8
第一列:1,2,4,8,全正,无符号变化 → 右半平面无极点。
但辅助方程给出 s=±j2,极点在虚轴上 → 临界稳定。
自测题 3:不稳定系统
特征方程:
s4+s3+3s2+5s+10=0
判断稳定性。
答案
系数:1,1,3,5,10,全部为正。
Routh 表:
| 行 | 列1 | 列2 | 列3 |
|---|
| s4 | 1 | 3 | 10 |
| s3 | 1 | 5 | 0 |
| s2 | 11×3−1×5=−2 | 11×10−1×0=10 | |
| s1 | −2−2×5−1×10=−2−20=10 | | |
| s0 | 1010×10−(−2)×0=10 | | |
第一列:1,1,−2,10,10
符号变化:1→−2(1次),−2→10(1次),共 2 次符号变化 → 2 个极点在右半平面 → 系统不稳定。
注意:这个例子说明”系数全正”只是必要条件,不是充分条件——四个系数全正,系统照样不稳定。
自测题 4:简答
为什么 Routh 判据只能判断”有几个极点在右半平面”,而不能给出极点的具体数值?
答案
Routh 判据利用了多项式系数与根的位置之间的拓扑关系(符号变化次数等于右半平面根的个数),但这个关系不携带根的具体位置信息。要知道极点的确切数值,只能直接求解特征方程。Routh 判据的优势在于”不解方程就能判断稳定性”,代价是丢失了根的精确位置。
学习路线
- 第一步: 彻底搞懂传递函数的定义和求法——从微分方程到 s 域。把例题 1 这种基础题练到闭眼能做。
- 第二步: 练习求极点和零点——二阶系统先做熟(求根公式),再做三阶(试根 + 多项式除法)。
- 第三步: 理解”极点位置 ↔ 系统行为”的对应关系——在 s 平面上画图,标出稳定区和不稳定区。
- 第四步: 掌握 Routh 表的构造方法——从简单(三阶)到复杂(五阶、六阶)。先做没有特殊情况的标准题。
- 第五步: 处理特殊情况——全零行和辅助方程、第一列出现零的 ϵ 替代法。这类题在考试中出现频率不高但一旦出现就是区分度所在。
建议做题顺序:先做二阶系统的极点判断(不需要 Routh),再做三阶系统的 Routh 判据(最常考),最后练带辅助方程的高阶问题。时间紧张的话,三阶 Routh 判据 + 极点位置与稳定性判断必须优先掌握。
和后续章节的关系
- 根轨迹(Root Locus): 根轨迹研究”当系统某个参数变化时,极点在 s 平面上怎么移动”。本章的极点和稳定性概念是根轨迹分析的基础——你得先知道”极点在哪边是稳定的”,才能看懂根轨迹图。
- 频率响应(Bode Plot / Nyquist): 频率响应分析把极点信息转化为”系统对不同频率输入的反应”。稳定性是频率响应分析的前提——不稳定的系统谈频率响应没有工程意义。
- PID 控制器设计: PID 控制器本质上是改变闭环系统的极点位置。设计 PID 时,第一件事就是确保闭环系统的极点全部落在左半平面。如果你不知道什么是”极点在左半平面”,PID 参数整定就是瞎调。
- 状态空间方法: 后续如果学到状态空间,稳定性分析对应的是系统矩阵 A 的特征值——本质上和传递函数的极点是同一件事的不同表述。
一句话:稳定性是控制理论的入场券。 不会判断稳定性,后面所有章节都没法学。