Routh-Hurwitz Criterion

Trigonometry

Trigonometry Logo

Routh-Hurwitz Criterion

Suppose we are given an nth order homogeneous system of differential equations with constant coefficients:

X(t)=AX(t),X(t)=[x1(t)x2(t)xn(t)],A=[a11a12a1na21a22a2nan1an2ann],

where X(t) is an n-dimensional vector containing the unknown functions, A is a square matrix of size n×n.

A nonlinear autonomous system can be reduced to the linear system by performing a linearization around an equilibrium point. Then without loss of generality, we may assume that the equilibrium point is at the origin. It is always possible to reach by choosing a suitable coordinate system.

The stability or instability of the equilibrium state is determined by the signs of the real parts of the eigenvalues of A. To find the eigenvalues λ, it is necessary to solve the auxiliary equation

det(AλI)=0,

which is reduced to an algebraic equation of the nth degree

a0λn+a1λn1+a2λn2++an1λ+an=0.

The roots of this equation can be easily calculated in the case n=2, and in some cases when n3. In other cases, solving the auxiliary equation can be a difficult problem. Moreover, Norwegian mathematician Niels Henrik Abel (18021829) proved a theorem according to which the general algebraic equation of degree n5 cannot be solved using four basic arithmetical operations, that is there is no formula expressing the roots of the equation through its coefficients in the case n5.

In such a situation, methods allowing to determine whether all roots have negative real parts and establish the stability of the system without solving the auxiliary equation itself, are of great importance. One of these methods is the Routh-Hurwitz criterion, which contains the necessary and sufficient conditions for the stability of the system.

Consider again the auxiliary equation

a0λn+a1λn1+a2λn2++an1λ+an=0,

describing the dynamic system. Note that the necessary condition for the stability is satisfied if all the coefficients ai>0. Therefore, we assume that the coefficient a0>0. We write the so-called Hurwitz matrix. It is composed as follows.

The main diagonal of the matrix contains elements a1, a2,, an. The first column contains numbers with odd indices a1, a3, a5, In each row, the index of each following number (counting from left to right) is 1 less than the index of its predecessor. All other coefficients ai with indices greater than n or less than 0 are replaced by zeros. The result is a matrix of kind

The principal diagonal minors Δi of the Hurwitz matrix are given by the formulas

Δ1=a1,Δ2=|a1a0a3a2|,Δ3=|a1a00a3a2a1a5a4a3|,Δn=|a1a000a3a2a10a5a4a30000an|.

We now formulate the Routh-Hurwitz stability criterion: The roots of the auxiliary equation have negative real parts if and only if all the principal diagonal minors of the Hurwitz matrix are positive provided that a0>0: Δ1>0, Δ2>0,, Δn>0. As Δn=anΔn1, the last inequality can be written as an>0.

For the most common systems of the 2nd, 3rd and 4th order, we obtain the following stability criteria:

For a second order system, the condition of the stability is given by

a0>0,Δ1=a1>0,Δ2=|a1a0a3a2|=a1a2>0

or

a0>0,a1>0,a2>0,

that is, all coefficients of the quadratic characteristic equation must be positive. In other words, for a system of 2nd order, the necessary condition of the stability is also the sufficient one. We emphasize that we consider here the asymptotic stability of the zero solution.

For a 3rd order system, the stability criterion is defined by the inequalities

a0>0,Δ1=a1>0,Δ2=|a1a0a3a2|=a1a2a0a3>0,Δ3=a3>0.

or

a0>0,a1>0,a2>0,a3>0,a1a2a0a3>0.

Similarly, for a 4th order system, we obtain the following set of inequalities:

a0>0,Δ1=a1>0,Δ2=|a1a0a3a2|=a1a2a0a3>0,Δ3=|a1a00a3a2a10a4a3|=a1a2a3a12a4a0a32>0,Δ4=a4>0

or

ai>0(i=0,,4),a1a2a0a3>0,a1a2a3a12a4a0a32>0.

If all the n1 principal minors of the Hurwitz matrix are positive and the nth order minor is zero: Δn=0, the system is at the boundary of stability. As Δn=anΔn1, then there are two options:

  1. The coefficient an=0. This corresponds to the case when one of the roots of the auxiliary equation is zero. The system is on the boundary of the aperiodic stability.
  2. The determinant Δn1=0. In this case, there are two complex conjugate imaginary roots. The system is on the boundary of the oscillatory stability.

The Routh-Hurwitz stability criterion belongs to the family of algebraic criteria. It can be conveniently used to analyze the stability of low order systems. The computational complexity grows significantly with the increase of the order. In such cases, it may be preferable to use other criteria such as the Lienard-Shipart theorem or the Nyquist frequency criterion.