查看“︁劳斯–赫尔维茨稳定性判据”︁的源代码

'''劳斯–赫尔维茨稳定性判据'''（{{lang-en|Routh–Hurwitz stability criterion}}）是[[控制理论]]中的一個數學判据，是[[線性系統|線性]][[时不变系统]]（LTI）[[穩定多項式|穩定]]的充分必要條件。劳斯測試是由英國數學家{{le|愛德華·勞斯|Edward Routh}}在1876年提出的快速演算法，可以判斷一[[線性系統]]其[[特徵方程式]]的[[根 (数学)|根]]是否都有負的實部<ref>{{cite book |last=Routh |first=E. J. |year=1877 |title=A Treatise on the Stability of a Given State of Motion: Particularly Steady Motion |url=https://archive.org/details/treatiseonstabil00routuoft |publisher=Macmillan |isbn= }}</ref>。德國數學家[[阿道夫·赫維茲]]在1895年獨立的提出將多項式的係數放到一個方陣中（此方陣稱為[[赫維茲矩陣]]），證明多項式穩定若且唯若赫維茲矩陣的主要子矩陣其行列式形成的數列均為正值<ref>{{cite journal |last=Hurwitz |first=A. |year=1895 |title=Ueber die Bedingungen, unter welchen eine Gleichung nur Wurzeln mit negativen reellen Theilen besitzt |journal=[[Mathematische Annalen|Math. Ann.]] |volume=46 |issue=2 |pages=273–284 |doi=10.1007/BF01446812 }} (English translation “On the conditions under which an equation has only roots with negative real parts” by H. G. Bergmann in ''Selected Papers on Mathematical Trends in Control Theory'' R. Bellman and R. Kalaba Eds. New York: Dover, 1964 pp. 70–82.)</ref>。二個程序是等價的，而劳斯測試提供一個有效計算赫維茲行列式的方法。滿足劳斯–赫尔维茨稳定性判据的多項式稱為[[赫爾維茨多項式]]。

此稳定性判据之所以重要，是因為若[[線性系統]]之特徵方程式的根''p''均有負的實部，表示其解''e<sup>pt</sup>''為穩定的（[[有界輸入有界輸出穩定性|BIBO穩定]]）。因此稳定性判据提供了方式，可以在不求解[[線性系統]]的[[运动方程]]的情形下，判斷其是否只有穩定解。對於離散系統，對應穩定性的測試可以由Schur–Cohn判据、[[Jury稳定性判据]]及{{le|Bistritz稳定性判据|Bistritz stability criterion}}來判斷。隨著電腦的進步，此稳定性判据變的較少使用，另一種判斷的方式則是用數值方法直接求解多項式，得到其解的近似值。

劳斯測試可以由[[輾轉相除法]]以及在計算{{le|柯西指標|cauchy index}}時用[[施图姆定理]]來[[勞斯陣列的推導|推導]]。赫尔维茨利用另一種方式來推導其稳定性判据<ref name="Gopal">{{cite book
  | last = Gopal
  | first = M.
  | title = Control Systems: Principles and Design, 2nd Ed.
  | publisher = Tata McGraw-Hill Education
  | date = 2002
  | location = 
  | pages = 14
  | url = http://books.google.com/books?id=FZak6CkrVLQC&pg=PA14&dq=%22flyball+governor
  | doi = 
  | id = 
  | isbn = 0070482896
  | access-date = 2015-11-03
  | archive-date = 2019-12-16
  | archive-url = https://web.archive.org/web/20191216000311/https://books.google.com/books?id=FZak6CkrVLQC&pg=PA14&dq=%22flyball+governor
  | dead-url = no
  }}</ref>。

==利用輾轉相除法求解==

劳斯–赫尔维茨稳定性判据和{{le|劳斯–赫尔维茨定理|Routh–Hurwitz theorem}}有關。由定理的陳述，可得<math>p-q=w(+\infty)-w(-\infty)</math>其中：
* ''p''為多項式''ƒ''(''z'')的根中實部為負值的個數。
* ''q''為多項式''ƒ''(''z'')的根中實部為正值的個數。（此假設''ƒ''(''z'')的根都不在虛軸上）
* ''w''(''x'')為由 <math>P_0(y)</math>及<math>P_1(y)</math>由[[施图姆定理]]得到的變號數（中間利用連續的[[輾轉相除法]]），其中<math>f(iy)=P_0(y)+iP_1(y)</math>，''y''為實數。
 
根據[[代数基本定理]]，每個''n''次的多項式在複數平面上會有''n'' 個根（也就是，對於根都不在虛軸上的''ƒ''，''p''&nbsp;+&nbsp;''q''&nbsp;=&nbsp;''n''）。因此可得到''ƒ''為（穩定的）[[赫爾維茨多項式]]若且唯若''p''&nbsp;−&nbsp;''q''&nbsp;=&nbsp;''n''<!--(the [[Routh-Hurwitz stability criterion#Appendix A|proof]] is given below).-->。利用劳斯–赫尔维茨定理，可以將''p''和''q''的條件改為以廣義施图姆鏈組成的條件，也就是以''ƒ''的係數組合而成的條件。

==使用矩阵==

令 ''f''(''z'') 为一个复多项式。过程如下：
# 计算多项式 <math>P_0(y)</math> 和 <math>P_1(y)</math> 使得 <math>f(iy)=P_0(y)+iP_1(y)</math> 其中 ''y'' 为一实数。
# 计算与 <math>P_0(y)</math> 和 <math>P_1(y)</math> 相关的[[西尔维斯特矩阵]]。
# 重新排列每一行，让奇数行和其下一行起始地零数量相同。
# 计算该矩阵的每个[[子式和余子式|主子式]]。
# 如果这些子式中至少一个为负（或零），则多项式 ''f'' 是不稳定的。

===范例===

* 令 <math>f(z)=az^2+bz+c</math> (为了简单起见，此處取实系数），其中 <math> c\neq 0</math>（为了避免一根为零，可以利用劳斯–赫尔维茨定理）。首先，要计算实多项式 <math>P_0(y)</math> 与 <math>P_1(y)</math>：
:: <math>f(iy)=-ay^2+iby+c=P_0(y)+iP_1(y)=-ay^2+c+i(by).</math>
: 接下来，将多项式辗转相除来得到广义施图姆链：
:* <math>P_0(y)=((-a/b)y)P_1(y)+c,</math> 得到 <math> P_2(y)=-c,</math>
:* <math>P_1(y)=((-b/c)y)P_2(y),</math> 得到 <math>P_3(y)=0</math>，于是[[輾轉相除法]]结束。
注意在第一次除法中，必须假设 ''b'' 不为零。在此情形下，广义施图姆链为 <math>(P_0(y),P_1(y),P_2(y))=(c-ay^2,by,-c)</math>。令 <math>y=+\infty</math>，<math>c-ay^2</math> 的符号与 ''a'' 相反，而 ''by'' 的符号与 ''b'' 相同。当令 <math>y=-\infty</math> 时，该链的第一部分符号也与 ''a'' 相反，而 ''by'' 的符号与 ''b'' 相反。最终，-''c'' 的符号总与 ''c'' 相反。

现在假设 ''f'' 是赫尔维茨稳定的。这意味着 <math>w(+\infty)-w(-\infty)=2</math>（''f'' 的阶数）。由函数 ''w'' 的性质，这与 <math>w(+\infty)=2</math> 及 <math>w(-\infty)=0</math> 相同。因此，''a'', ''b'' 和 ''c'' 必须符号相同。因此找到了二阶多项式稳定的必要条件.

===二阶、三阶、四阶多项式的劳斯–赫尔维茨判据===

在下面，假设最高阶的系数（例如二阶多项式中的 <math>a_2</math>）为正。如果有必要，总可以将该多项式乘以 <math>-1</math> 得到。

* 对于二阶多项式 <math> P(s) = a_2s^2 + a_1s + a_0 = 0 </math>，如果所有系数都满足  a_1，2.....n 符号相同，则所有根都在左半平面（特征方程为 <math> P(s) </math> 的系统稳定）.
* 对于三阶多项式 <math> P(s) = a_3s^3 + a_2s^2 + a_1s + a_0 = 0</math>，所有系数都满足 <math> a_n > 0 </math>，并且 <math> a_2a_1 > a_3a_0 </math>
* 对于四阶多项式 <math> P(s) = a_4s^4 + a_3s^3 + a_2s^2 + a_1s + a_0 = 0 </math>，所有系数都满足 <math> a_n > 0 </math>，并且 <math> a_3a_2 > a_4a_1</math> 且 <math> a_3a_2a_1 > a_4a_1^2 + a_3^2a_0 </math>
* 一般地，劳斯稳定性判据要求劳斯表的第一列所有元素符号相同。
满足上述判据的系统称为闭环稳定系统，否则会由于第一列元素符号改变而不稳定。

===高阶的例子===

当难以求得高阶特征多项式的根的时候，可以用表格的方法来确定稳定性。对一个 ''n'' 阶多项式
* <math>D(s)=a_ns^n+a_{n-1}s^{n-1}+\cdots+a_1s+a_0</math>
该表格由 ''n''&nbsp;+&nbsp;1 行，结构如下：
{| class="wikitable" width="200px"
 | <math>a_n</math> || <math>a_{n-2}</math> || <math>a_{n-4}</math> || <math>\dots</math>
 |-
 | <math>a_{n-1}</math> || <math>a_{n-3}</math> || <math>a_{n-5}</math> || <math>\dots</math>
 |-
 | <math>b_1</math> || <math>b_2</math> || <math>b_3</math> || <math>\dots</math>
 |-
 | <math>c_1</math> || <math>c_2</math> || <math>c_3</math> || <math>\dots</math>
 |-
 | <math>\vdots</math> || <math>\vdots</math> || <math>\vdots</math> || <math>\ddots</math>
 |}
其中元素 <math>b_i</math> 和 <math>c_i</math> 可以计算如下：
* <math>b_i=\frac{a_{n-1}\times{a_{n-2i}}-a_n\times{a_{n-2i-1}}}{a_{n-1}}.</math>
* <math>c_i=\frac{b_1\times{a_{n-2i-1}}-a_{n-1}\times{b_{i+1}}}{b_1}.</math>
算完之后，第一列中的符号数的变化将是非负极点的数目。
{| class="wikitable"
 |-
 | width="40px" |0.75
 | width="40px" |1.5
 | width="40px" |0
 | width="40px" |0
 |-
 |<nowiki>-3</nowiki>||6||0||0
 |-
 |3||0||0||0
 |-
 |6||0||0||0
|}
在第一列中，有2个符号的变化（0.75&nbsp;→&nbsp;−3，以及 −3&nbsp;→&nbsp;3），因此，有2个非负的根，系统是不稳定的。

有时虚轴上的极点会造成临界稳定的情形。在那种情形中，“劳斯表”的系数一整行都会变为零，因而不能进一步求解出符号的改变了。然后另一种方法可以发挥作用。含有零的这一行的上面一行叫做“辅助多项式”。

* <math>s^6+2s^5+8s^4+12s^3+20s^2+16s+16=0.\,</math>
可得以下的表格：
{| class="wikitable"
 |-
 | width="40px"|1|| width="40px"|8|| width="40px"|20|| width="40px"|16
 |-
 | 2    || 12  || 16 || 0
 |-
 | 2    || 12  || 16 || 0
 |-
 | 0    || 0   ||  0 || 0
 |}
在这个例子中，辅助多项式为 <math>A(s)=2s^4+12s^2+16.\,</math> 仍为零。下一步是对上面的方程求导，得到下面的多项式。<math>B(s)=8s^3+24s^1</math>。包含零的行现在变为 "8" 和 "24"。使用这些值继续建立劳斯表，就会得出虚轴上的两个点。这两个虚轴上的点是边缘稳定性的主要原因。<ref>{{cite book|last=Saeed|first=Syed Hasan|title=Automatic Control Systems|year=2008|publisher=Katson Publishers|location=Delhi|isbn=978-81-906919-2-5|pages=206, 207}}</ref>

==参见==
* [[控制工程]]
* [[勞斯陣列的推導]]
* [[奈奎斯特稳定判据]]
* {{le|劳斯–赫尔维茨定理|Routh–Hurwitz theorem}}
* [[根軌跡圖]]
* [[传递函数]]
* [[Jury稳定性判据]]
* {{le|Bistritz稳定性判据|Bistritz stability criterion}}
* [[哈利托诺夫定理]]
* [[林纳德–奇帕特判据]]

==参考文献==

{{reflist}}
*Felix Gantmacher (J.L. Brenner translator) (1959) ''Applications of the Theory of Matrices'', pp 177–80, New York: Interscience.
* {{cite journal
 |author      = Pippard, A. B.
 |author2     = Dicke, R. H.
 |year        = 1986
 |title       = Response and Stability, An Introduction to the Physical Theory
 |journal     = American Journal of Physics
 |volume      = 54
 |issue       = 11
 |pages       = 1052
 |url         = http://link.aip.org/link/?AJPIAS/54/1052/1
 |accessdate  = 2008-05-07
 |doi         = 10.1119/1.14826
 |bibcode     = 1986AmJPh..54.1052P
 |deadurl     = yes
 |archiveurl  = http://arquivo.pt/wayback/20160514133902/http://link.aip.org/link/?AJPIAS/54/1052/1
 |archivedate = 2016-05-14
}}
* {{cite book
 | author = Richard C. Dorf, Robert H. Bishop
 | year = 2001
 | title =Modern Control Systems | url = https://archive.org/details/moderncontrolsys0000dorf
 |edition=9th
 | publisher = Prentice Hall
 | isbn = 0-13-030660-6
}}
* {{cite book | last1=Rahman | first1=Q. I. | last2=Schmeisser | first2=G. | title=Analytic theory of polynomials | series=London Mathematical Society Monographs. New Series | volume=26 | location=Oxford | publisher=[[牛津大學出版社|Oxford University Press]] | year=2002 | isbn=0-19-853493-0 | zbl=1072.30006 }}
* {{cite journal
  |author = Weisstein, Eric W
  |title = Routh-Hurwitz Theorem.
  |url = http://mathworld.wolfram.com/Routh-HurwitzTheorem.html
  |journal = MathWorld--A Wolfram Web Resource
  |access-date = 2015-11-03
  |archive-date = 2022-07-20
  |archive-url = https://web.archive.org/web/20220720132429/https://mathworld.wolfram.com/Routh-HurwitzTheorem.html
  |dead-url = no
  }}

==外部链接==
* [http://www.mathworks.com/matlabcentral/fileexchange/25956-routh-hurwitz-stability-test A MATLAB script implementing the Routh-Hurwitz test] {{Wayback|url=http://www.mathworks.com/matlabcentral/fileexchange/25956-routh-hurwitz-stability-test |date=20210616213012 }}

{{DEFAULTSORT:劳斯–赫尔维茨稳定性判据}}
[[Category:稳定性理论]]
[[Category:电子反馈]]
[[Category:放大器電路]]
[[Category:信号处理]]
[[Category:多項式]]