Λογισμός των μεταβολών: Διαφορά μεταξύ των αναθεωρήσεων
Περιεχόμενο που διαγράφηκε Περιεχόμενο που προστέθηκε
μ Αντικατάσταση παρωχημένου προτύπου με references tag |
μ συναρτησιακό Ετικέτες: Οπτική επεξεργασία Επεξεργασία από κινητό Διαδικτυακή επεξεργασία από κινητό |
||
Γραμμή 1:
Ο '''λογισμός των μεταβολών''' ή '''μεταβολικός λογισμός''' είναι κλάδος της [[Μαθηματική ανάλυση|μαθηματικής ανάλυσης]] που ασχολείται με τη μεγιστοποίηση ή την ελαχιστοποίηση συναρτησιακών ή [[Συναρτησοειδές|συναρτησοειδών,]]
== Εισαγωγικά ==
Γραμμή 12:
== Ακρότατα ==
Ο λογισμός των μεταβολών ασχολείται με μέγιστα ή ελάχιστα
Τόσο τα ισχυρά όσο και τα ασθενή ακρότατα ενός
== Εξίσωση Euler-Lagrange ==
{{main|Εξίσωση Όιλερ-Λαγκράνζ}}
Η εύρεση των ακρότατων των συναρτησοειδών είναι παρόμοια με την εύρεση των μέγιστων και ελάχιστων μιας συναρτήσης. Για να ορίσουμε τα μέγιστα και τα ελάχιστα μίας συνάρτησης χρειάζεται να βρούμε τα σημεία που μηδενίζουν την παράγωγό της. Ανάλογα, το ακρότατο του συναρτησοειδούς μπορεί να βρεθεί με την κατασκευή μίας συνάρτησης όπου η [[συναρτησιακή παράγωγος]] είναι ίση με 0. Αυτό μας οδηγεί στην επίλυση της αντίστοιχης εξίσωσης Euler-Lagrange.<ref group="Σημείωση">Η ακόλουθη απόδειξη της εξίσωσης Euler–Lagrange αντιστοιχεί σε αυτή των σελ. 184–5 του:<br>{{cite book | authors = [[Richard Courant|Courant, R.]], [[David Hilbert|Hilbert, D.]] | title = Methods of Mathematical Physics | volume = Vol. I | edition = πρώτη αγγλική | publisher = Interscience Publishers, Inc | year = 1953 | location = New York | page = | accessdate = | isbn = 978-0471504474}}</ref>▼
▲Η εύρεση των ακρότατων των
Ας θεωρήσουμε το συναρτησοειδές▼
: <math> J[y] = \int_{x_1}^{x_2} L[x,y(x),y'(x)]\, dx \, ,</math>
όπου
Γραμμή 27 ⟶ 28 :
: {{math|''y'' ′(''x'') {{=}} ''dy / dx'' ,}}
: {{math|''L''[''x'', ''y'' (''x''), ''y'' ′(''x'')]}} είναι δύο φορές συνεχώς (ολικά) παραγωγίσημη ως προς τα {{math|''x'', ''y'', ''y'' ′ .}}
: Αν το
: <math>J[f] \le J[f + \varepsilon \eta] \, .</math>
: Ο όρος {{math|''εη''}} ονομάζεται μεταβολή της συνάρτησης {{math|''f''}} και συμβολίζεται με {{math|''δf'' .}}<ref name="CourHilb1953P1844">{{harvnb|Courant|Hilbert|1953|p = 184}}</ref>
Αντικαθιστώντας όπου {{math|''y''}} το {{math|''f'' + ''εη''}} στο
:: <math> \Phi(\varepsilon) = J[f+\varepsilon\eta] \, . </math>
Δεδομένου ότι το
:: <math> \Phi'(0) \equiv \left.\frac{d\Phi}{d\varepsilon}\right|_{\varepsilon = 0} = \int_{x_1}^{x_2} \left.\frac{dL}{d\varepsilon}\right|_{\varepsilon = 0} dx = 0 \, . </math>
Λαμβάνοντας την [[Ολική παράγωγος|ολική παράγωγο]] της {{math|''L''[''x'', ''y'', ''y'' ′] ,}} όπου οι {{math|''y'' {{=}} ''f'' + ''ε η''}} και {{math|''y'' ′ {{=}} ''f'' ′ + ''ε η''′}} είναι συναρτήσεις των {{math|''ε''}} αλλά το {{math|''x''}} δεν είναι,
Γραμμή 55 ⟶ 58 :
Σύμφωνα με το [[θεμελειώδες λήμμα του λογισμού των μεταβολών]], ο όρος στην παρένθεση είναι ίσος με το μηδέν, δηλαδή
::: <div style="{{divstylewhite}}; width:15em; margin:.3em"><center><math> \frac{\part L}{\part f} -\frac{d}{dx} \frac{\part L}{\part f'}=0 </math></center></div>
το οποίο λέγεται '''εξίσωση Euler–Lagrange.''' Το αριστερό μέλος της εξίσωσης ονομάζεται [[συναρτησιακή παράγωγος|συναρτησιακή]] ή [[συναρτησιακή παράγωγος|μεταβολική παράγωγος]] της {{math|''J''[''f'']}} και συμβολίζεται με {{math|''δJ''/''δf''(''x'') .}}
Σε γενικές γραμμές αυτό δίνει μία δεύτερης τάξης [[Συνήθης διαφορική εξίσωση|συνήθη διαφορική εξίσωση]] η οποία μπορεί να λυθεί για να πάρουμε την ακρότατη συνάρτηση {{math| ''f''(''x'').}} Η εξίσωση των Euler–Lagrange είναι αναγκαία, αλλά όχι ικανή συνθήκη για να είναι το {{math|''J''[''f'']}} ακρότατο. Μία ικανή συνθήκη για το ελάχιστο δίνεται στην ενότητα [[Λογισμός των μεταβολών#Μεταβολές και επαρκείς συνθήκες για ελάχιστο|Μεταβολές και επαρκείς συνθήκες για ελάχιστο.]]
Γραμμή 64 ⟶ 67 :
με
: <math> y\,'(x) = \frac{dy}{dx} \, , \ \ y_1=f(x_1) \, , \ \ y_2=f(x_2) \, . </math>
Η εξίσωση των Euler–Lagrange θα χρησιμοποιηθεί τώρα για την εύρεση της ακρότατης συνάρτησης {{math|''f'' (''x'')}} η οποία ελαχιστοποιεί
: <math> \frac{\part L}{\part f} -\frac{d}{dx} \frac{\part L}{\part f'}=0 </math>
με
Γραμμή 73 ⟶ 76 :
Αντικαθιστώντας το {{math|''L''}} και παίρνοντας τη μερική παράγωγο έχουμε,
: <math> \frac{d}{dx} \ \frac{ f'(x) } {\sqrt{1 + [ f'(x) ]^2}} \ = 0 \, . </math>
: <math>\frac{d^2 f}{dx^2}\ \cdot\ \frac{1}{\left[\sqrt{1+[f'(x)]^2}\ \right]^3} = 0 \, , </math>
αλλά αφού το {{math|1+[''f ′''(''x'')]<sup>2</sup>}} είναι μη-μηδενικό,
Γραμμή 79 ⟶ 83 :
απ' όπου συνάγουμε ότι η βραχύτερη καμπύλη που συνδέει δύο σημεία {{math|(''x''<sub>1</sub>, ''y''<sub>1</sub>)}} και {{math|(''x''<sub>2</sub>, ''y''<sub>2</sub>) }} είναι η
::: <math> f(x) = m x + b, \qquad \mu \varepsilon \ \ m = \frac{y_2 - y_1}{x_2 - x_1} \quad \kappa \alpha \iota \quad b = \frac{x_2 y_1 - x_1 y_2}{x_2 - x_1} ,</math>
και έχουμε τότε βρει την ακρότατη συνάρτηση {{math|''f''(''x'')}} που ελαχιστοποιεί το
== Ταυτότητα Beltrami ==
Γραμμή 94 ⟶ 98 :
== Φαινόμενο Lavrentiev ==
O πρώτος που έδωσε καλές συνθήκες για να έχουν οι εξισώσεις Euler-Lagrange στάσιμη λύση ήταν ο [[Ντάβιντ Χίλμπερτ|Χίλμπερτ.]] Μέσα σε μια κυρτή επιφάνεια και μία θετική τρεις φορές παραγωγίσιμη
Ωστόσο [[Μιχαήλ Λαβρέντιεφ|Lavrentiev]] το 1926 έδειξε ότι υπάρχουν περιπτώσεις όπου δεν υπάρχει βέλτιστη λύση αλλά μια από αυτές μπορεί να προσεγγιστεί με αυθαίρετη ακρίβεια αυξάνοντας τον αριθμό των τομών. Για παράδειγμα η ακόλουθη:<dl>
Γραμμή 111 ⟶ 115 :
Συνήθως φτάνει να θεωρήσουμε μικρές μετατοπίσεις της μεμβράνης, της οποίας η ενεργειακή διαφορά από την αρχική της θέση δίνεται προσεγγιστικά από τον τύπο
:<math> V[\varphi] = \frac{1}{2}\iint_D \nabla \varphi \cdot \nabla \varphi \, dx\, dy.\,</math>
Ζητείται η ελαχιστοποίηση του
:<math> \frac{d}{d\varepsilon} V[u + \varepsilon v]|_{\varepsilon=0} = \iint_D \nabla u \cdot \nabla v \, dx\,dy = 0.\,</math>
Δεδομένου ότι η <math display="inline">u</math> είναι δύο φορές παραγωγίσιμη, μπορούμε να εφαρμόσουμε το θεώρημα της απόκλισης απ' όπου παίρνουμε την
Γραμμή 141 ⟶ 146 :
όπου <math display="inline">\lambda</math> είναι ο λόγος
:<math> \lambda = \frac{Q[u]}{R[u]}. \,</math>
Αποδεικνύεται (δείτε [http://www.amazon.com/Calculus-Variations-Dover-Books-Mathematics/dp/0486414485 Gelfand και Fomin 1963]) ότι η
Η αμέσως μικρότερη ιδιοτιμή και ιδιοσυνάρτηση μπορεί να βρεθεί ελαχιστοποιώντας το <math display="inline">Q</math> με τον επιπλέον περιορισμό
Γραμμή 147 ⟶ 152 :
Η παρούσα διαδικασία μπορεί να επαναληφθεί για να πάρουμε μία πλήρη ακολουθία ιδιοτιμών και ιδιοσυναρτήσεων για το πρόβλημα.
Στο ίδιο πρόβλημα μεταβολών καταλήγουμε επίσης και όταν έχουμε πιο γενικές οριακές συνθήκες. Αντί να απαιτούμε το μηδενισμό της <math display="inline">\varphi</math> στα όρια ολοκλήρωσης μπορούμε να μην
:<math>Q[\varphi] = \int_{x_1}^{x_2} \left[ p(x) \varphi'(x)^2 + q(x)\varphi(x)^2 \right] \, dx + a_1 \varphi(x_1)^2 + a_2 \varphi(x_2)^2, \,</math>
όπου τα <math>a_1</math> και <math>a_2</math> είναι επιλεγμένα αυθαίρετα. Θέτοντας <math display="inline"> \varphi = u + \varepsilon v </math> η πρώτη μεταβολή του πηλίκου <math display="inline"> Q/R</math> είναι
Γραμμή 169 ⟶ 174 :
όπου
:<math> \lambda = \frac{Q[u]}{R[u]}.\,</math>
Η
:<math> p(S) \frac{\part u}{\part n} + \sigma(S) u =0</math>
στο σύνορο <math display="inline"> B.</math> Αυτό το αποτέλεσμα βασίζεται στη θεωρία κανονικότητας για τις ελλειπτικές μερικές διαφορικές εξισώσεις· δείτε Jost και Li-Jost (1998) για λεπτομέρειες. Πολλές επεκτάσεις, συμπεριλαμβανομένων αποτελεσμάτων πληρότητας, ασυμπτωτικών ιδιοτήτων των ιδιοτιμών και αποτελεσμάτων που αφορούν τους κόμβους ιδιοσυναρτήσεων, βρίσκονται στο σύγγραμμα των Courant και Hilbert (1953).
Γραμμή 194 ⟶ 200 :
<math> -\frac{d}{dx} \left[\frac{ n(x,f_0) f_0'}{\sqrt{1 + f_0'^2}} \right] + n_y (x,f_0) \sqrt{1 + f_0'(x)^2} =0. \,</math>
Ολοκληρώνοντας την παραπάνω εξίσωση μπορούν να υπολογιστούν οι ακτίνες φωτός. Αυτός ο φορμαλισμός χρησιμοποιείται στο πλαίσιο της [[Λαγκατζιανής οπτικής|Λαγρανζιανής οπτικής]], καθώς και της [[Χαμιλτονιανής οπτικής]].
==== Ο νόμος του Σνελ ====
Γραμμή 229 ⟶ 235 :
: <math> A[C] = \int_{t=t_0}^{t_1} P \cdot \dot X \, dt.\,</math>
Από την τελευταία μορφή έπεται ότι αν μπορέσουμε να βρούμε συνάρτηση <math display="inline">\psi (X)</math> τέτοια ώστε <math>\nabla \psi = P,</math> τότε το ολοκλήρωμα {{math|''A''}} δίνεται από τη διαφορά <math>\psi (X(t_1)) - \psi (X(t_0)).</math> Επομένως, το πρόβλημα της μελέτης των καμπυλών για τις οποίες το ολοκλήρωμα είναι στάσιμο συνδεέται με τη μελέτη των [[Σταθμικό σύνολο|σταθμικών επιφανειών]] (level surfaces) της <math display="inline">\psi .</math> Για την εύρεση μίας τέτοιας συνάρτησης θα χρησιμοποιήσουμε την κυματική εξίσωση, που περιγράφει τη διάδοση του φωτός. Αυτός ο φορμαλισμός χρησιμοποιείται ως γενικό πλαίσιο στη [[Λαγκατζιανή οπτική|Λαγρανζιανή οπτική]] και στην [[Χαμιλτονιανή οπτική]].
===== Συσχέτιση της <math>\psi</math> με την κυματική εξίσωση =====
Γραμμή 265 ⟶ 271 :
Κυρίως Άρθρο: [[Δράση (φυσική)]]
Στην κλασσική μηχανική, η δράση, ''S'', ορίζεται ως η χρονική ολοκλήρωση της
<math> L = T - U, \,</math>
όπου ''T'' είναι η [[κινητική ενέργεια]] ενός μηχανικού συστήματος και ''U'' η [[δυναμική του ενέργεια]]. Η [[αρχή του Hamilton]] (ή η αρχη της ελάχιστης δράσης) υποστηρίζει
:<math> S = \int_{t=t_0}^{t_1} L(x, \dot x, t) dt \,</math>
είναι στάσιμο σε σχέση με τις μεταβολές της διαδρομής ''x(t)''. Οι εξισώσεις Euler–Lagrange αυτού του συστήματος είναι γνωστές ως εξισώσεις Lagrange:
Γραμμή 283 ⟶ 289 :
<math> p = m \dot x. \,</math>
Η [[μηχανική του Hamilton]] προκύπτει εαν θέσουμε τις συζυγείς ορμές στη θέση του <math>\dot x</math>, και η
:<math> H(x, p, t) = p \,\dot x - L(x,\dot x, t).\,</math>
Η Χαμιλτονιανή είναι η συνολική ενέργεια του συστήματος: ''H'' = ''T'' + ''U''. Η αναλογία με τις αρχές του Fermat ορίζει ότι οι λύσεις των
:<math> \frac{\part \psi}{\part t} + H\left(x,\frac{\part \psi}{\part x},t\right) =0.\,</math>
== Μεταβολές και επαρκείς συνθήκες για ελάχιστο ==
Ο Λογισμός των μεταβολών ασχολείται με τις μεταβολές των
Για παράδειγμα, αν ''J''[''y''] είναι ένα
<math> \Delta J[h] = J[y+h] - J[y] </math>.<ref group="Σημείωση">Σημειωτέον ότι το Δ ''J''[''h''] και οι μεταταβολές παρακάτω,εξαρτώνται από τα ''y'' και ''h''. To ''y'' έχει εξαιρεθεί για να απλοποιηθεί ο συμβολισμός. Για παράδειγμα, Δ ''J''[''h''] θα μπορούσε να είχε γραφεί Δ ''J''[''y'' ; ''h''].</ref>
Το
: <math> \Delta J[h] = \phi [h] + \epsilon \|h\| </math>,
όπου ''φ''[''h''] είναι ένα γραμμικό
: <math>\delta J[h] = \phi(h) </math>.
Το
: <math> \Delta J[h] = \phi_1 [h] + \phi_2 [h] + \epsilon \|h\|^2 </math>,
όπου ''φ<sub>1</sub>''[''h''] είναι ένα γραμμικό
: <math>\delta^2 J[h] = \phi_2(h) </math>.
Η δεύτερη μεταβολή ''δ''<sup>2</sup>''J''[''h''] λέγεται '''αυστηρά θετική''' αν
Γραμμή 309 ⟶ 316 :
Επαρκής συνθήκη για ένα ελάχιστο''':'''
Το
== Δείτε επίσης ==
|