Routh Criterion: Master Stability Analysis! (60 Char)

The Routh-Hurwitz stability criterion, a powerful technique in control systems engineering, determines system stability by analyzing the coefficients of its characteristic equation. Nyquist plots offer a graphical method to assess stability, while the Laplace transform converts differential equations into algebraic ones, facilitating analysis with tools like the routh criterion. Organizations such as the IEEE frequently publish research and standards that utilize and advance understanding of the routh criterion for stability analysis. This article provides a comprehensive guide to mastering the routh criterion for robust stability analysis.

Mastering Stability Analysis with the Routh Criterion

This article provides a comprehensive guide to understanding and applying the Routh Criterion, a powerful tool for determining the stability of linear time-invariant (LTI) systems. We will cover the fundamental principles, the step-by-step procedure for constructing the Routh array, and techniques for handling special cases.

Introduction to Stability and the Routh Criterion

What is System Stability?

In control systems, stability refers to the ability of a system to return to its equilibrium state after being subjected to a disturbance. A stable system’s output remains bounded for any bounded input. Conversely, an unstable system’s output grows without bound, even for a small disturbance.

Why is Stability Important?

Ensuring system stability is crucial for predictable and safe operation. An unstable control system can lead to oscillations, uncontrolled movements, and potentially catastrophic failures.

Introducing the Routh Criterion

The Routh criterion is a mathematical method used to determine the stability of a linear time-invariant (LTI) system by analyzing the coefficients of its characteristic equation. It provides a straightforward way to assess stability without explicitly solving for the roots of the characteristic equation, which can be computationally expensive for higher-order systems.

Constructing the Routh Array

The core of the Routh criterion is the construction of a Routh array (also known as a Routh table). Here’s the step-by-step procedure:

1. Obtain the Characteristic Equation: The starting point is the characteristic equation of the system, typically represented as a polynomial in the Laplace variable ‘s’:

   P(s) = a_nsⁿ + a_n-1s^n-1 + … + a₁s + a₀ = 0

   Where a_n, a_n-1, …, a₁, a₀ are the coefficients of the polynomial and a_n > 0. It’s important that a_n is positive, otherwise multiply the whole equation by -1 to ensure this.

2. Form the First Two Rows: The first two rows of the Routh array are formed using the coefficients of the characteristic equation.

   • Row 1 (sⁿ): The coefficients of the even powers of ‘s’ (a_n, a_n-2, a_n-4, …) are placed in the first row.
   • Row 2 (s^n-1): The coefficients of the odd powers of ‘s’ (a_n-1, a_n-3, a_n-5, …) are placed in the second row.

   Example: For P(s) = s⁴ + 3s³ + 5s² + 4s + 2 = 0

   The first two rows would be:

   s4	1	5	2
   s3	3	4	0

3. Calculate the Remaining Rows: Each subsequent row is calculated using the elements of the two rows immediately above it. The general formula for calculating an element b_i,j in row i is:

   b_i,j = – (determinant of [[a_i-2,1, a_i-2,j+1], [a_i-1,1, a_i-1,j+1]]) / a_i-1,1

   Where:

   • a_i-2,1 is the first element of the row two rows above the current row.
   • a_i-2,j+1 is the (j+1)^th element of the row two rows above the current row.
   • a_i-1,1 is the first element of the row immediately above the current row.
   • a_i-1,j+1 is the (j+1)^th element of the row immediately above the current row.

   Continue calculating rows until you reach the s⁰ row.

4. Example Construction: Continuing with P(s) = s⁴ + 3s³ + 5s² + 4s + 2 = 0

   Calculating the s² row:

   • b_1,1 (element under ‘1’): – (determinant of [[1, 5], [3, 4]]) / 3 = -(4-15)/3 = 11/3
   • b_1,2 (element under ‘5’): – (determinant of [[1, 2], [3, 0]]) / 3 = -(0-6)/3 = 2

   Calculating the s¹ row:

   • b_2,1 (element under ‘3’): – (determinant of [[3, 4], [11/3, 2]]) / (11/3) = -(6-44/3) / (11/3) = 26/11

   Calculating the s⁰ row:

   • b_3,1 (element under ’11/3′): – (determinant of [[11/3, 2], [26/11, 0]]) / (26/11) = -(0 – 52/11) / (26/11) = 2

   The complete Routh array is:

   s4	1	5	2
   s3	3	4	0
   s2	11/3	2	0
   s1	26/11	0	0
   s0	2	0	0

Applying the Routh Criterion for Stability Determination

Routh’s Stability Theorem

The Routh criterion states that the number of roots of the characteristic equation with positive real parts (i.e., located in the right-half of the s-plane) is equal to the number of sign changes in the first column of the Routh array.

Stability Condition

• For a system to be stable, all the elements in the first column of the Routh array must be positive. In other words, there should be no sign changes in the first column.

Applying the Criterion to the Example

In our example, the first column of the Routh array contains the elements 1, 3, 11/3, 26/11, and 2. All these values are positive, so there are no sign changes. Therefore, according to the Routh criterion, the system is stable because all the roots of the characteristic equation have negative real parts.

Special Cases in Routh Array Construction

The Routh array construction can encounter two main special cases:

Case 1: Zero in the First Column

If a zero appears in the first column of the Routh array, but the rest of the row is not entirely zero, it can cause division by zero in subsequent calculations. To resolve this:

1. Replace the zero with a small positive number, ε (epsilon).
2. Complete the Routh array using ε.
3. Analyze the stability by taking the limit as ε approaches zero. Observe the sign changes in the first column as ε approaches zero.

Case 2: Entire Row of Zeros

If an entire row of the Routh array consists of zeros, it indicates the presence of roots that are either:

• Located symmetrically about the origin (e.g., complex conjugate pairs on the imaginary axis).
• Or are real roots with equal magnitudes but opposite signs.

To resolve this:

1. Form the auxiliary equation: Construct an auxiliary equation, A(s), using the coefficients of the row immediately above the row of zeros. The auxiliary equation will have the form:

   A(s) = c_ks^k + c_k-2s^k-2 + c_k-4s^k-4 + … = 0

   Where the coefficients c_k, c_k-2, c_k-4, … are the elements of the row above the row of zeros. The power of ‘s’ in the first term, s^k, is the same as the power of ‘s’ for the row where the coefficients are taken.

2. Differentiate the auxiliary equation: Calculate the derivative of A(s) with respect to s, denoted as dA(s)/ds.

3. Replace the row of zeros with the coefficients of dA(s)/ds. Complete the Routh array using these new coefficients.

4. Analyze the Stability: Analyze the sign changes in the first column of the completed Routh array. Additionally, the roots of the auxiliary equation A(s) = 0 are also roots of the original characteristic equation, and their location determines partial stability.

Alright, hopefully, you’ve got a solid grasp on the routh criterion now! Go ahead and test it out on some systems. Let me know if you run into any snags!