版权归原作者所有,如有侵权,请联系我们

[科普中国]-劳斯–赫尔维茨稳定性判据

科学百科
原创
科学百科为用户提供权威科普内容,打造知识科普阵地
收藏

劳斯–赫尔维茨稳定性判据(英语:Routh–Hurwitz stability criterion)是控制理论中的一个数学测试,是线性时不变系统(LTI)稳定的充份必要条件。劳斯测试是由英国数学家爱德华·劳斯在1876年提出的快速算法,可以判断一线性系统其特征多项式的根是否都有负的实部。德国数学家阿道夫·赫维兹在1895年独立的提出将多项式的系数放到一个方阵中(此方阵称为赫维兹矩阵),证明多项式稳定当且仅当赫维兹矩阵的主要子矩阵其行列式形成的数列均为正值。二个程序是等价的,而劳斯测试提供一个有效计算赫维兹行列式的方法。满足劳斯–赫尔维茨稳定性判据的多项式称为赫尔维茨多项式。

详解此稳定性判据之所以重要,是因为若线性系统之特征方程式的根p均有负的实部,表示其解e为稳定的(BIBO稳定)。1因此稳定性判据提供了方式,可以在不求解线性系统的运动方程的情形下,判断其是否只有稳定解。对于离散系统,对应稳定性的测试可以由Schur–Cohn判据、Jury稳定性判据及Bistritz稳定性判据来判断。随着电脑的进步,此稳定性判据变的较少使用,另一种判断的方式则是用数值方法直接求解多项式,得到其解的近似值。

劳斯测试可以由辗转相除法以及在计算柯西指标时用施图姆定理来推导。赫尔维茨利用另一种方式来推导其稳定性判据。

利用辗转相除法求解劳斯–赫尔维茨稳定性判据和劳斯–赫尔维茨定理有关。由定理的陈述,可得 其中:

1)p为多项式ƒ(z)的根中实部为负值的个数。

2)q为多项式ƒ(z)的根中实部为正值的个数。(此假设ƒ(z)的根都不在虚轴上)

3)w(x)为 由施图姆定理得到的变号数(中间利用连续的辗转相除法),其中 ,y为实数。

根据代数基本定理,每个n次的多项式在复数平面上会有n个根(也就是,对于根都不在虚轴上的ƒ,p+q=n)。因此可得到ƒ为(稳定的)赫尔维茨多项式当且仅当p−q=n。利用劳斯–赫尔维茨定理,可以将p和q的条件改为以广义施图姆链组成的条件,也就是以ƒ的系数组合而成的条件。

使用矩阵令f(z) 为一个复多项式。过程如下:

1)计算多项式和使得

其中y为一实数。

2)计算与相关的西尔维斯特矩阵。

3)重新排列每一行,让奇数行和其下一行起始地零数量相同。

4)计算该矩阵的每个主子式。

5)如果这些子式中至少一个为负(或零),则多项式f是不稳定的。

范例1)令(为了简单起见,我们取实系数),其中

(为了避免一根为零,我们可以利用劳斯–赫尔维茨定理)。首先,我们要计算实多项式

2)接下来,我们将多项式辗转相除来得到广义施图姆链:

得到 ,于是辗转相除法结束。

注意我们在第一次除法中,必须假设b不为零。在此情形下,广义施图姆链为

的符号与a相反,而by的符号与b相同。当我们令 时,该链的第一部分符号也与a相反,而by的符号与b相反。最终,-c的符号总与c相反。

现在假设f是赫尔维茨稳定的。这意味着(f的阶数)。由函数w的性质,这与相同。因此,a,b和c必须符号相同。我们找到了二阶多项式稳定的必要条件。

高阶多项式的劳斯–赫尔维茨判据在下面,我们假设最高阶的系数(例如二阶多项式中的 )为正。如果有必要,总可以将该多项式乘以 得到。

1)对于二阶多项式

如果所有系数都满足:

则所有根都在左半平面(特征方程为 的系统稳定).

2)对于三阶多项式

所有系数都满足

并且有

3)对于四阶多项式

所有系数都满足

并且

4)一般地,劳斯稳定性判据要求劳斯表的第一列所有元素符号相同。

满足上述判据的系统称为闭环稳定系统,否则会由于第一列元素符号改变而不稳定。

本词条内容贡献者为:

刘军 - 副研究员 - 中国科学院工程热物理研究所