Table of Contents
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 systématique pour évaluer la stabilité des systèmes linéaires invariants du temps (LTI) sans calculer explicitement les racines de l'équation caractéristique. Développé indépendamment par Edward John Routh en 1874 et Adolf Hurwitz en 1895, cette technique algébrique reste largement utilisée en génie du contrôle car elle offre une façon rapide et fiable de déterminer si un système réagira d'une manière stable, prévisible ou présente des oscillations et instabilités. La stabilité est fondamentale pour la conception du système – un système instable peut conduire à des débordements dangereux, des oscillations persistantes, voire à des défaillances catastrophiques, en particulier dans des applications telles que l'aérospatiale, les contrôles automobiles, les procédés chimiques et la robotique.
En analysant les coefficients d'un polynôme caractéristique du système, le critère de Routh-Hurwitz identifie le nombre de pôles (les racines de l'équation caractéristique) qui se trouvent dans la moitié droite du plan complexe. Les pôles avec des parties réelles positives causent un comportement instable, tandis que les pôles avec des parties réelles négatives contribuent à des réponses stables et convergentes. La méthode n'exige pas de résoudre des polynômes de haut degré, ce qui le rend pratique pour l'analyse manuelle et la conception en début de phase.
Cet article expliquera d'abord la base mathématique du critère de Routh-Hurwitz, puis vous guidera dans la construction du tableau de Routh avec un exemple étape par étape. Après cela, nous explorerons comment le critère prédit la réponse du système et le comportement transitoire, y compris les caractéristiques clés comme le dépassement et le temps de règlement. Un exemple numérique pratique illustrera le processus, et nous discuterons des limitations de la méthode et des outils complémentaires.
C'est quoi le critère de Routh-Hurwitz ?
Le critère de Routh-Hurwitz est une condition nécessaire et suffisante pour la stabilité des systèmes LTI. Il fonctionne sur l'équation caractéristique du système, qui est dérivée du dénominateur de la fonction de transfert en boucle fermée après l'avoir réglée à zéro. Pour un système décrit par la fonction de transfert , l'équation caractéristique est , où est un polynôme en degrés n:
Le critère indique que pour que le système soit stable (c.-à-d. que tous les pôles ont des parties réelles négatives), tous les coefficients ai doivent être non nuls et avoir le même signe (généralement positif), et le tableau Routh construit à partir de ces coefficients ne doit pas avoir de changement de signe dans la première colonne.
Il est important de noter que le critère de Routh-Hurwitz est purement algébrique, il ne donne pas les emplacements réels des pôles. Il fournit plutôt un compte de pôles de demi-plan droit (RHP) qui le rend inestimable lorsque le polynôme caractéristique est d'ordre élevé (par exemple, quatrième, cinquième ou plus) où la recherche de racine est intensive par calcul ou invraisemblable par analyse.
Le critère est directement lié à la matrice Hurwitz, une matrice carrée formée à partir des coefficients du polynôme caractéristique. Tous les principaux mineurs de la matrice Hurwitz doivent être positifs pour la stabilité, c'est la condition Hurwitz. La méthode tabulaire de Routh , simplifie ce processus de vérification déterminant dans le tableau plus simple de Routh.
Construction de la galerie Routh: étape par étape
Construire le réseau de Routh est une procédure systématique. Nous allons l'illustrer en utilisant un exemple concret polynôme:
Ce polynôme a des racines connues à s = -1, -2, -3, -4 (tous négatifs), donc le système est stable. Nous l'utiliserons pour vérifier que le tableau de Routh ne donne aucun changement de signe.
Étape 1: Écrire le polynôme sous forme standard
Pour notre exemple, les coefficients sont :[
a4 = 1 ([4), a[3 = 10 ([3), a[2] = 35 (2]), a[]1 = 50 [1]), a[0] = 24 (terme constant).
Étape 2: Initialiser les deux premières lignes
La première ligne du tableau de Routh contient les coefficients des termes indexés de façon uniforme (à partir de sn): ligne s[4: [1, 35, 24]. La deuxième ligne contient les coefficients des termes indexés de façon impaire: ligne s3: [10, 50, 0]. (Remplir les termes manquants avec zéros.)
rangées 4: 1, 35, 24
] rangées 3: 10, 50, 0
Étape 3: Calculer les lignes suivantes en utilisant la règle Rhombus
Chaque élément de la ligne suivante est calculé en utilisant les deux lignes au-dessus. Pour un élément de la ligne s2, première colonne, la formule est:
Mais il est plus facile d'appliquer la formule standard : pour la ligne k, colonne j (à partir de la colonne 1), l'élément est :
Faisons-le manuellement pour plus de clarté.
2:[
]Premier élément: Utiliser des lignes s4 et s3
R(s[,1) = [[ (1×50) - (10×35) ] / 10 = [50 - 350]/10 = [[300]/10 = 30
] Deuxième élément: colonne suivante (j=2). R(s4,2) = 35, R(s4
[F][LT:02][F][LT:02][2[F=10]][[[[[FLT:
Ainsi, ligne s2: 30, 24, 0.
Row s1[:[
Utiliser les lignes s[3 et s[2[[
R(s[]1,1) = -[ (10×24) - (30×50) ] / 30 = [240 - 1500]/30 = [[-1260]/30 = 42
] R(s[]]1[,2) = -[ (10×0) - (30×0) ] / 30 = 0.
Ainsi, ligne s1: 42, 0.
Row s0:[
Utiliser les lignes s[2 et s[1[
R(s[]0,1) = -[ (30×0) - (42×0) ] / 42 = 0. Mais la dernière ligne est toujours le terme constant du polynôme, qui devrait être 24. En fait, il y a une règle connue : le dernier élément de la dernière ligne est égal à a[]0 [le terme constant].
Donc le tableau complet de Routh est :
| s4 | 1 | 35 | 24 |
|---|---|---|---|
| s3 | 10 | 50 | 0 |
| s2 | 30 | 24 | 0 |
| s1 | 42 | 0 | 0 |
| s0 | 24 | 0 | 0 |
Étape 4: Examiner la première colonne
La première colonne contient: 1, 10, 30, 42, 24. Tous sont positifs, et il n'y a pas de changement de signe.
Cas particuliers
Deux cas particuliers peuvent survenir lors de la construction de la chaîne de Routh:
- Zero dans la première colonne : Si le premier élément d'une ligne est zéro mais que la ligne n'est pas tous zéros, remplacez le zéro par un petit épsilon (ε), puis procédez. Le signe de ε détermine les changements de signe. Ceci indique une stabilité marginale ou une instabilité selon le reste du tableau.
- Entrez la ligne de zéros: Cela se produit lorsqu'il y a des racines symétriquement placées (par exemple, des paires de racines imaginaires ou des paires symétriques autour de l'axe réel). Dans ce cas, formez un polynôme auxiliaire à partir de la ligne au-dessus de la ligne zéro et différenciez-le pour obtenir des coefficients pour la ligne manquante. Cela indique la présence d'une paire de racines purement imaginaires ou d'une paire de racines réelles avec des signes opposés (c.-à-d., instabilité).
Ces cas particuliers étendent l'applicabilité du critère aux systèmes à pôles sur l'axe imaginaire ou les motifs symétriques.
Prédiction de la réponse du système et comportement transitoire
Le critère de Routh-Hurwitz ne se limite pas à un verdict binaire stable/non stable. En révélant le nombre de pôles avec des parties réelles positives, il informe directement le comportement transitoire du système. Chaque pôle du demi-plan droit contribue à un mode de croissance exponentielle, provoquant la réponse transitoire à diverger. Inversement, tous les pôles du demi-plan gauche conduisent à des décroissances exponentielles, atteignant éventuellement un état stable.
Cependant, le critère ne donne pas les emplacements exacts des pôles, il ne peut donc pas quantifier précisément le dépassement, le temps de tassement ou le temps de montée.
- Si la première colonne n'a pas de changement de signe:[ Tous les pôles ont des parties réelles négatives. La réponse transitoire sera stable, se résiliera à un état stable. Mais il pourrait y avoir une décomposition lente ou oscillatoire (prochaine des pôles imaginaires) que le critère ne détecte pas – seul le signe de la partie réelle est garanti.
- S'il y a des changements de signe: Le nombre de changements de signe est égal au nombre de pôles RHP. Le système est instable, et la réponse transitoire va croître sans limite. Le nombre de pôles RHP indique la nature de l'instabilité : un seul pôle positif conduit à une croissance exponentielle monotonique, tandis qu'une paire de pôles RHP complexes produit des oscillations croissantes.
- Stabilisation marginale:[ Si l'équation caractéristique a des racines sur l'axe imaginaire (purement imaginaire ou plusieurs pôles sur l'axe), le tableau de Routh produira une rangée de zéros. Cela indique que le système est à la limite de la stabilité, et la réponse transitoire peut présenter des oscillations soutenues (si de simples pôles imaginaires) ou des oscillations en montée en rampe (si des pôles imaginaires répétés).
En combinant le résultat de Routh-Hurwitz avec la connaissance de l'ordre système, les ingénieurs peuvent faire des prédictions éclairées. Pour un système de deuxième ordre, le critère réduit à exiger que tous les coefficients soient positifs (ce qui est à la fois nécessaire et suffisant pour la stabilité).Pour les systèmes de plus haut ordre, le tableau fournit un indicateur de marge de stabilité – plus les éléments de première colonne sont proches de zéro, plus le système est proche de l'instabilité.
Relation avec les caractéristiques transitoires
La réponse transitoire d'un système stable est principalement déterminée par les pôles dominants, les plus proches de l'axe imaginaire. Le critère de Routh-Hurwitz ne localise pas directement ces pôles, mais il peut aider à l'accordage. Par exemple, si le critère montre qu'augmenter un gain de contrôleur entraînera éventuellement un changement de signe dans la première colonne, il marque la marge de gain. Les ingénieurs utilisent souvent cela pour fixer des gains proportionnels dans les contrôleurs P, PI ou PID sans avoir à simuler l'ensemble du système à chaque fois.
Les principales mesures transitoires qui sont influencées par la stabilité comprennent :
- Dépassement – La valeur maximale de pointe de la réponse par rapport à la valeur finale. Le dépassement est généralement causé par des pôles conjugués complexes avec un faible rapport d'amortissement. Bien que le critère Routh-Hurwitz ne calcule pas l'amortissement, un système avec une rangée de zéros (stabilité marginale) correspond souvent à un rapport d'amortissement proche de zéro, ce qui implique un dépassement élevé.
- Temps de règlement – Le temps nécessaire pour que la réponse reste dans un certain pourcentage (par exemple, 2%) de la valeur finale. Le temps de règlement est inversement proportionnel à la partie réelle du pôle dominant. Le tableau de Routh ne peut pas donner cette partie réelle, mais si le tableau montre que tous les pôles sont dans le demi-plan gauche, le temps de règlement est fini; si un pôle est très proche de l'axe imaginaire, le temps de règlement sera long.
- Heure de montée – Le temps de passer de 10 % à 90 % de la valeur finale. Ceci est lié à la fréquence naturelle des pôles complexes dominants. Là encore, le critère de Routh-Hurwitz ne fournit aucune valeur numérique directe, mais indique si le système est suffisamment stable pour avoir un temps de montée bien défini.
- Oscillations – La fréquence des oscillations transitoires est déterminée par les parties imaginaires des pôles complexes. Le critère ne révèle pas les parties imaginaires, mais la présence d'une rangée de zéros indique des pôles purement imaginaires, qui produisent des oscillations soutenues dans la réponse transitoire.
Ainsi, le critère de Routh-Hurwitz agit comme gardien de porte : il assure la stabilité comme condition préalable à une analyse transitoire significative à l'aide d'autres outils comme les locus racinaires, les tracés de Bode ou les simulations de réponse par étapes.
Exemple pratique : Système de troisième commande
Appliquez le critère de Routh-Hurwitz à un exemple plus complexe avec des paramètres inconnus. Considérons la fonction de transfert en boucle fermée avec retour d'unité. L'équation caractéristique est :
Nous voulons déterminer la plage de K pour laquelle le système est stable. Construire le tableau Routh pour polynôme .
s3[: 1, 10[
s[2[: 7, K
s[1[
[
s[]0: En fait, plus simple: l'élément de la dernière ligne est le terme constant K (à condition que la ligne ci-dessus n'est pas tous les zéros). Plus systématique: la ligne s[]0] est dérivée des lignes s2 et s1.
| s3 | 1 | 10 |
|---|---|---|
| s2 | 7 | K |
| s1 | (K-70)/7 | 0 |
| s0 | K | 0 |
Pour la stabilité, tous les éléments de la première colonne doivent être positifs : 1 > 0, 7 > 0, (K-70)/7 > 0, K > 0. Cela donne K > 70 et K > 0, donc K > 70. La limite inférieure est K > 70. Si K = 70, la ligne s[1 devient zéro, ce qui est un cas particulier indiquant une stabilité marginale (racines purement imaginaires). Ainsi le système est stable pour K > 70. Ceci indique à l'ingénieur que les gains au-dessus de 70 donnent une réponse transitoire stable, mais que K approche 70, les oscillations transitoires deviennent non amassées. La réponse transitoire pour K = 200, par exemple, aura deux pôles conjugués complexes dans le demi-plan gauche et un pôle réel, conduisant à une réponse de dépassement mais éventuellement de tassement.
Réaction transitoire Incidences pour l'exemple
Pour K juste au-dessus de 70, le système est stable mais les pôles dominants sont très proches de l'axe imaginaire. La réponse transitoire montrera une très lente désintégration et un grand dépassement. Le temps de décantation sera long. K augmente encore, les pôles se déplacent vers la gauche (généralement), améliorant ainsi le temps d'amortissement et de décantation. Cependant, si K est trop élevé, le tableau de Routh pour ce système particulier ne montrera aucun changement de signe supplémentaire – la stabilité persiste.
Ainsi, le critère donne des informations critiques : la plage de gain pour la stabilité et la proximité de l'instabilité. En évaluant les éléments de première colonne, les ingénieurs peuvent calculer les marges de gain et les contrôleurs de conception qui évitent la stabilité marginale.
Limitations du critère de la rugosité
Malgré sa puissance, le critère de Routh-Hurwitz a d'importantes limites que tout ingénieur de contrôle doit comprendre :
- Applicable uniquement aux systèmes LTI:[ Il ne peut pas être utilisé directement pour les systèmes de données de temps variable, non linéaires ou échantillonnés (bien que des extensions existent pour les systèmes à temps discret par des transformations bilinéaires).
- Ne quantifie pas les performances transitoires:[ Le critère ne donne aucune information sur le dépassement, le temps de pointe, le temps de montée ou le temps de tassement exact. Pour cela, le locus racine, la réponse de fréquence ou les simulations de domaine temporel sont nécessaires.
- Nécessite un polynôme caractéristique connu: Les coefficients doivent être connus exactement. En pratique, les incertitudes dans les paramètres peuvent nécessiter une analyse de stabilité robuste (p. ex., en utilisant le théorème de Kharitonov).
- Devient lourd pour les systèmes à très haut ordre: Les calculs manuels pour les polynômes à dixième ordre sont sujets à des erreurs, bien que le logiciel (MATLAB, Python) automatise la construction du tableau.
- Ne peut pas gérer les délais directement :[ Les systèmes avec des délais ont des équations caractéristiques transcendantales – le tableau de Routh ne peut être appliqué directement que si l'on utilise des approximations de Padé.
Pour surmonter ces limites, les ingénieurs utilisent couramment le critère de Routh-Hurwitz en conjonction avec d'autres méthodes. Par exemple, les tracés de locus racinaires montrent comment les pôles en boucle fermée se déplacent avec gain, fournissant des informations transitoires exactes. Le critère de Nyquist gère les retards et les non-linéarités.
Outils complémentaires
Le critère de Routh-Hurwitz est nécessaire et suffisant pour la stabilité des polynômes LTI, mais pour l'analyse des performances, les ingénieurs s'appuient sur :
- Root Locus: Attribue la trajectoire des pôles comme paramètre (p. ex. gain) varie. Il montre directement le rapport d'amortissement et la fréquence naturelle.
- Plots de la tige et diagrammes de Nyquist: Fournit une marge de gain et une marge de phase, qui sont des indicateurs de stabilité relative et de caractéristiques de réponse transitoire.
- Simulation de la réponse d'étape:[ Donne les mesures transitoires réelles (dépassement, temps de réglage) pour un système donné.
- Méthodes d'état-espace: Les valeurs propres de la matrice d'état (qui sont les pôles) peuvent être calculées numériquement, et la méthode Lyapunovs étend l'analyse de stabilité aux systèmes non linéaires.
Le critère Routh-Hurwitz occupe donc un point d'intérêt historique et pratique : il est assez mathématique pour être rigoureux, mais assez simple pour être réalisé sur un tableau blanc lors d'une revue de conception. Il reste une compétence essentielle pour tout ingénieur de contrôle.
Conclusion
Le critère de Routh-Hurwitz offre une façon claire et algébrique d'évaluer la stabilité du système de contrôle en analysant le polynôme caractéristique sans résolution pour les racines. Grâce à la construction du tableau de Routh, les ingénieurs peuvent déterminer le nombre de pôles instables, comprendre la sensibilité du système aux changements de paramètres et identifier les limites de stabilité qui affectent directement le comportement transitoire.
La maîtrise du critère de Routh-Hurwitz vous permet d'évaluer rapidement si un système de contrôle se comportera de façon fiable dans des conditions transitoires. Que vous conçoyiez un simple contrôleur proportionnel ou que vous accordiez un système multi-boucle complexe, le tableau de Routh est un point de départ fiable. Pour plus de détails, envisagez d'explorer le critère de stabilité Routh-Hurwitz sur Wikipedia pour une profondeur mathématique supplémentaire, ou consultez les Tutoriels de contrôle de l'Université du Michigan sur Routh Array pour des exemples plus travaillés.