Bootstrap

根轨迹的基本概念与绘制

根轨迹的基本概念与绘制

根轨迹概念

根轨迹:系统的某一参数由 0 → ∞ 0\rightarrow \infty 0变化时,闭环极点 λ \lambda λ(特征方程的根)在S平面相应变化所描绘出来的轨迹。(那么根轨迹里的根,指的就是特征方程的根。)
下面给出一个例题来讲解一些常用的参数定义;

系统结构图如下,分析 λ \lambda λ随开环增益K变化的趋势
在这里插入图片描述
写出系统开环传递函数:
G ( s ) = K s ( 0.5 s + 1 ) . . . . . . ( 1 ) = K ∗ s ( s + 2 ) . . . . . . . . . . . ( 2 ) \begin{aligned} G(s)&=\frac{K}{s(0.5s+1)} ......(1) \\ &=\frac{K^{*}}{s(s+2)}...........(2) \end{aligned} G(s)=s(0.5s+1)K......(1)=s(s+2)K...........(2)
其中, K ∗ = 2 K K^*=2K K=2K
{ 式 ( 1 ) 为 尾 一 标 准 型 , K 定 义 为 开 环 增 益 式 ( 2 ) 为 首 一 标 准 型 , K ∗ 定 义 为 根 轨 迹 增 益 \left\{ \begin{array}{c}式(1)为尾一标准型,K定义为开环增益\\ 式(2)为首一标准型,K^*定义为根轨迹增益\\ \end{array} \right. {(1)K(2)K
写出闭环传递函数为;
Φ ( s ) = C ( s ) R ( s ) = K ∗ s 2 + 2 s + K ∗ \Phi(s)=\frac{C(s)}{R(s)}=\frac{K^*}{s^2+2s+K^*} Φ(s)=R(s)C(s)=s2+2s+KK
由于研究 K K K 0 → ∞ 0\rightarrow \infty 0变化,与 K ∗ K^* K 0 → ∞ 0\rightarrow \infty 0变化相同。因此,我们可以直接对 K ∗ K^* K进行讨论。
那么,求闭环极点,特征方程的根,即解如下方程:
s 2 + 2 S + K ∗ = 0 s^2+2S+K^*=0 s2+2S+K=0
设, λ 1 \lambda_1 λ1 λ 2 \lambda_2 λ2为方程的两个根,解得:
λ 1 = − 1 + 1 − K ∗ λ 2 = − 1 − 1 − K ∗ \begin{array}{c} \lambda_1=-1+\sqrt{1-K^*}\\ \lambda_2=-1-\sqrt{1-K^*}\\ \end{array} λ1=1+1K λ2=11K
从而,当 K ∗ K^* K 0 → ∞ 0\rightarrow \infty 0变化时,每一个 K ∗ K^* K的值,都对应着有两个根。如下表。

Alt
由此,根据 λ 1 \lambda_1 λ1 λ 2 \lambda_2 λ2不同的数值,可以绘制得到两条曲线。如下图。

Alt
那么,如上图所示曲线,即为根轨迹。

利用根轨迹分析系统性能

既然画出了根轨迹,那么我们要清楚为什么画根轨迹,画根轨迹的意义是什么。由于根轨迹的根,指的是闭环传递函数的特征根,而通过根在复平面所在位置,我们可以分析并判断系统是否稳定。因此,通过根轨迹我们可以得知系统稳定的根出现在什么位置,并由根的值来反推传递函数中 K ∗ K^* K的值。
因此,接着上面这个例题,接下来分析系统稳定性的问题。

由特征根数值,我们可以知道系统是否稳定。
由于此例题中,根轨迹都在复域的左半平面上,即没有正实根的情况,因此,随着 K ∗ K^* K的值的变化,系统都将保持稳定。

根轨迹的基本方程

我们写出通用的闭环传递函数方程;
Φ ( s ) = G ( s ) H ( s ) 1 + G ( s ) H ( s ) (1) \Phi (s)=\frac{G(s)H(s)}{1+G(s)H(s)} \tag{1} Φ(s)=1+G(s)H(s)G(s)H(s)(1)
由闭环传递函数,我们可得闭环系统特征方程为:
1 + G ( s ) H ( s ) = 0 1+G(s)H(s)=0 1+G(s)H(s)=0
换个式子表示:
G ( s ) H ( s ) = − 1 (2) G(s)H(s)=-1 \tag{2} G(s)H(s)=1(2)
此时,我们假设有m个开环零点,有n个开环极点。则得以下公式;
G ( s ) H ( s ) = K ∗ ∏ j = 1 m ( s − z j ) ∏ i = 1 n ( s − p i ) (3) G(s)H(s)=K^*\frac{\prod_{j=1}^m{\left( s-z_j \right)}}{\prod_{i=1}^n{\left( s-p_i \right)}} \tag{3} G(s)H(s)=Ki=1n(spi)j=1m(szj)(3)
由,(2),(3)两式联立可得:
K ∗ ∏ j = 1 m ( s − z j ) ∏ i = 1 n ( s − p i ) = − 1 (4) K^*\frac{\prod_{j=1}^m{\left( s-z_j \right)}}{\prod_{i=1}^n{\left( s-p_i \right)}}=-1 \tag{4} Ki=1n(spi)j=1m(szj)=1(4)
因为,对于一个高阶系统来说,要求数值解会相当得麻烦,因此,我们会考虑换一种方法来获得根轨迹。而坐标轴除了直角坐标系以外,还有极坐标系。那么,我们这里可以考虑使用开环传递函数的模长和角度满足式(2)来获得根轨迹。
则,我们可以把特征方程分解为角度方程和模长方程两个方程。

复数向量角度及模长计算

在这里要讲解一个数学推导公式,以获得角度方程和模长方程。
在复数中,我们有一个最基本且最重要的公式叫做欧拉公式,复数的建立基本都来源于此公式,如下:
e i φ = cos ⁡ φ + i sin ⁡ φ (5) e^{i\varphi}=\cos \varphi +i\sin \varphi \tag{5} eiφ=cosφ+isinφ(5)
因此,我们将复数的一个形式依据欧拉公式做如下变换:
z = a + b i = a 2 + b 2 ( a a 2 + b 2 + b a 2 + b 2 i ) = a 2 + b 2 ( c o s θ + i s i n θ ) = a 2 + b 2 e i θ = ρ e i θ \begin{aligned} z&=a+bi\\ &=\sqrt{a^2+b^2}\left(\frac{a}{\sqrt{a^2+b^2}}+\frac{b}{\sqrt{a^2+b^2}}i\right)\\ &=\sqrt{a^2+b^2}\left(cos\theta+isin\theta \right)\\ &=\sqrt{a^2+b^2}e^{i\theta}\\ &=\rho e^{i\theta} \end{aligned} z=a+bi=a2+b2 (a2+b2 a+a2+b2 bi)=a2+b2 (cosθ+isinθ)=a2+b2 eiθ=ρeiθ
由上式,我们推导了将复数转化为了一种特殊的形式。其中, ρ \rho ρ为复数的模长,而 θ \theta θ表示复数向量的相角。
那么,我们重新推导复数的乘法及除法公式。首先假设两个复数如下:
z 1 = ρ 1 e i θ 1 z 2 = ρ 2 e i θ 2 z_1 = \rho _1e^{i\theta_1}\\ z_2 = \rho_2 e^{i\theta_2} z1=ρ1eiθ1z2=ρ2eiθ2
将两复数相乘,可得下式:
z ′ = z 1 z 2 = ρ 1 e i θ 1 ⋅ ρ 2 e i θ 2 = ρ 1 ρ 2 e i ( θ 1 + θ 2 ) = ρ ′ e i θ ′ (6) \begin{aligned} z' =z_1 z_2&= \rho _1e^{i\theta_1}\cdot\rho_2 e^{i\theta_2}\\ &=\rho_1\rho2e^{i(\theta_1+\theta_2)}\\ &=\rho 'e^{i\theta'} \end{aligned}\tag{6} z=z1z2=ρ1eiθ1ρ2eiθ2=ρ1ρ2ei(θ1+θ2)=ρeiθ(6)
那么,从上式,我们可以得到如下公式:
{ θ ′ = θ 1 + θ 2 ρ ′ = ρ 1 ρ 2 (7) \left\{ \begin{aligned} \theta'&=\theta _1+\theta _2\\ \rho'&=\rho _1\rho _2\\ \end{aligned} \right. \tag{7} {θρ=θ1+θ2=ρ1ρ2(7)
由上式,说明两复数相乘后的相角等于两复数相角之和,而两复数相乘后模长等于两复数模的积
同样的,将两复数做除法可得下式;
z ′ = z 1 z 2 = ρ 1 e i θ 1 ρ 2 e i θ 2 = ρ 1 ρ 2 e i ( θ 1 − θ 2 ) = ρ ′ e i θ ′ (8) \begin{aligned} z'=\frac{z_1}{z_2} &=\frac{\rho_1e^{i\theta _1}}{\rho_2 e^{i\theta 2}}\\ &=\frac {\rho_1}{\rho_2}e^{i(\theta_1-\theta _2)}\\ &=\rho'e^{i\theta '} \end{aligned} \tag{8} z=z2z1=ρ2eiθ2ρ1eiθ1=ρ2ρ1ei(θ1θ2)=ρeiθ(8)
于是,同理,可得
{ θ ′ = θ 1 − θ 2 ρ ′ = ρ 1 ρ 2 (7) \left\{ \begin{aligned} \theta'&=\theta _1-\theta _2\\ \rho'&=\frac {\rho_1}{\rho_2}\\ \end{aligned} \right. \tag{7} θρ=θ1θ2=ρ2ρ1(7)
由上式,说明两复数相除后的相角等于两复数相角之差,而复数相除后的模等于两复数模的商。

相角条件及模长条件

在有了上述复数运算基础公式(7)后,我们可以通过特征方程获得到相角条件和模长条件,也就是相角方程和模长方程。

K ∗ ∏ j = 1 m ∣ s − z j ∣ ∏ i = 1 n ∣ s − p i ∣ = 1 (8) K^*\frac{\prod_{j=1}^m{\left| s-z_j \right|}}{\prod_{i=1}^n{\left| s-p_i \right|}}=1 \tag{8} Ki=1nspij=1mszj=1(8)
∑ m j = 1 ∠ ( s − z j ) − ∑ n i = 1 ∠ ( s − p i ) = ( 2 k + 1 ) π (9) \sum_{m}^{j=1}\angle {\left(s-z_j \right)}-\sum_{n}^{i=1}\angle {\left( s-p_i \right)}=(2k+1)\pi \tag{9} mj=1(szj)ni=1(spi)=(2k+1)π(9)
那么由(8)和(9)式,我们就可以绘制出系统的根轨迹了。根轨迹8大法则见下篇。

;