Routh–Hurwitz stability criterion

From Wikipedia, the free encyclopedia

In control system theory, the Routh–Hurwitz stability criterion is a mathematical test that is a necessary and sufficient condition for the stability of a linear time invariant (LTI) control system. The Routh test is an efficient recursive algorithm that English mathematician Edward John Routh proposed in 1876 to determine whether all the roots of the characteristic polynomial of a linear system have negative real parts. German mathematician Adolf Hurwitz independently proposed in 1895 to arrange the coefficients of the polynomial into a square matrix, called the Hurwitz matrix, and showed that the polynomial is stable if and only if the sequence of determinants of its principal submatrices are all positive. The two procedures are equivalent, with the Routh test providing a more efficient way to compute the Hurwitz determinants than computing them directly. A polynomial satisfying the Routh-Hurwitz criterion is called a Hurwitz polynomial.

The importance of the criterion is that the roots p of the characteristic equation of a linear system with negative real parts represent solutions e^pt of the system that are stable (bounded). Thus the criterion provides a way to determine if the equations of motion of a linear system have only stable solutions, without solving the system directly. For discrete systems, the corresponding stability test can be handled by the Schur-Cohn criterion, the Jury test and the Bistritz test.

The Routh test can be derived through the use of the Euclidean algorithm and Sturm's theorem in evaluating Cauchy indices. Hurwitz derived his conditions differently.^[1]

Using Euclid's algorithm

The criterion is related to Routh–Hurwitz theorem. Indeed, from the statement of that theorem, we have $p-q=w(+\infty )-w(-\infty )$ where:

p is the number of roots of the polynomial ƒ(z) with negative Real Part;
q is the number of roots of the polynomial ƒ(z) with positive Real Part (let us remind ourselves that ƒ is supposed to have no roots lying on the imaginary line);
w(x) is the number of variations of the generalized Sturm chain obtained from $P_{0}(y)$ and $P_{1}(y)$ (by successive Euclidean divisions) where $f(iy)=P_{0}(y)+iP_{1}(y)$ for a real y.

By the fundamental theorem of algebra, each polynomial of degree n must have n roots in the complex plane (i.e., for an ƒ with no roots on the imaginary line, p + q = n). Thus, we have the condition that ƒ is a (Hurwitz) stable polynomial if and only if p − q = n (the proof is given below). Using the Routh–Hurwitz theorem, we can replace the condition on p and q by a condition on the generalized Sturm chain, which will give in turn a condition on the coefficients of ƒ.

Using matrices

Let f(z) be a complex polynomial. The process is as follows:

Compute the polynomials $P_{0}(y)$ and $P_{1}(y)$ such that $f(iy)=P_{0}(y)+iP_{1}(y)$ where y is a real number.
Compute the Sylvester matrix associated to $P_{0}(y)$ and $P_{1}(y)$ .
Rearrange each row in such a way that an odd row and the following one have the same number of leading zeros.
Compute each principal minor of that matrix.
If at least one of the minors is negative (or zero), then the polynomial f is not stable.

Example

Let $f(z)=az^{2}+bz+c$ (for the sake of simplicity we take real coefficients) where $c\neq 0$ (to avoid a root in zero so that we can use the Routh–Hurwitz theorem). First, we have to calculate the real polynomials $P_{0}(y)$ and $P_{1}(y)$ :

$f(iy)=-ay^{2}+iby+c=P_{0}(y)+iP_{1}(y)=-ay^{2}+c+i(by).$

Next, we divide those polynomials to obtain the generalized Sturm chain:

$P_{0}(y)=((-a/b)y)P_{1}(y)+c,$ yields $P_{2}(y)=-c,$
$P_{1}(y)=((-b/c)y)P_{2}(y),$ yields $P_{3}(y)=0$ and the Euclidean division stops.

Notice that we had to suppose b different from zero in the first division. The generalized Sturm chain is in this case $(P_{0}(y),P_{1}(y),P_{2}(y))=(c-ay^{2},by,-c)$ . Putting $y=+\infty$ , the sign of $c-ay^{2}$ is the opposite sign of a and the sign of by is the sign of b. When we put $y=-\infty$ , the sign of the first element of the chain is again the opposite sign of a and the sign of by is the opposite sign of b. Finally, -c has always the opposite sign of c.

Suppose now that f is Hurwitz-stable. This means that $w(+\infty )-w(-\infty )=2$ (the degree of f). By the properties of the function w, this is the same as $w(+\infty )=2$ and $w(-\infty )=0$ . Thus, a, b and c must have the same sign. We have thus found the necessary condition of stability for polynomials of degree 2.

Routh–Hurwitz criterion for second, third, and fourth-order polynomials

In the following, we assume the coefficient of the highest order (e.g. $a_{2}$ in a second order polynomial) to be positive. If necessary, this can always be achieved by multiplication of the polynomial with $-1$ .

For a second-order polynomial, $P(s)=a_{2}s^{2}+a_{1}s+a_{0}=0$ , all the roots are in the left half plane (and the system with characteristic equation $P(s)$ is stable) if all the coefficients satisfy $a_{n}>0$ .
For a third-order polynomial $P(s)=a_{3}s^{3}+a_{2}s^{2}+a_{1}s+a_{0}=0$ , all the coefficients must satisfy $a_{n}>0$ , and $a_{2}a_{1}>a_{3}a_{0}$
For a fourth-order polynomial $P(s)=a_{4}s^{4}+a_{3}s^{3}+a_{2}s^{2}+a_{1}s+a_{0}=0$ , all the coefficients must satisfy $a_{n}>0$ , and $a_{3}a_{2}>a_{4}a_{1}$ and $a_{3}a_{2}a_{1}>a_{4}a_{1}^{2}+a_{3}^{2}a_{0}$

Higher-order example

A tabular method can be used to determine the stability when the roots of a higher order characteristic polynomial are difficult to obtain. For an nth-degree polynomial

$D(s)=a_{n}s^{n}+a_{{n-1}}s^{{n-1}}+\cdots +a_{1}s+a_{0}$

the table has n + 1 rows and the following structure:

$a_{n}$	$a_{{n-2}}$	$a_{{n-4}}$	$\dots$
$a_{{n-1}}$	$a_{{n-3}}$	$a_{{n-5}}$	$\dots$
$b_{1}$	$b_{2}$	$b_{3}$	$\dots$
$c_{1}$	$c_{2}$	$c_{3}$	$\dots$
$\vdots$	$\vdots$	$\vdots$	$\ddots$

where the elements $b_{i}$ and $c_{i}$ can be computed as follows:

$b_{i}={\frac {a_{{n-1}}\times {a_{{n-2i}}}-a_{n}\times {a_{{n-2i-1}}}}{a_{{n-1}}}}.$
$c_{i}={\frac {b_{1}\times {a_{{n-2i-1}}}-a_{{n-1}}\times {b_{{i+1}}}}{b_{1}}}.$

When completed, the number of sign changes in the first column will be the number of non-negative poles.

Consider a system with a characteristic polynomial

$D(s)=s^{5}+4s^{4}+2s^{3}+5s^{2}+3s+6.\,$

We have the following table:

1	2	3	0
4	5	6	0
0.75	1.5	0	0
−3	6	0
3	0
6	0

In the first column, there are two sign changes (0.75 → −3, and −3 → 3), thus there are two non-negative roots where the system is unstable. " Sometimes the presence of poles on the imaginary axis creates a situation of marginal stability. In that case the coefficients of the "Routh Array" become zero and thus further solution of the polynomial for finding changes in sign is not possible. Then another approach comes into play. The row of polynomial which is just above the row containing the zeroes is called "Auxiliary Polynomial".

$s^{6}+2s^{5}+8s^{4}+12s^{3}+20s^{2}+16s+16=0.\,$

We have the following table:

1	8	20	16
2	12	16	0
2	12	16	0
0	0	0	0

In such a case the Auxiliary polynomial is $A(s)=2s^{4}+12s^{2}+16.\,$ which is again equal to zero. The next step is to differentiate the above equation which yields the following polynomial. $B(s)=8s^{3}+24s^{1}.\,$ . The coefficients of the row containing zero now become "8" and "24". The process of Routh array is proceeded using these values which yield two points on the imaginary axis. These two points on the imaginary axis are the prime cause of marginal stability.^[2]

References

↑ Gopal, M. (2002). Control Systems: Principles and Design, 2nd Ed.. Tata McGraw-Hill Education. p. 14. ISBN 0070482896.
↑ Saeed, Syed Hasan (2008). Automatic Control Systems. Delhi: Katson Publishers. pp. 206, 207. ISBN 978-81-906919-2-5.

Hurwitz, A. (1895). "Über die Bedingungen unter welchen eine Gleichung nur Wurzeln mit negativen reellen Teilen besitzt". Math. Ann. 46, 273-284 (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.).

Routh, E. J. (1877). A Treatise on the Stability of a Given State of Motion: Particularly Steady Motion. Macmillan and co.

Gantmacher, F. R. (1959). "Applications of the Theory of Matrices". Interscience, New York 641 (9): 1–8.

Pippard, A. B.; Dicke, R. H. (1986). "Response and Stability, An Introduction to the Physical Theory". American Journal of Physics 54: 1052. doi:10.1119/1.14826. Retrieved 2008-05-07. Cite uses deprecated parameters (help)
Richard C. Dorf, Robert H. Bishop (2001). Modern Control Systems, 9th Edition. Prentice Hall. ISBN 0-13-030660-6.
Rahman, Q. I.; Schmeisser, G. (2002). Analytic theory of polynomials. London Mathematical Society Monographs. New Series 26. Oxford: Oxford University Press. ISBN 0-19-853493-0. Zbl 1072.30006.

Weisstein, Eric W. "Routh-Hurwitz Theorem.". MathWorld--A Wolfram Web Resource.

External links

A MATLAB script implementing the Routh-Hurwitz test

This article is issued from Wikipedia. The text is available under the Creative Commons Attribution/Share Alike; additional terms may apply for the media files.