Endelig forskel - Finite difference

En endelig forskel er et matematisk udtryk for formen f  ( x + b ) - f  ( x + a ) . Hvis en begrænset forskel divideres med b - a , får man en differenskvotient . Tilnærmelse af derivater med begrænsede forskelle spiller en central rolle i endelige differensmetoder til numerisk løsning af differentialligninger , især grænseværdiproblemer .

Certain recurrence relations can be written as difference equations by replacing iteration notation with finite differences.

I dag betragtes udtrykket "begrænset forskel" ofte som synonymt med tilnærmelser til afgrænsede forskelle mellem derivater , især i forbindelse med numeriske metoder . Endelige forskelstilnærmelser er begrænsede forskelskvotienter i den terminologi, der er anvendt ovenfor.

Endelige forskelle blev introduceret af Brook Taylor i 1715 og er også blevet undersøgt som abstrakte selvstændige matematiske objekter i værker af George Boole (1860), LM Milne-Thomson (1933) og Károly Jordan (1939). Endelige forskelle sporer deres oprindelse tilbage til en af Jost Bürgis algoritmer ( ca.  1592 ) og værker af andre, herunder Isaac Newton . Den formelle beregning af begrænsede forskelle kan ses som et alternativ til beregningen af uendelige tal .

Grundlæggende typer

Image
De tre typer af de begrænsede forskelle. Den centrale forskel på x giver den bedste tilnærmelse af funktionens derivat ved x.

Tre grundtyper betragtes almindeligvis: fremad , bagud og centrale begrænsede forskelle.

En forskel fremad , betegnet med en funktion f er en funktion defineret som

Afhængigt af applikationen kan afstanden h være variabel eller konstant. Når den udelades, anses h for at være 1; det er,

En bagudrettet forskel bruger funktionsværdierne ved x og x - h , i stedet for værdierne ved x + h og  x :

Endelig er den centrale forskel givet ved

Forholdet til derivater

Endelig forskel bruges ofte som en tilnærmelse til derivatet, typisk i numerisk differentiering .

Den derivat af en funktion f ved et punkt x er defineret af grænsen .

Hvis h har en fast (ikke-nul) værdi i stedet for at nærme sig nul, vil højre side af ovenstående ligning blive skrevet

Derfor er terminsforskellen divideret med h tilnærmede derivatet, når h er lille. Fejlen i denne tilnærmelse kan udledes af Taylors sætning . Forudsat at f er to gange differentierbar, har vi

Den samme formel gælder for den tilbagestående forskel:

Den centrale (også kaldet centreret) forskel giver imidlertid en mere præcis tilnærmelse. Hvis f er tre gange differentierbar,

Hovedproblemet med den centrale differensmetode er imidlertid, at oscillerende funktioner kan give nul -derivater. Hvis f  ( nh ) = 1 for n ulige, og f  ( nh ) = 2 for n lige, så er f  ′ ( nh ) = 0, hvis det beregnes med det centrale forskelsskema . Dette er især besværligt, hvis domænet f er diskret. Se også Symmetrisk derivat

Forfattere, for hvem begrænsede forskelle betyder endelige forskelle, tilnærmelser definerer forskellene fremad/bagud/centralt som kvotienterne i dette afsnit (i stedet for at anvende definitionerne i det foregående afsnit).

Højere orden forskelle

På en analog måde kan man opnå endelige forskelstilnærmelser til derivater af højere orden og differentialoperatorer. For eksempel ved at bruge ovenstående centrale forskelsformel for f  ′ ( x + h/2) og f  ′ ( x -h/2) og ved at anvende en central differensformel for derivatet af f  ′ ved x , opnår vi den centrale forskels tilnærmelse til det andet derivat af f :

Andenordens central

På samme måde kan vi anvende andre differentieringsformler på en rekursiv måde.

Anden ordre frem
Anden ordre baglæns

Mere generelt er den n. Rækkefølge fremad, bagud og centrale forskelle givet ved henholdsvis

Frem

eller for h = 1 ,

Baglæns
Central

Disse ligninger bruger binomiske koefficienter efter summeringstegnet vist som (n
i
)
. Hver række iPascals trekantgiver koefficienten for hver værdi påi.

Bemærk, at den centrale forskel for ulige n har h ganget med ikke-heltal. Dette er ofte et problem, fordi det svarer til at ændre diskretiseringsintervallet. Problemet kan afhjælpes ved at tage gennemsnittet af δ n [  f  ] ( x -h/2) og δ n [  f  ] ( x +h/2) .

Fremadrettede forskelle anvendt på en sekvens kaldes undertiden sekvensens binomiske transformation og har en række interessante kombinatoriske egenskaber. Fremadrettede forskelle kan evalueres ved hjælp af Nörlund -Rice -integralet . Den integrerede repræsentation for disse serietyper er interessant, fordi integralet ofte kan evalueres ved hjælp af asymptotisk ekspansion eller sadelpunktsteknikker ; derimod kan fremadforskel -serien være ekstremt svær at evaluere numerisk, fordi binomiske koefficienter vokser hurtigt for store n .

Forholdet mellem disse højere ordensforskelle med de respektive derivater er ligetil,

Højere ordensforskelle kan også bruges til at konstruere bedre tilnærmelser. Som nævnt ovenfor tilnærmer førsteordensforskellen det første ordens derivat op til en ordreperiode h . Dog kombinationen

tilnærmer f  ′ ( x ) op til en ordensperiode h 2 . Dette kan bevises ved at udvide ovenstående udtryk i Taylor -serien eller ved at bruge beregningen af ​​begrænsede forskelle, forklaret nedenfor.

Om nødvendigt kan den begrænsede forskel centreres om et hvilket som helst punkt ved at blande fremad, bagud og centrale forskelle.

Kerner af vilkårlig størrelse

Ved hjælp af lineær algebra kan man konstruere endelige forskelstilnærmelser, der anvender et vilkårligt antal punkter til venstre og et (muligvis forskelligt) antal punkter til højre for evalueringspunktet, for ethvert ordrederivat. Dette indebærer løsning af et lineært system, så Taylor -udvidelsen af summen af ​​disse punkter omkring evalueringspunktet bedst tilnærmer Taylor -ekspansionen af ​​det ønskede derivat. Sådanne formler kan gengives grafisk på et sekskantet eller diamantformet gitter.

Dette er nyttigt for at differentiere en funktion på et gitter, hvor man, når man nærmer sig kanten af ​​gitteret, skal prøve færre og færre punkter på den ene side.

Detaljerne er beskrevet i disse noter .

The Finite Difference Koefficienter Calculator konstruerer endelig differens approksimationer for ikke-standard (og endog ikke-heltallige) stencils givet en vilkårlig stencil og en ønsket afledes orden.

Ejendomme

  • For alle positive k og n

I differentialligninger

En vigtig anvendelse af begrænsede forskelle er i numerisk analyse , især i numeriske differentialligninger , der sigter mod den numeriske løsning af almindelige og partielle differentialligninger . Ideen er at erstatte de derivater, der vises i differentialligningen, med begrænsede forskelle, der tilnærmer dem. De resulterende metoder kaldes endelige forskelsmetoder .

Almindelige anvendelser af metoden med begrænset forskel er inden for beregningsvidenskab og ingeniørdiscipliner, såsom termisk teknik , væskemekanik osv.

Newtons serie

Den Newton-serien består af vilkårene i Newton fremad forskel ligning , opkaldt efter Isaac Newton ; i det væsentlige er det Newton -interpoleringsformlen , der først blev offentliggjort i hans Principia Mathematica i 1687, nemlig den diskrete analog til den kontinuerlige Taylor -ekspansion,

som holder for ethvert polynomium funktion f og for mange (men ikke alle) analytiske funktioner (Det betyder ikke holde, når f er eksponentiel form Dette ses let, som den sinusfunktion forsvinder ved hele multipla af. fremstilles det tilsvarende Newton serien er identisk nul , da alle begrænsede forskelle er nul i dette tilfælde. Men klart er sinusfunktionen ikke nul.). Her udtrykket

er den binomiske koefficient , og

er " faldende faktor " eller "lavere faktor", mens det tomme produkt ( x ) 0 er defineret til at være 1. I dette særlige tilfælde er der en antagelse om enhedstrin for ændringerne i værdierne på x , h = 1 af generaliseringen herunder.

Bemærk den formelle korrespondance af dette resultat til Taylors sætning . Historisk set har dette såvel som Chu – Vandermonde -identiteten ,

(følger heraf og svarende til binomial sætning ), er inkluderet i de observationer, der modnet til systemet med paraplyregning .

For at illustrere, hvordan man kan bruge Newtons formel i praksis, overveje de første par termer for at fordoble Fibonacci -sekvensen f = 2, 2, 4, ... Man kan finde et polynom, der gengiver disse værdier, ved først at beregne en differenstabel, og derefter erstatte de forskelle, der svarer til x 0 (understreget) i formlen som følger,

For ikke -ensartede trin i værdierne x beregner Newton de opdelte forskelle ,

serien af ​​produkter,

og det resulterende polynom er skalarproduktet ,

.

I analysen med p -adic tal , Mahlers sætning hedder det, at den antagelse, at f er et polynomium funktion kan svækkes hele vejen til den antagelse, at f er blot kontinuerlig.

Carlsons sætning giver nødvendige og tilstrækkelige betingelser for, at en Newton -serie kan være unik, hvis den findes. Imidlertid findes der generelt ikke en Newton -serie.

Newton -serien er sammen med Stirling -serien og Selberg -serien et specialtilfælde af den generelle differensserie , som alle er defineret i form af passende skalerede forskelle fremad.

I en komprimeret og lidt mere generel form og ækvidistante noder læser formlen

Beregning af begrænsede forskelle

Forwardforskellen kan betragtes som en operatør , kaldet differensoperator , som kortlægger funktionenftilΔ h [  f  ]. Denne operatør beløber sig til

hvor T h er skifteoperatoren med trin h , defineret af T h [  f  ] ( x ) = f  ( x + h ) , og I er identitetsoperatoren .

Den begrænsede forskel mellem højere ordrer kan defineres på rekursiv måde som Δn
h
≡ Δ hn - 1
time
)
. En anden tilsvarende definition er Δn
h
= [ T h - I ] n
.

Forskellen operatør Δ h er en lineær operator , som sådan opfylder Δ h [ αf + βg ] ( x ) = a Δ h [  f  ] ( x ) + p Δ h [ g ] ( x ) .

Det opfylder også en særlig Leibniz -regel angivet ovenfor, Δ h ( f  ( x ) g ( x )) = (Δ h f  ( x )) g ( x + h ) + f  ( x ) (Δ h g ( x )) . Lignende udsagn gælder for tilbagestående og centrale forskelle.

Formelt at anvende Taylor -serien med hensyn til h , giver formlen

hvor D betegner kontinuumderivatoperatoren, kortlægning f til dens derivat f  ′ . Udvidelsen er gyldig, når begge sider virker på analytiske funktioner i tilstrækkeligt små timer . Således er T h = e hD og formelt invertering af de eksponentielle udbytter

Denne formel holder i den forstand, at begge operatorer giver det samme resultat, når de anvendes på et polynom.

Selv for analytiske funktioner er serien til højre ikke garanteret at konvergere; det kan være en asymptotisk serie . Den kan dog bruges til at opnå mere præcise tilnærmelser til derivatet. For eksempel giver bevarelse af de to første termer i serien andenordens tilnærmelse til f  ′ ( x ) nævnt i slutningen af afsnittet Højere ordenforskelle .

De analoge formler for tilbagestående og centrale differensoperatorer er

Beregningen af ​​begrænsede forskelle er relateret til paraplyregningen for kombinatorik. Denne bemærkelsesværdige systematiske korrespondance skyldes identiteten af kommutatorerne for paraplymængderne til deres kontinuumanaloger ( h → 0 grænser),

Et stort antal formelle differentialforhold i standardregning, der involverer funktioner f  ( x ), kortlægger således systematisk til paraply-endelige forskelsanaloger, der involverer f  ( xT−1
t
)
.

For eksempel er paraplyanalogen for et monomial x n en generalisering af ovenstående faldende faktor ( Pochhammer k-symbol ),

så det

derfor den ovennævnte Newton -interpoleringsformel (ved at matche koefficienter i udvidelsen af ​​en vilkårlig funktion f  ( x ) i sådanne symboler) og så videre.

For eksempel er paraply -sinus

Som i kontinuumgrænsen er egenfunktionen af Δ h/h tilfældigvis også er en eksponentiel,

og derfor kortlægges Fouriers summer af kontinuumfunktioner let til paraply Fourier -summer trofast , dvs. involverer de samme Fourier -koefficienter, der multiplicerer disse paraplybaserede eksponentielle værdier. Denne paraplyeksponentiel udgør således den eksponentielle genererende funktion af Pochhammer -symbolerne .

Således kortlægger f.eks. Dirac delta -funktionen til sin paraplykorrespondent, kardinal sinusfunktionen ,

og så videre. Differenceligninger kan ofte løses med teknikker, der meget ligner dem til løsning af differentialligninger .

Den inverse operatør af den fremadgående differensoperator, så paraplyintegralet, er den ubestemte sum eller antidifferensoperator.

Regler for beregning af endelige forskelsoperatorer

Analogt med reglerne for at finde derivatet har vi:

  • Konstant regel : Hvis c er en konstant , så

Alle ovennævnte regler gælder lige så godt for enhver differensoperator, herunder som Δ .

eller

Se referencer.

Generaliseringer

  • En generaliseret begrænset forskel defineres normalt som
    hvor μ = ( μ 0 ,…, μ N ) er dens koefficientvektor. En uendelig forskel er en yderligere generalisering, hvor den endelige sum ovenfor er erstattet af en uendelig række . En anden måde at generalisere på er at gøre koefficienter μ k afhængige af punkt x : μ k = μ k ( x ) , og dermed overveje en vægtet begrænset forskel . Man kan også få trin h til at afhænge af punkt x : h = h ( x ) . Sådanne generaliseringer er nyttige til konstruktion af forskellige kontinuitetsmoduler .
  • Den generaliserede forskel kan ses som polynomringene R [ T h ] . Det fører til forskel algebraer.
  • Differensoperatør generaliserer til Möbius -inversion over et delvist bestilt sæt .
  • Som en konvolutionsoperatør: Via formalismen af forekomstalgebraer kan differensoperatorer og anden Möbius -inversion repræsenteres ved konvolvering med en funktion på poset, kaldet Möbius -funktionen μ ; for differensoperatoren er μ sekvensen (1, -1, 0, 0, 0,…) .

Multivariate begrænsede forskelle

Endelige forskelle kan betragtes i mere end én variabel. De er analoge med partielle derivater i flere variabler.

Nogle delvise afledte tilnærmelser er:

Alternativt, for applikationer, hvor beregningen af f er det mest omkostningstunge trin, og både første og anden derivat skal beregnes, er en mere effektiv formel for det sidste tilfælde

da de eneste værdier til beregning, der ikke allerede er nødvendige for de foregående fire ligninger, er f  ( x + h , y + k ) og f  ( x - h , y - k ) .

Se også

Referencer

  1. ^ a b c Paul Wilmott; Sam Howison; Jeff Dewynne (1995). The Mathematics of Financial Derivatives: A Student Introduction . Cambridge University Press. s. 137 . ISBN 978-0-521-49789-3.
  2. ^ a b c Peter Olver (2013). Introduktion til delvise differentialligninger . Springer Science & Business Media. s. 182. ISBN 978-3-319-02099-0.
  3. ^ a b c M Hanif Chaudhry (2007). Åben kanalflow . Springer. s. 369. ISBN 978-0-387-68648-6.
  4. ^ Jordán, op. cit., s. 1 og Milne-Thomson, s. xxi. Milne-Thomson, Louis Melville (2000): Beregningen af ​​endelige forskelle (Chelsea Pub Co, 2000) ISBN  978-0821821077
  5. ^ Fraser, Duncan C. (1. januar 1909). "Om den grafiske afgrænsning af interpolationsformulæ" . Journal of the Institute of Actuaries . 43 (2): 235–241. doi : 10.1017/S002026810002494X . Hentet 17. april 2017 .
  6. ^ Newton, Isaac, (1687). Principia , bog III, Lemma V, sag 1
  7. ^ Richtmeyer, D. og Morton, KW, (1967). Difference Methods for Initial Value Problems , 2. udgave, Wiley, New York.
  8. ^ Boole, George , (1872). En afhandling om beregningen af ​​endelige forskelle , 2. udgave, Macmillan og Company. Online . Også [Dover udgave 1960]
  9. ^ Jordan, Charles, (1939/1965). "Beregning af endelige forskelle", Chelsea Publishing. Online: [1]
  10. ^ Zachos, C. (2008). "Paraplydeformationer på diskret rumtid". International Journal of Modern Physics A . 23 (13): 2005–2014. arXiv : 0710.2306 . Bibcode : 2008IJMPA..23.2005Z . doi : 10.1142/S0217751X08040548 . S2CID  16797959 .
  11. ^ Curtright, TL; Zachos, CK (2013). "Paraply Vade Mecum" . Grænser i fysik . 1 : 15. arXiv : 1304.0429 . Bibcode : 2013FrP ..... 1 ... 15C . doi : 10.3389/fphy.2013.00015 . S2CID  14106142 .
  12. ^ Levy, H .; Lessman, F. (1992). Endelige forskelsligninger . Dover. ISBN 0-486-67260-3.
  13. ^ Ames, WF, (1977). Numeriske metoder til partielle differentialligninger , afsnit 1.6. Academic Press, New York. ISBN  0-12-056760-1 .
  14. ^ Hildebrand, FB , (1968). Finite-Difference Equations and Simulations , afsnit 2.2, Prentice-Hall, Englewood Cliffs, New Jersey.
  15. ^ Flajolet, Philippe; Sedgewick, Robert (1995). "Mellin -transformationer og asymptotik: Endelige forskelle og Rices integraler" (PDF) . Teoretisk datalogi . 144 (1–2): 101–124. doi : 10.1016/0304-3975 (94) 00281-M ..
  • Richardson, CH (1954): En introduktion til beregningen af ​​endelige forskelle (Van Nostrand (1954) online kopi
  • Mickens, RE (1991): Difference Equations: Theory and Applications (Chapman og Hall/CRC) ISBN  978-0442001360

eksterne links