Πεπερασμένη διαφορά: Διαφορά μεταξύ των αναθεωρήσεων
Χωρίς σύνοψη επεξεργασίας |
μ Διόρθωση συνδέσμων |
||
(36 ενδιάμεσες εκδόσεις από 6 χρήστες δεν εμφανίζονται) | |||
Γραμμή 1: | Γραμμή 1: | ||
Η '''πεπερασμένη διαφορά''' είναι μία μαθηματική έκφραση της μορφής {{math|''f''(''x'' + ''b'') − ''f''(''x'' + ''a'')}}. Αν η πεπερασμένη διαφορά διαιρείται από το {{math|''b'' − ''a''}}, τότε αυτή λαμβάνει μόνο ένα [[πεπερασμένο πηλίκο]].Ο υπολογισμός των [[παραγώγων]]από τις πεπερασμένες διαφορές παίζει βασικό ρόλο στις [[μεθόδους πεπερασμένης διαφοράς]] για την [[αριθμητική]] λύση των [[διαφορικών εξισώσεων]], ειδικότερα στα [[προβλήματα συνοριακών τιμών]]. |
Η '''πεπερασμένη διαφορά''' είναι μία μαθηματική έκφραση της μορφής {{math|''f''(''x'' + ''b'') − ''f''(''x'' + ''a'')}}. Αν η πεπερασμένη διαφορά διαιρείται από το {{math|''b'' − ''a''}}, τότε αυτή λαμβάνει μόνο ένα [[πεπερασμένο πηλίκο]].Ο υπολογισμός των [[παραγώγων]] από τις πεπερασμένες διαφορές παίζει βασικό ρόλο στις [[μεθόδους πεπερασμένης διαφοράς]] για την [[αριθμητική]] λύση των [[διαφορικών εξισώσεων]], ειδικότερα στα [[προβλήματα συνοριακών τιμών]]. |
||
[[Ορισμένες επαναληπτικές σχέσεις]] μπορούν να γραφούν ως '''εξισώσεις διαφοράς''' αντικαθιστώντας την προσέγγιση της σημειογραφείας με τις πεπερασμένες διαφορές. |
|||
Σήμερα, ο όρος "πεπερασμένη διαφορά" συχνά λαμβάνεται ως συνώνυμο με [[τους υπολογισμούς της πεπερασμένης διαφοράς]] των παραγώγων, ειδικότερα στο πλαίσιο των [[αριθμητικών μεθόδων]].<ref name="WilmottHowison1995">{{cite book|author1=Paul Wilmott|author2=Sam Howison|author3=Jeff Dewynne|title=The Mathematics of Financial Derivatives: A Student Introduction|year=1995|publisher=Cambridge University Press|isbn=978-0-521-49789-3|page=[https://archive.org/details/mathematicsoffin00wilm/page/137 137]|url=https://archive.org/details/mathematicsoffin00wilm/page/137}}</ref><ref name="Olver2013">{{cite book|author=Peter Olver|author-link=Peter J. Olver|title=Introduction to Partial Differential Equations|year=2013|publisher=Springer Science & Business Media|isbn=978-3-319-02099-0|page=182}}</ref><ref name="Chaudhry2007">{{cite book|author=M Hanif Chaudhry|title=Open-Channel Flow|year=2007|publisher=Springer|isbn=978-0-387-68648-6|pages=369}}</ref> Οι υπολογισμοί των πεπερασμένων διαφορών είναι πηλίκα της πεπερασμένης διαφοράς στην ορολογία που χρησιμοποιήθηκε ανωτέρω. |
|||
Οι πεπερασμένες διαφορές,επίσης,έχουν γίνει το θέμα μελέτης ως αφηρημένα αυτοτελή μαθηματικά αντικείμενα, π.χ. σε έργα των [[Τζορτζ Μπουλ]] (1860), [[L. M. Milne-Thomson]] (1933), και [[Károly Jordan]] (1939), εντοπίζοντας τις προελέυσεις των αλγορίθμων του [[Jost Bürgi]] (ca. 1592) και άλλων συμπεριλαμβανομένου του [[Ισαάκ Νεύτων]]. Σε αυτό το σημείο, ο λογισμός (που ασχολείται με τα όρια) των πεπερασμένων τιμών είναι εναλλακτικός του λογισμού των [[απειροελάχιστων]].<ref>Jordán, op. cit., p. 1 and Milne-Thomson, p. xxi. |
|||
Milne-Thomson, Louis Melville (2000): ''The Calculus of Finite Differences'' (Chelsea Pub Co, 2000) ISBN 978-0821821077</ref> |
|||
==Προς τα εμπρός, προς τα πίσω, και κεντρική διαφορά== |
|||
Τρεις μορφές μελετώνται ευρέως: οι προς τα εμπρός,οι προς τα πίσω και οι κεντρικές διαφορές.<ref name="WilmottHowison1995"/><ref name="Olver2013"/><ref name="Chaudhry2007"/> |
|||
Η '''προς τα εμπρός διαφορά (προς τα πάνω)''' είναι μια έκφραση της μορφής |
|||
:<math> \Delta_h[f](x) = f(x + h) - f(x). \ </math> |
|||
Σύμφωνα με την εφαρμογή, η τάξη ''h'' μπορεί να είναι μεταβλητή ή σταθερή. Όταν παραλείπεται, η ''h'' θεωρείται ότι είναι 1: <math>\Delta[f](x) = \Delta_1[f](x)</math>. |
|||
Η '''προς τα πίσω διαφορά (προς τα κάτω)''' χρησιμοποιεί τις τιμές της συνάρτησης στο ''x'' και ''x'' − ''h'', αντί των τιμών στο ''x'' + ''h'' and ''x'': |
|||
:<math> \nabla_h[f](x) = f(x) - f(x-h). \ </math> |
|||
{{anchor|central difference}} |
|||
Τελικά, η '''κεντρική διαφορά (κεντρώα)''' δίνεται από |
|||
:<math> \delta_h[f](x) = f(x+\tfrac12h)-f(x-\tfrac12h)~.</math> |
|||
== Σχέση με τις παραγώγους == |
|||
{{anchor|finite difference approximation}} |
|||
Η [[παράγωγος]] της συνάρτησης {{mvar|f}} στο σημείο {{mvar| x}} ορίζεται από το [[όριο]] |
|||
:<math> f'(x) = \lim_{h\to0} \frac{f(x+h) - f(x)}{h}. </math> |
|||
Αν η {{mvar|h}} λαμβάνει μια σταθερή (μη μηδενική) τιμή αντί να πλησιάζει το μηδέν, τότε η δεξιά πλευρά της παραπάνω εξίσωσης γράφεται |
|||
:<math> \frac{f(x + h) - f(x)}{h} = \frac{\Delta_h[f](x)}{h}. </math> |
|||
Συνεπώς, η προς τα πάνω διαφορά που διαιρείται από το {{mvar|h}} πλησιάζει την παράγωγο όταν το {{mvar|h}} είναι μικρό. Το σφάλμα σε αυτή την προσέγγιση μπορεί να βρεθεί από τον [[Τύπο του Taylor]]. Λαμβάνοντας υπόψη, ότι η {{mvar|f}} είναι διαφορίσιμη , έχουμε |
|||
:<math> \frac{\Delta_h[f](x)}{h} - f'(x) = O(h)\to 0 \quad \text{as }(h \to 0). </math> |
|||
Ο ίδιος τύπος αναμένεται για την προς τα πίσω διαφορά: |
|||
:<math> \frac{\nabla_h[f](x)}{h} - f'(x) = O(h)\to 0 \quad \text{as }(h \to 0). </math> |
|||
Ωστόσο, η κεντρική (ή αλλιώς κεντρώα) διαφορά αποδίδει μια πιο ακριβή προσέγγιση. Αν η {{mvar|f}} είναι δύο φορές διαφορίσιμη, |
|||
:<math> \frac{\delta_h[f](x)}{h} - f'(x) = O(h^2) . \!</math> |
|||
Το κύριο πρόβλημα με τη μέθοδο της κεντρικής διαφοράς, ωστόσο , είναι ότι οι συναρτήσεις που ταλαντεύονται αποδίδουν μηδενική παράγωγο. Αν {{math|''f''(''nh''){{=}}1}} για {{mvar|n}} περιττό, και {{math|''f''(''nh''){{=}}2}} για {{mvar|n}} άρτιο, αλλιώς {{math|''f ' '' (''nh''){{=}}0}} αν υπολογίζεται με το συνδυασμό της κεντρικής διαφοράς. Αυτό είναι ιδιαίτερα ενοχλητικό αν ο κλάδος της {{mvar|f}} είναι διακεκριμένος. |
|||
Συγγραφείς, για τους οποίους οι πεπερασμένες διαφορές σημαίνουν προσεγγίσεις πεπερασμένων διαφορών, ορίζουν τις προς τα πάνω/προς τα κάτω/κεντρικές διαφορές ως πηλίκα που δίνονται σε αυτή την ενότητα. (αντί να χρησιμοποιούν τους ορισμούς που δόθηκαν στην προηγούμενη ενότητα.)<ref name="WilmottHowison1995"/><ref name="Olver2013"/><ref name="Chaudhry2007"/> |
|||
{{See also|Συμμετρική παράγωγος}} |
|||
==Μεγαλύτερης τάξης διαφορές== <!-- Το τμήμα αυτό συνδέεται με τα περαιτέρω που γράφονται στη συνέχεια του άρθρου --> |
|||
Με ανάλογο τρόπο,μπορεί κανείς να ανάγει τις προσεγγίσεις πεπερασμένης διαφοράς σε μεγαλύτερης τάξης παραγώγους και διαφορικούς φορείς. Για παράδειγμα, χρησιμοποιώντας τον παραπάνω τύπο κεντρικής διαφοράς για {{math|''f ' ''(''x''+''h''/2)}} και {{math|''f ' ''(''x''−''h''/2)}} και εφαρμόζοντας τον τύπο της κεντρικής διαφοράς για την παράγωγο της {{mvar| f '}} στο {{mvar|x}}, παίρνουμε την προσέγγιση της κεντρικής διαφοράς της δεύτερης παραγώγου της {{mvar|f}}: |
|||
'''Κεντρική 2ης τάξης''' |
|||
:<math> f''(x) \approx \frac{\delta_h^2[f](x)}{h^2} = \frac{f(x+h) - 2 f(x) + f(x-h)}{h^{2}} . </math> |
|||
Με παρόμοιο τρόπο μπορούμε να εφαρμόσουμε άλλους τύπους διαφοράς με έναν αναδρομικό τρόπο. |
|||
'''Προς τα εμπρός 2ης τάξης''' |
|||
:<math> f''(x) \approx \frac{\Delta_h^2[f](x)}{h^2} = \frac{f(x+2h) - 2 f(x+h) + f(x)}{h^{2}} . </math> |
|||
Γενικότερα , οι '''{{mvar|n}}-ης τάξης προς τα εμπρός, προς τα πίσω, και κεντρικές''' διαφορές δίνονται,αντίστοιχα, |
|||
'''Προς τα εμπρός''' |
|||
:<math>\Delta^n_h[f](x) = |
|||
\sum_{i = 0}^{n} (-1)^i \binom{n}{i} f(x + (n - i) h), |
|||
</math> |
|||
ή για h=1, |
|||
:<math>\Delta^n [f](x)= \sum_{k=0}^n\binom nk(-1)^{n-k}f(x + k)</math> |
|||
'''Προς τα πίσω''' |
|||
:<math>\nabla^n_h[f](x) = |
|||
\sum_{i = 0}^{n} (-1)^i \binom{n}{i} f(x - ih), |
|||
</math> |
|||
'''Κεντρική''' |
|||
:<math>\delta^n_h[f](x) = |
|||
\sum_{i = 0}^{n} (-1)^i \binom{n}{i} f\left(x + \left(\frac{n}{2} - i\right) h\right). |
|||
</math> |
|||
Αυτές οι εξισώσεις χρησιμοποιούν [[διωνυμικούς συντελεστές]] ύστερα από το σημείο άθροισης απεικονίζονται ως <math>\ \binom{n}{i}</math>. Κάθε σειρά από το [[Τρίγωνο του Pascal]] παρέχει τον συντελεστή για κάθε τιμή του i. |
|||
Παρατήρησε ότι η κεντρική διαφορά θα, για περιττό {{mvar|n}}, έχει {{mvar|h}} πολλαπλασιαζόμενο με μη ακέραιο. Αυτό είναι συχνά πρόβλημα επειδή ανέρχεται σε αλλαγή του διαστήματος διακριτοποίησης. Το πρόβλημα μπορεί να αντιμετωπιστεί παίρνοντας το μέσο όρο του <math>\delta^n[f](x - h/2)</math> και <math>\delta^n[f](x + h/2)</math>. |
|||
Οι προς τα εμπρός διαφορές που εφαρμόζονται σε μία [[ακολουθία]] ονομάζονται μερικές φορές [[διωνυμικός μετασχηματισμός]] της ακολουθίας, και έχουν μια σειρά από ενδιαφέρουσες συνδυαστικές ιδιότητες. |
|||
Οι προς τα εμπρός διαφορές μπορούν να εκτιμηθούν χρησιμοποιώντας το [[Nörlund–Rice ολοκλήρωμα]]. Η αναπαράσταση του ολοκληρώματος για αυτούς τους τύπους των σειρών είναι ενδιαφέρουσα, επειδή το ολοκλήρωμα μπορεί συχνά να εκτιμηθεί χρησιμοποιώντας [[ασυμπτωτική επέκταση]] ή τεχνικές [[απότομης καθόδου]] : αντίθετα, η προς τα εμπρός διαφορά των σειρών μπορεί να είναι εξαιρετικά δύσκολο να εκτιμηθεί αριθμητικά, επειδή οι διωνυμικοί συντελεστές αυξάνονται με ταχείς ρυθμούς για μεγάλα {{mvar|n}}. |
|||
Η σχέση των διαφορών μεγαλύτερης τάξης με τις αντίστοιχες παραγώγους είναι απλή, |
|||
:<math>\frac{d^n f}{d x^n}(x) = \frac{\Delta_h^n[f](x)}{h^n}+O(h) = \frac{\nabla_h^n[f](x)}{h^n}+O(h) = \frac{\delta_h^n[f](x)}{h^n} + O(h^2).</math> |
|||
Οι μεγαλύτερης τάξης διαφορές μπορούν να χρησιμοποιηθούν για να φτιάξουμε καλύτερες προσεγγίσεις. Όπως αναφέρθηκε παραπάνω, η πρώτης τάξης διαφορά προσεγγίζει την παράγωγο πρώτης τάξης σε έναν όρο της τάξης {{mvar|h}}. Ωστόσο, ο συνδυασμός |
|||
:<math> \frac{\Delta_h[f](x) - \frac12 \Delta_h^2[f](x)}{h} = - \frac{f(x+2h)-4f(x+h)+3f(x)}{2h} </math> |
|||
προσεγγίζει την ''f'''(''x'') σε έναν όρο της τάξης {{math|''h''<sup>2</sup>}}. Αυτό αποδεικνύεται αναπτύσσοντας την παραπάνω έκφραση στη [[Σειρά Τέιλορ]], ή χρησιμοποιώντας το λογισμό των πεπερασμένων διαφορών,που εξηγείται παρακάτω. |
|||
Αν χρειάζεται,η πεπερασμένη διαφορά μπορεί να κεντράρεται γύρω από οποιοδήποτε σημείο συνδυάζοντας τις προς τα πάνω,προς τα κάτω και κεντρικές διαφορές. |
|||
===Πυρήνες αυθαίρετου μεγέθους=== |
|||
Χρησιμοποιώντας λίγο την γραμμική άλγεβρα, μπορεί κανείς αρκετά εύκολα να φτιάξει προσεγγίσεις, της οποίας παίρνω δείγμα ένα αυθαίρετο αριθμό σημείων από τα αριστερά και έναν (ενδεχομένως διαφορετικό) αριθμό σημείων από τα δεξιά του κεντρικού σημείου,για κάθε τάξη της παραγώγου. Αυτό συνεπάγεται την επίλυση ενός γραμμικού συστήματος έτσι ώστε η επέκταση Τέιλορ του αθροίσματος των εν λόγω σημείων ,γύρω από το κεντρικό σημείο,να προσεγγίζει καλά την επέκταση Τέιλορ της επιθυμητής παραγώγου. |
|||
Αυτό είναι χρήσιμο για την διαφοροποίηση μιας συνάρτησης σε μια περιοχή,όπου,καθώς αυτή θα πλησιάζει την άκρη της περιοχής, πρέπει να λαμβάνεται δείγμα όλο και λιγότερων σημείων στη μία πλευρά. |
|||
Οι λεπτομέρειες συνοψίζονται σε αυτές τις [http://commons.wikimedia.org/wiki/File:FDnotes.djvu σημειώσεις]. |
|||
===Ιδιότητες=== |
|||
* Για όλα τα θετικά ''k'' and ''n'' |
|||
:<math>\Delta^n_{kh} (f, x) = \sum\limits_{i_1=0}^{k-1} \sum\limits_{i_2=0}^{k-1} \cdots \sum\limits_{i_n=0}^{k-1} \Delta^n_h (f, x+i_1h+i_2h+\cdots+i_nh).</math> |
|||
* [[Κανόνας Leibniz (generalized product rule)|Κανόνας Leibniz]]: |
|||
:<math>\Delta^n_h (fg, x) = \sum\limits_{k=0}^n \binom{n}{k} \Delta^k_h (f, x) \Delta^{n-k}_h(g, x+kh).</math> |
|||
==Μέθοδοι πεπερασμένης διαφοράς== |
|||
Μία σημαντική εφαρμογή των πεπερασμένων διαφορών είναι στην [[αριθμητική ανάλυση]], ειδικότερα στις [[αριθμητικές διαφορικές εξισώσεις]], οι οποίες στοχεύουν στην αριθμητική επίλυση των [[συνήθων]] και [[μερικών διαφορικών εξισώσεων]] αντίστοιχα. Η ιδέα είναι να αντικατασταθούν οι παράγωγοι που εμφανίζονται στην διαφορική εξίσωση από τις πεπερασμένες διαφορές που τις προσεγγίζουν. Οι μέθοδοι που προκύπτουν ονομάζονται [[μέθοδοι πεπερασμένης διαφοράς]]. |
|||
Κοινές εφαρμογές της μεθόδου πεπερασμένης τιμής συναντώνται στην υπολογιστική επιστήμη και στην εφαρμοσμένη μηχανική, όπως [[θερμική μηχανική]], [[μηχανική των ρευστών]], κτλ. |
|||
==Σειρές του Νεύτων== |
|||
Οι '''[[Σειρές του Νεύτων]]''' ,αποτελούνται από τους όρους της '''εξίσωσης του Νεύτων της προς τα εμπρός διαφοράς''', πήραν το όνομα τους από τον [[Ισαάκ Νεύτων]]; στην ουσία, είναι ο '''τύπος παρεμβολής του Νεύτων''', που δημοσιεύτηκε για πρώτη φορά στο ''[[Philosophiæ Naturalis Principia Mathematica|Principia Mathematica]]'' του το 1687,<ref>Newton, Isaac, (1687). [https://books.google.com/books?id=KaAIAAAAIAAJ&dq=sir%20isaac%20newton%20principia%20mathematica&as_brr=1&pg=PA466#v=onepage&q&f=false ''Principia'', Book III, Lemma V, Case 1]</ref> και συγκεκριμένα το διακριτό ανάλογο της συνεχούς επέκτασης Τέιλορ, |
|||
{{Equation box 1 |
|||
|indent =:: |
|||
|equation = <math>f(x)=\sum_{k=0}^\infty\frac{\Delta^k [f](a)}{k!} ~(x-a)_k |
|||
= \sum_{k=0}^\infty {x-a \choose k}~ \Delta^k [f](a) ~, |
|||
</math> |
|||
|cellpadding= 6 |
|||
|border |
|||
|border colour = #0073CF |
|||
|background colour=#F9FFF7}} |
|||
η οποία ισχύει για κάθε [[πολυωνυμική]] συνάρτηση ''f'' και για τις περισσότερες (αλλά όχι όλες) τις [[αναλυτικές συναρτήσεις]]. Εδώ,η έκφραση |
|||
:<math>{x \choose k} = \frac{(x)_k}{k!}</math> |
|||
είναι ο [[διωνυμικός συντελεστής]], και |
|||
:<math>(x)_k=x(x-1)(x-2)\cdots(x-k+1)</math> |
|||
είναι το "[[φθίνων παραγοντικό]]" ή το "κατώτερο παραγοντικό", ενώ το [[κενό γινόμενο]] (''x'')<sub>0</sub> ορίζεται να είναι είναι 1. Σε αυτή τη συγκεκριμένη περίπτωση, υπάρχει μια υπόθεση των βημάτων μονάδα για τις αλλαγές στις τιμές των ''x'', ''h'' = 1 της παρακάτω γενίκευσης. |
|||
Σημειώστε την συνήθη αντιστοιχία αυτού του αποτελέσματος στο [[Θεώρημα του Τέιλορ]]. Ιστορικά,αυτό, καθώς και η [[ταυτότητα Βαντερμόντ]], |
|||
:<math>(x+y)_n=\sum_{k=0}^n {n \choose k} (x)_{n-k} ~(y)_k ~,</math> |
|||
(απορρέουσες από αυτή,και αντιστοιχίζοντας την [[διωνυμική θεωρία]]), συμπεριλαμβάνονται στις παρατηρήσεις που έληξαν στο σύστημα του λογισμού που αναφέρεται στην ομοιότητα μεταξύ των φαινομενικά άσχετων πολυωνυμικών συναρτήσεων. |
|||
Για να φανεί το πώς κανείς μπορεί να χρησιμοποιήσει τον τύπο του Νεύτων στην πράξη, εξετάστε τους πρώτους όρους διπλασιασμού της [[ακολουθία Φιμπονάτσι|ακολουθίας Φιμπονάτσι]] {{mvar|f}} = 2, 2, 4, ... Μπορεί κανείς να βρει ένα [[πολυώνυμο]] που αναπαράγει αυτές τις τιμές, με πρώτο τον υπολογισμό ενός πίνακα διαφοράς, και στη συνέχεια αντικαθιστώντας τις διαφορές που αντιστοιχούν στο ''x''<sub>0</sub> (υπογραμμισμένο) στο τύπο όπως ακολουθεί, |
|||
:<math> |
|||
\begin{matrix} |
|||
\begin{array}{|c||c|c|c|} |
|||
\hline |
|||
x & f=\Delta^0 & \Delta^1 & \Delta^2 \\ |
|||
\hline |
|||
1&\underline{2}& & \\ |
|||
& &\underline{0}& \\ |
|||
2&2& &\underline{2} \\ |
|||
& &2& \\ |
|||
3&4& & \\ |
|||
\hline |
|||
\end{array} |
|||
& |
|||
\quad \begin{align} |
|||
f(x) & =\Delta^0 \cdot 1 +\Delta^1 \cdot \dfrac{(x-x_0)_1}{1!} + \Delta^2 \cdot \dfrac{(x-x_0)_2}{2!} \quad (x_0=1)\\ |
|||
\\ |
|||
& =2 \cdot 1 + 0 \cdot \dfrac{x-1}{1} + 2 \cdot \dfrac{(x-1)(x-2)}{2} \\ |
|||
\\ |
|||
& =2 + (x-1)(x-2) \\ |
|||
\end{align} |
|||
\end{matrix} |
|||
</math> |
|||
Στην περίπτωση των μη-ομοιόμορφων μέτρων στις τιμές του ''x'', ο Νεύτων υπολογίζει τις [[διαχωρίσιμες διαφορές]], |
|||
:<math>\Delta _{j,0}=y_j,\quad \quad \Delta _{j,k}=\frac{\Delta _{j+1,k-1}-\Delta _{j,k-1}}{x_{j+k}-x_j}\quad \ni \quad \left\{ k>0,\ \ j\le \max \left( j \right)-k \right\},\quad \quad \Delta 0_k=\Delta _{0,k}</math> |
|||
τις σειρές των γινομένων, |
|||
:<math>{P_0}=1,\quad \quad P_{k+1}=P_k\cdot \left( \xi -x_k \right) ~,</math> |
|||
και το πολυώνυμο που προκύπτει είναι το [[βαθμωτό γινόμενο]], <math>f(\xi ) = \Delta 0 \cdot P\left( \xi \right)</math> .<ref>[[Robert D. Richtmyer|Richtmeyer, D.]] and Morton, K.W., (1967). ''Difference Methods for Initial Value Problems'', 2nd ed., Wiley, New York.</ref> |
|||
Σε ανάλυση με τους αριθμούς που έχουν την ιδιότητα να είναι κοντά όταν η διαφορά τους διαιρείται από αριθμό μεγαλύτερης δύναμης του p, διαιρούνται με έναν, το [[Θεώρημα του Mahler]] αναφέρει ότι η υπόθεση ότι η ''f'' είναι μια πολυωνυμική συνάρτηση μπορεί να αποδυναμώσει όλη τη διαδρομή προς την υπόθεση ότι η ''f'' είναι απλώς συνεχής. |
|||
Το [[Θεώρημα του Carlson]] παρέχει σημαντικές και επαρκείς συνθήκες για τη μοναδικότητα των σειρών του Νεύτων,αν υπάρχει.Ωστόσο, οι σειρές του Νεύτων,γενικά,δεν υπάρχουν. |
|||
Οι σειρές του Νεύτων,μαζί με τις [[σειρές του Stirling]] και τις [[σειρές του Selberg]], είναι μια ειδική περίπτωση της γενικής [[διαφοράς των σειρών]], οι οποίες ορίζονται στο πλαίσιο της κατάλληλης κλιμάκωσης των προς προς τα εμπρός διαφορών. |
|||
Σε μία σύντομη και ελαφρώς πιο γενική μορφή και σε ίση απόσταση κόμβων ο τύπος διαβάζει |
|||
:<math>f(x)=\sum_{k=0}{\frac{x-a}h \choose k} \sum_{j=0}^k (-1)^{k-j}{k\choose j}f(a+j h).</math> |
|||
==Λογισμός των πεπερασμένων διαφορών== |
|||
Η προς τα εμπρός διαφορά μπορεί να θεωρηθεί ως [[τελεστής]] διαφοράς,<ref>[[George Boole|Boole, George]], (1872). ''A Treatise On The Calculus of Finite Differences'', 2nd ed., Macmillan and Company. [https://archive.org/details/cu31924031240934 On line]. Also, [Dover edition 1960]</ref><ref>Jordan, Charles, (1939/1965). "Calculus of Finite Differences", Chelsea Publishing. On-line: [https://books.google.com/books?hl=en&lr=&id=3RfZOsDAyQsC&oi=fnd&pg=PA1&ots=AqSuAgOKs3&sig=fzPpAdvnzp7sG6PorqIe5qFjD2Q#v=onepage]</ref> ο οποίος καθορίζει την συνάρτηση {{mvar|''f''}} σε {{math|Δ<sub>''h''</sub>[''f'' ]}}. Αυτός ο τελεστής ισοδυναμεί με |
|||
::<math>\Delta_h = T_h-I, \,</math> |
|||
όπου {{math|''T''<sub>''h''</sub>}} είναι ο [[μεταβαλλόμενος τελεστής]] με δείκτη ''h'', ορίζεται από {{math|''T''<sub>''h''</sub>[''f'' ](''x'') {{=}} ''f''(''x''+''h'')}}, και {{mvar|''I''}} είναι ο [[τελεστής ταυτότητας]]. |
|||
Η πεπερασμένη διαφορά των μεγαλύτερων τάξεων μπορεί να οριστεί με αναδρομικό τρόπο όπως {{math| Δ<sub>''h''</sub><sup>''n''</sup> ≡ Δ<sub>''h''</sub> (Δ<sub>''h''</sub><sup>''n''−1</sup>}}). Ένας άλλος ισοδύναμος ορισμός είναι {{math| Δ<sub>''h''</sub><sup>''n''</sup> {{=}} [''T''<sub>''h''</sub> −''I'']<sup>''n''</sup>}}. |
|||
Ο τελεστής διαφοράς {{math|Δ<sub>''h''</sub>}} είναι ένας [[γραμμικός τελεστής]] και ικανοποιεί έναν ειδικό [[κανόνα του Leibniz]] που αναφέρθηκε παραπάνω, |
|||
{{math|Δ<sub>''h''</sub>(''f''(''x'')''g''(''x'')) {{=}} (Δ<sub>''h''</sub>''f''(''x'')) ''g''(''x''+''h'') + ''f''(''x'') (Δ<sub>''h''</sub>''g''(''x''))}}. Όμοιες δηλώσεις ισχύουν και για τις προς τα πίσω και κεντρικές διαφορές. |
|||
Τυπικά,εφαρμόζωντας τη [[Σειρά Τέιλορ]] σε σχέση με το ''h'', δίνει τον τύπο |
|||
:<math> \Delta_h = hD + \frac{1}{2} h^2D^2 + \frac{1}{3!} h^3D^3 + \cdots = \mathrm{e}^{hD} - I ~, </math> |
|||
όπου το ''D'' δηλώνει την συνεχή παράγωγο του τελεστή, καθορίζοντας την ''f'' στην παράγωγο της ''f'''. Η επέκταση είναι έγκυρη όταν και οι δυο πλευρές δρουν ως [[αναλυτικές συναρτήσεις]], για αρκετά μικρό ''h''. Έτσι, {{math|''T''<sub>''h''</sub>{{=}}e<sup>''hD''</sup>}}, και τυπικά αναστρέφοντας τις εκθετικές αποδόσεις |
|||
:<math> hD = \log(1+\Delta_h) = \Delta_h - \tfrac{1}{2} \Delta_h^2 + \tfrac{1}{3} \Delta_h^3 + \cdots. \, </math> |
|||
Ο τύπος αυτός οδηγεί στην έννοια ότι και οι δυο τύποι δίνουν το ίδιο αποτέλεσμα όταν εφαρμόζεται σε ένα πολυώνυμο. |
|||
Ακόμα και για τις αναλυτικές συναρτήσεις,οι σειρές στα δεξιά δεν εγγυάται ότι συγκλίνουν: μπορεί να είναι [[ασυμπτωτική σειρά]]. Ωστόσο, μπορούν να χρησιμοποιηθούν για να βρεθούν πιο ακριβείς μετασχηματισμοί για την παράγωγο. Για παράδειγμα, διατηρώντας τους πρώτους δύο όρους της σειράς παράγει η προσέγγιση της δεύτερης τάξης την {{math|''f’''(''x'')}} όπως αναφέρθηκε στο τέλος της mentioned at the end of the [[ενότητας των μεγαλύτερης τάξης διαφορών]]. |
|||
Οι ανάλογοι τύποι για τους τελεστές της προς τα πίσω και κεντρικής διαφοράς είναι |
|||
:<math> hD = -\log(1-\nabla_h) \quad\text{and}\quad hD = 2 \, \operatorname{arsinh}(\tfrac12\delta_h). </math> |
|||
Ο λογισμός των πεπερασμένων διαφορών σχετίζεται με το λογισμό (που αναφέρεται στην ομοιότητα μεταξύ των φαινομενικά άσχετων πολυωνυμικών εξισώσεων) της συνδυαστικής. Αυτή η εξαιρετικά συστηματική αντιστοιχία οφείλεται στην ταυτότητα των [[μετατροπέων]] των ποσοτήτων του παραπάνω λογισμού στο αναλογικό συνεχές τους ({{math|''h''→0}} limits), |
|||
{{Equation box 1 |
|||
|indent =:: |
|||
|equation = |
|||
<math> \Bigl[ \frac{\Delta_h}{h} ~,~ x\, T^{-1}_h \Bigr] = [ D ~,~ x ] = I ~.</math> |
|||
|cellpadding= 6 |
|||
|border |
|||
|border colour = #0073CF |
|||
|background colour=#F9FFF7}} |
|||
Ένας μεγάλος αριθμός των διαφορικών σχέσεων του λογισμού με νόρμα αφορά συναρτήσεις {{math|''f''(''x'')}} άρα ''καθορίζει συστηματικά την αναλογική πεπερασμένη διαφορά'' που αφορά την {{math|''f''(''xT''<sub>h</sub><sup>−1</sup>)}}. |
|||
Για παράδειγμα, η αναλογική ενός μονώνυμου ''x''<sup>n</sup> είναι μια γενίκευση του πτωτικού παραγοντικού ([[Pochhammer k-symbol]]), |
|||
:<math>~(x)_n\equiv (xT_h^{-1})^n=x (x-h) (x-2h) \cdots (x-(n-1)h)</math> , |
|||
έτσι ώστε |
|||
::<math>\frac{\Delta_h}{h} ~(x)_n=n ~(x)_{n-1} ~,</math> |
|||
άρα ο παραπάνω τύπος παρεμβολής Τέιλορ (αντιστοιχίζοντας συντελεστές στην επέκταση μιας αυθαίρετης συνάρτησης ''f''(''x'') σε τέτοια σύμβολα), και ούτω καθεξής. |
|||
Για παράδειγμα, το ημίτονο είναι |
|||
:<math>\sin (x\,T_h^{-1}) = x -\frac{(x)_3}{3!} + \frac{(x)_5}{5!} - \frac{(x)_7}{7!} + \cdots .</math> |
|||
Όπως στο συνεχές όριο, η ιδιοσυνάρτηση της {{math| Δ<sub>''h''</sub> /''h''}} επίσης συμβαίνει ''να είναι μια εκθετική'', |
|||
::<math>\frac{\Delta_h}{h}~(1+\lambda h)^{x/h} =\frac{\Delta_h}{h} ~e^{\ln (1+\lambda h) ~x/h}= \lambda ~e^{\ln (1+\lambda h) ~x/h} ~,</math> |
|||
και άρα ''τα αθροίσματα Φουριέρ των συνεχών συναρτήσεων αντιστοιχίζονται εύκολα σε αθροίσματα Φουριέρ που σχετίζονται με την ομοιότητα των φαινομενικά άσχετων πολυωνυμικών εξισώσεων'', που αφορούν τους ίδιους συντελεστές Φουριέρ πολλαπλασιάζοντας αυτές τις εκθετικές βάσεις.<ref>{{cite journal |last =Zachos|first =C.| authorlink =Cosmas Zachos| year =2008| title =Umbral Deformations on Discrete Space-Time | journal =International Journal of Modern Physics A| volume =23 | issue=13| pages =2005–2014 | doi = 10.1142/S0217751X08040548 }}</ref> Αυτή η εκθετική (που έχει σχέση με την ομοιότητα των πολυωνυμικών εξισώσεων)ισοδυναμεί με την εκθετική [[συνάρτηση παραγωγής]] των [[Pochhammer συμβόλων]]. |
|||
Έτσι,για παράδειγμα, η [[συνάρτηση δέλτα του Ντιράκ]] καθορίζει τους αντιπροσώπους, την [[απόλυτη ημιτονοειδή συνάρτηση]], |
|||
:<math>\delta (x) \mapsto \frac{\sin \bigl[ \frac{\pi}{2}(1+x/h) \bigr]}{ \pi (x+h) }~,</math> |
|||
και ούτω καθεξής.<ref>{{Cite journal | last1 = Curtright | first1 = T. L. | last2 = Zachos | first2 = C. K. | doi = 10.3389/fphy.2013.00015 | title = Umbral Vade Mecum | journal = Frontiers in Physics | volume = 1 | year = 2013 | pmid = | pmc = }}</ref> Οι [[πεπερασμένες εξισώσεις]] μπορούν συχνά να λυθούν με τεχνικές πολύ όμοιες με εκείνες για την επίλυση των [[διαφορικών εξισώσεων]]. |
|||
Ο αντίστροφος τελεστής του τελεστή της προς τα εμπρός διαφοράς, ώστε στη συνέχεια να έχουμε το ολοκλήρωμα, είναι το [[αόριστο άθροισμα]] ή ο τελεστής της αντί διαφοράς. |
|||
==Κανόνες λογισμού για τους τελεστές της πεπερασμένης διαφοράς== |
|||
Ανάλογο με τους [[κανόνες για την εύρεση της παραγώγου]], έχουμε: |
|||
* '''Συνεχής κανόνας''': Αν ''c'' είναι μια [[σταθερά]], τότε |
|||
:<math>\Delta c = 0{\,}</math> |
|||
* '''[[Γραμμικότητα]]''': Αν ''a'' και ''b'' είναι [[σταθερές]], |
|||
:<math>\Delta (a f + b g) = a \,\Delta f + b \,\Delta g</math> |
|||
Όλοι οι παραπάνω κανόνες ισχύουν εξίσου καλά σε κάθε τελεστή διαφοράς, συμπεριλαμβανομένου <math>\nabla</math> ως προς <math>\Delta</math>. |
|||
* '''[[Κανόνας γινομένου]]''': |
|||
:<math> \Delta (f g) = f \,\Delta g + g \,\Delta f + \Delta f \,\Delta g </math> |
|||
:<math> \nabla (f g) = f \,\nabla g + g \,\nabla f - \nabla f \,\nabla g </math> |
|||
* '''[[Κανόνας πηλίκου]]''': |
|||
:<math>\nabla \left( \frac{f}{g} \right) = \frac{1}{g} \det \begin{bmatrix} \nabla f & \nabla g \\ f & g \end{bmatrix} |
|||
\left( \det {\begin{bmatrix} g & \nabla g \\ 1 & 1 \end{bmatrix}}\right)^{-1} </math> |
|||
::ή |
|||
:<math>\nabla\left( \frac{f}{g} \right)= \frac {g \,\nabla f - f \,\nabla g}{g \cdot (g - \nabla g)}</math> |
|||
* '''[[Κανόνες αθροίσματος]]''': |
|||
:<math>\sum_{n=a}^{b} \Delta f(n) = f(b+1)-f(a)</math> |
|||
:<math>\sum_{n=a}^{b} \nabla f(n) = f(b)-f(a-1)</math> |
|||
<br />Δες τις παραπομπές <ref>{{cite book|last=Levy|first=H.|author2=Lessman, F.|title=Finite Difference Equations|url=https://archive.org/details/finitedifference0000levy|year=1992|publisher=Dover|isbn=0-486-67260-3}}</ref><ref>Ames, W. F., (1977). ''Numerical Methods for Partial Differential Equations'', Section 1.6. Academic Press, New York. ISBN 0-12-056760-1.</ref><ref>[[Francis B. Hildebrand|Hildebrand, F. B.]], (1968). ''Finite-Difference Equations and Simulations'', Section 2.2, Prentice-Hall, Englewood Cliffs, New Jersey.</ref><ref>{{Cite journal |
|||
| first1 = Philippe | last1 = Flajolet |
|||
| authorlink2 = Robert Sedgewick (computer scientist) | first2 = Robert | last2 = Sedgewick |
|||
| url = http://algo.inria.fr/flajolet/Publications/FlSe95.pdf |
|||
| title = Mellin transforms and asymptotics: Finite differences and Rice's integrals |
|||
| journal=Theoretical Computer Science |
|||
| volume = 144 | issue = 1–2 | year = 1995 | pages = 101–124 |
|||
| doi = 10.1016/0304-3975(94)00281-M |
|||
| postscript = <!-- Bot inserted parameter. Either remove it; or change its value to "." for the cite to end in a ".", as necessary. -->{{inconsistent citations}}}}.</ref> |
|||
==Γενικεύσεις== |
|||
*Μία '''γενικευμένη πεπερασμένη διαφορά''' συνήθως ορίζεται ως |
|||
:<math>\Delta_h^\mu[f](x) = \sum_{k=0}^N \mu_k f(x+kh),</math> |
|||
όπου <math>\mu = (\mu_0,\ldots,\mu_N)</math> είναι το διάνυσμα των συντελεστών του. Μία '''άπειρη διαφορά''' είναι μια περαιτέρω γενίκευση,όπου το πεπερασμένο άθροισμα παραπάνω αντικαθίσταται από μία [[άπειρη σειρά]]. Ένας άλλος τρόπος γενίκευσης είναι η κατασκευή συντελεστών <math>\mu_k</math> που εξαρτώνται από το σημείο <math>x</math> : <math>\mu_k=\mu_k(x)</math>, θεωρώντας έτσι '''σταθμισμένη πεπερασμένη διαφορά'''. Επίσης,θα μπορούσε κανείς να φτιάξει έναν δείκτη <math>h</math> που να εξαρτάται από το σημείο <math>x</math> : <math>h=h(x)</math>. Τέτοιες γενικεύσεις είναι χρήσιμες για την κατασκευή διαφορετικού [[μέτρου συνέχειας]]. |
|||
*Η γενικευμένη διαφορά μπορεί να θεωρηθεί ως πολυώνυμο δακτυλίων <math>R[T_h]</math> . Αυτό οδηγεί σε αλγεβρική διαφορά. |
|||
*Ο φορέας διαφοράς γενικεύει την [[Möbius μετάθεση]] πάνω από ένα [[μερικώς διατεταγμένο σύνολο]]. |
|||
*Ως τελεστής ελίκωσης: Μέσω του φορμαλισμού των [[αλγεβρικών συχνοτήτων]], τελεστές διαφοράς και άλλη Möbius αναστροφή μπορεί να εκπροσωπείται από [[συνεξέλιξη]] με μία συνάρτηση on the poset, γνωστή ως [[Möbius συνάρτηση]] μ; για τον τελεστή διαφοράς, μ είναι η ακολουθία (1, −1, 0, 0, 0, ...). |
|||
==Πεπερασμένη διαφορά πολλών μεταβλητών== |
|||
Οι πεπερασμένες διαφορές μπορούν να θεωρηθούν περισσότερες από μία μεταβλητές. Είναι ανάλογες με τις [[μερικές παραγώγους]] σε αρκετές μεταβλητές. |
|||
Μερικές προσεγγίσεις των μερικών παραγώγων είναι (χρησιμοποιώντας την κεντρική μέθοδο του δείκτη): |
|||
:<math> f_{x}(x,y) \approx \frac{f(x+h ,y) - f(x-h,y)}{2h} \ </math> |
|||
:<math> f_{y}(x,y) \approx \frac{f(x,y+k ) - f(x,y-k)}{2k} \ </math> |
|||
:<math> f_{xx}(x,y) \approx \frac{f(x+h ,y) - 2 f(x,y) + f(x-h,y)}{h^2} \ </math> |
|||
:<math> f_{yy}(x,y) \approx \frac{f(x,y+k) - 2 f(x,y) + f(x,y-k)}{k^2} \ </math> |
|||
:<math> f_{xy}(x,y) \approx \frac{f(x+h,y+k) - f(x+h,y-k) - f(x-h,y+k) + f(x-h,y-k)}{4hk} ~. </math> |
|||
Εναλλακτικά,σε εφαρμογές στις οποίες ο υπολογισμός της {{mvar|f}} είναι πιο δύσκολο βήμα, και πρέπει να υπολογιστούν και οι δύο,πρώτοι και δεύτεροι παράγωγοι, ένας πιο αποτελεσματικός τύπος για την τελευταία περίπτωση είναι |
|||
:<math> f_{xy}(x,y) \approx \frac{f(x+h, y+k) - f(x+h, y) - f(x, y+k) + 2 f(x,y) - f(x-h, y) - f(x, y-k) + f(x-h, y-k)}{2hk} ~,</math> |
|||
δεδομένου ότι οι μόνες τιμές για υπολογισμό που δεν είναι απαραίτητες για τις προηγούμενες τέσσερις εξισώσεις είναι {{math|''f''(''x''+''h'', ''y''+''k'')}} και {{math|''f''(''x''−''h'', ''y''−''k'')}}. |
|||
== Παραπομπές == |
|||
<references /> |
|||
== Πηγές == |
|||
* Richardson, C. H. (1954): ''An Introduction to the Calculus of Finite Differences'' (Van Nostrand (1954)'' [http://babel.hathitrust.org/cgi/pt?id=mdp.39015000982945;view=1up;seq=5 online copy] |
|||
* Mickens, R. E. (1991): ''Difference Equations: Theory and Applications'' (Chapman and Hall/CRC) ISBN 978-0442001360 |
|||
== Εξωτερικοί σύνδεσμοι == |
|||
* {{springer|title=Finite-difference calculus|id=p/f040230}} |
|||
* [http://reference.wolfram.com/mathematica/tutorial/NDSolvePDE.html#c:4 Table of useful finite difference formula generated using Mathematica ] |
|||
* D. Gleich (2005), [http://www.stanford.edu/~dgleich/publications/finite-calculus.pdf ''Finite Calculus: A Tutorial for Solving Nasty Sums''] {{Webarchive|url=https://web.archive.org/web/20090419132601/http://www.stanford.edu/~dgleich/publications/finite-calculus.pdf |date=2009-04-19 }} |
|||
* [http://mathformeremortals.wordpress.com/2013/01/12/a-numerical-second-derivative-from-three-points/ Discrete Second Derivative from Unevenly Spaced Points] |
|||
{{Portal bar|Μαθηματικά}} |
|||
[[Κατηγορία:Πεπερασμένες διαφορές]] |
Τρέχουσα έκδοση από την 14:07, 11 Φεβρουαρίου 2024
Η πεπερασμένη διαφορά είναι μία μαθηματική έκφραση της μορφής f(x + b) − f(x + a). Αν η πεπερασμένη διαφορά διαιρείται από το b − a, τότε αυτή λαμβάνει μόνο ένα πεπερασμένο πηλίκο.Ο υπολογισμός των παραγώγων από τις πεπερασμένες διαφορές παίζει βασικό ρόλο στις μεθόδους πεπερασμένης διαφοράς για την αριθμητική λύση των διαφορικών εξισώσεων, ειδικότερα στα προβλήματα συνοριακών τιμών.
Ορισμένες επαναληπτικές σχέσεις μπορούν να γραφούν ως εξισώσεις διαφοράς αντικαθιστώντας την προσέγγιση της σημειογραφείας με τις πεπερασμένες διαφορές.
Σήμερα, ο όρος "πεπερασμένη διαφορά" συχνά λαμβάνεται ως συνώνυμο με τους υπολογισμούς της πεπερασμένης διαφοράς των παραγώγων, ειδικότερα στο πλαίσιο των αριθμητικών μεθόδων.[1][2][3] Οι υπολογισμοί των πεπερασμένων διαφορών είναι πηλίκα της πεπερασμένης διαφοράς στην ορολογία που χρησιμοποιήθηκε ανωτέρω.
Οι πεπερασμένες διαφορές,επίσης,έχουν γίνει το θέμα μελέτης ως αφηρημένα αυτοτελή μαθηματικά αντικείμενα, π.χ. σε έργα των Τζορτζ Μπουλ (1860), L. M. Milne-Thomson (1933), και Károly Jordan (1939), εντοπίζοντας τις προελέυσεις των αλγορίθμων του Jost Bürgi (ca. 1592) και άλλων συμπεριλαμβανομένου του Ισαάκ Νεύτων. Σε αυτό το σημείο, ο λογισμός (που ασχολείται με τα όρια) των πεπερασμένων τιμών είναι εναλλακτικός του λογισμού των απειροελάχιστων.[4]
Προς τα εμπρός, προς τα πίσω, και κεντρική διαφορά
[Επεξεργασία | επεξεργασία κώδικα]Τρεις μορφές μελετώνται ευρέως: οι προς τα εμπρός,οι προς τα πίσω και οι κεντρικές διαφορές.[1][2][3]
Η προς τα εμπρός διαφορά (προς τα πάνω) είναι μια έκφραση της μορφής
Σύμφωνα με την εφαρμογή, η τάξη h μπορεί να είναι μεταβλητή ή σταθερή. Όταν παραλείπεται, η h θεωρείται ότι είναι 1: .
Η προς τα πίσω διαφορά (προς τα κάτω) χρησιμοποιεί τις τιμές της συνάρτησης στο x και x − h, αντί των τιμών στο x + h and x:
Τελικά, η κεντρική διαφορά (κεντρώα) δίνεται από
Σχέση με τις παραγώγους
[Επεξεργασία | επεξεργασία κώδικα]Η παράγωγος της συνάρτησης f στο σημείο x ορίζεται από το όριο
Αν η h λαμβάνει μια σταθερή (μη μηδενική) τιμή αντί να πλησιάζει το μηδέν, τότε η δεξιά πλευρά της παραπάνω εξίσωσης γράφεται
Συνεπώς, η προς τα πάνω διαφορά που διαιρείται από το h πλησιάζει την παράγωγο όταν το h είναι μικρό. Το σφάλμα σε αυτή την προσέγγιση μπορεί να βρεθεί από τον Τύπο του Taylor. Λαμβάνοντας υπόψη, ότι η f είναι διαφορίσιμη , έχουμε
Ο ίδιος τύπος αναμένεται για την προς τα πίσω διαφορά:
Ωστόσο, η κεντρική (ή αλλιώς κεντρώα) διαφορά αποδίδει μια πιο ακριβή προσέγγιση. Αν η f είναι δύο φορές διαφορίσιμη,
Το κύριο πρόβλημα με τη μέθοδο της κεντρικής διαφοράς, ωστόσο , είναι ότι οι συναρτήσεις που ταλαντεύονται αποδίδουν μηδενική παράγωγο. Αν f(nh)=1 για n περιττό, και f(nh)=2 για n άρτιο, αλλιώς f ' (nh)=0 αν υπολογίζεται με το συνδυασμό της κεντρικής διαφοράς. Αυτό είναι ιδιαίτερα ενοχλητικό αν ο κλάδος της f είναι διακεκριμένος.
Συγγραφείς, για τους οποίους οι πεπερασμένες διαφορές σημαίνουν προσεγγίσεις πεπερασμένων διαφορών, ορίζουν τις προς τα πάνω/προς τα κάτω/κεντρικές διαφορές ως πηλίκα που δίνονται σε αυτή την ενότητα. (αντί να χρησιμοποιούν τους ορισμούς που δόθηκαν στην προηγούμενη ενότητα.)[1][2][3]
Μεγαλύτερης τάξης διαφορές
[Επεξεργασία | επεξεργασία κώδικα]Με ανάλογο τρόπο,μπορεί κανείς να ανάγει τις προσεγγίσεις πεπερασμένης διαφοράς σε μεγαλύτερης τάξης παραγώγους και διαφορικούς φορείς. Για παράδειγμα, χρησιμοποιώντας τον παραπάνω τύπο κεντρικής διαφοράς για f ' (x+h/2) και f ' (x−h/2) και εφαρμόζοντας τον τύπο της κεντρικής διαφοράς για την παράγωγο της f ' στο x, παίρνουμε την προσέγγιση της κεντρικής διαφοράς της δεύτερης παραγώγου της f:
Κεντρική 2ης τάξης
Με παρόμοιο τρόπο μπορούμε να εφαρμόσουμε άλλους τύπους διαφοράς με έναν αναδρομικό τρόπο.
Προς τα εμπρός 2ης τάξης
Γενικότερα , οι n-ης τάξης προς τα εμπρός, προς τα πίσω, και κεντρικές διαφορές δίνονται,αντίστοιχα,
Προς τα εμπρός
ή για h=1,
Προς τα πίσω
Κεντρική
Αυτές οι εξισώσεις χρησιμοποιούν διωνυμικούς συντελεστές ύστερα από το σημείο άθροισης απεικονίζονται ως . Κάθε σειρά από το Τρίγωνο του Pascal παρέχει τον συντελεστή για κάθε τιμή του i.
Παρατήρησε ότι η κεντρική διαφορά θα, για περιττό n, έχει h πολλαπλασιαζόμενο με μη ακέραιο. Αυτό είναι συχνά πρόβλημα επειδή ανέρχεται σε αλλαγή του διαστήματος διακριτοποίησης. Το πρόβλημα μπορεί να αντιμετωπιστεί παίρνοντας το μέσο όρο του και .
Οι προς τα εμπρός διαφορές που εφαρμόζονται σε μία ακολουθία ονομάζονται μερικές φορές διωνυμικός μετασχηματισμός της ακολουθίας, και έχουν μια σειρά από ενδιαφέρουσες συνδυαστικές ιδιότητες. Οι προς τα εμπρός διαφορές μπορούν να εκτιμηθούν χρησιμοποιώντας το Nörlund–Rice ολοκλήρωμα. Η αναπαράσταση του ολοκληρώματος για αυτούς τους τύπους των σειρών είναι ενδιαφέρουσα, επειδή το ολοκλήρωμα μπορεί συχνά να εκτιμηθεί χρησιμοποιώντας ασυμπτωτική επέκταση ή τεχνικές απότομης καθόδου : αντίθετα, η προς τα εμπρός διαφορά των σειρών μπορεί να είναι εξαιρετικά δύσκολο να εκτιμηθεί αριθμητικά, επειδή οι διωνυμικοί συντελεστές αυξάνονται με ταχείς ρυθμούς για μεγάλα n.
Η σχέση των διαφορών μεγαλύτερης τάξης με τις αντίστοιχες παραγώγους είναι απλή,
Οι μεγαλύτερης τάξης διαφορές μπορούν να χρησιμοποιηθούν για να φτιάξουμε καλύτερες προσεγγίσεις. Όπως αναφέρθηκε παραπάνω, η πρώτης τάξης διαφορά προσεγγίζει την παράγωγο πρώτης τάξης σε έναν όρο της τάξης h. Ωστόσο, ο συνδυασμός
προσεγγίζει την f'(x) σε έναν όρο της τάξης h2. Αυτό αποδεικνύεται αναπτύσσοντας την παραπάνω έκφραση στη Σειρά Τέιλορ, ή χρησιμοποιώντας το λογισμό των πεπερασμένων διαφορών,που εξηγείται παρακάτω.
Αν χρειάζεται,η πεπερασμένη διαφορά μπορεί να κεντράρεται γύρω από οποιοδήποτε σημείο συνδυάζοντας τις προς τα πάνω,προς τα κάτω και κεντρικές διαφορές.
Πυρήνες αυθαίρετου μεγέθους
[Επεξεργασία | επεξεργασία κώδικα]Χρησιμοποιώντας λίγο την γραμμική άλγεβρα, μπορεί κανείς αρκετά εύκολα να φτιάξει προσεγγίσεις, της οποίας παίρνω δείγμα ένα αυθαίρετο αριθμό σημείων από τα αριστερά και έναν (ενδεχομένως διαφορετικό) αριθμό σημείων από τα δεξιά του κεντρικού σημείου,για κάθε τάξη της παραγώγου. Αυτό συνεπάγεται την επίλυση ενός γραμμικού συστήματος έτσι ώστε η επέκταση Τέιλορ του αθροίσματος των εν λόγω σημείων ,γύρω από το κεντρικό σημείο,να προσεγγίζει καλά την επέκταση Τέιλορ της επιθυμητής παραγώγου.
Αυτό είναι χρήσιμο για την διαφοροποίηση μιας συνάρτησης σε μια περιοχή,όπου,καθώς αυτή θα πλησιάζει την άκρη της περιοχής, πρέπει να λαμβάνεται δείγμα όλο και λιγότερων σημείων στη μία πλευρά.
Οι λεπτομέρειες συνοψίζονται σε αυτές τις σημειώσεις.
Ιδιότητες
[Επεξεργασία | επεξεργασία κώδικα]- Για όλα τα θετικά k and n
Μέθοδοι πεπερασμένης διαφοράς
[Επεξεργασία | επεξεργασία κώδικα]Μία σημαντική εφαρμογή των πεπερασμένων διαφορών είναι στην αριθμητική ανάλυση, ειδικότερα στις αριθμητικές διαφορικές εξισώσεις, οι οποίες στοχεύουν στην αριθμητική επίλυση των συνήθων και μερικών διαφορικών εξισώσεων αντίστοιχα. Η ιδέα είναι να αντικατασταθούν οι παράγωγοι που εμφανίζονται στην διαφορική εξίσωση από τις πεπερασμένες διαφορές που τις προσεγγίζουν. Οι μέθοδοι που προκύπτουν ονομάζονται μέθοδοι πεπερασμένης διαφοράς.
Κοινές εφαρμογές της μεθόδου πεπερασμένης τιμής συναντώνται στην υπολογιστική επιστήμη και στην εφαρμοσμένη μηχανική, όπως θερμική μηχανική, μηχανική των ρευστών, κτλ.
Σειρές του Νεύτων
[Επεξεργασία | επεξεργασία κώδικα]Οι Σειρές του Νεύτων ,αποτελούνται από τους όρους της εξίσωσης του Νεύτων της προς τα εμπρός διαφοράς, πήραν το όνομα τους από τον Ισαάκ Νεύτων; στην ουσία, είναι ο τύπος παρεμβολής του Νεύτων, που δημοσιεύτηκε για πρώτη φορά στο Principia Mathematica του το 1687,[5] και συγκεκριμένα το διακριτό ανάλογο της συνεχούς επέκτασης Τέιλορ,
Πρότυπο:Equation box 1 η οποία ισχύει για κάθε πολυωνυμική συνάρτηση f και για τις περισσότερες (αλλά όχι όλες) τις αναλυτικές συναρτήσεις. Εδώ,η έκφραση
είναι ο διωνυμικός συντελεστής, και
είναι το "φθίνων παραγοντικό" ή το "κατώτερο παραγοντικό", ενώ το κενό γινόμενο (x)0 ορίζεται να είναι είναι 1. Σε αυτή τη συγκεκριμένη περίπτωση, υπάρχει μια υπόθεση των βημάτων μονάδα για τις αλλαγές στις τιμές των x, h = 1 της παρακάτω γενίκευσης.
Σημειώστε την συνήθη αντιστοιχία αυτού του αποτελέσματος στο Θεώρημα του Τέιλορ. Ιστορικά,αυτό, καθώς και η ταυτότητα Βαντερμόντ,
(απορρέουσες από αυτή,και αντιστοιχίζοντας την διωνυμική θεωρία), συμπεριλαμβάνονται στις παρατηρήσεις που έληξαν στο σύστημα του λογισμού που αναφέρεται στην ομοιότητα μεταξύ των φαινομενικά άσχετων πολυωνυμικών συναρτήσεων.
Για να φανεί το πώς κανείς μπορεί να χρησιμοποιήσει τον τύπο του Νεύτων στην πράξη, εξετάστε τους πρώτους όρους διπλασιασμού της ακολουθίας Φιμπονάτσι f = 2, 2, 4, ... Μπορεί κανείς να βρει ένα πολυώνυμο που αναπαράγει αυτές τις τιμές, με πρώτο τον υπολογισμό ενός πίνακα διαφοράς, και στη συνέχεια αντικαθιστώντας τις διαφορές που αντιστοιχούν στο x0 (υπογραμμισμένο) στο τύπο όπως ακολουθεί,
Στην περίπτωση των μη-ομοιόμορφων μέτρων στις τιμές του x, ο Νεύτων υπολογίζει τις διαχωρίσιμες διαφορές,
τις σειρές των γινομένων,
και το πολυώνυμο που προκύπτει είναι το βαθμωτό γινόμενο, .[6]
Σε ανάλυση με τους αριθμούς που έχουν την ιδιότητα να είναι κοντά όταν η διαφορά τους διαιρείται από αριθμό μεγαλύτερης δύναμης του p, διαιρούνται με έναν, το Θεώρημα του Mahler αναφέρει ότι η υπόθεση ότι η f είναι μια πολυωνυμική συνάρτηση μπορεί να αποδυναμώσει όλη τη διαδρομή προς την υπόθεση ότι η f είναι απλώς συνεχής.
Το Θεώρημα του Carlson παρέχει σημαντικές και επαρκείς συνθήκες για τη μοναδικότητα των σειρών του Νεύτων,αν υπάρχει.Ωστόσο, οι σειρές του Νεύτων,γενικά,δεν υπάρχουν.
Οι σειρές του Νεύτων,μαζί με τις σειρές του Stirling και τις σειρές του Selberg, είναι μια ειδική περίπτωση της γενικής διαφοράς των σειρών, οι οποίες ορίζονται στο πλαίσιο της κατάλληλης κλιμάκωσης των προς προς τα εμπρός διαφορών.
Σε μία σύντομη και ελαφρώς πιο γενική μορφή και σε ίση απόσταση κόμβων ο τύπος διαβάζει
Λογισμός των πεπερασμένων διαφορών
[Επεξεργασία | επεξεργασία κώδικα]Η προς τα εμπρός διαφορά μπορεί να θεωρηθεί ως τελεστής διαφοράς,[7][8] ο οποίος καθορίζει την συνάρτηση f σε Δh[f ]. Αυτός ο τελεστής ισοδυναμεί με
όπου Th είναι ο μεταβαλλόμενος τελεστής με δείκτη h, ορίζεται από Th[f ](x) = f(x+h), και I είναι ο τελεστής ταυτότητας.
Η πεπερασμένη διαφορά των μεγαλύτερων τάξεων μπορεί να οριστεί με αναδρομικό τρόπο όπως Δhn ≡ Δh (Δhn−1). Ένας άλλος ισοδύναμος ορισμός είναι Δhn = [Th −I]n.
Ο τελεστής διαφοράς Δh είναι ένας γραμμικός τελεστής και ικανοποιεί έναν ειδικό κανόνα του Leibniz που αναφέρθηκε παραπάνω, Δh(f(x)g(x)) = (Δhf(x)) g(x+h) + f(x) (Δhg(x)). Όμοιες δηλώσεις ισχύουν και για τις προς τα πίσω και κεντρικές διαφορές.
Τυπικά,εφαρμόζωντας τη Σειρά Τέιλορ σε σχέση με το h, δίνει τον τύπο
όπου το D δηλώνει την συνεχή παράγωγο του τελεστή, καθορίζοντας την f στην παράγωγο της f'. Η επέκταση είναι έγκυρη όταν και οι δυο πλευρές δρουν ως αναλυτικές συναρτήσεις, για αρκετά μικρό h. Έτσι, Th=ehD, και τυπικά αναστρέφοντας τις εκθετικές αποδόσεις
Ο τύπος αυτός οδηγεί στην έννοια ότι και οι δυο τύποι δίνουν το ίδιο αποτέλεσμα όταν εφαρμόζεται σε ένα πολυώνυμο.
Ακόμα και για τις αναλυτικές συναρτήσεις,οι σειρές στα δεξιά δεν εγγυάται ότι συγκλίνουν: μπορεί να είναι ασυμπτωτική σειρά. Ωστόσο, μπορούν να χρησιμοποιηθούν για να βρεθούν πιο ακριβείς μετασχηματισμοί για την παράγωγο. Για παράδειγμα, διατηρώντας τους πρώτους δύο όρους της σειράς παράγει η προσέγγιση της δεύτερης τάξης την f’(x) όπως αναφέρθηκε στο τέλος της mentioned at the end of the ενότητας των μεγαλύτερης τάξης διαφορών. Οι ανάλογοι τύποι για τους τελεστές της προς τα πίσω και κεντρικής διαφοράς είναι
Ο λογισμός των πεπερασμένων διαφορών σχετίζεται με το λογισμό (που αναφέρεται στην ομοιότητα μεταξύ των φαινομενικά άσχετων πολυωνυμικών εξισώσεων) της συνδυαστικής. Αυτή η εξαιρετικά συστηματική αντιστοιχία οφείλεται στην ταυτότητα των μετατροπέων των ποσοτήτων του παραπάνω λογισμού στο αναλογικό συνεχές τους (h→0 limits),
Ένας μεγάλος αριθμός των διαφορικών σχέσεων του λογισμού με νόρμα αφορά συναρτήσεις f(x) άρα καθορίζει συστηματικά την αναλογική πεπερασμένη διαφορά που αφορά την f(xTh−1).
Για παράδειγμα, η αναλογική ενός μονώνυμου xn είναι μια γενίκευση του πτωτικού παραγοντικού (Pochhammer k-symbol),
- ,
έτσι ώστε
άρα ο παραπάνω τύπος παρεμβολής Τέιλορ (αντιστοιχίζοντας συντελεστές στην επέκταση μιας αυθαίρετης συνάρτησης f(x) σε τέτοια σύμβολα), και ούτω καθεξής.
Για παράδειγμα, το ημίτονο είναι
Όπως στο συνεχές όριο, η ιδιοσυνάρτηση της Δh /h επίσης συμβαίνει να είναι μια εκθετική,
και άρα τα αθροίσματα Φουριέρ των συνεχών συναρτήσεων αντιστοιχίζονται εύκολα σε αθροίσματα Φουριέρ που σχετίζονται με την ομοιότητα των φαινομενικά άσχετων πολυωνυμικών εξισώσεων, που αφορούν τους ίδιους συντελεστές Φουριέρ πολλαπλασιάζοντας αυτές τις εκθετικές βάσεις.[9] Αυτή η εκθετική (που έχει σχέση με την ομοιότητα των πολυωνυμικών εξισώσεων)ισοδυναμεί με την εκθετική συνάρτηση παραγωγής των Pochhammer συμβόλων.
Έτσι,για παράδειγμα, η συνάρτηση δέλτα του Ντιράκ καθορίζει τους αντιπροσώπους, την απόλυτη ημιτονοειδή συνάρτηση,
και ούτω καθεξής.[10] Οι πεπερασμένες εξισώσεις μπορούν συχνά να λυθούν με τεχνικές πολύ όμοιες με εκείνες για την επίλυση των διαφορικών εξισώσεων.
Ο αντίστροφος τελεστής του τελεστή της προς τα εμπρός διαφοράς, ώστε στη συνέχεια να έχουμε το ολοκλήρωμα, είναι το αόριστο άθροισμα ή ο τελεστής της αντί διαφοράς.
Κανόνες λογισμού για τους τελεστές της πεπερασμένης διαφοράς
[Επεξεργασία | επεξεργασία κώδικα]Ανάλογο με τους κανόνες για την εύρεση της παραγώγου, έχουμε:
- Συνεχής κανόνας: Αν c είναι μια σταθερά, τότε
- Γραμμικότητα: Αν a και b είναι σταθερές,
Όλοι οι παραπάνω κανόνες ισχύουν εξίσου καλά σε κάθε τελεστή διαφοράς, συμπεριλαμβανομένου ως προς .
-
- ή
Δες τις παραπομπές [11][12][13][14]
Γενικεύσεις
[Επεξεργασία | επεξεργασία κώδικα]- Μία γενικευμένη πεπερασμένη διαφορά συνήθως ορίζεται ως
όπου είναι το διάνυσμα των συντελεστών του. Μία άπειρη διαφορά είναι μια περαιτέρω γενίκευση,όπου το πεπερασμένο άθροισμα παραπάνω αντικαθίσταται από μία άπειρη σειρά. Ένας άλλος τρόπος γενίκευσης είναι η κατασκευή συντελεστών που εξαρτώνται από το σημείο : , θεωρώντας έτσι σταθμισμένη πεπερασμένη διαφορά. Επίσης,θα μπορούσε κανείς να φτιάξει έναν δείκτη που να εξαρτάται από το σημείο : . Τέτοιες γενικεύσεις είναι χρήσιμες για την κατασκευή διαφορετικού μέτρου συνέχειας.
- Η γενικευμένη διαφορά μπορεί να θεωρηθεί ως πολυώνυμο δακτυλίων . Αυτό οδηγεί σε αλγεβρική διαφορά.
- Ο φορέας διαφοράς γενικεύει την Möbius μετάθεση πάνω από ένα μερικώς διατεταγμένο σύνολο.
- Ως τελεστής ελίκωσης: Μέσω του φορμαλισμού των αλγεβρικών συχνοτήτων, τελεστές διαφοράς και άλλη Möbius αναστροφή μπορεί να εκπροσωπείται από συνεξέλιξη με μία συνάρτηση on the poset, γνωστή ως Möbius συνάρτηση μ; για τον τελεστή διαφοράς, μ είναι η ακολουθία (1, −1, 0, 0, 0, ...).
Πεπερασμένη διαφορά πολλών μεταβλητών
[Επεξεργασία | επεξεργασία κώδικα]Οι πεπερασμένες διαφορές μπορούν να θεωρηθούν περισσότερες από μία μεταβλητές. Είναι ανάλογες με τις μερικές παραγώγους σε αρκετές μεταβλητές.
Μερικές προσεγγίσεις των μερικών παραγώγων είναι (χρησιμοποιώντας την κεντρική μέθοδο του δείκτη):
Εναλλακτικά,σε εφαρμογές στις οποίες ο υπολογισμός της f είναι πιο δύσκολο βήμα, και πρέπει να υπολογιστούν και οι δύο,πρώτοι και δεύτεροι παράγωγοι, ένας πιο αποτελεσματικός τύπος για την τελευταία περίπτωση είναι
δεδομένου ότι οι μόνες τιμές για υπολογισμό που δεν είναι απαραίτητες για τις προηγούμενες τέσσερις εξισώσεις είναι f(x+h, y+k) και f(x−h, y−k).
Παραπομπές
[Επεξεργασία | επεξεργασία κώδικα]- ↑ 1,0 1,1 1,2 Paul Wilmott· Sam Howison· Jeff Dewynne (1995). The Mathematics of Financial Derivatives: A Student Introduction. Cambridge University Press. σελ. 137. ISBN 978-0-521-49789-3.
- ↑ 2,0 2,1 2,2 Peter Olver (2013). Introduction to Partial Differential Equations. Springer Science & Business Media. σελ. 182. ISBN 978-3-319-02099-0.
- ↑ 3,0 3,1 3,2 M Hanif Chaudhry (2007). Open-Channel Flow. Springer. σελ. 369. ISBN 978-0-387-68648-6.
- ↑ Jordán, op. cit., p. 1 and Milne-Thomson, p. xxi. Milne-Thomson, Louis Melville (2000): The Calculus of Finite Differences (Chelsea Pub Co, 2000) ISBN 978-0821821077
- ↑ Newton, Isaac, (1687). Principia, Book III, Lemma V, Case 1
- ↑ Richtmeyer, D. and Morton, K.W., (1967). Difference Methods for Initial Value Problems, 2nd ed., Wiley, New York.
- ↑ Boole, George, (1872). A Treatise On The Calculus of Finite Differences, 2nd ed., Macmillan and Company. On line. Also, [Dover edition 1960]
- ↑ Jordan, Charles, (1939/1965). "Calculus of Finite Differences", Chelsea Publishing. On-line: [1]
- ↑ Zachos, C. (2008). «Umbral Deformations on Discrete Space-Time». International Journal of Modern Physics A 23 (13): 2005–2014. doi: .
- ↑ Curtright, T. L.; Zachos, C. K. (2013). «Umbral Vade Mecum». Frontiers in Physics 1. doi: .
- ↑ Levy, H.· Lessman, F. (1992). Finite Difference Equations. Dover. ISBN 0-486-67260-3.
- ↑ Ames, W. F., (1977). Numerical Methods for Partial Differential Equations, Section 1.6. Academic Press, New York. ISBN 0-12-056760-1.
- ↑ Hildebrand, F. B., (1968). Finite-Difference Equations and Simulations, Section 2.2, Prentice-Hall, Englewood Cliffs, New Jersey.
- ↑ Flajolet, Philippe; Sedgewick, Robert (1995). «Mellin transforms and asymptotics: Finite differences and Rice's integrals». Theoretical Computer Science 144 (1–2): 101–124. doi:. http://algo.inria.fr/flajolet/Publications/FlSe95.pdf.
Πηγές
[Επεξεργασία | επεξεργασία κώδικα]- Richardson, C. H. (1954): An Introduction to the Calculus of Finite Differences (Van Nostrand (1954) online copy
- Mickens, R. E. (1991): Difference Equations: Theory and Applications (Chapman and Hall/CRC) ISBN 978-0442001360
Εξωτερικοί σύνδεσμοι
[Επεξεργασία | επεξεργασία κώδικα]- Hazewinkel, Michiel, επιμ.. (2001), «Finite-difference calculus», Encyclopedia of Mathematics, Springer, ISBN 978-1-55608-010-4, http://www.encyclopediaofmath.org/index.php?title=p/f040230
- Table of useful finite difference formula generated using Mathematica
- D. Gleich (2005), Finite Calculus: A Tutorial for Solving Nasty Sums Αρχειοθετήθηκε 2009-04-19 στο Wayback Machine.
- Discrete Second Derivative from Unevenly Spaced Points