劳斯阵列的推导

劳斯阵列是劳斯–赫尔维茨稳定性判据中，用来判断系统是否稳定的方式，是透过系统的特征多项式系数所建立的阵列。劳斯阵列和劳斯–赫尔维茨理论是古典控制理论的核心，结合了欧几里得算法和施图姆定理来计算柯西指标。

## 柯西指标

给定系统

$$\begin{align} <pre> f(x) & {} = a_0x^n+a_1x^{n-1}+\cdots+a_n & {} \quad (1) \\ & {} = (x-r_1)(x-r_2)\cdots(x-r_n) & {} \quad (2) \\ </pre> \end{align}$$

假设$f(x) = 0$的根都不在虚轴上，并且令

$N$ = 是$f(x) = 0$的根的实部为负数的个数。

$P$ = 是$f(x) = 0$的根的实部为正数的个数。

因此可得

$$N+P=n \quad (3)$$

将$f(x)$以极座标型式表示，可得

$$f(x) = \rho(x)e^{j\theta(x)} \quad (4)$$

其中

$$\rho(x) = \sqrt{\mathfrak{Re}^2[f(x)]+\mathfrak{Im}^2[f(x)]} \quad (5)$$

且

$$\theta(x) = \tan^{-1}\big(\mathfrak{Im}[f(x)]/\mathfrak{Re}[f(x)]\big) \quad (6)$$

根据(2)会发现

$$\theta(x) = \theta_{r_1}(x)+\theta_{r_2}(x)+\cdots+\theta_{r_n}(x) \quad (7)$$

其中

$$\theta_{r_i}(x) = \angle(x-r_i) \quad (8)$$

若$f(x) = 0$的第i个根的实部为正，则（用y=(REy,IMy)的表示法）

$$\begin{align} \theta_{r_i}(x)\big|_{x=-j\infty} & = \angle(x-r_i)\big|_{x=-j\infty} \\ <pre> & = \angle(0-\mathfrak{Re}[r_i],-\infty-\mathfrak{Im}[r_i]) \\ & = \angle(-|\mathfrak{Re}[r_i]|,-\infty) \\ & = \pi + \lim_{\phi \to \infty}\tan^{-1}\phi=\frac{3\pi}{2} \quad (9)\\ </pre> \end{align}$$

且

$$\theta_{r_i}(x)\big|_{x=j0} = \angle(-|\mathfrak{Re}[r_i]|,0) = \pi - \tan^{-1}0=\pi \quad (10)$$

且

$$\theta_{r_i}(x)\big|_{x=j\infty} = \angle(-|\mathfrak{Re}[r_i]|,\infty) = \pi - \lim_{\phi \to \infty}\tan^{-1}\phi=\frac{\pi}{2} \quad (11)$$

同样地，若$f(x) = 0$的第i个根的实部为负。

$$\begin{align} \theta_{r_i}(x)\big|_{x=-j\infty} & = \angle(x-r_i)\big|_{x=-j\infty} \\ <pre> & = \angle(0- \mathfrak{Re}[r_i],-\infty-\mathfrak{Im}[r_i]) \\ & = \angle(|\mathfrak{Re}[r_i]|,-\infty) \\ & = 0 - \lim_{\phi \to \infty}\tan^{1}\phi=-\frac{\pi}{2} \quad (2)\\ </pre> \end{align}$$

且

$$\theta_{r_i}(x)\big|_{x=j0} = \angle(|\mathfrak{Re}[r_i]|,0) = \tan^{-1}0=0\, \quad (13)$$

且

$$\theta_{r_i}(x)\big|_{x=j\infty} = \angle(|\mathfrak{Re}[r_i]|,\infty) = \lim_{\phi \to \infty}\tan^{-1}\phi=\frac{\pi}{2}\, \quad (14)$$

由(9)至(11)式可知，若$f(x)$的第i个根实部为正，则$\theta_{r_i}(x)\Big|_{x=-j\infty}^{x=j\infty} = -\pi$，由(12)至(14)式可知，若$f(x)$的第i个根实部为负，则$\theta_{r_i}(x)\Big|_{x=-j\infty}^{x=j\infty} = \pi$。因此

$$\theta_{r_i}(x)\Big|_{x=-j\infty}^{x=j\infty} = \angle(x-r_1)\Big|_{x=-j\infty}^{x=j\infty}+\angle(x-r_2)\Big|_{x=-j\infty}^{x=j\infty}+\cdots+\angle(x-r_n)\Big|_{x=-j\infty}^{x=j\infty} = \pi N - \pi P \quad (15)$$

若定义

$$\Delta=\frac{1}{\pi}\theta(x)\Big|_{-j\infty}^{j\infty} \quad (16)$$

则可以得到以下的关系

$$N - P = \Delta \quad (17)$$

结合(3)式及(17)式可得

$N = \frac{n+\Delta}{2}$且$P = \frac{n-\Delta}{2} \quad (18)$

因此，给定$n$次的方程$f(x)$，只需要计算$\Delta$，就可以得到根的实部为负的个数$N$，以及根的实部为正的个数$P$。

<image>

图1 $\tan(\theta)$相对$\theta$的图

配合(6)式及图1，$\tan(\theta)$相对$\theta$的图，将$x$在区间(a,b)之间变化，其中$\theta_a=\theta(x)|_{x=ja}$，而$\theta_b=\theta(x)|_{x=jb}$，都是$\pi$的整数倍，若此变化会使函数$\theta(x)$增加$\pi$，表示在从点a到点b的过程中，$\tan \theta (x) = \mathfrak{Im}[f(x)]/\mathfrak{Re}[f(x)]$从$+\infty$“跳到”$-\infty$的次数比从$-\infty$“跳到”$+\infty$的次数多一次。相反的，此变化会使函数$\theta(x)$减少$\pi$，表示在从点a到点b的过程中，$\tan(\theta)$从$+\infty$“跳到”$-\infty$的次数比从$-\infty$“跳到”$+\infty$的次数少一次。

因此，$\theta(x)\Big|_{-j\infty}^{j\infty}$是$\mathfrak{Im}[f(x)]/\mathfrak{Re}[f(x)]$从$-\infty$跳到$+\infty$的次数，减掉同函数从$+\infty$跳到$-\infty$的次数，两者差的$\pi$倍。假设在$x=\pm j\infty$处，$\tan[\theta(x)]$有定义

<image>

图2 $-\cot(\theta)$相对$\theta$的图

若起始点是在不连续点（$\theta_a=\pi/2 \pm i\pi$, i = 0, 1, 2, ...），则因为公式(17)（$N$和$P$都是整数，因此$\Delta$也是整数），其结束点也会在不连续点。此时可以调整指标函数（正跳跃和负跳跃的差值）的计算方式，将正切函数的X轴移动$\pi/2$，也就是在$\theta$上加$\pi/2$。此时的指标函数在各种$f(x)$的系数组合下都有定义，就是在起始点（及结束点）连续的区间(a,b) = $(+j\infty, -j\infty)$内计算$\tan[\theta]=\mathfrak{Im}[f(x)]/\mathfrak{Re}[f(x)]$，再在起始点连续的区间，计算

$$\tan[\theta'(x)]=\tan[\theta + \pi/2] = -\cot[\theta(x)] = -\mathfrak{Re}[f(x)]/\mathfrak{Im}[f(x)] \quad (19)$$

差值$\Delta$是$x$从正跳跃和负跳跃的差值，若计算从$-j\infty$到$+j\infty$所产生的差值，即为相角正切的柯西指标，其相角为$\theta(x)$或$\theta'(x)$，视$\theta_a$是否是$\pi$的整数倍而定。

## 劳斯准则

为了要推导劳斯准则，会将$f(x)$的奇次方项和偶次方项分开来列：

$$f(x) = a_0x^n + b_0x^{n-1} + a_1x^{n-2} + b_1x^{n-3} + \cdots \quad (20)$$

因此可得到

$$\begin{align} <pre> f(j\omega) & = a_0(j\omega)^n+b_0(j\omega)^{n-1}+a_1(j\omega)^{n-2}+b_1(j\omega)^{n-3}+\cdots & {} \quad (21)\\ & = a_0(j\omega)^n+a_1(j\omega)^{n-2}+a_2(j\omega)^{n-4}+\cdots & {} \quad (22)\\ & + b_0(j\omega)^{n-1}+b_1(j\omega)^{n-3}+b_2(j\omega)^{n-5}+\cdots \\ </pre> \end{align}$$

若$n$为偶数：

$$\begin{align} <pre> f(j\omega) & = (-1)^{n/2}\big[a_0\omega^n-a_1\omega^{n-2}+a_2\omega^{n-4}-\cdots \big] & {} \quad (23)\\ & + j(-1)^{(n/2)-1}\big[b_0\omega^{n-1}-b_1\omega^{n-3}+b_2\omega^{n-5}-\cdots \big] & {} \\ </pre> \end{align}$$

若$n$为奇数：

$$\begin{align} <pre> f(j\omega) & = j(-1)^{(n-1)/2}\big[a_0\omega^n-a_1\omega^{n-2}+a_2\omega^{n-4}-\cdots \big] & {} \quad (24)\\ & + (-1)^{(n-1)/2}\big[b_0\omega^{n-1}-b_1\omega^{n-3}+b_2\omega^{n-5}-\cdots \big] & {}\\ </pre> \end{align}$$

可以看出若$n$为奇数，根据(3)式，$N+P$为奇数。若$N+P$为奇数，$N-P$也是奇数。同样的，若 $n$是偶数，$N-P$也是偶数。(15)式可以看出若$N-P$是偶数，$\theta$是$\pi$的整数倍。因此在$n$为偶数时，$\tan(\theta)$有定义，是n为偶数时使用的正确指标，在而在$n$为奇数时，$\tan(\theta') = \tan(\theta+\pi) = -\cot(\theta)$有定义，也是n为奇数时使用的正确指标。

因此，根据(6)式及(23)式，$n$为偶数时：

$$\Delta = I_{-\infty}^{+\infty}\frac{-\mathfrak{Im}[f(x)]}{\mathfrak{Re}[f(x)]}= I_{-\infty}^{+\infty}\frac{b_0\omega^{n-1}-b_1\omega^{n-3}+\cdots}{a_0\omega^n-a_1\omega^{n-2}+\ldots} \quad (25)$$

因此，根据(19)式及(24)式，$n$为奇数时：

$$\Delta = I_{-\infty}^{+\infty}\frac{\mathfrak{Re}[f(x)]}{\mathfrak{Im}[f(x)]}= I_{-\infty}^{+\infty}\frac{b_0\omega^{n-1}-b_1\omega^{n-3}+\ldots}{a_0\omega^n-a_1\omega^{n-2}+\ldots} \quad (26)$$

因此可以计算相同的柯西指标：

$$\Delta = I_{-\infty}^{+\infty}\frac{b_0\omega^{n-1}-b_1\omega^{n-3}+\ldots}{a_0\omega^n-a_1\omega^{n-2}+\ldots} \quad (27)$$

## 参考资料

- Hurwitz, A., "On the Conditions under which an Equation has only Roots with Negative Real Parts", Rpt. in Selected Papers on Mathematical Trends in Control Theory, Ed. R. T. Ballman et al. New York: Dover 1964
- Routh, E. J., A Treatise on the Stability of a Given State of Motion. London: Macmillan, 1877. Rpt. in Stability of Motion, Ed. A. T. Fuller. London: Taylor & Francis, 1975
- Felix Gantmacher (J.L. Brenner translator) (1959) Applications of the Theory of Matrices, pp 177–80, New York: Interscience.