Σε αυτή την εργασία παρουσιάζουμε μεθόδους για την κατασκευή συναρτήσεων Lyapunov για δυναμικά συστήματα αλλά και για τον καθορισμό του ελκτικού συνόλου ενός σημείου ισορροπίας.
Η μελέτη των διαφορικών εξισώσεων έχει ως κίνητρο τις πολλαπλές εφαρμογές τους στη Φυσική, τη Χημεία, τα Οικονομικά, τη Βιολογία, κ.λ.π.. Εστιάζουμε στις αυτόνομες διαφορικές εξισώσεις της μορφής οι οποίες ορίζουν ένα δυναμικό σύστημα. Οι πιο απλές λύσεις μίας τέτοιας εξίσωσης καλούνται σημεία ισορροπίας. Πολύ σημαντικός είναι επίσης και ο καθορισμός του ελκτικού συνόλου. Ο καθορισμός του ελκτικού συνόλου επιτυγχάνεται μέσω υποεπίπεδων συνόλων μίας συνάρτησης Lyapunov, δηλαδή μίας συνάρτησης με αρνητική παράγωγο κατά μήκος των τροχιών στη περιοχή ισορροπίας.
Σε αυτή την εργασία παρουσιάζουμε μεθόδους κατασκευής συναρτήσεων Lyapunov για ένα σημείο ισορροπίας. Υπάρχει πλούσια βιβλιογραφία πάνω στις συναρτήσεις Lyapunov. Το 1893, ο Lyapunov εισήγαγε την άμεση ή δεύτερη μέθοδό του, όπου κατάφερε να εξασφαλίσει αποτελέσματα για την ευστάθεια ενός σημείου ισορροπίας χωρίς να γνωρίζει τη λύση της διαφορικής εξίσωσης, αλλά χρησιμοποιώντας μόνο την ίδια τη διαφορική εξίσωση. Από τότε έχει δοθεί πλήθος αντίστροφων θεωρημάτων, που εξασφαλίζουν την ύπαρξη μίας συνάρτησης Lyapunov, από διάφορους συγγραφείς. Το πρώτο κύριο θεώρημα για ασυμπτωτική ευστάθεια δόθηκε από τον Massera το 1949 και από τότε έχει βελτιωθεί από πολλούς συγγραφείς προς διάφορες κατευθύνσεις. Ωστόσο, κανένα από τα θεωρήματα ύπαρξης δεν παρέχει μία μέθοδο σαφούς κατασκευής μίας συνάρτησης Lyapunov.
Για γραμμικά συστήματα μπορεί κάποιος να κατασκευάσει μία τετραγωνικής μορφής συνάρτηση Lyapunov της μορφής με ένα συμμετρικό, θετικά ορισμένο πίνακα , όπου συμβολίζει το σημείο ισορροπίας. Ο Hahn περιγράφει πως μπορεί κάποιος, ξεκινώντας από ένα μη-γραμμικό σύστημα, να χρησιμοποιήσει την τετραγωνικής μορφής συνάρτηση Lyapunov του γραμμικοποιημένου συστήματος σα μία συνάρτηση Lyapunov για το μη-γραμμικό σύστημα.
Πολλές προσεγγίσεις θεωρούν ειδικές συναρτήσεις Lyapunov, όπως τετραγωνικής μορφής, πολυωνυμικές, κατά τμήματα γραμμικές, ή κατά τμήματα τετραγωνικής μορφής. Οι μέθοδοι όμως αυτές μπορούν να χρησιμοποιηθούν μόνο σε συγκεκριμμένες διαφορικές εξισώσεις.
Σε αυτή την εργασία θα ασχοληθούμε με δύο μεθόδους κατασκευής συναρτήσεων Lyapunov. Για τη πρώτη μέθοδο κατασκευής συναρτήσεων Lyapunov για ένα σημείο ισορροπίας, ξεκινούμε με ένα θεώρημα που εξασφαλίζει την ύπαρξη μίας συνάρτησης Lyapunov η οποία ικανοποιεί την ισότητα , όπου είναι μία γνωστή σταθερά. Βασικός στόχος της μεθόδου είναι να προσεγγίσει τη λύση αυτής της μερικής διαφορικής εξίσωσης με τη χρήση συναρτήσεων ακτινωτής βάσης. Τότε και η προσέγγιση είναι μία συνάρτηση Lyapunov και έτσι, μπορούμε να τη χρησιμοποιήσουμε για να καθορίσουμε το ελκτικό σύνολο. Επειδή η συνάρτηση δεν ορίζεται στο , μελετούμε και μία δεύτερη κλάση συναρτήσεων Lyapunov , οι οποίες ορίζονται και είναι ομαλές στο . Αυτές ικανοποιούν την ισότητα , όπου είναι μία δοθείσα συνάρτηση με συγεκριμμένες ιδιότητες, μία εκ των οποίων είναι ότι . Για την προσέγγιση χρησιμοποιούμε συναρτήσεις ακτινωτής βάσης.
Στη δεύτερη μέθοδο κατασκευάζουμε μια κατά τμήματα γραμμική συνάρτηση Lyapunov για το αρχικό μη-γραμμικό σύστημα χρησιμοποιώντας γραμμικό προγραμματισμό. / In this diploma work we present methods for the construction of Lyapunov functions for dynamical systems but also we determine the basin of attraction of
an equilibrium.
The study of differential equations is motivated from numerous applications in physics, chemistry, economics, biology, etc.
We focus on autonomous differential equations x’ = f(x), x ∈ Rn which define a dynamical system. The simplest solutions x(t) of such an equation are equilibria, i.e. solutions x(t) = x0 which remain constant. An important and non-trivial task is thedetermination of their basin of attraction.
The determination of the basin of attraction is achieved through sublevel sets of a Lyapunov function, i.e. a function with negative orbital derivative. The orbital derivative V ‘(x) of a function V (x) is the derivative along solutions of the differential equation.
In this book we present a method to construct Lyapunov functions for an equilibrium. There is a rich literature on the functions of Lyapunov. In 1893, Lyapunov introduced the direct method, where he managed to secure results for the stability of an equilibrium point without knowing the solution of the differential equation, but using only the same differential equation.
Since then many inverse theorems have been given that ensure the existence of a function Lyapunov, by various authors. The first main theorem on asymptotic stability given by Massera in 1949 and since then has been improved by many authors in different directions. However, none of the theorem of existence does not provide a clear method of manufacturing a Lyapunov function.
For linear systems, one can construct a quadratic form of a Lyapunov function with a symmetric positive definite table. The Hahn describes how people, starting from a non-linear system, use the like it were a Lyapunov function for the nonlinear system.
Many approaches consider special functions Lyapunov, such as quadratic form, polynomial. These methods can be used only in specific differential equations.
In this book we present a method to construct Lyapunov functions for an equilibrium. We start from a theorem which ensures the existence of a Lyapunov function T which satisfies the equation T’(x) = −c, where -c > 0
is a given constant. This equation is a linear first-order partial differential equation. The main goal of this method is to approximate the solution T of this partial differential equation using radial basis functions. Then the approximation itself is a Lyapunov function, and thus can be used to determine the basin of attraction.
Since the function T is not defined at x0, we also study a second class of Lyapunov functions V which are defined and smooth at x0. They satisfy the equation V ‘(x) = −p(x), where p(x) is a given function with certain properties, in particular p(x0) = 0.
For the approximation we use radial basis functions, a powerful meshless approximation method.
In the second method we construct a linear Lyapunov function for the original non-linear system using linear programming.
Identifer | oai:union.ndltd.org:upatras.gr/oai:nemertes:10889/6445 |
Date | 14 October 2013 |
Creators | Αλωνιάτη, Μαρία |
Contributors | Κούσουλας, ΝΙκόλαος, Aloniati, Maria, Κούσουλας, Νικόλαος, Καζάκος, Δημοσθένης |
Source Sets | University of Patras |
Language | gr |
Detected Language | Greek |
Type | Thesis |
Rights | 0 |
Page generated in 0.0135 seconds