Introduction au critère de Routh-Hurwitz

Le critère de Routh-Hurwitz est une pierre angulaire de la théorie classique du contrôle, fournissant une méthode algébrique rapide pour évaluer la stabilité d'un système linéaire invariant dans le temps (LTI) sans calculer explicitement les racines de son polynôme caractéristique. Développé indépendamment par Edward John Routh en 1876 et Adolf Hurwitz en 1895, la méthode reste indispensable pour les ingénieurs qui doivent garantir un comportement stable et prévisible dans les systèmes de contrôle de rétroaction, les réseaux électriques et les structures mécaniques.

Cet article fournit un guide détaillé et étape par étape pour construire et interpréter la table Routh-Hurwitz, en mettant l'accent sur l'application réelle. Vous apprendrez à détecter l'instabilité à partir des changements de signe, gérer des cas spéciaux tels que zéros dans la première colonne ou des lignes entières de zéros, et appliquer le critère à une variété de domaines d'ingénierie.

Construction de la table Routh-Hurwitz

P(s) = an[snn[n-1s[n-1[n+1[[FLT:]n[s+[n[n1[FLT:]s+[0[FLT:]]0[FLT:]0[FLT:]] = 0], où tous les coefficients sont des nombres réels et [[FLT:]an] > 0]] par convention (si a][F[

Étape 1: Disposer les deux premières rangées

Énumérer les coefficients du polynôme en rangées alternées:

  • sn: contient des coefficients an, an-2, an-4, ... (coefficients de termes indexés de façon uniforme, commençant par le plus haut degré).
  • sn-1: contient des coefficients an-1, an-3, a]n-5, ... (coefficients de termes indexés impairs).

Continuer jusqu'à ce qu'il ne reste plus de coefficients; remplir les entrées vides avec des zéros.

Étape 2: Calculer les lignes suivantes

Chaque ligne suivante (pour sn-2, s[n-3, ...) est calculée à l'aide d'une formule déterminante. Pour un élément de la ligne i et de la colonne j, la formule est:

bj = -(det[[a[1,1[, a1,j+1[], [a[2,1, a[2,j+1]])/a[2,1]]

où a1,1[ et a2,1 sont les premiers éléments des deux lignes précédentes. En pratique, le motif est:

  • Prenez le premier élément de la ligne précédente (qui devient le dénominateur).
  • Créer une matrice 2×2 en utilisant la première colonne des deux lignes précédentes et les colonnes suivantes de ces lignes.
  • Calculez le déterminant, changez son signe et divisez par le dénominateur.

Par exemple, pour calculer le premier élément de la ligne sn-2ssnssn-1:

b1 = (a[n-1 * an-2 - a[n * an-3]) / a[n-1]]

Répétez pour chaque colonne jusqu'à ce que la ligne soit complète (le nombre d'éléments dans chaque ligne diminue au fur et à mesure que vous descendez). Continuez jusqu'à ce que vous atteignez la ligne s0, qui ne contient qu'un seul élément.

Interprétation de la stabilité à partir de la première colonne

Une fois le tableau terminé, vous devez vous concentrer exclusivement sur la première colonne (la colonne la plus à gauche de chaque ligne).

  • Toutes les entrées de la première colonne doivent être positives (ou toutes négatives) pour que le système soit stable. Si une entrée est nulle ou qu'un changement de signe se produit, le système a des racines dans la moitié droite s-plan et est instable.
  • Le nombre de changements de signe dans la première colonne (de positif à négatif ou vice versa) est égal au nombre de racines avec des parties réelles positives. Une racine sur l'axe imaginaire (stabilité marginale) apparaît comme un zéro dans la première colonne ou une rangée entière de zéros.

Exemple : Système stable

Considérez le polynôme P(s) = s4 + 2s3 + 3s2 + 4s + 5]. La table Routh-Hurwitz est construite comme suit:

s rowCoefficients
s41, 3, 5
s32, 4, 0
s2(2*3 - 1*4)/2 = 1, (2*5 - 1*0)/2 = 5, 0
s1(1*4 - 2*5)/1 = -6, (1*0 - 2*0)/1 = 0
s0(-6*5 - 1*0)/(-6) = 5

Les premières entrées de colonnes sont : 1, 2, 1, -6, 5. Il y a un changement de signe (de 1 à -6, puis de retour à 5 compte comme deux changements de signe? En fait, un seul changement de signe: positif à négatif à la ligne [. Puis négatif à positif à 0[? Attendre, attention: Les valeurs de première colonne: ligne s^4: +1, s^3: +2, s^2: +1, s^1: -6, s^0: +5. Sign va +,+,+,-,+ => il y a deux changements de signe: + à - (entre s^2 et s^1), et - à + (entre s^1 et s^0).

Pour un système vraiment stable, par exemple P(s) = s2 + 2s + 1, la première colonne serait [1, 2, (2*1 - 1*0)/2 = 1] => tous positifs => stables.

Erreur fréquente : Changements de signe de comptage

Chaque transition de positif à négatif ou négatif à positif compte comme changement de signe. Ne sautez pas les lignes; incluez la ligne s0.

Applications pratiques en génie

Le critère de Routh-Hurwitz n'est pas seulement un exercice théorique, il est largement utilisé dans la conception et l'analyse dans plusieurs disciplines d'ingénierie.

Conception du système de contrôle

Dans le contrôle de la rétroaction, l'équation caractéristique provient de la fonction de transfert en boucle fermée. Les ingénieurs règlent les gains du contrôleur (p. ex., paramètres PID) et utilisent la table Routh-Hurwitz pour trouver des plages de gains qui maintiennent tous les pôles dans le plan de la moitié gauche. Ceci est particulièrement utile pendant la phase de conception précoce avant la simulation, car il fournit des limites de stabilité immédiates sans résoudre les polynômes à haut ordre.

Génie aérospatial

Aircraft and spacecraft flight control systems must be stable under all flight conditions. The Routh-Hurwitz criterion is used to assess stability margins for autopilots, stability augmentation systems, and structural mode suppression. Given the complexity of aircraft dynamics (often represented by state-space models of order 10 or higher), the criterion offers a quick check after linearization. It also helps in analyzing the effect of parameter variations (e.g., fuel consumption) on stability.

Circuits électriques et systèmes d'alimentation

La stabilité des réseaux électriques — des circuits simples d'opamp aux réseaux électriques à grande échelle — peut être évaluée à l'aide du polynôme caractéristique de la fonction de transfert du système. Par exemple, le lieu de base d'un stabilisateur du système d'alimentation peut être informé par l'analyse de Routh-Hurwitz pour empêcher les oscillations de basse fréquence.

Systèmes mécaniques et robotique

Les manipulateurs robotiques, les suspensions de véhicules et les structures flexibles ont souvent des équations caractéristiques dérivées des équations Newton-Euler ou Lagrange. Le critère de Routh-Hurwitz aide à déterminer si un ensemble donné de constantes de ressort/d'éboulement conduit à un mouvement stable.

Contrôle des procédés chimiques

Dans les industries de process, les régulateurs de température, de pression et de débit doivent être réglés pour éviter les oscillations ou les conditions de fuite. La table Routh-Hurwitz est appliquée à l'équation caractéristique du modèle linéaire d'un réacteur chimique, aidant à fixer des gains de PID qui assurent la stabilité malgré les retards de temps (approximation via Padé polynomials).

Traitement des cas spéciaux dans la table Routh-Hurwitz

Deux anomalies communes se produisent lors de la construction du tableau : un zéro dans la première colonne (mais pas une rangée entière), et une rangée entière de zéros. Les deux nécessitent un traitement spécial pour poursuivre l'analyse.

Zéro dans la première colonne (mais pas tous les zéros)

Si le premier élément d'une ligne est zéro alors que le reste de la ligne est non zéro, la table ne peut pas être poursuivie avec la formule normale parce que la division par zéro est nécessaire. La solution standard est de remplacer le zéro par un petit nombre positif ε[ (épsilon) et de procéder aux calculs. Après avoir complété la table, analyser les signes de la première colonne comme ε → 0+. Le nombre de changements de signe comme ε approche zéro du côté positif indique le nombre de pôles de demi-plan droit. Cette méthode est valable parce qu'une première colonne zéro correspond souvent à une paire de racines d'axe imaginaire ou à un système légèrement stable.

Exemple:] Polynominal P(s) = s3+2[fLT:]2[FLT:]2[F]2[f]2[fLT:]2[fLT:]2[f]2[fLT:]2[f]2[fLT:]2[f]2[fLT:]2[fLT:]2[f]2[fLT:]2[fLT:]2[f]2[fLT:]22[fLT:]2[fLT:]2[fLT:]2[f]2[fLT:]2[fLT:]2[fLT:]2[f]2[fLT:]2[f]2[fLT:]2[f]2[fLT:]2[fLT:]2[f]2[fLT:]2[fLT:]2[f]2[f]2[f.]2[f.

Toute la rangée de zéros

[FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][FLT:][F:][F][F.][F

+ + + + + + + + + + + + + + + + + +

Limitations du critère de la rugosité

Bien que puissant, le critère de Routh-Hurwitz comporte plusieurs limites que chaque ingénieur doit comprendre :

  • Seulement applicable aux systèmes linéaires invariants du temps (LTI) avec des équations caractéristiques polynômes. Il ne peut pas gérer directement les retards, les non-linéarités ou les paramètres de variation du temps.
  • Aucune information directe sur la réponse transitoire (telle que le temps de tassement, le dépassement, le rapport d'amortissement).Le critère indique seulement la stabilité ou l'instabilité, et non la performance.
  • Les erreurs de coefficient sont sensibles — de petites erreurs de mesure ou d'arrondi peuvent conduire à des évaluations incorrectes des signes, en particulier près des limites de stabilité.
  • Ne quantifie pas les marges de stabilité (gain, marge de phase), pour cela, des méthodes de domaine de fréquence comme les placettes de Nyquist ou de Bode sont nécessaires.
  • Il est impossible de manipuler facilement les polynômes avec des coefficients symboliques au-delà de paramétrages simples (p. ex., un gain unique).

Malgré ces limites, la table Routh-Hurwitz reste un outil de première ligne rapide, intuitif et mathématiquement rigoureux pour l'analyse de la stabilité.

Conseils pratiques pour l'utilisation de la table Routh-Hurwitz

  1. Pour la stabilité, tous les coefficients du polynôme caractéristique doivent être positifs. Tout coefficient zéro ou négatif (surtout pour les termes à haut degré) signale immédiatement une instabilité possible — mais cela ne suffit pas; le tableau est toujours nécessaire pour confirmer.
  2. Utiliser la précision numérique:[ Lors de la construction de grandes tables à la main, utiliser des fractions symboliques ou l'arithmétique de haute précision pour éviter les erreurs de signe.
  3. Combine avec le locus racine: Pour la conception, utilisez le critère de Routh-Hurwitz pour trouver la plage d'un paramètre gain, puis utilisez le locus racine pour examiner la migration du pôle boucle fermée pour déterminer le comportement transitoire.
  4. En pratique, vérifier avec simulation:[ Après avoir appliqué le critère, toujours simuler le système (p. ex., réponse par étapes) dans des logiciels comme MATLAB, Simulink ou Python Control Library pour confirmer.
  5. Pour les systèmes à ordre élevé (order > 10), envisager d'utiliser des méthodes numériques (calcul de la valeur propre) pour déterminer la précision, mais utiliser la table Routh-Hurwitz pour comprendre la sensibilité des paramètres.

Ressources externes et lectures complémentaires

Pour une compréhension plus approfondie du critère de Routh-Hurwitz et de sa place dans la théorie du contrôle, les ressources suivantes sont recommandées :

Conclusion

La maîtrise de la table Routh-Hurwitz permet aux ingénieurs d'évaluer rapidement la stabilité du système pendant la phase de conception initiale, sans recourir à des algorithmes de résolution de racine complexes. En comprenant comment construire la table, interpréter les changements de signe et gérer des cas spéciaux (méthode d'épsilon et polynôme auxiliaire), vous obtenez un outil fiable et intuitif qui fonctionne pour les polynômes de n'importe quel ordre. La connexion directe du critère aux modèles de coefficients le rend particulièrement utile pour les études paramétriques, l'ajustement et l'estimation de la marge de stabilité dans un large éventail de disciplines d'ingénierie, de la commande aérospatiale et automobile à l'électronique de puissance et à l'automatisation des processus.