Ordinary Differential Equations

First Order D.O.:


F(t,y,y)=0


Λυμένη (ή Κανονική) μορφή F.O.D.O.:


y=f(t,y)

Εξισώσεις Χωριζομένων Μεταβλητών:


\( \frac{dy}{dt} = g(t) \; h(y) \)
\(\Rightarrow \frac{dy}{h(y)} = g(t) \; dt \; , \; \forall \; h(y) \neq 0 \)


Note: For \( h(y) =0 \) there's still a (constant) solution.

Ορισμοί


Γραμμικές Εξισώσεις σε κανονική μορφή (Γ.Δ.Ε.):


\( y'=-p(t) y + g(t) \)


Μέθοδος του Ολοκληρωτικού Παράγοντα


Κατάλληλο \( μ(t) \ne 0 \):
\( μ(t) \frac{dy}{dt} + μ(t) p(t) y = μ(t) g(t) \Leftrightarrow\frac{d}{dt} \left( μ(t) y(t) \right) = μ(t) g(t) \).


Ισχύει για: \( \Rightarrow μ(t) = \exp \left( \int dt \; p(t) \right) \) ή \( \Rightarrow μ(t) = \exp \left( \int^t_{t_0} dτ \; p(τ) \right) \)

Μέθοδος Μεταβολής των Παραμέτρων


Η λύση της ομογενούς εξίσωσης \( g(t)=0 , \forall t \) είναι: \( y(t) = c \; \exp\left( \int dt \; p(t) \right) \).


Αναζητούμε λύση της μορφής: \( y(t) = c(t) \; \exp\left( \int dt \; p(t) \right) \).


Ισχύει γιά: \( \frac{d c(t)}{dt} = g(t) \exp \left( \int dt \; p(t) \right) \)

Γενική λύση (Τύπος της μεταβολής των παραμέτρων):


\( \Rightarrow y(t) = \frac{1}{μ(t)} \left[ \int dt \; μ(t) g(t) + c \right] = e^{ - \int dt \; p(t) } \left[ c + \int dt \; g(t) e^{ \int dt \; p(t) } \right] \)


Για Π.Α.Τ.: \( \Rightarrow y(t) = y_0 e^{ - \int^t_{t_0} p(s) \; ds } + \int_{t_0}^t dτ \; g(τ) \; e^{ - \int^t_{τ} p(s) \; ds } \)


Ίδια γενική λύση.


Note: Η γενική λύση είναι το άθροισμα της γενικής λύσης της ομογενούς εξίσωσης και μίας λύσης της μη ομογενούς εξόισωσης.


Περιοδικές Εξισώσεις


Της μορφής: \( y'=-p(t) y + g(t) \) για τις οποίες: \( p(t) \) και \( g(t) \) είναι περιοδικές με περίοδο Τ.


Note: Η περιοδική συνάρτηση \( g(t) \) ονομάζεται και συνάρτηση εξαναγκασμού.

Τετριμμένη λύση

Χώρος λύσεων

Theorem 1.1: Fredholm Theorem


Fredholm_Theorem

Ο χώρος των λύσεων είναι μονοδιάστατος

Theorem 1.2


img_20180224_155549.8281

Απεικόνιση Poincaré:


07 03 2018 15 27 Office Lens


Note:

  • Ορισμός σταθερών σημείων: \( Π(y_0^*) = y_0^* \).
  • Δεν συνεπάγει περιοδική λύση κάθε αρχική συνθήκη.
  • Υπάρχουν Δ.Ε.: \( \lim_{n\to\infty} Π^n(y_0) = y_0^* \).

Βασικές Εξισώσεις

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

Οι αναδρομικές σχέσεις είναι εύκολες!

Εξισώσεις ομογενείς και αναγόμενες σε ομογενείς


Εξισώσεις της μορφής: \( \frac{dy}{dt} = f(at+by) \; , \; a,b \in \mathcal{R} \) ανάγονται σε εξισώσεις χωριζομένων μεταβλητών μέσω του μετασχηματισμού: \( z = at + by \).

n-Order Homogeneous Equations:


Homogeneous_Equations


n-Order_Homogeneous_Equations

Εξισώσεις που ανάγονται σε ομογενείς


02 03 2018 16 09 Office Lens

SecondOrder D.O.:


\( y''=f(y,y') \)

Μετασχηματισμός: \( V = \frac{dy}{dt} \) μετατρέπει σε 1ου βαθμού: \( v \frac{dV}{dt} = h(y,V) \)
Note: Συνεπώς η αρχική Δ.Ε. 2ης τάξης είναι ισοδύναμη με δύο εξισώσεις πρώτης τάξης.

Ακριβείς Εξισώσεις

Ορισμός 1.6 (Ακριβής διαφορική μορφή:):


28 02 2018 16 08 Office Lens #

Θεώρημα 1.4 (Ικανή και αναγκαία συνθήκη για ακριβή διαφορική μορφή):


28 02 2018 16 20 Office Lens

Επίλυση με επικαμπύλια ολοκλήρωση:


28 02 2018 17 29 Office Lens
Note: Δεδομένης της απλής συνεκτικότητας του πεδίου ορισμού της \( F(t,y) \) συνάρτησης.

Ολοκληρωτικοί παράγοντες:


Μη ακριβής διαφορικές εξισώσεις: \( M(t,y) dt + N(t,y) dy = 0 \) μπορούν να μετατραπούν σε ακριβείς με κατάλληλους ολοκληρωτικούς παράγοντες \( μ(t,y) \): \( μ(t,y) M(t,y) dt + μ(t,y) N(t,y) dy = 0 \).


Η ύπαρξη ενός τέτοιου ολοκληρωτικού παράγοντα προϋποθέτει την ύπαρξη γενικής λύσης της αρχικής διαφορικής εξίσωσης.

Γενική διαδικασία κατασκευής ενός ολοκληρωτικού παράγοντα:


\( \frac{\partial(μM)}{\partial y} = \frac{\partial(μN)}{\partial t} \)


\( \Rightarrow μ\frac{\partial M}{\partial y} + M\frac{\partial μ}{\partial y} = μ\frac{\partial N}{\partial t} + N\frac{\partial μ}{\partial t} \)


\( \Rightarrow \frac{1}{μ} \left( N \frac{\partial μ}{\partial t} - M \frac{\partial μ}{\partial y} \right) = \frac{\partial M}{\partial y} - \frac{\partial N}{\partial t} \)

Περιπτώσεις:

Συνάρτηση μόνο του \( t \):


\( \Rightarrow \frac{1}{μ} \frac{dμ}{dt} = \frac{1}{N} \left( \frac{\partial M}{\partial y} - \frac{\partial N}{\partial t} \right) \equiv g(t) \)


\( \Rightarrow μ(t) = e^{\int g(t) dt} \)

Συνάρτηση μόνο του \( y \):


\( \Rightarrow \frac{1}{μ} \frac{dμ}{dy} = -\frac{1}{M} \left( \frac{\partial M}{\partial y} - \frac{\partial N}{\partial t} \right) \equiv g(y) \)


\( \Rightarrow μ(y) = e^{\int g(y) dy} \)

Συνάρτηση του \( ty \):


\( \Rightarrow \frac{1}{μ} \frac{dμ}{d(ty)} = -\frac{1}{Ny-Mt} \left( \frac{\partial M}{\partial y} - \frac{\partial N}{\partial t} \right) \equiv g(ty) \)


Θέτουμε \( s \equiv ty \):
\( \Rightarrow \frac{1}{μ} \frac{dμ}{ds} = g(s) \Rightarrow μ(ty) = e^{G(s)} \), όπου \( G'(s) = g(s) \)

Αυτόνομη εξίσωση: Δεν εξαρτάται ευθέως από το \( t \).

Συστήματα

Σύστημα Lotka-Volterra:


\( \frac{dx}{dt} = (α - β y) x \)
\( \frac{dy}{dt} = (δ x - γ) y \)


Γενική Λύση:


Απαλοιφή του \( t \):


\( \Rightarrow - (δ x - γ) y \; dx + (α - β y) x \; dy = 0 \)
η οποία εξίσωση μετατρέπεται σε ακριβή εξίσωση με τον ολοκληρωτικό παράγοντα \( \frac{1}{xy} \):


\( \Rightarrow - \bigg( δ - \frac{γ}{x} \bigg) dx + \bigg(-β + \frac{α}{y} \bigg) dy = 0 \)


\( \Rightarrow F(x,y) = - δ x + γ \ln(x) - β y + α \ln(y) = c \)

Χαμιλτονιανά Συστήματα:


5a772795daf6ee0001851c20-7666561c-d4ea-49c9-9505-8816984814aa


Η πιο σημανική τους ιδιότητα:
02 03 2018 16 15 Office Lens

Θεώρημα 1.3:


02 03 2018 16 02 Office Lens

Τα Βασικά Θεωρήματα

Απόλυτη σύκλιση συνεπάγει σύγκλιση.

Προσεγγιστικές Μέθοδοι

Μέθοδος Euler:


05 03 2018 13 31 Office Lens

Μέθοδος Picard:


Μετατρέπει το Π.Α.Τ. σε ισοδύναμη ολοκληρωτική εξίσωση., την οποία στην συνέχεια και επιλύει:

\( y(t) = y_0 + \int^{t}_{t_0} f(s,y(s)) \; ds \)


05 03 2018 13 13 Office Lens

Σύγκλιση Προσεγγίσεων Picard – Θεώρημα Τοπικής Ύπαρξης Picard–Lindelöf (Μονοσήμαντο):


05 03 2018 19 16 Office Lens

Συνθήκη Lipschit

(πιο χαλαρή συνθήκη για μονοσήμαντο):


08 03 2018 12 47 Office Lens


08 03 2018 12 49 Office Lens


Note:

  • Υπάρχει ενδεχόμενο η μοναδικότητα να χαθεί αν η \( f \) δεν ικανοποιεί τη συνθήκη του Lipschitz, παράδειγμα: \( y'=|y|^{1/2} \).
  • Η απαίτηση της κανονικότητας επηρεάζει και το πεδίο ορισμού της λύσης για συγκεκριμένη Α.Σ..

Γραμμικές Εξισώσεις με Σταθερούς Συντελεστές

Ορισμός:


16 03 2018 17 07 Office Lens_1

Περιεχόμενα.

Κεφάλαιο 1: Εξισώσεις Πρώτης Τάξης #

Κεφάλαιο 2: Τα Βασικά Θεωρήματα #

Κεφάλαιο 3: Ποιοτική Θεωρία #

Κεφάλαιο 4: Γραμμικές Εξισώσεις με Σταθερούς Συντελεστές #

Κεφάλαιο 4:

Ομογενής και μη ομογενής εξισώσεις:


16 03 2018 17 07 Office Lens_2

Κεφάλαιο 1:

Κεφάλαιο 2:

Η Ομογενής Εξίσωση Τάξης 2:


16 03 2018 17 47 Office Lens


Note: Λύση επίσης είναι και η \( φ(t) = 0, \; \forall t \).

Χαρακτηριστικό Πολυώνυμο:


\( L(e^{r t}) = (r^2 + a_1 r + a_2) e^{r t} = 0 \; , \)


\( e^{r t} \ne 0 \Rightarrow r^2 + a_1 r + a_2 = 0 \Rightarrow \)


Θέτουμε τότε ως χαρακτηριστικό πολυώνυμο:\( p(r) \equiv r^2 + a_1 r + a_2 \)


\( \Rightarrow L(e^{r t}) = p(r) e^{r t} \)

Ορισμοί

Εξισώσεις Χωριζομένων Μεταβλητών #

Γραμμικές Εξισώσεις #

Περιοδικές Εξισώσεις #

Εξισώσεις Bernoullί Και Rίccatί #

Εξισώσεις Ομογενείς και Αναγόμενες σε Ομογενείς #

Εξισώσεις Δεύτερης Τάξης #

Ακριβείς Εξισώσεις #

Συστήματα #

Εξίσωση Bernoulli:


\( y' + a(t) y = b(t) y^r \), όπου \( a(t), b(t) \) πραγματικές συναρτήσεις.

  • Για \( r = 0,1 \): γραμμική.
  • Για \( r > 0 \): \( y(t) = 0 \) λύση.
  • Για \( r \ne 0,1 \): μη γραμμική. Με τον μετασχηματισμό όμως \( u = y^{1-r} \) και πολλαπλασιάζοντας με \( (1-r) y^{-r} \) μπορεί να λυθεί ως γραμμική.

Εξίσωση Riccati:


\( y' + p(t) y + q(t) y^2 = f(t) \), όπου \( p(t), q(t), g(t) \) πραγματικές συναρτήσεις.


Με αλλαγή μεταβλητής : \( y(t) = y_1(t) + u(t) \), όπου \( y_1(t) \) ειδική λύση,
μετατρέπεται στην εξίσωση Bernoulli.


Note: \( \Rightarrow u(t) = y(t) - y_1(t) \), σκοπός είναι η απαλοιφή του όρου \( y(t) \).

Η Άλλη Αντίληψη #

Μέθοδοι Picard και Euler #

Ανισότητα του Gronwall #

Συνεχής Εξάρτηση

Το Θεώρημα της Πεπλεγμένης Συνάρτηση

Σύγκλιση Προσεγγίσεων Picard - Θεώρημα Τοπικής Ύπαρξης Picard–Lindelöf #

Η άλλη αντίληψη


Επίλυση ενός Π.Α.Τ.:


\( \frac{dy}{dt} = f(t,y) \; , \; y(t_0) = y_0 \)


Η μέθοδος των διαδοχικών προσεγγίσεων: η λύση ως όριο γνωστών συναρτήσεων,


\( y(t) = \lim_{n \to \infty} y_n(t) \)

Ποιοτική Θεωρία

Κεφάλαιο 3:

Πληθυσμιακά Μοντέλα

Διαγράμματα Φάσης

Γραμμικοποιήση

Δυναμικά Συστήματα

Το Διανυσματικό Πεδίο

Διακλαδώσεις

Η Ομογενής Εξίσωση Τάξης 2 #

Γραμμική Ανεξαρτησία #

Η Μη Ομογενής Εξίσωση Τάξης 2 #

Μέθοδος Lagrange #

Η Ομογενής Εξίσωση Τάξης n #

Η Μη Ομογενής Εξίσωση Τάξης n #

Μέθοδος Απροσδιόριστων Συντελεστών

Μηχανικές Ταλαντώσεις

Εισαγωγή #

Πεδίο Ορισμού Π.Α.Τ.:


19 03 2018 14 11 Office Lens

Βασικό Παράδειγμα:


19 03 2018 14 58 Office Lens

Ανισότητα του Gronwall:


19 03 2018 15 01 Office Lens

Γραμμική Ανεξαρτησία:


19 03 2018 16 41 Office Lens


19 03 2018 16 52 Office Lens

Απόδειξη λύσης για \( r_1 = r_2 \):


19 03 2018 16 29 Office Lens

Κατασκευή του Χώρου των Λύσεων:


19 03 2018 17 18 Office Lens


19 03 2018 17 22 Office Lens_1

15 03 2018 14 28 Office Lens

15 03 2018 14 42 Office Lens

15 03 2018 14 45 Office Lens

Ιδιότητα της Ορίζουσας Wronski 2 λύσεων:


19 03 2018 17 22 Office Lens

Η Ομογενής Εξίσωση Τάξης n:


\( L(y) = y^{(n)} + a_1 y^{(n-1)} + a_2 y^{(n-2)} + ... + a_n y = 0 \)


Για εκθετική λύση της μορφής: \( e^{r t} \) \( \Rightarrow L(e^{r t}) = p(r) e^{r t} \), όπου
\( p(r) \) το χαρακτηριστικό πολυώνυμο της εξίσωσης \( L(y) = 0 \).

Λύσεις:


20 03 2018 12 33 Office Lens


Proof:

Γραμμική Ανεξαρτησία:


20 03 2018 12 33 Office Lens

Πρόβλημα Αρχικών Τιμών:


20 03 2018 13 40 Office Lens

Η Μη Ομογενής Εξίσωση Τάξης 2:


Μέθοδος Lagrange:


20 03 2018 15 08 Office Lens

Cauchy–Euler Equation:


Η Μη Ομογενής Εξίσωση Τάξης n

Κεφάλαιο 5: Μέθοδος Δυναμοσειρών #

Κεφάλαιο 6: Γραμμικά Συστήματα #

Κεφάλαιο 7: Εξισώσεις Διαφορών

Κεφάλαιο 8: Μετασχηματισμός Laplace

Κεφάλαιο 9: Αρχή Μεγίστου και Θεωρία Sturm

Κεφάλαιο 10: Επίπεδο Φάσης

Μέθοδος Δυναμοσειρών

Κεφάλαιο 5

5.1 Δυναμοσειρές #

5.2 Ομαλά σημεία #

5.3 Εξίσωση Legendre #

5.4 Κανονικό Ιδιάζον Σημείο #

5.5 Εξίσωση Bessel #

Δυναμοσειρές

Ομαλά σημεία

Εξίσωση Legendre

Κανονικό Ιδιάζον Σημείο

Εξίσωση Bessel

Πρόβλημα Αρχικών Τιμών (Π.Α.Τ.): Η διαφορική εξίσωση μαζί με μία αρχική συνθήκη.

Εξίσωση Legendre τάξη α:


\( (1 - t^2) y'' - 2ty'+ α(α+1) y = 0\), όπου α μία πραγματική σταθερά


3 25 2018 18 00 Office Lens

Check: Σύγκλιση σειράς!

Λύση:


Το \( t=0 \) είναι ομαλό σημείο της εξίσωσης Legendre και συνεπώς επιδιώκοντας λύση της μορφής:
\( \sum_{m=0}^\infty a_m t^m\),
καταλήγουμε:
3 25 2018 18 07 Office Lens
Οι σειρές συγκλίνουν για \( |t|<1 \).

Γραμμικά Συστήματα

Κεφάλαιο 6

6.1 Ορισμοί, Ύπαρξη-Μονοσήμαντο #

6.2 Ομογενείς γραμμικές εξισώσεις #

6.3 Ο τύπος της μεταβολής των παραμέτρων #

6.4 Συναρτήσεις Green #

6.5 Γραμμικές εξισώσεις τάξης n #

6.6 Σταθεροί συντελεστές Ι: Πίνακες διαγωνοποιήσιμοι #

6.7 Σταθεροί συντελεστές ΙΙ: Πίνακες μη απλής δομής

6.8 Περιοδικά γραμμικά συστήματα

Ορισμοί, Ύπαρξη-Μονοσήμαντο [Start]

Σύστημα Διαφορικών Εξισώσεων Πρώτης Τάξης:


28 03 2018 21 33 Office Lens

Μετατροπή διαφορικής εξίσωσης τάξης n λυμένης μορφής σε ένα σύστημα διαφορικών εξισώσεων πρώτης τάξης:


28 03 2018 21 33 Office Lens_1


28 03 2018 21 36 Office Lens


Μη-ομογενές και ομογενές σύστημα:


28 03 2018 23 24 Office Lens

Ορισμός διαφορίσιμης και ολοκληρώσιμης διανυσματικής συνάρτησης και πίνακα-συνάρτησης και ιδιότητες:


28 03 2018 21 36 Office Lens 1


28 03 2018 21 37 Office Lens

Ευκλείδια νόρμα της διανυσματικής συνάρτησης και του πίνακα-συνάρτησης:


28 03 2018 21 37 Office Lens_1


28 03 2018 23 04 Office Lens

Πρόβλημα Αρχικών Τιμών:


28 03 2018 23 04 Office Lens_1

Θεώρημα Ύπαρξης:


28 03 2018 23 10 Office Lens

Θεώρημα Μονοσήμαντου:


28 03 2018 23 17 Office Lens

Ομογενείς Γραμμικές Εξισώσεις:


3 29 2018 12 33 Office Lens

Το σύνολο των λύσεων της ομογενούς:


3 29 2018 12 33 Office Lens_1


Δομή διανυσματικού χώρου:
29 03 2018 12 38 Office Lens

Θεώρημα Διάστασης Χώρου Λύσεων:


29 03 2018 12 38 Office Lens_1

Θεμελιώδες Σύστημα Λύσεων:


29 03 2018 12 41 Office Lens

Θεμελιώδης Πίνακας Λύσεων:


29 03 2018 12 41 Office Lens_1

Ιδιότητες Πίνακα Λύσεων:


29 03 2018 12 41 Office Lens_2

Ορίζουσα Wronski του Πίνακα Λύσεων:


29 03 2018 12 41 Office Lens_3

Θεώρημα Γραμμικής Ανεξαρτησίας Λύσεων:


29 03 2018 12 41 Office Lens_4

Πόρισμα: Μηδενική Ορίζουσα Wronski:


29 03 2018 23 27 Office Lens

Μεταγραφή Λύσεως Π.Α.Τ.


30 03 2018 00 40 Office Lens

Θεώρημα Μετασχηματισμού Πίνακα Λύσεων


30 03 2018 00 45 Office Lens

Ιδιότητες Κύριου Πίνακα:


31 03 2018 22 27 Office Lens

Ορισμός Κύριου Πίνακα:


31 03 2018 22 27 Office Lens_1

Γενική Λύση του Π.Α.Τ.:


31 03 2018 22 52 Office Lens

Τύπος του Liouville:


31 03 2018 22 53 Office Lens


31 03 2018 23 32 Office Lens #

Γεωμετρική Ερμηνεία του Τύπου του Liouville:


31 03 2018 23 32 Office Lens_1

Ο Τύπος του Liouville για Θεμελιώδη Πίνακα


31 03 2018 23 32 Office Lens_2 #

Ο Τύπος της Μεταβολής των Παραμέτρων

Η Γενική Περίπτωση

Η Εκθετική Συνάρτηση \( e^{At} \)

Το Συζυγές Πρόβλημα, Εφαρμογές στα Περιοδικά Συστήματα

Π.Α.Τ. για Μη Ομογενές Σύστημα:


01 04 2018 19 04 Office Lens

Ο Τύπος της Μεταβολής των Παραμέτρων


01 04 2018 19 04 Office Lens_1


01 04 2018 19 06 Office Lens
#

Ο Χώρος Λύσεων της Ομογενούς και Μη Ομογενούς


01 04 2018 19 27 Office Lens


01 04 2018 19 07 Office Lens

Λύση Π.Α.Τ. για Α Πίνακα Σταθερό (Με Μέθοδο Μεταβολής των Παραμέτρων):


02 04 2018 16 05 Office Lens_2


02 04 2018 16 03 Office Lens_1

Ορισμός Πίνακα \( e^{A t} \):


02 04 2018 16 03 Office Lens

Ιδιότητες Πίνακα \( e^{A t} \):


02 04 2018 16 02 Office Lens

Ορισμός Συζυγούς Εξίσωσης:


02 04 2018 16 31 Office Lens


02 04 2018 16 47 Office Lens

Λύση Π.Α.Τ. για Α Πίνακα Σταθερό (Με Μέθοδο του Ολοκληρωτικού Παράγοντα):


02 04 2018 17 53 Office Lens


02 04 2018 16 31 Office Lens_1

Το Κριτήριο Εναλλαγής του Fredholm:


02 04 2018 16 47 Office Lens_1

Συναρτήσεις Green

Γενική Μέθοδος


03 04 2018 00 13 Office Lens


03 04 2018 00 14 Office Lens

Παράδειγμα Προβλήματος Συνοριακών Συνθηκών:


03 04 2018 00 16 Office Lens_1


03 04 2018 00 17 Office Lens

Ενδεικτικό Παράδειγμα:


03 04 2018 00 14 Office Lens_1


03 04 2018 00 16 Office Lens_2

Γραμμικές Εξισώσεις Τάξης n

Σταθεροί συντελεστές Ι: Πίνακες διαγωνοποιήσιμοι