Introduction à l'analyse de stabilité dans les systèmes discrets

L'analyse de stabilité constitue l'épine dorsale de l'ingénierie des systèmes de contrôle, en veillant à ce qu'un système réponde de façon prévisible et demeure limité dans des conditions normales d'exploitation. Pour les systèmes à temps discret et no 8212;communs dans les applications de contrôle numérique, de traitement des signaux et d'intégration et no 8212;le critère de stabilité passe du plan s continu à la moitié gauche du plan à l'unité du plan z. Les ingénieurs comptent sur des outils analytiques robustes pour déterminer si tous les pôles d'une fonction de transfert discrète résident à l'intérieur du cercle de l'unité sans résoudre explicitement les polynômes à haut ordre.

Cet article présente une comparaison complète de l'essai de stabilité de Routh-Hurwitz et Jury’, en explorant leurs fondements mathématiques, les étapes de procédure, les adaptations pratiques et les forces relatives.

Stabilité du système discret : le critère du cercle unitaire

Avant de plonger dans les méthodes spécifiques, il est essentiel de rappeler la condition de stabilité fondamentale des systèmes à temps discret. Un système discret linéaire invariant dans le temps décrit par la fonction de transfert H(z) = N(z) / D(z) est stable si et seulement si tous les pôles et #8212;les racines du polynôme caractéristique D(z)[ = 0—s'insèrent dans le cercle de l'unité dans le plan en z, c'est-à-dire:=======================================================================================

Les algorithmes de recherche de racines peuvent calculer directement les emplacements des pôles, mais ils deviennent encombrants pour les polynômes à ordre élevé. Par conséquent, les critères algébriques comme Roth-Hurwitz et Jury’s test offrent des vérifications efficaces et non-itératives. Le défi clé est que Jury’s test fonctionne directement sur les coefficients de D(z), alors que Roth-Hurwitz travaille généralement sur les polynômes dans le domaine s. Une transformation est donc nécessaire pour appliquer Roth-Hurwitz à des systèmes discrets.

Critère de Routh-Hurwitz et son adaptation aux systèmes discrets

Routh-Hurwitz classique pour systèmes continus

Le critère de Routh-Hurwitz fournit des conditions nécessaires et suffisantes pour la stabilité des systèmes linéaires continus.

a0sn +a]1s[n–1+]an] = 0

La méthode construit un tableau de Routh à partir des coefficients. Le système est stable si tous les éléments de la première colonne du tableau sont positifs et non nuls. Un changement de signe indique les racines dans le plan de droite et donc l'instabilité. Cette approche évite le calcul explicite de la racine et est simple pour les polynômes jusqu'à des ordres modérés.

Transformation bilinéaire : Conversion du plan z en plan s

Pour appliquer Routh-Hurwitz aux systèmes discrets, les ingénieurs utilisent la transformation bilinéaire (également appelée transformation Tustin’s). Cette cartographie relie la variable z à une nouvelle variable complexe w (souvent désignée comme s pour plus de commodité après la transformation):

z = (1 + wT[/2) / (1 – wT[/2)

ou, pour une période d'échantillonnage normalisée T[ = 2, la forme la plus simple z = (1 + w / (1 – w[). Sur cette carte bilinéaire, l'intérieur du cercle de l'unité dans z-plan plan de la moitié gauche du plan w-plan. Par conséquent, le polynôme transformé dans w peut être analysé à l'aide du tableau standard Routh-Hurwitz.

Étapes pour appliquer Routh-Hurwitz aux systèmes discrets

  1. Observer le polynôme caractéristique discret D(z) du système et de la fonction de transfert #8217;s.
  2. Soustituer la transformation bilinéaire en D(z) pour obtenir un polynôme dans w[. Cette étape peut impliquer une expansion et une simplification algébriques, en particulier pour les systèmes à ordre supérieur.
  3. Appliquer le critère de Routh-Hurwitz au polynôme transformé dans w[.Construire le tableau de Routh et examiner les signes des entrées de première colonne.
  4. Conclusion de stabilité :[ Si toutes les entrées de première colonne sont positives, le système discret est stable (tous les pôles à l'intérieur du cercle de l'unité).

Bien que cette approche fonctionne, elle introduit une complexité algébrique supplémentaire. La transformation bilinéaire peut augmenter significativement l'ordre polynôme si elle n'est pas manipulée avec soin. De plus, la transformation elle-même peut devenir une source d'erreur numérique, en particulier pour les systèmes à haut ordre ou lorsque les périodes d'échantillonnage ne sont pas normalisées correctement.

Limitations de l'adaptation Routh-Hurwitz

  • Fonctionnement algébrique: La transformation bilinéaire élargit le polynôme, introduisant potentiellement des erreurs d'arrondi et rendant les calculs manuels fastidieux.
  • Sensibilité numérique : Pour les polynômes à ordre élevé, les grandeurs des coefficients peuvent varier considérablement, ce qui conduit à des tableaux mal conditionnés.
  • Cas particuliers: La méthode exige la manipulation de polynômes auxiliaires pour les racines sur l'axe imaginaire (dans le w[-plan) qui correspondent à des pôles sur le cercle unitaire – une condition de stabilité marginale souvent non désirée.
  • Inefficacité pour les systèmes purement discrets:[ Comme la méthode n'est pas native du domaine discret, de nombreux manuels de contrôle recommandent Jury’s test comme l'alternative plus directe et moins sujette aux erreurs.

Néanmoins, le critère de Routh-Hurwitz est parfois enseigné dans les cours de contrôle comme un pont entre l'analyse continue et discrète, et il peut être appliqué à l'aide d'outils de calcul symboliques pour atténuer les difficultés algébriques.

Test de stabilité du jury: Conçu pour le domaine discret

Origines et principe

Le jury a introduit dans les années 1960 un test de stabilité du jury, qui est la contrepartie discrète du critère de Routh-Hurwitz. Il examine directement les coefficients du polynôme caractéristique D(z) sans exiger de transformation de domaine. Au lieu de vérifier le signe des entrées de première colonne, le jury et le jury ont construit une table (la table du jury) à partir des coefficients polynômes et applique un ensemble de conditions qui garantissent que toutes les racines ont une magnitude inférieure à un.

Considérez le polynôme caractéristique de l'ordre des moniques n:

D(z) = a0zn + a1zn–1 + … + an–1z + an = 0, with a0 > 0.

Le test Jury’s fournit une procédure systématique pour déterminer si toutes les racines se trouvent à l'intérieur du cercle unitaire en formant une table avec n + 1 lignes. Les lignes sont générées séquentiellement à l'aide de calculs déterminant impliquant les coefficients.

Procédure de construction du tableau du jury

  1. Listez les coefficients[ dans la première rangée (ligne 0) dans l'ordre naturel: a0, a1, ..., a]n].
  2. Former la deuxième ligne (ligne 1)[ en inversant l'ordre des coefficients de la ligne 0: a[n, an–1, ..., a]0.
  3. [f][f][f.[f.[
  4. Continuer jusqu'à ce qu'un seul élément non zéro reste dans la première colonne, ou jusqu'à ce que la table s'effondre. La dernière ligne fournit la dernière condition.
  5. Vérifier les conditions:[ Pour la stabilité, les éléments suivants doivent être maintenus:
    • D(1) > 0
    • D[–1) > 0
    • Pour le tableau du jury, le premier élément de chaque ligne impair (lignes 1, 3, 5, ...) doit être supérieur à zéro. Les conditions équivalentes impliquent les entrées principales des polynômes réduits.

Ces conditions sont nécessaires et suffisantes pour que toutes les racines puissent satisfaire à la condition suivante : 1. En pratique, le test est mis en œuvre dans les systèmes d'algèbres informatiques, mais le processus manuel est gérable pour les commandes jusqu'à environ 4 ou 5.

Exemple : Système de deuxième commande

D(z) = z[2 + a1z] + a]2] = 0. Appliquant le jury et #8217;s test:

  • Row 0: 1, a1, a2
  • Tige 1: a2, a1, 1
  • Row 2:[ La ligne suivante comporte des entrées b0 = 1 – a22, b1] = [a]1[ – a]1a]2[]][[[[

Conditions de stabilité pour le second ordre:

  1. D(1) = 1 + a1 + a2 > 0
  2. D(–1) = 1 – a1 + a2 > 0
  3. a2

Ce sont les conditions familières pour un système discret de deuxième ordre stable. Notez que le jury et #8217;s test produit directement ces inégalités simples sans aucune transformation.

Avantages du jury et du test no 8217;

  • Application directe: Fonctionne sur le z-domaine polynôme; aucune cartographie n'est requise.
  • L'ensemble complet des conditions:[ Fournit des critères de stabilité à la fois nécessaires et suffisants.
  • Poignées cas spéciaux:[ Comprend les vérifications des racines sur le cercle unitaire (stabilité marginale) par le biais des D[(1) et D[(–1) conditions.
  • Convient pour l'implémentation numérique:[ La construction de la table est algorithmique et facilement programmée.

Comparaison entre les deux essais de stabilité de Routh-Hurwitz et Jury

AspectRouth-Hurwitz (Adapted)Jury’s Stability Test
DomainContinuous (s-plane) originally; adapted via bilinear transform to w-planeDiscrete (z-plane) directly
Stability regionLeft-half plane (after transformation, maps to inside unit circle)Inside unit circle (|z|<1)
ProcedureConstruct Routh array from w-polynomial; check sign of first columnConstruct Jury table from z-polynomial; check series of inequalities
ComplexityModerate after transformation; transformation adds extra algebraic stepsLow to moderate; table construction straightforward for low to medium orders
Numerical stabilityCan suffer from coefficient scaling after bilinear transformGenerally good; determinants may cause swelling for high orders
Handling of marginal stabilityAuxiliary polynomial required; roots on imaginary axis in w-planeDirect checks through D(1) and D(-1) conditions
Educational accessibilityStudents often learn continuous version first; adaptation extends existing knowledgeRequires understanding of unit circle concept; but method is self-contained
Software implementationEasily coded via array construction; bilinear transform may require symbolic toolboxEasily coded via recursive table; native to z-domain
Typical use caseWhen system is already in s-domain or hybrid continuous/discreteWhen characteristic polynomial is directly from a discrete transfer function

Points de comparaison détaillés

1. Applicabilité et Fit Domaine

Le test de Jury’s est le choix naturel pour les systèmes discrets purs, en particulier dans les applications de commande numérique où la sortie du contrôleur et l'entrée de l'usine sont échantillonnées. Routh-Hurwitz nécessite une transformation qui ajoute de la complexité et peut masquer l'interprétation physique des pôles par rapport au cercle unitaire.

2. Complexité algébrique

Cependant, le test Jury’s évite l'expansion bilinéaire, qui peut produire des termes avec des puissances mixtes et de grands coefficients. Par exemple, un D(z)ssssssss[F.][F.[F.[F

3. Robustesse numérique

Le domaine de Routh-Hurwitz dans w- peut souffrir d'un mauvais conditionnement si la transformation bilinéaire crée des coefficients avec des grandeurs très différentes. Le test Jury’s est également confronté à des problèmes numériques pour les polynômes à haut ordre (p. ex., ordre > 10) en raison de soustractions répétées et de calculs déterminants, mais ceux-ci sont moins sévères que ceux introduits par la transformation.

4. Capacité de déterminer le nombre de pôles instables

Routh-Hurwitz, en vertu du tableau de Routh et des modifications de signe du numéro 8217; peut indiquer directement le nombre de pôles dans la moitié droite w-plan, qui correspond au nombre de pôles à l'extérieur du cercle unitaire. Jury et #8217;s test, sous sa forme de base, fournit une réponse oui/pas de stabilité. Cependant, les extensions du test du jury et #8217;s (par exemple, l'algorithme Jury-Gantmacher) peuvent aussi compter le nombre de racines à l'extérieur du cercle unitaire, bien que cela soit moins enseigné couramment.

5. Manipulation des coefficients polynômes spéciaux

Lorsqu'un système discret a des coefficients zéro (pouvoirs manquants), les deux méthodes exigent une prudence. Routh-Hurwitz traite des entrées zéro dans le tableau en utilisant de petites approximations d'épsilon. Jury’s test doit également gérer les zéros pendant la construction de la table, mais l'algorithme reste valide tant que le coefficient de tête a0 est non nul. Si une ligne de la table du jury devient entièrement zéro, le polynôme a des racines symétriques indiquant une stabilité marginale ou une antisymétrie.

Applications pratiques et quand choisir chaque méthode

Utiliser Jury , test quand

  • Vous travaillez exclusivement dans le domaine discret (z-plan).
  • Vous voulez un test en une seule étape sans transformation de domaine.
  • Vous analysez un contrôleur numérique ou un système de données échantillonnées avec une période d'échantillonnage connue.
  • Il faut des conditions claires pour la stabilité marginale (racines sur le cercle unitaire).
  • L'ordre polynôme est modéré (≤ 8 pour les vérifications manuelles; plus élevé pour les logiciels).

Utiliser Routh-Hurwitz (via la transformation bilinéaire)

  • Votre conception du système de contrôle est déjà dans le domaine continu et vous êtes discrétisant plus tard.
  • Vous préférez un cadre d'analyse unifié pour les systèmes continus et discrets.
  • Vous avez un logiciel (par exemple MATLAB, Python avec SymPy) pour gérer automatiquement la transformation bilinéaire.
  • Vous devez compter le nombre exact de pôles instables (par des changements de signe dans le tableau de Routh).
  • Vous enseignez le lien entre les critères de stabilité continus et discrets.

Exemples réels mondiaux

Exemple : Contrôleur numérique de vitesse du moteur Un microcontrôleur met en œuvre un contrôleur PI pour un moteur DC. Le polynôme caractéristique de boucle fermée est dérivé dans le domaine z. Un ingénieur utiliserait Jury’s test pour vérifier la stabilité pour divers taux d'échantillonnage et gains proportionnels. Le test révèle rapidement si les pôles se déplacent à l'extérieur du cercle de l'unité lorsque le gain augmente.

Exemple : Conception de filtre adaptatif Dans le traitement des signaux, les filtres IIR nécessitent des pôles à l'intérieur du cercle de l'unité pour éviter l'auto-oscillation. Jury’s test est généralement intégré dans les bibliothèques de conception de filtre pour valider la stabilité après quantification des coefficients.

Exemple : Système hybride à diffusion continue. Une installation continue contrôlée par un compensateur numérique nécessite une discrétisation via une cale à ordre zéro. L'équation caractéristique globale du domaine z peut être analysée soit par Jury’s test directement, soit par l'utilisation de la transformation bilinéaire pour appliquer Routh-Hurwitz. De nombreux ingénieurs de contrôle préfèrent Jury’s test parce qu'ils peuvent rester dans le domaine z et éviter d'autres transformations.

Mise en œuvre et automatisation des logiciels

Les deux méthodes sont disponibles dans les boîtes à outils du système de contrôle.

  • MATLAB: La fonction (ou scripts personnalisés) peut construire le tableau Routh, tandis que est disponible dans la Boîte à outils robotique ou peut être programmée. La fonction effectue la discrétisation et la transformation bilinéaire si nécessaire.
  • Python (SciPy et bibliothèques de contrôle): La bibliothèque comprend pour le test Jury’s. Pour Routh-Hurwitz, la fonction ou la construction manuelle de tableaux est utilisée.
  • Outils symboliques: SymPy et Mathematica permettent la manipulation algébrique pour la transformation bilinéaire et la génération de réseaux subséquente.

Lors de l'utilisation du logiciel, l'étape de transformation de Routh-Hurwitz est automatisée, réduisant ainsi les erreurs manuelles. Cependant, pour le prototypage rapide des systèmes discrets, le test Jury’s reste le choix le plus direct.

Limitations et réserves

Les deux méthodes ne sont pas toutes deux une panacée, mais elles deviennent peu fiables pour les polynômes à ordre élevé (ordre > 10) en raison de l'accumulation d'erreurs numériques. Dans de tels cas, la recherche de racine (par exemple, les méthodes de valeur propre) est souvent plus fiable.

Autre mise en garde importante : l'adaptation de Routh-Hurwitz par transformation bilinéaire n'est valable que pour les systèmes linéaires invariants dans le temps. Elle ne s'étend pas directement aux systèmes discrets non linéaires ou à temps variable.

Enfin, il peut arriver que les conditions soient mal interprétées.Pour le test du jury et du numéro 8217;s, la condition D[(1) > 0 et (–1)nD[[–1) > 0 sont nécessaires mais ne suffisent pas seuls; ils doivent être combinés aux conditions du tableau.

Conclusion

Le choix entre Routh-Hurwitz et Jury’s test de stabilité pour les systèmes discrets dépend du contexte, de la familiarité et des ressources informatiques. Jury’s test est la méthode plus directe, conçue pour les polynômes à temps discret, offrant des conditions claires et évitant la transformation bilinéaire. Le critère adapté de Routh-Hurwitz fournit un pont pour les ingénieurs formés dans les systèmes continus et peut être utile lors de l'intégration d'analyses continues et discrètes.

Pour la plupart des systèmes modernes de contrôle numérique et d'analyse de systèmes discrets, le test de stabilité Jury & #8217;s est l'approche recommandée. Cependant, la compétence des deux techniques permet à l'ingénieur de s'adapter à divers énoncés de problèmes et de tirer parti du logiciel disponible pour vérifier les résultats.

Lecture supplémentaire et ressources externes