Ανάλυση Σ.Α.Ε στο χώρο κατάστασης

1. ΚΕΣ 01 – Αυτόματος Έλεγχος Ανάλυση Σ.Α.Ε στο χώρο κατάστασης

2.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Βιβλιογραφία Ενότητας • Παρασκευόπουλος [2004]: Κεφάλαιο 5: Ενότητες 5.1-5.2 • Παρασκευόπουλος [2005]: Εφαρμογές, Κεφάλαιο 5 • DiStefano [1995]: Chapter 3: Section 3.15, Chapter 20: Sections 20.1 & 20.2 • Tewari [2005]: Chapters3 & 4

3.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Ορισμός Καταστάσεων και Χώρου Κατάστασης • Οι εξισώσεις κατάστασης είναι μια περιγραφή στο πεδίο του χρόνου η οποία μπορεί να χρησιμοποιηθεί για μια μεγάλη γκάμα συστημάτων όπως γραμμικά, μη γραμμικά, χρονικά αναλλοίωτα ή μη, με ή χωρίς αρχικές συνθήκες • Κατάσταση ονομάζουμε ένα σύνολο εσωτερικών μεταβλητών του συστήματος η παρακολούθηση των οποίων στον χρόνο μας περιγράφει το σύστημα. • Οι παραπάνω μεταβλητές ονομάζονται μεταβλητές κατάστασης • Χώρος κατάστασης ονομάζεται ο Ευκλείδιος χώρος ο οποίος δημιουργείται από τις μεταβλητές κατάστασης • Ορισμός: • Οι μεταβλητές κατάστασης x1(t), x2(t), …, xn(t)ενός συστήματος ορίζονται ως ένας (ελάχιστος) αριθμός μεταβλητών τέτοιων ώστε αν γνωρίζουμε τις τιμές τους για οποιαδήποτε χρονική στιγμή t0, τη συνάρτηση εισόδου που εφαρμόζεται στο σύστημα για t≥ t0, και το μαθηματικό νόμο που συνδέει την είσοδο, τις μεταβλητές κατάστασης και το σύστημα, να καθίσταται δυνατός ο προσδιορισμός της κατάστασης του συστήματος για οποιαδήποτε χρονική στιγμή t≥ t0.

4.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Πλήθος και επιλογή μεταβλητών κατάστασης • Ο ελάχιστος αριθμός των μεταβλητών κατάστασης είναι ίσος με την τάξη του συστήματος: • Αυτό είναι απαραίτητο διότι για τον πλήρη προσδιορισμό της εξόδου ενός συστήματος τάξης nχρειάζονται nαρχικές συνθήκες. Εφόσον οι μεταβλητές κατάστασης μπορούν να προδιαγράψουν πλήρως το σύστημα για οποιαδήποτε χρονική στιγμή είναι φανερό ότι πρέπει να είναι ίσες σε πλήθος με το πλήθος των αρχικών συνθηκών. • Οι μεταβλητές κατάστασης για να μπορούν να περιγράψουν πλήρως το σύστημα πρέπει να είναι γραμμικά ανεξάρτητες. • Για την περιγραφή ενός συστήματος μπορούν να επιλεγούν διάφορα σύνολα μεταβλητών κατάστασης φτάνει να έχουν πλήθος n(όσο η τάξη του συστήματος) και να είναι γραμμικά ανεξάρτητες. • Σε ηλεκτρικά και ηλεκτρονικά κυκλώματα οι μεταβλητές κατάστασης είναι συνήθως γραμμικές συναρτήσεις των (α) φορτίων των πυκνωτών, (β) ρευμάτων στα πηνία. Τα στοιχεία αυτά μπορούν να έχουν αρχικές συνθήκες οι οποίες επηρεάζουν τον προσδιορισμό της εξόδου του συστήματος.

5.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Περιγραφή Σ.Α.Ε στο χώρο κατάστασης • Έστω το σύστημα πολλών εισόδων – πολλών εξόδων του σχήματος. Μπορούμε να εκφράσουμε τις mεισόδους, pεξόδους και nμεταβλητές κατάστασης ως διανύσματα: • Οι εξισώσεις κατάστασης ενός συστήματος είναι ένα σύστημα n διαφορικών εξισώσεων πρώτης τάξης που συνδέει το διάνυσμα εισόδου u(t) με το διάνυσμα κατάστασης x(t) και έχει τη μορφή: όπου fείναι μια στήλη με nστοιχεία. Η συνάρτηση fείναι γενικά μια πεπλεγμένη μη γραμμική συνάρτηση των x(t) και u(t) • Το διάνυσμα εξόδουy(t) συνδέεται με τα διανύσματα εισόδου u(t) και κατάστασης x(t) με την εξίσωση εξόδου:

6.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Εξισώσεις Κατάστασης(ΙΙΙ) όπου gείναι μια στήλη με pστοιχεία. Η συνάρτηση gείναι γενικά μια πεπλεγμένη μη γραμμική συνάρτηση των x(t) και u(t) • Οι αρχικές συνθήκες των εξισώσεων κατάστασης είναι οι τιμές του διανύσματος κατάστασης x(t) για t=t0(t0 ισούται συνήθως με 0) και συμβολίζονται ως εξής: • Οι εξισώσεις κατάστασης, η εξίσωση εισόδου και οι αρχικές συνθήκες συνθέτουν την περιγραφή ενός δυναμικού συστήματος στο χώρο κατάστασης:

7.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Εξισώσεις κατάστασης • Αν ένα γραμμικό μη χρονικά μεταβαλλόμενο σύστημα μπορεί να περιγραφεί από ένα σύστημα συνήθων διαφορικών εξισώσεων, τότε οι εξισώσεις κατάστασης παίρνουν την ειδική μορφή: • Ο πίνακας Α έχει διαστάσεις nxn και ονομάζεται πίνακας του συστήματος, ο πίνακας Β έχει διαστάσεις nxm και ονομάζεται πίνακας εισόδου, ο πίνακας C έχει διαστάσεις pxn και ονομάζεται πίνακας εξόδου, ο πίνακας D έχει διαστάσεις pxm και ονομάζεται απευθείας πίνακας.

8.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Περιγραφή γραμμικών χρονικά μεταβαλλόμενωνσυστημάτων • Αν ένα γραμμικό χρονικά μεταβαλλόμενο σύστημα μπορεί να περιγραφεί από ένα σύστημα συνήθων διαφορικών εξισώσεων, τότε οι εξισώσεις κατάστασης παίρνουν τη μορφή:

9.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Παράδειγμα • Το ηλεκτρικό κύκλωμα του σχήματος περιγράφεται από την Ο.Δ.Ε (έξοδος η τάση στα άκρα της αντίστασης): • Θεωρώντας ως μεταβλητές κατάστασης • το ρεύμα στο πηνίο, x1(t)=iL(t) • τo φορτίο του πυκνωτή • τότε ισχύει

10.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Λύση Εξισώσεων Κατάστασης • Η λύση των εξισώσεων κατάστασης στοχεύει στον προσδιορισμό του διανύσματος κατάστασης x(t)για κάθε χρονική στιγμή t≥ t0 (συνήθως το t0 λαμβάνεται ίσο με μηδέν). • Η λύση των εξισώσεων κατάστασης με αρχικές συνθήκες: περιλαμβάνει την εύρεση της λύσης της ομογενούς εξίσωσης: η οποία ονομάζεται ελεύθερη απόκριση του συστήματος, καθώς και την εύρεση της απόκρισης του συστήματος στη διέγερση u(t)η οποία ονομάζεται διεγερμένη απόκριση. • Για την εύρεση της ελεύθερης απόκρισης μετασχηματίζουμε κατά Laplace την ομογενή εξίσωση: => =>

11.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Μεταβατικός πίνακας κατάστασης • Ο πίνακας ονομάζεται μεταβατικός πίνακας κατάστασης διότι μας προσδιορίζει τη μετάβαση του διανύσματος κατάστασης x(t) από την αρχική κατάσταση (Ι είναι ο πίνακας με μοναδικά μη μηδενικά στοιχεία αυτά της κύριας διαγωνίου) σε οποιαδήποτε τελική κατάσταση x(t). • Ο πίνακας μπορεί να προσδιοριστεί από το ανάπτυγμα Taylor: για αυτό και συμβολίζεται με • Ιδιότητες μεταβατικού πίνακα κατάστασης:

12.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Υπολογισμός μεταβατικού πίνακα κατάστασης • Υπάρχουν διάφορες μέθοδοι για τον υπολογισμό του μεταβατικού πίνακα κατάστασης: Μέθοδος 1: Απευθείας υπολογισμός από τη σχέση Η μέθοδος αυτή είναι δύσκολή όταν ο πίνακας Α έχει διαστάσεις μεγαλύτερες από 3x3 εξαιτίας της δυσκολίας αντιστροφής του πίνακα: Παράδειγμα: Να ευρεθεί ο μεταβατικός πίνακας κατάστασης για το σύστημα με και να υπολογίσετε το διάνυσμα κατάστασης

13.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Παράδειγμα • Σχηματίζουμε τον πίνακα: οπότε ο μεταβατικός πίνακας κατάστασης θα είναι: και το διάνυσμα κατάστασης ισούται με:

14.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Παράδειγμα (συν.) • Η χρονική μορφή του διανύσματος κατάστασης φαίνεται στο σχήμα. • Είναι φανερό πως η ελεύθερη απόκριση του συστήματος στις συγκεκριμένες αρχικές συνθήκες προοδευτικά μηδενίζεται. • Για τον υπολογισμό του μεταβατικού πίνακα κατάστασης στη Matlab χρειάζεται η χρήση του symbolic math toolbox: • Εντολές: • syms t (για ορισμό της t ως συμβολικής μεταβλητής) • Phi = expm(A*t); (όπου ο Α έχει οριστεί σύμφωνα με τις τιμές που δόθηκαν στην εκφώνηση του παραδείγματος)

15.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Υπολογισμός μεταβατικού πίνακα κατάστασης (ΙΙ) Μέθοδος 2: Υπολογισμός από τη σχέση με τη βοήθεια του αλγορίθμου του Leverrier o οποίος χρησιμοποιείται για την αντιστροφή του πίνακα: όπου: και οι πίνακες Fiκαι οι συντελεστές aiυπολογίζονται επαναληπτικά από τις σχέσεις: Ο αλγόριθμος Leverrier αναπτύχθηκε για ευκολία υπολογισμού του αντίστροφου του πίνακα μέσω υπολογιστή.

16.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Υπολογισμός μεταβατικού πίνακα κατάστασης (ΙΙΙ) Μέθοδος 3: Με διαγωνοποίηση του πίνακα A (ισχύει εφόσον οι ιδιοτιμές του Aείναι διακριτές, δηλαδή δεν έχουμε ιδιοτιμές με πολλαπλότητα μεγαλύτερη από 1). • Οι ιδιοτιμές λi, i=1,…,n, του πίνακα A δίνονται αποτελούν λύσεις της εξίσωσης: δηλαδή είναι τιμές του λ για τις οποίες μηδενίζεται η ορίζουσα • Σημειώνεται ότι η ορίζουσα μας δίνει το Χαρακτηριστικό Πολυώνυμο του συστήματος. Επομένως οι ιδιοτιμές του πίνακα A αποτελούν τους πόλους του πίνακα συναρτήσεων μεταφοράς. Εφόσον ο πίνακας Α διαγωνοποιείται μπορεί να γραφεί ως: όπου Tο πίνακας με στήλες τα ιδιοδιανύσματα του Α και Λ ο πίνακας με διαγώνια στοιχεία τις ιδιοτιμές του Α.

17.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Υπολογισμός μεταβατικού πίνακα κατάστασης (ΙV) Μέθοδος 3 (συν): Με βάση τα προηγούμενα ο μεταβατικός πίνακας μπορεί να γραφεί: όπου

18.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Παράδειγμα: Να ευρεθεί ο μεταβατικός πίνακας κατάστασης (με διαγωνοποίηση του πίνακα Α) για το σύστημα με και να υπολογίσετε το διάνυσμα κατάστασης Λύση: => λ1=-1, λ2=-2. Οι ιδιοτιμές είναι διακριτές άρα ο πίνακας Α μπορεί να διαγωνοποιηθεί. Τα ιδιοδιανύσματα που αντιστοιχούν στις παραπάνω ιδιοτιμές ικανοποιούν τη σχέση: και είναι (βλέπε εντολή eigστη Matlab) οπότε ο πίνακας Τ είναι:

19.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Παράδειγμα (συν.) Επομένως ο μεταβατικός πίνακας κατάστασης θα είναι: τελικά και το διάνυσμα κατάστασης ισούται με:

20.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Παράδειγμα (συν.) • Η χρονική μορφή του διανύσματος κατάστασης φαίνεται στο σχήμα. • Είναι φανερό πως η ελεύθερη απόκριση του συστήματος στις συγκεκριμένες αρχικές συνθήκες προοδευτικά μηδενίζεται.

21.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Υπολογισμός μεταβατικού πίνακα κατάστασης (V) Μέθοδος 4: Η τελευταία μέθοδος για υπολογισμό του μεταβατικού πίνακα κατάστασης βασίζεται στη σχέση: και είναι καθαρά προγραμματιστική. Επειδή η παραπάνω σειρά έχει άπειρους όρους ο υπολογισμός του μεταβατικού πίνακα κατάστασης είναι προσεγγιστικός:

22.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Γενική λύση εξισώσεων κατάστασης • Η γενική λύση των εξισώσεων κατάστασης • στοχεύει υπολογίζει και τη διεγερμένη απόκριση η οποία βασίζεται και αυτή στον μεταβατικό πίνακα κατάστασης και δίνεται από το συνελικτικό ολοκλήρωμα: • Ο υπολογισμός του παραπάνω συνελικτικού ολοκληρώματος είναι δύσκολος για τις περισσότερες μορφές εισόδου (εξαίρεση αποτελούν η κρουστική και η βηματική συνάρτηση). • Ο απλούστερος τρόπος για την εύρεση του διανύσματος κατάστασης είναι η χρήση του μετασχηματισμού Laplace: => =>

23.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Παράδειγμα Παράδειγμα: Να υπολογίσετε το διάνυσμα κατάστασης για το σύστημα: όταν η είσοδος είναι η βηματική συνάρτηση us(t) Λύση:

24.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Παράδειγμα (συν.) Οπότε τελικά:

25.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Ελέγξιμο διανύσματος κατάστασης • Η έννοια του ελέγξιμου αναφέρεται στη δυνατότητα ελέγχου του διανύσματος κατάστασης από το διάνυσμα εισόδου. • Η δυνατότητα προσδιορισμού της ελεγξιμότητας είναι ένα από τα βασικά πλεονεκτήματα της περιγραφής συστημάτων μέσω των εξισώσεων κατάστασης. Ορισμός: • Το διάνυσμακατάστασης x(t)είναι ελέγξιμο αν υπάρχει κάποια τμηματικά συνεχής συνάρτηση εισόδου (ελέγχου) u(t) που μπορεί να οδηγήσει το x(t) από την αρχική συνθήκη x(t0) στη τελική του τιμή x(tf) σε ένα πεπερασμένο χρονικό διάστημα (tf-t0)≥0. Θεώρημα: • Το διάνυσμακατάστασης x(t)του συστήματος που περιγράφεται από τις παραπάνω εξισώσεις κατάστασης είναι ελέγξιμο τότε και μόνο τότε η τάξη του πίνακα S (διαστάσεων nxnm) είναι ίση με n(υπάρχουν δηλαδή nανεξάρτητες στήλες στον πίνακα S)

26.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Ελέγξιμο διανύσματος εξόδου • Η έννοια του ελέγξιμου της εξόδου αναφέρεται στη δυνατότητα ελέγχου του διανύσματος εξόδου από το διάνυσμα εισόδου. Ορισμός: • Το διάνυσμαεξόδου y(t)είναι ελέγξιμο αν υπάρχει κάποια τμηματικά συνεχής συνάρτηση εισόδου (ελέγχου) u(t) που μπορεί να οδηγήσει το y(t) από την αρχική συνθήκη y(t0) στη τελική του τιμή y(tf) σε ένα πεπερασμένο χρονικό διάστημα (tf-t0)≥0. Θεώρημα: • Το διάνυσμαεξόδου y(t)του συστήματος που περιγράφεται από τις παραπάνω εξισώσεις κατάστασης είναι ελέγξιμο τότε και μόνο τότε η τάξη του πίνακα Q (διαστάσεων px(m+1)n) είναι ίση με p(υπάρχουν δηλαδή pανεξάρτητες στήλες στον πίνακα Q)

27.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Παράδειγμα Να ελεγχθεί το ελέγξιμο του διανύσματος κατάστασης και του διανύσματος εξόδου για το σύστημα με Λύση: Κατασκευάζουμε τους πίνακεςSκαι Q(έχουμε n=m=p=2): η τάξη του Sείναι 2, άρα το διάνυσμα κατάστασης είναι ελέγξιμο. Η τάξη του Q είναι 2 άρα το διάνυσμα εξόδου είναι ελέγξιμο

28.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Παράδειγμα (ΙΙ) Να ελεγχθεί το ελέγξιμο του διανύσματος κατάστασης και του διανύσματος εξόδου για το σύστημα με Λύση: Κατασκευάζουμε τους πίνακεςSκαι Q(έχουμε n=2,m=p=1): η τάξη του Sείναι 1 (ορίζουσα |S |=0), άρα το διάνυσμα κατάστασης δεν είναι ελέγξιμο. Η τάξη του Q είναι 1 άρα το διάνυσμα εξόδου είναι ελέγξιμο

29.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Παρατηρήσιμο διανύσματος κατάστασης • Η έννοια του παρατηρήσιμου αναφέρεται στη δυνατότητα προσδιορισμού των αρχικών συνθηκών x(t0)(αρχική κατάσταση συστήματος) με βάση τα διανύσματα εισόδου u(t)και εξόδου y(t)τα οποία μετράμε σε ένα πεπερασμένο χρονικό διάστημα. • Σε περίπτωση που έστω καιμια μεταβλητή κατάστασης δεν είναι παρατηρήσιμη τότε το σύστημα συνολικά δεν είναι παρατηρήσιμο. • Η δυνατότητα προσδιορισμού του παρατηρήσιμου ενός συστήματος είναι ένα από τα βασικά πλεονεκτήματα της περιγραφής συστημάτων μέσω των εξισώσεων κατάστασης. Ορισμός: • Το διάνυσμακατάστασης x(t)είναι παρατηρήσιμο στο διάστημα [t0tf] όταν γνωρίζοντας τα τα διανύσματα εισόδου u(t)και εξόδου y(t) για tє[t0tf]μπορούμε να προσδιορίσουμε το διάνυσμα αρχικών συνθηκών x(t0)

30.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Παρατηρήσιμο διανύσματος κατάστασης (ΙΙ) Θεώρημα: • Το διάνυσμακατάστασης x(t)του συστήματος που περιγράφεται από τις παρακάτω εξισώσεις κατάστασης είναι παρατηρήσιμο τότε και μόνο τότε η τάξη του πίνακα R (διαστάσεων nxnp) είναι ίση με n(υπάρχουν δηλαδή nανεξάρτητες στήλες στον πίνακα R) Παράδειγμα: • Να βρεθεί αν το διάνυσμα κατάστασης του συστήματος με είναι παρατηρήσιμο: Λύση Ο πίνακας Rείναι τάξης 2, άρα το σύστημα δεν είναι παρατηρήσιμο

31.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Πίνακας Συναρτήσεων Μεταφοράς • Από τις εξισώσεις κατάστασης μπορούμε να μεταβούμε σε περιγραφή με πίνακα συναρτήσεων μεταφοράς H(s)(υπό την προϋπόθεση μηδενικών αρχικών συνθηκών) με χρήση το μετασχηματισμού Laplace: => => => οπότε άρα

32.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Συναρτήσεις Μεταφοράς και ελέγξιμο και παρατηρήσιμο Θεώρημα 1: • Αν ο πίνακας συναρτήσεων μεταφοράς H(s) ενός συστήματος παρουσιάζει απαλοιφή πόλων μηδενικών τότε το σύστημα είναι είτε μη ελέγξιμο ή μη παρατηρήσιμο ή και τα δύο. Αν δεν έχουμε απαλοιφή πόλων μηδενικών τότε το σύστημα με πίνακα συναρτήσεων μεταφοράς H(s) μπορεί να περιγραφεί με εξισώσεις κατάστασης ως ένα ελέγξιμο και παρατηρήσιμο σύστημα. Θεώρημα 2: • Ο πίνακας συναρτήσεων μεταφοράς H(s) περιέχει μόνο το ελέγξιμο και παρατηρήσιμο μέρος ενός συστήματος (εκτός ειδικών περιπτώσεων όπου το χαρακτηριστικό πολυώνυμο έχει ρίζες πολλαπλότητας μεγαλύτερης από 1).

33.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Συναρτήσεις Μεταφοράς και ελέγξιμο και παρατηρήσιμο (ΙΙ) Παράδειγμα: • Έστω το ηλεκτρονικό κύκλωμα του σχήματος: • Ο πίνακας συναρτήσεων μεταφοράς H(s) έχει τη μορφή: Παρατηρούμε ότι αν C1Rf1= C2Rf2 η H12(s) παρουσιάζει απαλοιφή πόλου μηδενικού και γίνεται οπότε το σύστημα εμφαίνεται ως τάξης 1 (ενώ είναι φανερό ότι είναι τάξης 2)

34.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Έλεγχος ευστάθειας στο χώρο κατάστασης • Συστήματα με περιγραφή στο χώρο κατάστασης μπορούν να ελεγχθούν ως προς την ευστάθεια τους με τη βοήθεια των πιο κάτω ορισμών: • Ορισμός 1: • Έστω το Γ.Χ.Α σύστημα Θεωρούμε ότι το σύστημα είναι μηδενικής διέγερσης (u(t)=0), εξετάζουμε δηλαδή την ελεύθερη απόκριση του συστήματος. Το σύστημα είναι ευσταθές αν για κάθε πεπερασμένη αρχική συνθήκη x(0) υπάρχει αριθμός Μ(x(0)) τέτοιος ώστε να ικανοποιούνται οι συνθήκες: όπου είναι το μέτρο του διανύσματος κατάστασης για ελεύθερη είσοδο. Με δεδομένο ότι η ελεύθερη απόκριση του συστήματος δίνεται από τη σχέση: και ότι το x(0) είναι πεπερασμένο, οι παραπάνω σχέσεις ισοδυναμούν με δηλαδή όλα τα στοιχεία του μεταβατικού πίνακα κατάστασης μηδενίζονται με την πάροδο του χρόνου

35.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Έλεγχος ευστάθειας στο χώρο κατάστασης (ΙΙ) • Ορισμός 2: • Έστω το Γ.Χ.Α σύστημα: για το οποίο ισχύει (εφόσον είναι ελέγξιμο και παρατηρήσιμο): Το σύστημα είναι ευσταθές αν οι ιδιοτιμές του πίνακα A (ισοδύναμα οι πόλοι του Χ.Π. ) βρίσκονται στο αριστερό μιγαδικό ημιεπίπεδο (έχουν πραγματικό μέρος αρνητικό).

36.  Ορισμός Κατάστασης  Περιγραφή Σ.Α.Ε στο χώρο κατάστασης  Λύση εξισώσεων κατάστασης  Ελέγξιμο Σ.Α.Ε  Παρατηρήσιμο Σ.Α.Ε  Ευστάθεια στο Χώρο Κατάστασης Παράδειγμα Παράδειγμα: • Να ελεγχθεί η ευστάθεια του συστήματος με: και να υπολογιστεί ο πίνακας συναρτήσεων μεταφοράς Λύση Οι ιδιοτιμές του πίνακα Αείναι λ1=-3, λ2=-2, λ3=-1. Όλες έχουν αρνητικό πραγματικό μέρος άρα το σύστημα είναι ευσταθές. Ο μεταβατικός πίνακας κατάστασης είναι: όλα τα στοιχεία του μηδενίζονται όταν t->∞, άρα και με αυτό το κριτήριο προκύπτει ότι το σύστημα είναι ευσταθές. Ο πίνακας συναρτήσεων μεταφοράς είναι: Είναι φανερό ότι έχουμε απαλοιφή πόλων μηδενικών διότι σύμφωνα με τον πίνακα συναρτήσεων μεταφοράς το σύστημα είναι πρώτης τάξης με Χ.Π. p(s)=(s+3)