Table of Contents
Introducción a la Criterio Routh-Hurwitz
El criterio de Routh-Hurwitz es una piedra angular de la teoría del control clásico, proporcionando un método sistemático para evaluar la estabilidad de los sistemas de tiempo-invariantes lineales (LTI) sin calcular explícitamente las raíces de la ecuación característica. Desarrollado independientemente por Edward John Routh en 1874 y Adolf Hurwitz en 1895, esta técnica algebraica sigue siendo ampliamente utilizada en la ingeniería de control porque ofrece una manera rápida y fiable para responder a un sistema estable
Analizando los coeficientes del polinomio característico de un sistema, el criterio de Routh-Hurwitz identifica el número de polos (raíz de la ecuación característica) que se encuentran en la mitad derecha del plano complejo. Los polos con partes reales positivas causan comportamiento inestable, mientras que los polos con partes reales negativas contribuyen a respuestas estables y convergentes. El método no requiere resolver polinomios de alto grado, lo que hace práctico
Este artículo explicará primero la base matemática del criterio Routh-Hurwitz, luego guiarlo mediante la construcción de la matriz de Routh con un ejemplo paso a paso. Después de eso, exploraremos cómo el criterio predice la respuesta del sistema y el comportamiento transitorio, incluyendo características clave como el sobresueldo y el tiempo de fijación. Un ejemplo numérico práctico ilustrará el proceso, y el autor discutiremos las limitaciones del método y las herramientas complementarias.
¿Cuál es el Criterión Routh-Hurwitz?
El criterio de Routh-Hurwitz es una condición necesaria y suficiente para la estabilidad en los sistemas LTI. Funciona en la ecuación característica del sistema, que se deriva del denominador de la función de transferencia de cierre cerrado después de establecerlo a cero. Para un sistema descrito por la función de transferencia , la ecuación característica es ], donde
El criterio indica que para que el sistema sea estable (es decir, todos los polos tienen partes reales negativas), todos los coeficientes ai] deben ser no cero y tener el mismo signo (normalmente positivo), y el array Routh construido a partir de estos coeficientes no debe tener cambios de signo en la primera columna. Si se producen cambios de signo, el número de cambios de signos equivale al número de polos positivos.
Es importante señalar que el criterio de Routh-Hurwitz es puramente algebraico, no da las ubicaciones de polos reales. En lugar de ello, proporciona un recuento de polos de derecha-half-plane (RHP). Esto lo hace invaluable cuando el polinomio característico es de alto orden (por ejemplo, cuarto, quinto o superior) donde la determinación de raíz es computacionalmente infesible o analíticamente infeasible.
El criterio está directamente ligado a la matriz Hurwitz, una matriz cuadrada formada por los coeficientes del polinomio característico. Todos los principales menores de la matriz Hurwitz deben ser positivos para la estabilidad, esta es la condición Hurwitz. El método tabular de Routh simplifica este proceso de determinación de verificación en la matriz Routh más directa.
Construyendo el Array de Routh: Paso a paso
Construir la matriz de Routh es un procedimiento sistemático. Lo ilustraremos usando un ejemplo de concreto polinomio:
Este polinomio tiene raíces conocidas en s = -1, -2, -3, -4 (todo negativo), por lo que el sistema es estable. Lo utilizaremos para verificar que el array Routh no da cambios de señal.
Paso 1: Escriba el polinomio en forma estándar
[LT] [LT] [LT] [L] [L]] [L] [L]] [L]] [L]] [L]] [L]] [L] [L]] [L] [L] [L] [L]] [L] [L] [L]] [L] [L]
Paso 2: Iniciar las dos primeras filas
La primera fila del array Routh contiene los coeficientes de términos uniformes (comenzando desde sn): fila s4: [1, 35, 24]. La segunda fila contiene los coeficientes de términos de índices impares: fila s3[FLT] [50]
Así:
Row s4: 1, 35, 24
Row s3: 10, 50, 0
Paso 3: Computar las filas posteriores utilizando la regla del Rhombus
Cada elemento en la siguiente fila se calcula utilizando las dos filas por encima de ella. Para un elemento en la fila s2, primera columna, la fórmula es:
Pero es más fácil aplicar la fórmula estándar: para la fila k, columna j (comenzando desde la columna 1), el elemento es:
Hagamos manualmente para la claridad.
[LT] [LT] [24] [FLT] [23] [4]
Así hilera s2: 30, 24, 0.
[LT:0] [FLT]
)] [(FLT:5]] [FLT] [30×2]] [FLT] [4]] [4]] [4]
Así hilera s1: 42, 0.
[LT:2]
[FLT: [FLT]] [FLT] [4]]] [FLT] [4]] [FLT] [4]] [4]] [4]]
Así que el conjunto completo de Routh es:
| s4 | 1 | 35 | 24 |
|---|---|---|---|
| s3 | 10 | 50 | 0 |
| s2 | 30 | 24 | 0 |
| s1 | 42 | 0 | 0 |
| s0 | 24 | 0 | 0 |
Paso 4: Examinar la Primera Columna
La primera columna contiene: 1, 10, 30, 42, 24. Todos son positivos, y no hay cambios de señal. Por lo tanto el sistema es estable, que coincide con nuestras raíces conocidas.
Casos especiales
Dos casos especiales] pueden surgir durante la construcción de la matriz de Routh:
- Zero en la primera columna: Si el primer elemento de una fila es cero, pero la fila no es todos los ceros, sustituya el cero con un pequeño epsilón (ε), luego proceda. El signo del ε determina los cambios de firma. Esto indica estabilidad o inestabilidad marginal dependiendo del resto de la matriz.
- Filera entera de ceros: Esto ocurre cuando hay raíces simétricamente colocadas (por ejemplo, pares de raíces imaginarias o pares simétricos sobre el eje real). En este caso, forman un polinomio auxiliar de la fila por encima de la fila cero y la diferencian para obtener coeficientes para la fila perdida. Esto indica la presencia de un par de raíces puramente imaginarias.
Estos casos especiales extienden la aplicabilidad del criterio a los sistemas con polos sobre el eje imaginario o patrones simétricos.
Respuesta del sistema de procesamiento y comportamiento transitorio
El criterio de Routh-Hurwitz no se limita a un veredicto binario estable/inestable. Al revelar el número de polos con partes reales positivas, informa directamente el comportamiento transitorio del sistema. Cada polo en el medio plano derecho contribuye a un modo de crecimiento exponencial, causando la respuesta transitoria a divergencia. Por el contrario, todos los polos en el medio plano izquierdo conducen a descomponer exponencialmente.
Sin embargo, el criterio no da los puntos de polo exactos, por lo que no puede cuantificar precisamente el sobresueldo, el tiempo de ajuste o el tiempo de ascenso. Sin embargo, proporciona pistas esenciales:
- Si la primera columna no tiene cambios de señal: Todos los polos tienen partes reales negativas. La respuesta transitoria será estable, eventualmente estable a un estado estable. Pero podría haber decaimiento lento o oscilatorio (polos casi imaginarios) que el criterio no detecta, sólo se garantiza la señal de la parte real.
- Si hay cambios de señal: El número de cambios de señal equivale al número de polos RHP. El sistema es inestable, y la respuesta transitoria crecerá sin límites. El número de polos RHP indica la naturaleza de la inestabilidad: un solo polo positivo conduce al crecimiento exponencial monotónico, mientras que un par de polos RHP complejos produce oscilaciones crecientes.
- Estabilidad marginal: Si la ecuación característica tiene raíces en el eje imaginario (puremente imaginarios o múltiples polos en el eje), la matriz de Routh producirá una fila de ceros. Esto indica que el sistema está en el límite de la estabilidad, y la respuesta transitoria puede exhibir oscilaciones sostenidas (si polos simples imaginarios) o oscilaciones de rampa.
Al combinar el resultado de Routh-Hurwitz con el conocimiento del orden del sistema, los ingenieros pueden hacer predicciones informadas. Para un sistema de segundo orden, el criterio reduce a exigir que todos los coeficientes sean positivos (lo cual es necesario y suficiente para la estabilidad).Para sistemas de mayor orden, el array proporciona un indicador de margen de estabilidad, cuanto más cerca los elementos de primer ciclo son a cero, más cerca el sistema es a la inestabilidad.
Relación con las características transitorias
La respuesta transitoria de un sistema estable está determinada principalmente por los polos dominantes —los más cercanos al eje imaginario. El criterio de Routh-Hurwitz no localiza directamente estos polos, pero puede ayudar en el ajuste. Por ejemplo, si el criterio muestra que aumentar un aumento de controlador eventualmente llevará a un cambio de signo en la primera columna, marca el margen de ganancia.Los ingenieros a menudo utilizan esto para establecer ganancias proporcionales en PID
Las métricas transitorias clave que están influenciadas por la estabilidad incluyen:
- Overshoot – El valor máximo de la respuesta relativa al valor final. Overshoot es típicamente causado por polos complejos conjugados con una baja relación de amortiguación. Mientras que el criterio de Routh-Hurwitz no compute damping, un sistema con una fila de ceros (estabilidad marginal) suele corresponder a una relación de amortiguación cercana a cero, lo cual implica un alto.
- Tiempo de instalación] – El tiempo necesario para la respuesta para permanecer dentro de un determinado porcentaje (por ejemplo, 2%) del valor final. El tiempo de asentamiento es inversamente proporcional a la parte real del polo dominante. El arreglo de la matriz de Routh no puede dar esta parte real, pero si el array muestra que todos los polos están en el medio plano izquierdo, el tiempo de asentamiento es finito muy cercano;
- Tiempo de ida] – El tiempo para pasar del 10% al 90% del valor final. Esto se relaciona con la frecuencia natural de los polos complejos dominantes. De nuevo, el criterio de Routh-Hurwitz no proporciona ningún valor numérico directo, pero indica si el sistema es lo suficientemente estable para tener un tiempo de ascenso bien definido.
- Oscilaciones] – La frecuencia de oscilaciones transitorias se determina por las partes imaginarias de los polos complejos. El criterio no revela partes imaginarias, pero la presencia de una fila de ceros indica polos puramente imaginarios, que producen oscilaciones sostenidas en la respuesta transitoria.
Así, el criterio de Routh-Hurwitz actúa como portero: garantiza la estabilidad como requisito para un análisis transitorio significativo utilizando otras herramientas como locus raíz, diagramas de Bode o simulaciones de respuesta paso.
Ejemplo práctico: Sistema de Tercer Orden
Aplicamos el criterio de Routh-Hurwitz a un ejemplo más complejo con parámetros desconocidos. Consideremos la función de transferencia de cierre con retroalimentación de unidad. La ecuación característica es:
Queremos determinar el rango de K para el cual el sistema es estable. Construir el array Routh para polinomial .
[LT] [LT] [FLT] [14]] s[FLT] [4]] [FLT] [4]] [FLT] [4]] [4]] [4]] [4] [4] [4]] [4]] [4]] [L] [L] [L] [L] [L] [L] [L] [L] [L] [L] [L]
| s3 | 1 | 10 |
|---|---|---|
| s2 | 7 | K |
| s1 | (K-70)/7 | 0 |
| s0 | K | 0 |
Para la estabilidad, todos los elementos de primer nivel deben ser positivos: 1 Ø 0, 7 √≥n 0, (K-70)/7 , K √≥ 0. Esto da K ¢70 y K √≥ 0, por lo que K √≥n 70. El límite inferior es K √≥ 70. Si K = 70, la ganancia de la fila s1 se convierte en cero, lo que indica una estabilidad de un caso estable
Implicaciones de respuesta transitoria para el ejemplo
Para K justo por encima de 70, el sistema es estable pero los polos dominantes están muy cerca del eje imaginario. La respuesta transitoria mostrará una decaimiento muy lento y una sobresuelción grande. El tiempo de ajuste será largo. Mientras K aumenta aún más, los polos se mueven hacia la izquierda (generalmente), mejorando el tiempo de amortiguación y de fijación.
Así, el criterio da información crítica: el rango de ganancia para la estabilidad, y la proximidad a la inestabilidad. Al evaluar los elementos de primer ciclo, los ingenieros pueden calcular los márgenes de ganancia y los controladores de diseño que evitan la estabilidad marginal.
Limitaciones de la Criterio Routh-Hurwitz
A pesar de su poder, el criterio de Routh-Hurwitz tiene importantes limitaciones que todo ingeniero de control debe entender:
- Aplicable únicamente a los sistemas LTI: No se puede utilizar directamente para sistemas de datos de tiempo, no lineales o de datos muestreados (aunque existen extensiones para sistemas discretos a tiempo mediante transformaciones bilineales).
- No cuantifica el rendimiento transitorio: El criterio no da información sobre el sobresueldo, el tiempo máximo, el tiempo de ascenso o el tiempo de fijación exacto. Para ello, se necesitan locus raíz, respuesta de frecuencia o simulaciones de tiempo-dominio.
- Requiere un polinomio característico conocido: Los coeficientes deben ser conocidos exactamente. En la práctica, las incertidumbres en los parámetros pueden requerir un análisis de estabilidad robusto (por ejemplo, usando el teorema de Kharitonov).
- Se vuelve engorroso para sistemas de muy alto orden: Los cálculos manuales para polinomios de décimo orden son propensas a errores, aunque el software (MATLAB, Python) automatiza la construcción de la matriz.
- No se puede manejar los retrasos del tiempo directamente: Los sistemas con retrasos en el tiempo tienen ecuaciones características trascendentales: la matriz de Routh no puede aplicarse directamente a menos que usen las aproximaciones de Padé.
Para superar estas limitaciones, los ingenieros utilizan comúnmente el criterio Routh-Hurwitz en conjunto con otros métodos. Por ejemplo, las tramas de langosta de raíz muestran cómo los polos de cierre se mueven con ganancia, proporcionando información de transito exacta. El criterio de Nyquist maneja retrasos de tiempo y no linearidades. Sin embargo, la matriz de Routh suele servir como un primer cheque rápido, especialmente durante el diseño preliminar o solución de problemas.
Herramientas complementarias
El criterio de Routh-Hurwitz es “necesario y suficiente” para la estabilidad de los polinomios LTI, pero para el análisis de rendimiento, los ingenieros dependen de:
- Locus de arranque:) Parcela la trayectoria de los polos como parámetro (por ejemplo, ganancia) varía. Muestra directamente la relación de amortiguación y la frecuencia natural.
- Parcelas de bodas y Diagramas de Nyquist: Proveer margen de ganancia y margen de fase, que son indicadores de estabilidad relativa y características de respuesta transitoria.
- Simulación de Respuesta de Paso: Da la métrica transitoria real (sobresueldo, tiempo de fijación) para un sistema dado.
- Métodos de Estado-Espacio: Los valores eigenados de la matriz estatal (que son los polos) pueden ser calculados numéricamente, y el método de Lyapunov extiende el análisis de estabilidad a sistemas no lineales.
El criterio de Routh-Hurwitz ocupa así un lugar dulce histórico y práctico: es lo suficientemente matemático para ser riguroso, pero lo suficientemente simple como para ser hecho en una pizarra durante una revisión de diseño.
Conclusión
El criterio Routh-Hurwitz ofrece una manera clara y algebraica de evaluar la estabilidad del sistema de control analizando el polinomio característico sin resolver para las raíces. Mediante la construcción de la matriz Routh, los ingenieros pueden determinar el número de polos inestables, entender la sensibilidad del sistema a los cambios de parámetro, e identificar los límites de estabilidad que afectan directamente el comportamiento transitorio.
El sistema de control de la universidad [LT] permite la implementación de los sistemas de control de la universidad, pero no es un punto de partida fiable.Si usted está diseñando un simple controlador proporcional o sintonizando un sistema complejo de multi-op, la matriz de la rutina es un punto de partida confiable.