I have been aware of the Routh-Hurwitz criterion for stability for a long time and I have applied it in three dimensions in my research and tried to apply it in four. Unfortunately I never felt that I really understood it completely. Here I want to finally clear this up. A source which I found more helpful than other things I have seen is https://www.math24.net/routh-hurwitz-criterion/. One problem I have had is that the Hurwitz matrices, which play a central role in this business, are often written in a form with lots of … and I was never sure that I completely understood the definition. I prefer to have a definite algorithm for constructing these matrices. The background is that we would like to understand the stability of steady states of a system of ODE. Suppose we have a system and a steady state , i.e. a solution of . It is well-known that this steady state is asymptotically stable if all eigenvalues of the linearization have negative real parts. This property of the eigenvalues is of course a property of the roots of the characteristic equation . It is always the case here that but I prefer to deal with a general polynomial with real coefficients and a criterion for the situation where all its roots have negative real parts. It is tempting to number the coefficients in the opposite direction, so that, for instance, becomes but I will stick to this convention. Note that it is permissible to replace by in any criterion of this type since if we multiply the polynomial by we get a polynomial in where the order of the coefficients has been reversed. Moreover, if the real part of is non-zero then it has the same sign as the real part of . I find it important to point this out since different authors use different conventions for this. It is convenient to formally extend the definition of the to the integers so that these coefficients are zero for and .

For a fixed value of the Hurwitz matrix is an by matrix defined as follows. The th diagonal element is , with . Starting from a diagonal element and proceeding to the left along a row the index increases by one in each step. Similarly, proceeding to the right along a row the index decreases by one. In the ranges where the index is negative or greater than the element can be replaced by zero. The leading principal minors of the Hurwitz matrix, in other words the determinants of the submatrices which are the upper left hand corner of the original matrix, are the Hurwitz determinants . The Hurwitz criterion says that the real parts of all roots of the polynomial are negative if and only if and for all . Note that a necessary condition for all roots to have negative real parts is that all are positive. Now and so the last condition can be replaced by . Note that the form of the does not depend on . For we get the conditions , and . For we get the conditions , , and . Note that the third condition is invariant under the replacement of by . When , and then the conditions and are equivalent to each other. In this way the invariance under reversal of the order of the coefficients becomes manifest. For we get the conditions , , , and .

Next we look at the issue of loss of stability. If is the region in matrix space where the Routh-Hurwitz criteria are satisfied, what happens on the boundary of ? One possibility is that at least one eigenvalue becomes zero. This is equivalent to the condition . Let us look at the situation where the boundary is approached while remains positive, in other words the determinant of the matrix remains non-zero. Now and so one of the quantities with must become zero. In terms of eigenvalues what happens is that a number of complex conjugate pairs reach the imaginary axis away from zero. The generic case is where it is just one pair. An interesting question is whether and how this kind of event can be detected using the alone. The condition for exactly one pair of roots to reach the imaginary axis is that while the remain positive for . In a paper of Liu (J. Math. Anal. Appl. 182, 250) it is shown that the condition for a Hopf bifurcation that the derivative of the real part of the eigenvalues with respect to a parameter is non-zero is equivalent to the condition that the derivative of with respect to the parameter is non-zero. In a paper with Juliette Hell (Math. Biosci. 282, 162), not knowing the paper of Liu, we proved a result of this kind in the case .

## Leave a Reply