Hurwitz matrix explained

In mathematics, a Hurwitz matrix, or Routh–Hurwitz matrix, in engineering stability matrix, is a structured real square matrix constructed with coefficients of a real polynomial.

Hurwitz matrix and the Hurwitz stability criterion

Namely, given a real polynomial

p(z)=a0

n+a
z
1

zn-1+ … +an-1z+an

the

n x n

square matrix

H= \begin{pmatrix} a1&a3&a5&...&...&...&0&0&0\\ a0&a2&a4&&&&\vdots&\vdots&\vdots\\ 0&a1&a3&&&&\vdots&\vdots&\vdots\\ \vdots&a0&a2&\ddots&&&0&\vdots&\vdots\\ \vdots&0&a1&&\ddots&&an&\vdots&\vdots\\ \vdots&\vdots&a0&&&\ddots&an-1&0&\vdots\\ \vdots&\vdots&0&&&&an-2&an&\vdots\\ \vdots&\vdots&\vdots&&&&an-3&an-1&0\\ 0&0&0&...&...&...&an-4&an-2&an \end{pmatrix}.

is called Hurwitz matrix corresponding to the polynomial

p

. It was established by Adolf Hurwitz in 1895 that a real polynomial with

a0>0

is stable(that is, all its roots have strictly negative real part) if and only if all the leading principal minors of the matrix

H(p)

are positive:

\begin{align} \Delta1(p)&=\begin{vmatrix}a1\end{vmatrix}&&=a1>0\\[2mm] \Delta2(p)&=\begin{vmatrix} a1&a3\\ a0&a2\\ \end{vmatrix}&&=a2a1-a0a3>0\\[2mm] \Delta3(p)&=\begin{vmatrix} a1&a3&a5\\ a0&a2&a4\\ 0&a1&a3\\ \end{vmatrix}&&=a3\Delta2-a1(a1a4-a0a5)>0 \end{align}

and so on. The minors

\Deltak(p)

are called the Hurwitz determinants. Similarly, if

a0<0

then the polynomial is stable if and only if the principal minors have alternating signs starting with a negative one.

Hurwitz stable matrices

A

is called a Hurwitz matrix if every eigenvalue of

A

has strictly negative real part, that is,

\operatorname{Re}[λi]<0

for each eigenvalue

λi

.

A

is also called a stable matrix, because then the differential equation
x

=Ax

is asymptotically stable, that is,

x(t)\to0

as

t\toinfty.

If

G(s)

is a (matrix-valued) transfer function, then

G

is called Hurwitz if the poles of all elements of

G

have negative real part. Note that it is not necessary that

G(s),

for a specific argument

s,

be a Hurwitz matrix — it need not even be square. The connection is that if

A

is a Hurwitz matrix, then the dynamical system
x(t)=A

x(t)+Bu(t)

y(t)=Cx(t)+Du(t)

has a Hurwitz transfer function.

Any hyperbolic fixed point (or equilibrium point) of a continuous dynamical system is locally asymptotically stable if and only if the Jacobian of the dynamical system is Hurwitz stable at the fixed point.

The Hurwitz stability matrix is a crucial part of control theory. A system is stable if its control matrix is a Hurwitz matrix. The negative real components of the eigenvalues of the matrix represent negative feedback. Similarly, a system is inherently unstable if any of the eigenvalues have positive real components, representing positive feedback.

See also

References