Simplifier les polynomiaux complexes en utilisant l'analyse de stabilité de Routh-hurwitz

Q: Travailler avec des Coefficients complexes?

Pour les polynômes avec des coefficients complexes ( a i -C), des modifications telles que le [Michailov critère peuvent être employées. Cependant, dans la plupart des problèmes pratiques d'ingénierie de contrôle, l'équation caractéristique provient de paramètres physiques évalués en valeur réelle, de sorte que les coefficients réels sont la norme.

Introduction au critère de stabilité de Routh-Hurwitz

Le critère de stabilité de Routh-Hurwitz est une pierre angulaire de la théorie du contrôle et de l'analyse polynôme. Il fournit une méthode systématique pour déterminer le nombre de racines d'un véritable polynôme qui se trouve dans la moitié droite du plan complexe, sans exiger un calcul de racine explicite.

Un polynôme complexe désigne ici un polynôme avec des coefficients réels (comme c'est le cas dans les applications de contrôle) mais dont les racines peuvent être complexes. La stabilité de nombreux systèmes physiques – des circuits électroniques aux liaisons mécaniques – dépend du signe des parties réelles des racines caractéristiques des polynômes. Le critère Routh-Hurwitz répond efficacement à cette question, même pour les polynômes à ordre élevé qui ne sont pas pratiques à prendre en compte.

Cet article s'étend sur le résumé original, fournissant une présentation détaillée de la méthode, des exemples, des cas spéciaux et des applications pratiques. D'ici la fin, les lecteurs pourront appliquer la technique Routh-Hurwitz pour simplifier et analyser des polynômes complexes avec confiance.

Qu'est-ce que les polynomes complexes?

Dans le domaine de l'ingénierie de contrôle et des mathématiques, un polynôme complexe est généralement une expression de la forme:

P(s) = a_n sⁿ + a_n−1] sⁿ⁻¹ + ... + a_]1 s + a₀]

où les coefficients a_i sont des nombres réels et s est une variable complexe. Le terme -complexe -complexe ne fait pas référence aux coefficients (qui restent réels) mais au fait que les racines de l'équation P(s)=0] sont généralement des nombres complexes. Ces racines déterminent le comportement dynamique d'un système décrit par le polynôme.

Par exemple, l'équation caractéristique d'un système de deuxième ordre comme un amortisseur de masse est ms² + bs + k = 0. Selon le rapport d'amortissement, les racines peuvent être réelles et négatives (surdémontées) ou une paire de conjugués complexes avec partie réelle négative (sous-déprimée). Le critère de Routh-Hurwitz permet de vérifier si all] racines ont des parties réelles négatives, sans résoudre l'équation.

Pourquoi le calcul direct de la racine est-il peu pratique?

Pour les polynômes de l'ordre cinq ou plus, il existe rarement des solutions racinaires analytiques exactes (théorème Abel-Ruffini). Les méthodes numériques peuvent approximativement approximativement les racines, mais elles sont sensibles à l'incertitude des coefficients et peuvent manquer de signes subtils d'instabilité. Le critère Routh-Hurwitz, en revanche, donne une réponse définitive oui/non au signe des parties réelles, en utilisant uniquement les coefficients polynômes.

Le critère de stabilité de Routh–Hurwitz

Développé indépendamment par Edward John Routh (1876) et Adolf Hurwitz (1895), le critère indique que pour un polynôme avec des coefficients réels, toutes les racines ont des parties réelles négatives si et seulement si tous les principaux mineurs de la matrice de Hurwitz (ou, équivalent, la première colonne du tableau de Routh) sont positifs. La méthode est habituellement enseignée en utilisant le tableau de Routh parce qu'il est plus simple à construire à la main.

Construction de la galerie Routh

Compte tenu du polynôme:

P(s) = a_n sⁿ + a_n−1 sⁿ⁻¹ + ... + a₁ s + a[₀] (avec a_]n]] > 0[]),

le tableau Routh est construit comme suit:

[[a][[a[[FLT:]][[n−2[[FLT:]]][_{_[n−1]}_^[[

[[FLT:][[FLT:][[[FLT:]][[[FLT:]][[[FLT:]][[[FLT:]][[FLT:]][, chaque élément est dérivé des deux lignes ci-dessus à l'aide d'une formule déterminante:
[b_{[a[[[FLT:]][n−1[FLT:] [[FLT:]] [n−2[FLT:][FLT:]][n[FLT:][[FLT:][FLT:][[FLT:]][[[FLT:]][[FLT:][[FLT:]][[FLT:]][FLT:][FLT:][[FLT:][FLT:][[FLT:][[[FLT:]][[[[}

Continuer jusqu'à ce que la seule ligne constante reste. La dernière ligne correspond à s⁰ et contient un seul élément, qui est le terme constant (typiquement a₀.

Examinez la première colonne. Si toutes les entrées de la première colonne du tableau sont positives (et aucune ne sont nulles), alors toutes les racines sont dans le demi-plan gauche (le système est stable).Tout changement de signe indique le nombre de racines avec des parties réelles positives – chaque changement de signe compte une telle racine.

Cas spéciaux dans le quartier de Routh

Deux situations particulières exigent une manipulation attentive :

Zero dans la première colonne. Si le premier élément d'une ligne est zéro, mais le reste de la ligne est non zéro, remplacer le zéro par un petit nombre positif ε et procéder. Le signe de ε détermine le résultat de la stabilité. On peut aussi utiliser la méthode [voir ci-dessous].

Entrez la ligne de zéros. Cela indique que le polynôme contient des racines symétriquement localisées (paires de racines qui sont négatives l'une de l'autre ou des conjugués complexes avec la partie réelle zéro). Dans ce cas, formez un polynôme auxiliaire à partir de la ligne au-dessus des zéros, prenez sa dérivée et utilisez les coefficients de cette dérivée pour remplacer la ligne zéro. Continuez ensuite le tableau.

Exemple : Un polynôme de quatrième ordre

P(s) = s⁴ + 2s³ + 3s² + 4s + 5.

s⁴: 1, 3, 5

s³: 2, 4, 0

s²: b₁ = (2·3 - 1·4)/2 = (6-4)/2 = 1, b₂ = (2·5 - 1·0)/2 = 10/2 = 5]

Ligne s¹: c₁ = (1·4 - 2·5)/1 = (4-10)/1 = -6], puis zéro

s⁰: d[₁ = (-6·5 - 1·0)/(-6) = (-30)/(-6) = 5

Première colonne : 1, 2, 1, -6, 5 → le signe change une fois (de 1 à -6, et de -6 à 5), indiquant une racine avec partie réelle positive. En effet, le polynôme a des racines : -1,65 ± 0,35i, 0,65 ± 1,72i — la dernière paire a partie réelle positive.

Simplifier les polynomes complexes en utilisant le critère

La méthode Routh-Hurwitz simplifie l'analyse de stabilité des polynômes complexes de plusieurs façons :

Pour les polynômes à ordre élevé, l'affacturage ou la recherche numérique de racines est coûteuse et parfois peu fiable. Le tableau ne nécessite que des opérations d'addition, de soustraction, de multiplication et de division, facilement effectuées à la main ou à la feuille de calcul.

Incertitude des paramètres de poignées En exprimant des coefficients en termes de paramètres du système (gains, constantes de temps), on peut déduire des conditions de stabilité en termes d'inégalités. Par exemple, la plage de gains d'un contrôleur K pour laquelle le système de boucle fermée reste stable peut être trouvée en appliquant le tableau de Routh et en exigeant que tous les éléments de première colonne soient positifs.

Identifie la stabilité marginale. Un zéro dans la première colonne (le cas ε) ou une rangée de zéros indique les racines sur l'axe imaginaire ou symétriquement placé sur l'origine.Ces conditions définissent souvent la limite entre le comportement stable et instable.

Travailler avec des Coefficients complexes?

Pour les polynômes avec des coefficients complexes (a_i -C), des modifications telles que le [Michailov critère peuvent être employées. Cependant, dans la plupart des problèmes pratiques d'ingénierie de contrôle, l'équation caractéristique provient de paramètres physiques évalués en valeur réelle, de sorte que les coefficients réels sont la norme.

Applications dans les systèmes de contrôle et au-delà

La principale application du critère de stabilité de Routh–Hurwitz consiste à déterminer la stabilité des systèmes linéaires invariants du temps (LTI) décrits par les fonctions de transfert. Le dénominateur de la fonction de transfert est le polynôme caractéristique, et la stabilité exige que toutes les racines aient des parties réelles négatives.

Conception des systèmes de contrôle

Pendant la conception du contrôleur, les ingénieurs ajustent souvent les gains (p. ex., un contrôleur PID) proportionnel, intégral et dérivé pour obtenir les performances souhaitées. Le critère de Routh permet de trouver rapidement les marges gain—la gamme de gains qui maintiennent le système stable. Par exemple, étant donné un système de rétroaction d'unité avec fonction de transfert vers l'avant G(s) = K / [s+1)(s+2)(s+3)]], le polynôme caractéristique est s³ + 6s[² + 11s + (6+K) = 0. L'application du tableau de Routh donne la condition ]K < 60 pour la stabilité.

Analyse de réseau

En génie électrique, le polynôme caractéristique d'un circuit RLC détermine si les oscillations transitoires se désintègrent. Le critère Routh-Hurwitz aide les ingénieurs à choisir les valeurs des composants pour éviter les oscillations durables (instabilité).

Contrôle des procédés industriels

Dans les industries chimiques et pétrolières, les régulateurs de processus régulent la température, la pression et le débit. Le réseau de Routh est utilisé lors de l'accordage des régulateurs proportionnels – intégraux – dérivés (PID) pour s'assurer que l'usine en boucle fermée demeure stable dans des conditions de fonctionnement variables.

Avantages et limites

Comme tout outil, le critère Routh-Hurwitz a des forces et des faiblesses.

Forces

Aucun calcul de racine n'est nécessaire Le tableau fonctionne uniquement sur les coefficients.
Évoluable. La méthode fonctionne pour les polynômes de tout ordre, bien que les commandes plus élevées deviennent fastidieuses à la main (facilement automatisée).
Capacité symbolique Les coefficients peuvent être des expressions algébriques, permettant une analyse de stabilité paramétrique.
Fournit le nombre de racines instables. Chaque changement de signe dans la première colonne correspond à une racine avec une partie réelle positive.

Limitations

Les coefficients réels seulement. La forme classique ne traite pas directement les polynômes avec des coefficients complexes, bien que des extensions existent.
Aucune information sur les emplacements des racines. Le critère indique seulement combien de racines instables existent, et non leurs valeurs exactes.
Sensibilité aux erreurs numériques. Lorsque les coefficients sont presque nuls ou contiennent des erreurs d'arrondi, la construction du tableau peut produire des ambiguïtés. L'utilisation de la méthode ε est un moyen de contourner les erreurs, mais nécessite une interprétation minutieuse.
Ne gère pas les systèmes de retard Les systèmes avec retard de temps ont des équations caractéristiques transcendantales, et non des polynômes. Le critère de Routh–Hurwitz n'est pas directement applicable; au contraire, des méthodes comme les approximations Padé ou le prédicteur Smith sont utilisées.

Contexte historique et lectures complémentaires

Le critère Routh–Hurwitz est nommé d'après Edward John Routh, mathématicien anglais qui a publié la méthode de tableau en 1876, et Adolf Hurwitz, mathématicien allemand qui a dérivé indépendamment une formulation de matrice similaire en 1895. Leur travail s'est basé sur des contributions antérieures de Cauchy et Hermite. Le critère reste un élément de base dans les manuels de contrôle modernes—voir l'article Wikipedia pour un aperçu. Pour un traitement mathématique plus profond, ScienceDirect=s topic page fournit des références aux sources originales.

Les ingénieurs qui étudient le contrôle avancé rencontreront également la matrice Hurwitz dans le contexte de la stabilité de Lyapunov. La connexion entre le tableau de Routh et le déterminant de la matrice Hurwitz est un résultat classique en algèbre linéaire et théorie du contrôle.

Exemple illustré étape par étape

Laissez-les passer par un exemple complet avec un polynôme complexe (de haut ordre) pour démontrer la puissance de la méthode. Considérez le polynôme caractéristique d'un système de contrôle :

P(s) = 2 s⁵ + 4 s⁴ + 6 s³ + 8 s² + 10 s + 12

Étape 1: Disposer les deux premières lignes

s⁵ (puissance 5): coefficients de s⁵, s³, s^]1[ → 2, 6, 10
s^4[(puissance 4): coefficients de s⁴, s², s^]0 → 4, 8, 12

Étape 2: Calculer la troisième ligne (s³

Utilisation de la formule déterminante:

Premier élément: (4·6 - 2·8)/4 = (24 - 16)/4 = 8/4 = 2
Deuxième élément: (4·10 - 2·12)/4 = (40 - 24)/4 = 16/4 = 4
Troisième élément: zéro (puisqu'il n'y a plus de coefficients dans la première ligne)

s³: 2, 4, 0

Étape 3: Calculer la quatrième ligne (s²)

s⁴[,8,12) et s³ (2,4,0):

Premier élément: (2·8 - 4·4)/2 = (16 - 16)/2 = 0/2 = 0
Deuxième élément: (2·12 - 4·0)/2 = (24 - 0)/2 = 12

s²: 0, 12

Ici, nous avons un zéro dans la première colonne. Le reste de la ligne est non zéro (12). Remplacer le zéro par ε (un petit nombre positif).

Étape 4: Calculer la cinquième rangée (s¹)

s³[2,4,0]ss² (ε,12):

Premier élément: (ε4 - 2·12)/ε = (4ε - 24)/ε = 4 - 24/ε
Deuxième élément: (ε0 - 2·0)/ε = 0

s¹: (4 - 24/ε), 0

Comme ε → 0⁺, le terme -24/ε domine, rendant cet élément négatif (−--).

Étape 5: Calculer la sixième ligne (s⁰)

En utilisant la ligne s² (ε,12) et la ligne s¹ (4 - 24/ε, 0):

Premier élément: [(4 - 24/ε)·12 - 0/ε] / (4 - 24/ε) = 12

s⁰: 12

Étape 6: Analyser la première colonne

Les premières entrées de colonne sont : rangées s[^5[: 2 (positif), rangées s⁴: 4 (positif), rangées s[³: 2 (positif), rangées s²: ε (positif, comme ε>0), rangées s[¹: 4 - 24/ε (négatif pour les petits ε), rangées s⁰: 12 (positif).

Il y a deux changements de signe : de positif (s²) à négatif (s[¹), et de négatif de retour à positif (s⁰. Par conséquent, le polynôme a deux racines avec des parties réelles positives.

Vérification:[ À l'aide d'outils numériques, les racines de 2 s[⁵ + ... + 12 = 0 sont approximativement: -1,09, -0,50 ± 0,87i, 0,79 ± 1,04i. En effet, les deux racines conjuguées complexes ont une partie réelle positive (0,79), confirmant la prédiction de Routh‐Hurwitz.

Approche alternative: le polynôme auxiliaire pour le cas de zéro ligne

Si nous avions rencontré une rangée complète de zéros au lieu de juste un zéro dans la première colonne, nous formerions un polynôme auxiliaire à partir de la rangée ci-dessus. Cette technique est essentielle pour des cas comme s⁴ + 5s² + 4 (qui a des racines purement imaginaires ±j, ±2j). La rangée zéro indique des racines symétriques sur l'origine, et le polynôme auxiliaire (découvert de la rangée ci-dessus) donne les racines réelles.

Conclusion

Le critère de stabilité de Routh-Hurwitz transforme l'analyse complexe des polynômes à haut ordre en un algorithme simple. En construisant le tableau de Routh et en inspectant les signes de la première colonne, les ingénieurs et mathématiciens peuvent déterminer la stabilité du système sans résoudre pour les racines. La méthode est indispensable dans la conception du système de contrôle, l'analyse des circuits et tout domaine où les racines d'un vrai polynôme dictent le comportement.

La maîtrise de cette technique permet de gérer avec confiance des polynômes complexes, que ce soit pour des études universitaires, des études industrielles ou des recherches. Combiné à des outils informatiques modernes (p. ex. MATLAB, Python), le critère Routh-Hurwitz demeure un élément fondamental de l'analyse dynamique du système. Pour ceux qui souhaitent explorer davantage, le blanc-papier Control Systems Principles fournit des exemples travaillés et des extensions avancées supplémentaires.