Routh Criterion Calculator

What is the Routh Criterion?

The Routh-Hurwitz stability criterion is a mathematical test used in control system theory to determine the stability of a linear time-invariant (LTI) system without having to solve for the system’s poles. A stable system is one where the output remains bounded (does not go to infinity) over time. This test is performed on the characteristic polynomial of the system’s transfer function. The routh criterion calculator automates this process, making stability analysis quick and efficient.

Engineers and students use this method to check if all the roots of the characteristic polynomial have negative real parts, a necessary and sufficient condition for system stability. A common misunderstanding is that all coefficients being positive guarantees stability; while it’s a necessary condition (for most cases), it is not sufficient, which is why the full Routh test is required.

The Routh Criterion Formula and Explanation

The core of the criterion is the construction of a Routh Array or Routh Table. This table is built from the coefficients of the characteristic polynomial:

a_n*s^n + a_{n-1}*s^{n-1} + ... + a_1*s + a_0 = 0

The first two rows are populated directly from the polynomial coefficients. Subsequent rows are calculated using determinants of the two rows above. For example, the elements of the third row (b_1, b_2, …) are calculated as:

b_1 = (a_{n-1}*a_{n-2} - a_n*a_{n-3}) / a_{n-1}

b_2 = (a_{n-1}*a_{n-4} - a_n*a_{n-5}) / a_{n-1}

The stability is then determined by examining the first column of the completed array. If there are no sign changes in the first column, the system is stable. The number of sign changes in the first column is equal to the number of roots with positive real parts (i.e., the number of poles in the right-half plane).

Variables Table

Variables used in the Routh Criterion analysis.

Variable	Meaning	Unit	Typical Range
a_n, a_{n-1},…	Coefficients of the characteristic polynomial	Unitless (context-dependent)	Real Numbers
s	Complex frequency variable (Laplace domain)	–	Complex Numbers
First Column Elements	Values in the first column of the Routh Array used for the stability test	Unitless	Real Numbers

Practical Examples

Example 1: A Stable System

Consider a system with the characteristic equation: s^3 + 6s^2 + 11s + 6 = 0.

• Inputs (Coefficients): 1, 6, 11, 6
• Routh Array: The first column of the generated array would be. All values are positive.
• Result: Since there are no sign changes, the system is Stable.

Example 2: An Unstable System

Consider a system with the characteristic equation: s^4 + 3s^3 + s^2 + 8s + 1 = 0.

• Inputs (Coefficients): 1, 3, 1, 8, 1
• Routh Array: The first column of the array would be [1, 3, -1.67, 10, 1].
• Result: There are two sign changes (from 3 to -1.67, and from -1.67 to 10). Therefore, the system is Unstable with two poles in the right-half plane. You can verify this with our Root Locus Plotter.

How to Use This Routh Criterion Calculator

1. Enter Coefficients: Type the coefficients of your system’s characteristic equation into the input field. Start from the highest power of ‘s’ and go down to the constant term. Separate each coefficient with a comma.
2. Calculate Stability: Click the “Calculate Stability” button.
3. Interpret Results:
   • The primary result will clearly state if the system is “Stable” or “Unstable”.
   • A summary will tell you the number of roots in the right-half plane (RHP), which cause instability.
   • The complete Routh Array is displayed, allowing you to see the calculations and verify the result. This is especially useful for learning, and can be compared with results from a Nyquist Plot Calculator.

Key Factors That Affect System Stability

• Coefficient Values: The magnitude and sign of the coefficients directly influence the calculations in the Routh array. A negative or zero coefficient can be an early indicator of instability.
• System Gain (K): In feedback systems, a gain ‘K’ is often part of the coefficients. Varying this gain can move a system from stable to unstable. The Routh criterion is excellent for finding the range of ‘K’ for which the system remains stable.
• System Order: Higher-order systems (with higher powers of ‘s’) have more complex characteristic equations and more rows in the Routh array, increasing the potential for instabilities.
• Pole Locations: The criterion is fundamentally a test for pole locations. Any pole in the right-half of the complex plane leads to an unstable response. Our Bode Plot Analyzer can provide more frequency-domain insights.
• Time Delays: Pure time delays in a system can introduce transcendental terms (like e^-sT), which must be approximated by a polynomial (e.g., using a Padé approximation) before the Routh criterion can be applied.
• Special Cases: A zero in the first column or an entire row of zeros requires special handling (the epsilon method or an auxiliary polynomial) and often indicates poles on the imaginary axis (marginal stability) or other specific root patterns.

Frequently Asked Questions (FAQ)

What does a sign change in the first column mean?

A sign change indicates the presence of a root of the polynomial in the right-half of the complex plane (a pole with a positive real part). Each sign change corresponds to one such root. The system is unstable if there is one or more sign changes.

What if a zero appears in the first column?

This is a special case. To proceed, the zero is replaced by a small positive number, epsilon (ε), and the calculation continues. The sign of the elements in the first column is then evaluated as ε approaches zero. This calculator handles this “epsilon method” automatically.

What if an entire row becomes zero?

This indicates that the polynomial has roots that are symmetric about the origin, such as roots on the imaginary axis (leading to marginal stability) or pairs of real roots with opposite signs. An auxiliary polynomial is formed from the row just above the row of zeros to continue the analysis.

Are the coefficients’ units important?

For the mathematical procedure, the coefficients are treated as dimensionless numbers. However, they are derived from the physical parameters of a system (like mass, resistance, etc.), and their units implicitly define the system’s behavior. The Routh test itself is a purely mathematical process.

Can this calculator handle any polynomial?

Yes, it can handle a polynomial of any practical order. Simply enter all the coefficients, and the array will be generated accordingly.

Is the Routh criterion always enough to guarantee stability?

Yes, for linear time-invariant (LTI) systems, the Routh criterion is a necessary and sufficient condition. If the test shows no right-half plane poles, the system is stable.

What is a “characteristic equation”?

In control systems, it is the denominator of the closed-loop transfer function. The roots of this equation are the poles of the system, which determine its stability. You can learn more with our guide to PID Controller Tuning.

What is the difference between stable, unstable, and marginally stable?

A stable system’s response to a bounded input will decay to zero or a finite value. An unstable system’s response will grow infinitely. A marginally stable system will oscillate with a constant amplitude indefinitely, which occurs when there are non-repeated poles on the imaginary axis (a case the Routh test can identify).

Related Tools and Internal Resources

Explore these other tools and resources to deepen your understanding of control systems and stability analysis:

• Nyquist Plot Calculator: Analyze stability using the Nyquist criterion, which is particularly useful for systems with delays.
• Bode Plot Analyzer: Visualize system gain and phase margins to understand frequency response and relative stability.
• Root Locus Plotter: See how system poles move as a parameter (like gain ‘K’) changes.
• PID Controller Basics: An introductory guide to the most common type of controller in industrial applications.