Polynôme de Hurwitz - Définition

Source: Wikipédia sous licence CC-BY-SA 3.0.
La liste des auteurs de cet article est disponible ici.

Introduction

Un polynôme à coefficients réels, de degré n, est appelé polynôme de Hurwitz si ses n racines complexes sont toutes à partie réelle négative.

On montre aisément que nécessairement, tous les coefficients du polynôme sont de même signe, et donc sans restreindre la généralité, on prend le terme constant T0 égal à 1 en physique.

Résultats pour n<5

n=2 : pas d'autre condition (évident pour un physicien habitué au circuit RLC).

n=3 : rajouter : ET T2.T1T3.T0 > 0

n=4 : rajouter : ET  T_2.T_1 - T_3.T_0 > T_4.T_1^2/T_3

Remarque : pour un physicien, les coefficientsTn sont des constantes de temps à la puissance n. Il pourra aisément contrôler l'homogénéité des formules.

Le cas n = 3 se retient aisément par des considérations énergétiques (voir l'article filtre électrique linéaire).

Le cas n=4 se retient mnémotechniquement aisément grâce au cas n=3 : effectivement si T4 est quasi-nul, on devait retrouver n=3 . Par contre, il faut connaître la situation de 2 oscillateurs couplés avec injection de puissance, pour intuiter aisément l'inéquation T2.T1 / T3.T0 > 1 + k.(T4 / T0).(T1 / T3)2, k=1 .

Au-delà, il vaut mieux traiter la question par l'algorithme des tableaux de Routh. Néanmoins, une remarque de pertinence : si P(x) est de Hurwitz, x^n .P(1/x) est de Hurwitz : les conditions sont donc « symétriques » sur les coefficients.

Par exemple, le cas n= 5 se retient assez aisément: il y a deux inégalités:

A: = (T4.T3T5.T2) > 0

évidemment et

B2A.C < 0

avec B: = (T4.T1T5.T0)

et

C: = (T2.T1T3.T0)

Cas général : les tableaux de Routh et critère de Routh-Hurwitz

Le tableau de Routh est un critère de stabilité utilisé en automatique. Il permet de déterminer, au regard du dénominateur de la fonction de transfert en boucle fermée, si le système bouclé, représenté par cette fonction, est stable ou non (cf stabilité EBSB).

Soit s la variable de la fonction de transfert. Le dénominateur s'écrit sous la forme d'un polynôme en s :

ansn + an − 1sn − 1 + ... + a0

  • Condition nécessaire :
    • Tous les coefficients ai doivent être non nuls
    • Tous les coefficients ai doivent être de même signe

Rem : autre manière de dire : ao = 1 et tous les coefficients strictement positifs.

Si ces 2 conditions ne sont pas réunies, le système est instable.

Pour déterminer en plus la condition suffisante, il nous faut dresser le tableau de Routh.

La 1e ligne du tableau, indexée par sn, comporte les coefficients an, an − 2,....

La 2e ligne du tableau, indexée par sn − 1, comporte les coefficients an − 1, an − 3,...

Puis pour les autres lignes, on complète de la manière suivante :

1e colonne de la 3e ligne : b_{n-2}=\frac{-1}{a_{n-1}}\begin{vmatrix}   a_{n} & a_{n-2} \\   a_{n-1} & a_{n-3}  \end{vmatrix}

2e colonne de la 3e ligne : b_{n-4}=\frac{-1}{a_{n-1}}\begin{vmatrix}   a_{n} & a_{n-4} \\   a_{n-1} & a_{n-5}  \end{vmatrix}

ect.... (si on ne peut pas compléter une case, on y met un zéro)

ce qui nous donne le tableau suivant :

Tableau de Routh
sn
an
an − 2
an − 4
sn − 1
an − 1
an − 3
an − 5
sn − 2 b_{n-2}=\frac{-1}{a_{n-1}}\begin{vmatrix}   a_{n} & a_{n-2} \\   a_{n-1} & a_{n-3}  \end{vmatrix} b_{n-4}=\frac{-1}{a_{n-1}}\begin{vmatrix}   a_{n} & a_{n-4} \\   a_{n-1} & a_{n-5}  \end{vmatrix}
...
sn − 3
\frac{-1}{b_{n-2}}\begin{vmatrix}   a_{n-1} & a_{n-3} \\   b_{n-2} & b_{n-4}  \end{vmatrix}
...
...
  • Condition nécessaire et suffisante :

Si tous les coefficients de la 1re colonne du tableau sont tous non nuls et de même signe, alors le système bouclé est stable.


Remarque : la condition nécessaire et suffisante de stabilité pour un système du 2e ordre est que tous les coefficents soient de même signe (on obtient simplement ce résultat en appliquant le critère de Routh).


Exemple : af''(t) + bf(t) = 0 avec a > 0


si b > 0 alors les coefficents sont de même signe donc le système est stable. En effet, la solution f(t) est :  f(t) = A\cos{(\sqrt{\frac{b}{a}}t+\phi)}

f(t) est alors bornée comme l'entrée (ici 0) donc le système est stable.


si b < 0, le système est instable. En effet, la solution est :  f(t) = A\exp{\sqrt{\frac{b}{a}}t} + B\exp{(-\sqrt{\frac{b}{a}}t)}

qui diverge lorsque t tend vers l'infini.

Page générée en 0.138 seconde(s) - site hébergé chez Contabo
Ce site fait l'objet d'une déclaration à la CNIL sous le numéro de dossier 1037632
A propos - Informations légales
Version anglaise | Version allemande | Version espagnole | Version portugaise