ConceptComplete

Stability of Linear Systems

For linear systems with constant coefficients, stability is completely determined by the eigenvalues of the coefficient matrix. This provides a complete and computable stability theory.


Eigenvalue Criterion

Theorem7.3Stability criterion for linear systems

For the system x=Ax\mathbf{x}' = A\mathbf{x} with constant matrix AA:

  1. The origin is asymptotically stable if and only if all eigenvalues have Re(λj)<0\text{Re}(\lambda_j) < 0.
  2. The origin is stable (not asymptotically) if and only if all eigenvalues satisfy Re(λj)0\text{Re}(\lambda_j) \leq 0 and those with Re(λj)=0\text{Re}(\lambda_j) = 0 are semisimple (geometric multiplicity equals algebraic multiplicity).
  3. The origin is unstable if any eigenvalue has Re(λj)>0\text{Re}(\lambda_j) > 0, or if an eigenvalue with Re(λj)=0\text{Re}(\lambda_j) = 0 is not semisimple.
ExampleStability analysis

A=(1102)A = \begin{pmatrix}-1 & 1\\0 & -2\end{pmatrix}: eigenvalues 1,2-1, -2, both negative. Asymptotically stable.

A=(0110)A = \begin{pmatrix}0 & 1\\-1 & 0\end{pmatrix}: eigenvalues ±i\pm i, pure imaginary, semisimple. Stable (center), not asymptotically stable.

A=(0100)A = \begin{pmatrix}0 & 1\\0 & 0\end{pmatrix}: eigenvalue 00 with algebraic multiplicity 22, geometric multiplicity 11. Unstable (x(t)=x0+y0tx(t) = x_0 + y_0 t grows).


The Routh-Hurwitz Criterion

Definition7.3Hurwitz matrix

For a polynomial p(λ)=λn+an1λn1++a0p(\lambda) = \lambda^n + a_{n-1}\lambda^{n-1} + \cdots + a_0, the Hurwitz matrix is constructed from the coefficients. All roots have negative real parts if and only if all leading principal minors of the Hurwitz matrix are positive.

Theorem7.4Routh-Hurwitz criterion

For a polynomial p(λ)=λn+an1λn1++a0p(\lambda) = \lambda^n + a_{n-1}\lambda^{n-1} + \cdots + a_0, all roots satisfy Re(λ)<0\text{Re}(\lambda) < 0 if and only if:

  • n=1n = 1: a0>0a_0 > 0.
  • n=2n = 2: a1>0a_1 > 0 and a0>0a_0 > 0.
  • n=3n = 3: a2>0a_2 > 0, a0>0a_0 > 0, and a2a1>a0a_2 a_1 > a_0.
ExampleApplying Routh-Hurwitz

Is λ3+3λ2+4λ+2=0\lambda^3 + 3\lambda^2 + 4\lambda + 2 = 0 stable? Check: a2=3>0a_2 = 3 > 0, a0=2>0a_0 = 2 > 0, and a2a1=12>2=a0a_2 a_1 = 12 > 2 = a_0. Yes, all roots have Re(λ)<0\text{Re}(\lambda) < 0.

Is λ3+λ2+2λ+8=0\lambda^3 + \lambda^2 + 2\lambda + 8 = 0 stable? a2a1=2<8=a0a_2 a_1 = 2 < 8 = a_0. No, the system is unstable.


Exponential Stability

Definition7.4Exponential stability

The equilibrium x0=0\mathbf{x}_0 = \mathbf{0} is exponentially stable if there exist C,α>0C, \alpha > 0 such that x(t)Ceαtx(0)\|\mathbf{x}(t)\| \leq Ce^{-\alpha t}\|\mathbf{x}(0)\| for all t0t \geq 0 and all initial conditions.

RemarkLinear systems and exponential stability

For constant-coefficient linear systems, asymptotic stability and exponential stability are equivalent. The decay rate α\alpha can be taken as any number less than minj(Re(λj))\min_j(-\text{Re}(\lambda_j)). For nonautonomous or nonlinear systems, the distinction matters.