Le critère de stabilité de Routh-Hurwitz est une pierre angulaire de la théorie classique du contrôle, fournissant une méthode algébrique directe pour déterminer si un système linéaire invariant dans le temps est stable sans résolution pour ses pôles. Pour les systèmes d'ordre entier – où le polynôme caractéristique n'a que des pouvoirs entiers de la variable Laplace ss]ssssssssssssssssssss[FLT]s]s[FLT]s]s[FLT:]s[FLT]s[F]]s[F][F][F]

Comprendre le critère classique de la Routh–Hurwitz

Avant d'explorer les extensions, il est essentiel de recapter la méthode classique.

Le critère de Routh–Hurwitz construit une table (la matrice de Routh) à partir des coefficients. La condition de stabilité est que chaque élément de la première colonne du tableau doit avoir le même signe (généralement positif). Si un changement de signe survient, le nombre de changements de signe est égal au nombre de pôles de demi-plan droit. La beauté du test réside dans sa nature purement algébrique : aucune recherche de racine, aucune analyse complexe.

Le critère s'étend également aux systèmes avec délais (par l'approximation Padé) et aux systèmes à temps discret (en utilisant la transformation bilinéaire), mais toutes ces extensions dépendent du cadre polynôme.

Systèmes de commande fractionnée : un bref aperçu

Un système à ordre fractionnaire est celui dont la dynamique est décrite par une équation différentielle fractionnelle de la forme

où les ordres de différenciation αi et βi[ ne sont pas nécessairement des entiers. La transformation Laplace d'une telle équation (dans des conditions initiales nulles) donne une fonction de transfert contenant des termes comme sγ avec γ fractionnel. Par exemple, un intégrateur à ordre fractionnel simple a la fonction de transfert 1/sα] avec 0 < α < 2. Ces systèmes présentent des propriétés non vues dans les systèmes d'ordre entier : représentation de l'espace à dimension infinie, relaxation non exponentielle et mémoire à longue portée.

Le calcul fractionnel a trouvé des applications dans le contrôle des structures flexibles, la modélisation des tissus biologiques, la caractérisation de l'impédance de la batterie, et même les séries temporelles financières. La nécessité de vérifier la stabilité de tels systèmes est aiguë, mais la table classique Routh-Hurwitz ne peut pas être construite parce que l'équation caractéristique n'est pas polynôme.

Pourquoi les Fails de test classiques Routh–Hurwitz pour les systèmes fractionnels

L'obstacle est fondamental: une équation caractéristique d'ordre fractionnaire, comme

Les pouvoirs de s ne sont pas des multiples entiers les uns des autres. Le tableau de Routh exige des pouvoirs entiers afin qu'une élimination systématique puisse être effectuée. De plus, les pouvoirs fractionnels conduisent à des coupes de branches dans le plan complexe; le concept de -pieu droit-demi-plan , devient plus subtil parce qu'une puissance fractionnelle peut avoir plusieurs feuilles Riemann. Un système fractionnel peut avoir un nombre infini de pôles (si l'ordre est irrationnel) ou un nombre fini si les ordres sont proportionnés.

Comme le critère de Routh-Hurwitz est essentiellement un test algébrique pour la propriété Hurwitz d'un polynôme, il ne peut pas être appliqué directement à une expression non polynôme. Il faut plutôt modifier l'esprit du test, en déterminant la stabilité en examinant une condition algébrique transformée, tout en tenant compte de la nature fractionnelle.

Le critère de stabilité Matignon : un remplacement algébrique direct

L'extension la plus connue est le critère de stabilité Matignon (également appelé théorème Matignon). Il fournit une condition nécessaire et suffisante pour une grande classe de systèmes à ordre fractionnel avec des ordres proportionnés. Un système fractionnel est dit avoir des ordres proportionnés si toutes les puissances fractionnelles de l'équation caractéristique sont des multiples entiers d'un ordre de base α (par exemple α, 2α, 3α, ...). De nombreux systèmes pratiques, en particulier ceux dérivés de contrôleurs PID fractionnels ou de modèles viscoélastiques fractionnels, sont proportionnés.

Critère de la matière:[ Pour un système à ordre fractionnaire avec équation caractéristique

où α est un nombre réel positif, le système est stable si et seulement si

απ/2 pour chaque racine i de P.

En d'autres termes, les racines doivent se trouver dans une région du plan complexe à l'extérieur d'un secteur d'angle απ/2 centré sur l'axe réel négatif. Cette condition réduit à l'état classique de gauche-demi-plan quand α = 1 (depuis lors απ/2 = π/2), mais pour α < 1 la région de stabilité est un coin.

Le critère de Matignon est puissant car il transforme le problème de stabilité en un problème de localisation racine pour le polynôme transformé Pw = sα. On peut appliquer le test Routh-Hurwitz au polynôme P(w)w = s[α], mais la condition sur l'argument de sn'est pas capturée par le tableau classique de Routh seulement. Cependant, pour des systèmes correspondants, la condition de stabilité peut être vérifiée en construisant d'abord le tableau Routh pour le polynôme [P(w]]][et en vérifiant ensuite l'état de l'argument.

Cela conduit à une procédure pratique:

  1. Identifiez l'ordre fractionnaire de base α de sorte que l'équation caractéristique devienne polynôme dans sα.
  2. Définit w = sα] et réécrit l'équation caractéristique comme P(w) = 0, un polynôme conventionnel.
  3. Appliquer le test classique Routh–Hurwitz à P(w). Si des changements de signe se produisent dans la première colonne, le système est instable sous la définition polynôme classique dans w.
  4. Vérifiez la condition de l'argument : Pour chaque racine wi, calculez les si = w[i[1/α][et vérifiez que]i]i]]]]i[FLT: απ/2. Cette étape peut nécessiter un calcul numérique parce que les racines w]i sont généralement complexes.

Le critère Matignon est souvent utilisé comme outil principal pour l'analyse de stabilité de systèmes à ordre fractionnel proportionnels. Pour les systèmes non-commensés (où les ordres fractionnels ne sont pas multiples d'une base commune), il faut utiliser d'autres méthodes comme l'approche par domaine de fréquence ou l'inversion numérique.

Autres approches analytiques

La fonction Mittag–Leffler et la représentation de l'espace public

Les systèmes d'ordre fractionnel peuvent être représentés sous une forme d'espace pseudo-étatique en utilisant la dérivée fractionnelle de l'ordre α. La condition de stabilité devient une contrainte sur les valeurs propres de la matrice du système.

le système est stable si et seulement si

où λ(A) sont les valeurs propres de la matrice A. C'est l'analogue multivariable direct du critère Matignon. Pour les systèmes correspondants, cette condition se réduit à la même région de stabilité en forme de coin. Les ingénieurs peuvent calculer les valeurs propres et vérifier leurs arguments en utilisant l'algèbre linéaire numérique standard.

Transformation en un système entier de commande par substitution variable

Une autre voie analytique consiste à introduire une nouvelle variable z = sα et ensuite à appliquer le test classique Routh–Hurwitz au polynôme transformé. Cependant, cette approche ne vérifie que la partie -polynomiale de la stabilité. Une condition nécessaire est que le polynôme dans z a toutes les racines dans le demi-plan gauche, mais la condition suffisante nécessite le test d'argument supplémentaire. En pratique, les ingénieurs effectuent souvent les deux : premier essai du tableau Routh sur P(z)[ pour détecter l'instabilité polynôme évidente; puis, pour les candidats stables restants, appliquer la condition d'argument.

Techniques de fréquence-domaine

La stabilité peut également être évaluée par les critères Nyquist ou Bode adaptés aux fonctions de transfert d'ordre fractionnel. La fonction de transfert de boucle ouverte G(s) d'un système d'ordre fractionnel est une fonction non-rationnelle, mais sa réponse de fréquence peut être évaluée numériquement. Le critère de stabilité Nyquist tient toujours : le nombre d'encerclements du point −1 détermine la stabilité de boucle fermée.

Outils numériques et informatiques

Compte tenu de la complexité des systèmes d'ordre fractionnaire, les pratiques d'ingénierie reposent fortement sur les outils logiciels.La boîte à outils FOMCON (Factional-Order Modeling and Control) pour MATLAB est un choix populaire. Elle fournit des fonctions pour calculer les contrôles de stabilité en utilisant le critère Matignon, pour construire des approximations (ex.: Oustaloup recursive approximation) des opérateurs fractionnels, et pour effectuer des simulations de temps-domaine.

Les autres outils sont les suivants :

  • CRONE toolbox (Commande Robuste d-Ordre Non Entre) : développé à l'Université de Bordeaux, il soutient la conception de contrôleurs de commande fractionnaire et l'analyse de stabilité.
  • Les bibliothèques de python: fractionnel[ (pour le calcul fractionnel) et le contrôle (pour l'analyse générale du système de contrôle) peuvent être utilisés ensemble pour mettre en œuvre le critère Matignon.
  • Le calcul symbolique dans Mathematica ou Maple peut gérer des manipulations algébriques de pouvoirs fractionnaires et de recherche de racines dans le plan complexe.

Pour un contrôle rapide, de nombreux ingénieurs utilisent l'algorithme numérique suivant:

  1. Définir l'équation caractéristique comme une expression symbolique dans ss avec des pouvoirs fractionnaires.
  2. Convertir en un polynôme dans w = sα si le système est proportionnel.
  3. Calculer les racines du polynôme dans w en utilisant un résolveur numérique.
  4. Pour chaque racine, calculer les s correspondants et leur argument.
  5. Vérifiez si απ/2arg(s) απ/2 pour toutes les racines.

Cet algorithme peut être mis en œuvre en quelques lignes de code et est robuste pour les systèmes les plus appropriés. Pour les systèmes non-commensés, il faut s'appuyer sur des méthodes de domaine de fréquence ou sur une évaluation numérique par grille de l'équation caractéristique dans le plan complexe.

Procédure pratique étape par étape pour les ingénieurs

Voici un workflow pragmatique qui combine l'esprit classique Routh-Hurwitz avec les modifications fractionnelles nécessaires.

  1. Obligation de l'équation caractéristique de la fonction de transfert du système et du numéro 8217;s. Pour un système à ordre fractionnaire, cela implique des termes comme sγ avec γ non entier.
  2. Choisir la mesurabilité Déterminer si tous les ordres fractionnels sont des multiples entiers d'une base commune α. Par exemple, les ordres 0,5, 1,5, 2,0 sont des multiples de 0,5; les ordres 0,7 et 1,2 ne sont pas des multiples d'une base unique (sauf si la base est 0,1, mais cela est rarement rationnel). Si le système est proportionnel, procéder à la transformation polynôme. Sinon, passer à l'étape 6 (méthode de fréquence-domaine).
  3. [P(w) = s[P(w) = 0. Par exemple, 2,5+3 s]1,2+2[ne sont pas des multiples d'une base commune], de sorte que cette étape échouerait.
  4. Appliquer le test classique Routh–Hurwitz à P(w) Construire le tableau Routh. Si un changement de signe survient dans la première colonne, le système est instable (sous l'interprétation polynomiale). Si toutes les entrées de première colonne ont le même signe, passer à l'étape suivante.
  5. Vérifier la condition de l'argument. Calculer les racines de P(w)[. Pour chaque racine w[isi]i[i[1/α]. La puissance complexe est multivalorisée; utiliser la valeur principale.
  6. Pour les systèmes non-commensés, utiliser une approche par domaine de fréquence. Évaluer l'équation caractéristique le long de la limite de la région de stabilité, qui pour les systèmes fractionnels est un secteur. La courbe Nyquist de la fonction de transfert en boucle ouverte peut indiquer la stabilité en boucle fermée. Ou utiliser la recherche numérique de la racine dans le s]‐plan en discrétisant l'opérateur fractionnaire, par exemple en utilisant l'approximation de Oustaloup. Ensuite, appliquer le test classique Routh–Hurwitz au système d'ordre entier approximatif. Cette méthode donne un résultat approximatif et doit être validée par simulation.

De nombreux ingénieurs estiment que le cas correspondant couvre une grande partie des modèles pratiques de contrôle fractionnel, en particulier les contrôleurs PID fractionnels (où les ordres dérivés et intégrés sont multiples d'un α commun, souvent 0,5 ou 0,1). La procédure ci-dessus est donc largement applicable.

Étude de cas illustrative: Contrôleur de l'IDP de commande fractionnelle

Considérer un contrôleur PID à ordre fractionnaire (PI[αD[β) contrôlant une usine de premier ordre.

Ce système n'est pas comparable parce que les ordres fractionnels 2.5, 1.2, 0.1 ne sont pas des multiples entiers d'une base commune. Nous ne pouvons pas appliquer directement le critère Matignon. Nous approchons plutôt les opérateurs fractionnels utilisant un filtre d'ordre Oustaloup, par exemple 5 sur une plage de fréquences [0.01, 100] rad/s. Cela donne une fonction de transfert d'ordre entier de haut niveau. Le polynôme caractéristique approximatif sera de l'ordre 5×5 = 25 pour chaque terme fractionnel, mais nous pouvons réduire l'ordre en utilisant une troncation équilibrée. Ensuite, le test classique Routh–Hurwitz peut être appliqué au polynôme réduit. Un contrôle de stabilité par simulation confirme le résultat. Cette méthode approximative est la norme de fait dans l'industrie lorsque les outils analytiques ne sont pas disponibles.

Par contre, si le contrôleur est conçu avec des ordres proportionnés (par exemple α = 0,5, donc les ordres sont de 0,5, 1,0, 1,5, 2,0), le système devient proportionné. Ensuite, nous pouvons utiliser le critère Matignon directement, obtenant une condition exacte nécessaire et suffisante (modulo l'interprétation de la coupe de branche).

Limites et recherche continue

L'extension de l'analyse de Routh-Hurwitz aux systèmes à ordre fractionnel n'est pas une formule unique, mais une collection de techniques, chacune ayant son propre domaine d'applicabilité. Le critère Matignon est exact uniquement pour les systèmes à proportion. Pour les systèmes à ordre non-commensuré, l'approche à ordre entier approximatif peut introduire des erreurs, surtout si le système a légèrement amorti les modes à haute fréquence. Les chercheurs explorent des critères algébriques directs pour les systèmes à ordre non-commensuré en utilisant le concept de réseaux à ordre fractionnel (p. ex., en utilisant la transformation de variables à ordre fractionnel s = (z)1/α et en appliquant un réseau modifié).

De plus, la région de stabilité des systèmes fractionnels avec des ordres irrationnels n'est pas un coin simple; elle peut devenir une région plus compliquée dans la surface de Riemann. Les ingénieurs pratiques s'appuient généralement sur l'évaluation numérique de l'équation caractéristique de magnitude et de phase le long de la limite de stabilité, qui est calculable intensive mais donne des résultats fiables pour des valeurs de paramètres spécifiques.

Conclusion

L'extension de l'analyse de Routh-Hurwitz aux systèmes à ordre fractionnel nécessite de dépasser le cadre polynôme classique. Le critère de stabilité Matignon pour les systèmes à proportion retransmet le problème en une condition de base qui peut être vérifiée à l'aide du tableau classique de Routh et de la recherche numérique de racine. Pour les systèmes non-commérés, les méthodes de domaine de fréquence et les approximations numériques restent les outils primaires.

Pour plus de détails, les lecteurs peuvent consulter l'entrée Wikipedia sur le critère Routh–Hurwitz et le résumé du calcul fractionnel. L'œuvre originale de Matignon est détaillée dans ce document (lien DOI).