Κριτήριο ευστάθειας Routh-Hurwitz

Το κριτήριο ευστάθειας Routh-Hurwitz δείχνει πόσες ρίζες ενός χαρακτηριστικού πολυωνύμου βρίσκονται στο δεξί ημιεπίπεδο χωρίς να χρειάζεται να λυθούν άμεσα οι ρίζες. Στα συστήματα ελέγχου, αυτό δίνει έναν γρήγορο έλεγχο για το αν ένα σύστημα συνεχούς χρόνου είναι ασυμπτωτικά ευσταθές.

Για ένα πολυώνυμο με πραγματικούς συντελεστές

$$p(s) = a_n s^n + a_{n-1} s^{n-1} + \cdots + a_1 s + a_0,$$

με $a_n > 0$, κατασκευάζεις τον πίνακα Routh και εξετάζεις την πρώτη του στήλη. Αφού αντιμετωπιστεί κάθε ειδική περίπτωση, όπως ένα μηδέν στην πρώτη στήλη ή μια ολόκληρη γραμμή από μηδενικά, ο αριθμός των αλλαγών προσήμου στην πρώτη στήλη ισούται με τον αριθμό των ριζών με θετικό πραγματικό μέρος. Αν δεν υπάρχουν αλλαγές προσήμου, όλες οι ρίζες βρίσκονται στο ανοιχτό αριστερό ημιεπίπεδο.

Τι ελέγχει το κριτήριο Routh-Hurwitz

Στα περισσότερα προβλήματα ελέγχου, ευστάθεια σημαίνει ότι κάθε χαρακτηριστική ρίζα ικανοποιεί $\operatorname{Re}(s) < 0$. Αυτή η συνθήκη είναι σημαντική: το συνηθισμένο κριτήριο Routh-Hurwitz αφορά συστήματα συνεχούς χρόνου στο επίπεδο $s$, όχι συστήματα διακριτού χρόνου στο επίπεδο $z$.

Η πρακτική του αξία είναι η ταχύτητα. Μπορείς να αποφασίσεις για την ευστάθεια μόνο από τους συντελεστές, κάτι που συχνά είναι ευκολότερο από το να βρεις τις ακριβείς ρίζες.

Πώς να κατασκευάσεις τον πίνακα Routh

Ξεκίνα γράφοντας τους συντελεστές κατά φθίνουσες δυνάμεις του $s$. Οι δύο πρώτες γραμμές εναλλάσσουν τους συντελεστές:

$$\begin{array}{c|ccc} s^n & a_n & a_{n-2} & a_{n-4} \\ s^{n-1} & a_{n-1} & a_{n-3} & a_{n-5} \end{array}$$

Έπειτα υπολόγισε τις κατώτερες γραμμές από τις δύο γραμμές που βρίσκονται ακριβώς από πάνω τους. Για ένα κυβικό πολυώνυμο,

$$p(s) = s^3 + a_1 s^2 + a_2 s + a_3,$$

ο πίνακας είναι:

$$\begin{array}{c|cc} s^3 & 1 & a_2 \\ s^2 & a_1 & a_3 \\ s^1 & \dfrac{a_1 a_2 - a_3}{a_1} & 0 \\ s^0 & a_3 & 0 \end{array}$$

Άρα, για ένα κυβικό πολυώνυμο με θετικό πρώτο συντελεστή, η ασυμπτωτική ευστάθεια απαιτεί

$$a_1 > 0,\quad a_2 > 0,\quad a_3 > 0,\quad a_1 a_2 > a_3.$$

Αυτή η τελευταία ανισότητα είναι το σημείο που συχνά παραβλέπεται. Οι θετικοί συντελεστές από μόνοι τους δεν αρκούν.

Λυμένο παράδειγμα: θετικοί συντελεστές αλλά ασταθές

Θεώρησε το

$$p(s) = s^3 + s^2 + 2s + 8.$$

Όλοι οι συντελεστές είναι θετικοί, οπότε με μια πρώτη ματιά μπορεί να φαίνεται ευσταθές. Ο πίνακας Routh δείχνει γιατί αυτό το συμπέρασμα είναι λάθος.

$$\begin{array}{c|cc} s^3 & 1 & 2 \\ s^2 & 1 & 8 \\ s^1 & \dfrac{(1)(2) - (1)(8)}{1} & 0 \\ s^0 & 8 & 0 \end{array}$$

Απλοποιώντας τη γραμμή $s^1$ παίρνουμε

$$\begin{array}{c|cc} s^3 & 1 & 2 \\ s^2 & 1 & 8 \\ s^1 & -6 & 0 \\ s^0 & 8 & 0 \end{array}$$

Τώρα κοίτα μόνο την πρώτη στήλη:

$$1,\quad 1,\quad -6,\quad 8.$$

Υπάρχουν δύο αλλαγές προσήμου, από $1$ σε $-6$ και από $-6$ σε $8$. Άρα το πολυώνυμο έχει $2$ ρίζες στο δεξί ημιεπίπεδο, πράγμα που σημαίνει ότι το σύστημα είναι ασταθές.

Αυτή είναι η βασική ιδέα που πρέπει να κρατήσεις: για βαθμό $3$ και πάνω, οι θετικοί συντελεστές δεν εγγυώνται ευστάθεια.

Συνηθισμένα λάθη με το κριτήριο Routh-Hurwitz

Έλεγχος μόνο των συντελεστών

Αν όλοι οι συντελεστές είναι θετικοί, αυτό δεν σημαίνει αυτόματα ότι το σύστημα είναι ευσταθές. Το λυμένο παράδειγμα παραπάνω είναι ακριβώς το αντιπαράδειγμα.

Παράβλεψη του πεδίου εφαρμογής του ελέγχου

Το συνηθισμένο κριτήριο Routh-Hurwitz εφαρμόζεται σε χαρακτηριστικά πολυώνυμα συνεχούς χρόνου στο $s$. Αν μελετάς σύστημα διακριτού χρόνου, χρειάζεσαι διαφορετικό έλεγχο.

Αγνόηση ειδικών περιπτώσεων

Ένα μηδέν στην πρώτη στήλη ή μια ολόκληρη γραμμή από μηδενικά είναι ειδική περίπτωση, όχι ένα κανονικό σημείο διακοπής. Αυτές οι περιπτώσεις χρειάζονται επιπλέον διαδικασία, συχνά με μια μικρή διαταραχή ή με ένα βοηθητικό πολυώνυμο.

Μη κανονικοποίηση του πολυωνύμου

Είναι πιο ασφαλές να γράφεις το πολυώνυμο κατά φθίνουσες δυνάμεις του $s$ και να κάνεις θετικό τον πρώτο συντελεστή πριν κατασκευάσεις τον πίνακα. Αλλιώς ο έλεγχος προσήμων μπορεί εύκολα να παρερμηνευτεί.

Πότε χρησιμοποιείται το κριτήριο Routh-Hurwitz

Χρησιμοποίησε το κριτήριο Routh-Hurwitz όταν ένα μοντέλο οδηγεί σε χαρακτηριστικό πολυώνυμο και χρειάζεσαι γρήγορα μια απάντηση για την ευστάθεια.

Στα συστήματα ελέγχου, ελέγχει αν ένα σύστημα κλειστού βρόχου είναι ευσταθές χωρίς να υπολογίζονται ρητά οι πόλοι. Σε ηλεκτρικά κυκλώματα και μηχανικά μοντέλα, βοηθά να ελεγχθεί αν οι επιλογές παραμέτρων οδηγούν σε αποκρινόμενες που φθίνουν ή αυξάνονται. Στον σχεδιασμό, είναι χρήσιμο για τον εντοπισμό περιοχών τιμών παραμέτρων που διατηρούν ένα σύστημα ευσταθές.

Είναι ιδιαίτερα χρήσιμο όταν η εύρεση των ακριβών ριζών θα ήταν αργή ή περιττή.

Δοκίμασε έναν παρόμοιο έλεγχο ευστάθειας

Δοκίμασε την ίδια διαδικασία στο

$$p(s) = s^3 + 4s^2 + 5s + 2.$$

Κατασκεύασε την πρώτη στήλη και έλεγξε αν εμφανίζονται αλλαγές προσήμου. Έπειτα σύγκρινε την απάντησή σου με τη συνθήκη για κυβικό πολυώνυμο $a_1 a_2 > a_3$ ώστε να δεις ότι και οι δύο μέθοδοι συμφωνούν.