根轨迹绘制八大法则推导
点击各节标题可跳转至法则速查表中的对应规则行。
根轨迹本质上画的是闭环特征根随开环增益 K 从 0 变到 ∞ 的变化轨迹。
设开环传递函数为
G(s)H(s)=KA(s)B(s)
其中 A(s)=0 是开环极点方程,B(s)=0 是开环零点方程。闭环特征方程为
1+G(s)H(s)=0
代入后整理得
A(s)+KB(s)=0
这就是根轨迹真正研究的方程。
当 K=0 时:
A(s)+0⋅B(s)=0⟹A(s)=0
闭环特征根刚好等于开环极点,所以根轨迹的起点是开环极点。
当 K→∞ 时:
两边除以 K,得
KA(s)+B(s)=0
因为 K→∞,所以 KA(s)→0,于是 B(s)=0,闭环特征根趋向开环零点,所以根轨迹的终点是开环零点。
但要注意:闭环特征方程的阶数由开环极点数 n 决定,所以根轨迹共有 n 条分支。如果开环零点数为 m 且 n>m,则只有 m 条分支终止于有限零点,剩下的 n−m 条分支没有有限零点可去,只能趋向无穷远处,终止于无限零点。
K=0 时闭环特征方程退化为开环分母,根轨迹从开环极点出发;K→∞ 时方程主要由开环分子决定,根轨迹终止于开环零点;不够的零点看作在无穷远处。
根轨迹画的是闭环特征根随增益 K 变化的轨迹。
设开环传递函数为
G(s)H(s)=K∏i=1n(s−pi)∏j=1m(s−zj)
闭环特征方程整理得
i=1∏n(s−pi)+Kj=1∏m(s−zj)=0
根轨迹分支数等于该方程的根的个数。
当 n>m 时,方程最高次由 ∏i=1n(s−pi) 决定,为 n 次,有 n 个根,因此
根轨迹分支数=n
当 m>n 时,方程最高次由 K∏j=1m(s−zj) 决定,为 m 次,有 m 个根,因此
根轨迹分支数=m
统一记为
根轨迹分支数=max(n,m)
实际控制系统中通常 n≥m,所以更常见的是根轨迹分支数等于开环极点数。
对称性: 根轨迹关于实轴对称,原因是闭环特征方程为实系数方程。若 s=σ+jω 是特征根,则其共轭 s=σ−jω 也必是特征根,二者关于实轴对称。
根轨迹分支数由闭环特征方程的次数决定,等于 max(n,m);由于方程是实系数方程,复根必成共轭对出现,根轨迹关于实轴对称。
这里默认 n>m。根轨迹共 n 条,其中 m 条终止于有限零点,剩下 n−m 条趋向无穷远,对应 n−m 条渐近线。
渐近线交角
根轨迹满足相角条件
j=1∑m∠(s−zj)−i=1∑n∠(s−pi)=(2k+1)π
当根轨迹上的点 s 趋向无穷远时,所有有限极点和零点看去方向几乎相同,每个角度都近似等于 φa,代入相角条件得
mφa−nφa=(2k+1)π
即
−(n−m)φa=(2k+1)π
整理后得
φa=n−m(2k+1)π,k=0,1,⋯,n−m−1
共取 n−m 个不同方向。
渐近线交点
考察 K→∞ 时无穷远分支的中心位置,研究
∏j=1m(s−zj)∏i=1n(s−pi)
当 s 很大时展开,最高次部分近似为
sn−m−(i=1∑npi−j=1∑mzj)sn−m−1+⋯
令 q=n−m,设平移 s=w+σa,为使渐近线围绕某中心点对称展开,需消去 wq−1 项,其系数为
qσa−(i=1∑npi−j=1∑mzj)=0
解得
σa=n−mi=1∑npi−j=1∑mzj
这个点称为渐近线交点,也叫渐近线重心。
n−m 条分支跑向无穷远,渐近线角度来自相角条件在无穷远处的简化,渐近线交点来自令展开多项式次高次项系数为零。
这条规则同样由根轨迹的相角条件推出。
根轨迹满足
j=1∑m∠(s−zj)−i=1∑n∠(s−pi)=(2k+1)π
即从所有开环零点指向试探点 s 的角度之和,减去从所有开环极点指向试探点 s 的角度之和,必须等于奇数倍 180∘。
现在只看实轴上的某一点 s:
- 若某个开环实零点或实极点在 s 的左边,从它指向 s 的方向向右,角度为 0∘
- 若某个开环实零点或实极点在 s 的右边,从它指向 s 的方向向左,角度为 180∘
设 s 右边有 Nz 个开环实零点、Np 个开环实极点,则总相角贡献为
Nz⋅π−Np⋅π=(Nz−Np)π
要满足相角条件,需要
(Nz−Np)π=(2k+1)π
即 Nz−Np 必须是奇数。而 Nz−Np 是奇数等价于 Nz+Np 是奇数,所以结论为:
实轴上某一区域,如果它右方的开环实数零点与实数极点总数为奇数,则该区域是根轨迹。
为什么不用管复数零极点?
实际系统的复数零极点成共轭对出现。对于实轴上的点,一对共轭复零点(或复极点)产生的两个角度大小相等、正负相反,相加为 360∘,对奇偶判断没有影响,因此只需数实轴上该区域右方的开环实零点和实极点。
左边的点贡献 0∘,右边的点贡献 180∘;右边有奇数个零极点时,总相角恰好是奇数倍 180∘,该实轴区域即为根轨迹。
分离点坐标的推导
设开环传递函数为
G(s)H(s)=K∏i=1n(s−pi)∏j=1m(s−zj)
由闭环特征方程 1+G(s)H(s)=0 整理得
K=−∏j=1m(s−zj)∏i=1n(s−pi)
在根轨迹上,增益 K 可以看成 s 的函数 K=K(s)。
分离点是几条根轨迹分支在实轴上相遇再离开的点,此处闭环根重合,相当于 K(s) 的极值点(转折点),因此分离点满足
dsdK=0
分式法推导(对数求导):
令 A(s)=∏i=1n(s−pi),B(s)=∏j=1m(s−zj),则 K=−A(s)/B(s),两边取对数后对 s 求导:
K1dsdK=i=1∑ns−pi1−j=1∑ms−zj1
令 dsdK=0,由于 K=0,得
i=1∑ns−pi1−j=1∑ms−zj1=0
即分离点坐标 d 满足
j=1∑md−zj1=i=1∑nd−pi1
导数法推导(直接求导):
直接对 K=−A(s)/B(s) 用商的求导法则:
dsdK=−[B(s)]2A′(s)B(s)−A(s)B′(s)
令分子为零,得
A(s)B′(s)−A′(s)B(s)=0
两种方法的关系:
两种方法本质相同,都来自 dsdK=0,只是求导方式不同:
- 导数法是直接对商式求导,结论为 A(s)B′(s)−A′(s)B(s)=0,适合 A(s)、B(s) 已展开为多项式的情形;
- 分式法是对对数形式求导,结论为各极点和零点到 d 的倒数之和相等,适合零极点形式直接给出的情形。
同一个 dsdK=0,两种等价的写法。
注意: 两种方法求出的 d 都只是候选分离点,还需验证该点在实轴根轨迹区域内,且对应的 K>0。
分离角的推导
设有 l 条根轨迹分支在分离点 d 相遇。在 d 附近,这 l 条分支从同一点向不同方向离开,局部可理解为 l 重根附近的分裂,即
(s−d)l≈C(某实数)
设离开方向的角度为 θ,令 s−d=rejθ,则
(s−d)l=rlejlθ
为满足根轨迹的相角条件,离开方向的总相角必须是奇数倍 π:
lθ=(2k+1)π
因此
θ=l(2k+1)π,k=0,1,⋯,l−1
最常见的情况是 l=2(两条分支相遇),分离角为 90∘ 和 270∘,两条分支垂直离开实轴。
分离点处 K(s) 取极值,所以用 dsdK=0 求;对数求导得分式法,直接求导得导数法,两者等价。分离角由相角条件在重根处的局部展开决定,l 条分支均匀分开,角度为 (2k+1)π/l。
这条规则同样来自根轨迹相角条件:
j=1∑m∠(s−zj)−i=1∑n∠(s−pi)=(2k+1)π
起始角
若根轨迹从复数极点 pi 出发,离开时的方向角称为起始角 θpi。
当根轨迹上的点 s 非常靠近 pi 时,从其他所有零点 zj 和其他极点 pj(j=i)指向 s 的角度,可以近似看作它们指向 pi 的角度。唯一不能近似的,是从 pi 指向 s 的那个角度——它正是未知的起始角 θpi。
将这些角度代入相角条件:
j=1∑mφzjpi−j=1j=i∑nθpjpi+θpi=(2k+1)π
整理,解出 θpi:
θpi=(2k+1)π+j=1∑mφzjpi−j=1j=i∑nθpjpi
根轨迹刚从某复极点出发时,其他零极点对它的角度几乎固定,只有离开的方向未知,代入相角条件即可求出起始角。
终止角
若根轨迹终止于复数零点 zi,进入时的方向角称为终止角 φzi。
当根轨迹上的点 s 非常靠近 zi 时,从其他零点 zj(j=i)和所有极点 pj 指向 s 的角度,可以近似看作它们指向 zi 的角度。唯一不能近似的,是从 zi 指向 s 的那个角度——它正是未知的终止角 φzi。
将这些角度代入相角条件:
j=1j=i∑mφzjzi+φzi−j=1∑nθpjzi=(2k+1)π
整理,解出 φzi:
φzi=(2k+1)π−j=1j=i∑mφzjzi−j=1∑nθpjzi
根轨迹快要到达某复零点时,其他零极点的角度几乎固定,只有进入的方向未知,代入相角条件即可求出终止角。
符号含义
| 符号 | 含义 |
|---|
| φzjpi | 从零点 zj 指向极点 pi 的角度 |
| θpjpi | 从极点 pj 指向极点 pi 的角度 |
| φzjzi | 从零点 zj 指向零点 zi 的角度 |
| θpjzi | 从极点 pj 指向零点 zi 的角度 |
起始角和终止角本质相同:在目标点附近,其他零极点的角度固定,只有一个未知角,代入相角条件解出即可。
根轨迹上的点本质是闭环特征根。虚轴上的点形式为 s=jω,所以根轨迹与虚轴相交,意味着闭环特征方程存在纯虚根 s=±jω,系统处于临界稳定状态——虚轴正好是稳定与不稳定的分界线。
方法一:直接代入 s=jω
将 s=jω 代入闭环特征方程
1+G(jω)H(jω)=0
这是一个复数方程,复数等于零要求实部和虚部同时为零:
{Re[1+G(jω)H(jω)]=0Im[1+G(jω)H(jω)]=0
方程中通常含有两个未知量 K 和 ω,联立求解即可得到穿越虚轴时的增益 K 及交点位置 s=±jω。
原理: 既然交点在虚轴上,直接令根等于 jω,代入特征方程求解。
方法二:劳斯判据
劳斯判据可判断特征根是否全在左半平面。根轨迹穿越虚轴时,系统恰好处于临界稳定,对应劳斯表中的临界情况。
步骤:
- 写出闭环特征方程 D(s)=0,其系数含 K
- 列劳斯表,找使劳斯表出现全零行的 K 值——此时特征方程存在关于原点对称的根,即 s=±jω
- 用全零行上一行构造辅助方程。例如上一行对应
as2+b=0
将 s=jω 代入:
a(jω)2+b=0⟹−aω2+b=0⟹ω=ab
从而求出虚轴交点 s=±jω。
两种方法的关系
两种方法本质相同,只是出发点不同:
- 方法一 从根的位置出发:根在虚轴上,令 s=jω 代入求解
- 方法二 从稳定性边界出发:根在虚轴上意味着系统临界稳定,用劳斯判据找临界 K,再由辅助方程求 ω
虚轴交点就是闭环特征根变成纯虚根的位置。直接法令 s=jω 代入求解;劳斯法找临界稳定的 K,再用辅助方程求 ω。两者等价。
这个结论分两层:第一层是韦达定理,第二层是根轨迹中的根之和规则。
闭环极点的和与积(韦达定理)
设闭环特征方程写成首一形式:
D(s)=sn+a1sn−1+a2sn−2+⋯+an−1s+an=0
它的 n 个闭环极点为 s1,s2,⋯,sn,因式分解形式为:
D(s)=(s−s1)(s−s2)⋯(s−sn)
展开后:
D(s)=sn−(s1+s2+⋯+sn)sn−1+⋯+(−s1)(−s2)⋯(−sn)
与 D(s)=sn+a1sn−1+⋯+an 逐项比较:
sn−1 项:−(s1+⋯+sn)=a1,所以
i=1∑nsi=−a1
常数项:(−s1)(−s2)⋯(−sn)=an,所以
i=1∏n(−si)=an
这就是韦达定理在闭环特征方程上的直接应用。
根之和规则的推导
设开环传递函数为:
G(s)H(s)=KA(s)B(s)=K∏i=1n(s−pi)∏j=1m(s−zj)
其中 pi 是开环极点,zj 是开环零点。闭环特征方程整理为:
A(s)+KB(s)=0
将 A(s)=∏i=1n(s−pi) 展开,最高两项为:
A(s)=sn−(i=1∑npi)sn−1+⋯
当 m≤n−2(即 n−m≥2)时,B(s) 的次数不超过 n−2,KB(s) 不影响 sn−1 项的系数。因此闭环特征方程 A(s)+KB(s)=0 中 sn−1 的系数始终为
−i=1∑npi
由韦达定理,闭环特征方程的 sn−1 系数 a1 满足 ∑si=−a1,所以
i=1∑nsi=i=1∑npi
也就是说,闭环极点虽然随 K 变化而移动,但它们的总和不变,始终等于开环极点之和。
注意: 若 n−m=1,则 KB(s) 会影响 sn−1 项系数,根之和规则不再成立。
闭环极点的和与积来自韦达定理;根之和规则来自 A(s)+KB(s)=0 中 sn−1 项系数与 K 无关(需 n−m≥2),故闭环极点之和恒等于开环极点之和。
Discussion
Comments
Share questions, corrections, or extra notes about this post.