Konečný rozdíl - Finite difference
Konečný rozdíl je matematické vyjádření ve tvaru f ( x + b ) - f ( x + ) . Je -li konečný rozdíl vydělen b - a , získá se rozdílový kvocient . Aproximace derivátů podle konečných rozdílů hraje ústřední roli v konečných rozdílných metod pro numerické řešení diferenciálních rovnic , zejména okrajových úloh .
Některé relace opakování lze zapsat jako diferenciální rovnice nahrazením iteračního zápisu konečnými rozdíly.
Dnes je termín „konečný rozdíl“ často brán jako synonymum pro aproximace konečných rozdílů derivátů , zejména v kontextu numerických metod . Aproximace konečných rozdílů jsou kvocienty konečných rozdílů ve výše použité terminologii.
Konečné rozdíly zavedl Brook Taylor v roce 1715 a byly také studovány jako abstraktní samostatně stojící matematické objekty v dílech George Boole (1860), LM Milne-Thomson (1933) a Károly Jordan (1939). Konečné rozdíly vysledují jejich původ zpět k jednomu z algoritmů Josta Bürgiho ( c. 1592 ) a k práci jiných, včetně Isaaca Newtona . Formální počet konečných rozdílů lze považovat za alternativu k počtu nekonečně malých .
Základní typy
Obvykle se uvažují tři základní typy: dopředné , zpětné a centrální konečné rozdíly.
Dopředu rozdíl , označený na funkce f je funkce definovaná jako
V závislosti na aplikaci může být rozteč h proměnná nebo konstantní. Pokud je vynecháno, h je považováno za 1; to znamená,
Dozadu rozdíl používá hodnoty funkce v x a x - h , místo hodnot na x + h a x :
Nakonec je centrální rozdíl dán vztahem
Vztah s deriváty
Konečný rozdíl se často používá jako aproximace derivátu, obvykle v numerické diferenciaci .
Derivát z funkce f v bodě x je definován limitem .
Pokud má h místo blížící se nuly pevnou (nenulovou) hodnotu, pak by byla napsána pravá strana výše uvedené rovnice
Dopředný rozdíl dělený h se tedy blíží derivaci, když h je malá. Chybu v této aproximaci lze odvodit z Taylorovy věty . Za předpokladu, že f je dvakrát diferencovatelné, máme
Stejný vzorec platí pro zpětný rozdíl:
Centrální (také nazývaný středový) rozdíl však poskytuje přesnější aproximaci. Pokud je f třikrát diferencovatelné,
Hlavním problémem metody centrální diference je však to, že oscilační funkce mohou poskytnout nulovou derivaci. Pokud f ( nh ) = 1 pro n liché, a f ( nh ) = 2 pro n sudé, pak f ′ ( nh ) = 0, pokud se vypočítá pomocí centrálního rozdílového schématu . To je obzvláště problematické, pokud je doména f diskrétní. Viz také symetrický derivát
Autoři, pro které konečné rozdíly znamenají aproximace konečných rozdílů, definují dopředné/zpětné/centrální rozdíly jako kvocienty uvedené v této části (namísto použití definic uvedených v předchozí části).
Rozdíly vyššího řádu
Analogickým způsobem lze získat aproximace konečných rozdílů k derivacím vyšších řádů a diferenciálním operátorům. Například pomocí výše uvedeného centrálního rozdílového vzorce pro f '( x + h/2) a f '( x -h/2) a použitím vzorce centrální diference pro derivaci f ' v x získáme aproximaci centrálního rozdílu druhé derivace f :
- Centrála druhého řádu
Podobně můžeme rekurzivně aplikovat další diferenční vzorce.
- Druhá objednávka vpřed
- Druhý řád zpět
Obecněji platí, že n -té pořadí vpřed, vzad a centrální rozdíly jsou dány, resp.
- Vpřed
nebo pro h = 1 ,
- Dozadu
- Centrální
Tyto rovnice používají binomické koeficienty za součtovým znaménkem zobrazeným jako (n
i) . Každý řádekPascalova trojúhelníkuposkytuje koeficient pro každou hodnotui.
Všimněte si, že střední rozdíl bude, pro liché n , mají h vynásobený non-celá čísla. To je často problém, protože to znamená změnu intervalu diskretizace. Problém lze vyřešit průměrem δ n [ f ] ( x -h/2) a δ n [ f ] ( x +h/2) .
Dopředné rozdíly aplikované na sekvenci se někdy nazývají binomická transformace sekvence a mají řadu zajímavých kombinatorických vlastností. Předběžné rozdíly lze vyhodnotit pomocí integrálu Nörlund – Rice . Integrální reprezentace pro tyto typy sérií je zajímavá, protože integrál lze často vyhodnotit pomocí asymptotické expanze nebo technik sedlového bodu ; naproti tomu dopředné diferenční řady může být extrémně těžké numericky vyhodnotit, protože binomické koeficienty rychle rostou pro velká n .
Vztah těchto rozdílů vyššího řádu s příslušnými deriváty je přímočarý,
Rozdíly vyššího řádu lze také použít ke konstrukci lepších aproximací. Jak bylo uvedeno výše, rozdíl prvního řádu se blíží derivaci prvního řádu až do termínu řádu h . Nicméně kombinace
aproximuje f '( x ) až do termínu řádu h 2 . To lze dokázat rozšířením výše uvedeného výrazu v Taylorových řadách nebo pomocí počtu konečných rozdílů, vysvětleného níže.
V případě potřeby lze konečný rozdíl vycentrovat na jakýkoli bod smícháním dopředných, zpětných a centrálních rozdílů.
Jádra libovolné velikosti
Pomocí lineární algebry lze konstruovat aproximace konečných rozdílů, které využívají libovolný počet bodů vlevo a (případně jiný) počet bodů vpravo od vyhodnocovacího bodu pro libovolnou řádovou derivaci. To zahrnuje řešení lineárního systému tak, aby Taylorova expanze součtu těchto bodů kolem vyhodnocovacího bodu nejlépe odpovídala Taylorově expanzi požadované derivace. Takové vzorce lze znázornit graficky na šestiúhelníkové nebo kosočtverečné mřížce.
To je užitečné pro odlišení funkce na mřížce, kde když se člověk přiblíží k okraji mřížky, musí na jedné straně vzorkovat méně a méně bodů.
Podrobnosti jsou uvedeny v těchto poznámkách .
Konečných diferencí Koeficienty kalkulačka konstrukce konečných rozdílných přiblížení pro nestandardní (a dokonce i non-celočíselné) šablony daných libovolné šablony a požadovaný derivát pořadí.
Vlastnosti
- Pro všechna kladná k a n
V diferenciálních rovnicích
Důležitá aplikace konečných rozdílů je v numerické analýze , zejména v numerických diferenciálních rovnicích , jejichž cílem je numerické řešení obyčejných a parciálních diferenciálních rovnic . Cílem je nahradit deriváty objevující se v diferenciální rovnici konečnými rozdíly, které je aproximují. Výsledné metody se nazývají metody konečných rozdílů .
Běžné aplikace metody konečných rozdílů jsou ve výpočetních vědách a technických oborech, jako je tepelné inženýrství , mechanika tekutin atd.
Newtonova série
Série Newton se skládá z výrazů Newtonovy dopředné diferenční rovnice pojmenované po Isaacovi Newtonovi ; v podstatě je to Newtonův interpolační vzorec , poprvé publikovaný v jeho Principia Mathematica v roce 1687, konkrétně diskrétní analog kontinuální Taylorovy expanze,
který platí pro jakoukoli polynomickou funkci f a pro mnoho (ale ne všechny) analytické funkce (neplatí, když f je exponenciální typ . To je snadno vidět, protože sinusová funkce mizí při celočíselných násobcích ; odpovídající Newtonova řada je identicky nulová , protože všechny konečné rozdíly jsou v tomto případě nulové. Přesto je jasné, že sinusová funkce není nulová.). Tady ten výraz
je binomický koeficient a
je " klesající faktoriál " nebo "nižší faktoriál", zatímco prázdný součin ( x ) 0 je definován jako 1. V tomto konkrétním případě existuje předpoklad jednotkových kroků pro změny hodnot x , h = 1 generalizace níže.
Všimněte si formální shody tohoto výsledku s Taylorovou větou . Historicky to, stejně jako identita Chu -Vandermonde ,
(z toho vyplývající a odpovídající binomické větě ) jsou zahrnuty v pozorováních, která dospěla do systému umbrálního počtu .
Abychom ilustrovali, jak lze ve skutečné praxi použít Newtonův vzorec, zvažte několik prvních podmínek zdvojnásobení Fibonacciho posloupnosti f = 2, 2, 4, ... Polynom, který tyto hodnoty reprodukuje, lze najít tak , že nejprve vypočítáme rozdílovou tabulku, a poté nahrazení rozdílů, které odpovídají x 0 (podtrženo) do vzorce následujícím způsobem,
Pro případ nejednotných kroků v hodnotách x vypočítá Newton dělené rozdíly ,
řada produktů,
a výsledný polynom je skalární součin ,
- .
Při analýze se p -adic čísly , Mahlerův teorém říká, že předpoklad, že f je polynom funkce může být oslabena až k předpokladu, že f je pouze spojité.
Carlsonova věta poskytuje nezbytné a dostatečné podmínky pro to, aby byla Newtonova řada jedinečná, pokud existuje. Newtonova řada však obecně neexistuje.
Série Newton, společně se sériemi Stirling a Selberg , je zvláštním případem řady obecných rozdílů , z nichž všechny jsou definovány z hlediska vhodně zmenšených dopředných rozdílů.
V komprimované a mírně obecnější formě a stejně vzdálených uzlech vzorec čte
Počet konečných rozdílů
Dopředný rozdíl lze považovat za operátor , nazývaný diferenční operátor , který mapuje funkcifnaΔ h [ f ]. Tento operátor činí
kde T h je operátor posunu s krokem h , definovaným T h [ f ] ( x ) = f ( x + h ) , a I je operátor identity .
Konečný rozdíl vyšších řádů lze definovat rekurzivně jako Δn
h≡ Δ h (Δn - 1
hodina) . Další ekvivalentní definice je Δn
h= [ T h - I ] n .
Rozdíl operátor Δ h je lineární operátor , jako takový, že splňuje delta h [ αf + βg ] ( x ) = alfa Δ h [ f ] ( x ) + beta Δ h [ g ] ( x ) .
Splňuje také speciální Leibnizovo pravidlo uvedené výše, Δ h ( f ( x ) g ( x )) = (Δ h f ( x )) g ( x + h ) + f ( x ) (Δ h g ( x )) . Podobná tvrzení platí pro zpětné a centrální rozdíly.
Formálně použitím Taylorovy řady s ohledem na h se získá vzorec
kde D označuje derivační operátor kontinua, mapování f na jeho derivát f ' . Rozšíření je platné, když obě strany působí na analytické funkce , po dostatečně malou h . Tak, T h = e hD , a formálně invertující exponenciální výtěžky
Tento vzorec platí v tom smyslu, že oba operátory dávají při použití na polynom stejný výsledek.
I u analytických funkcí není zaručeno, že se řady vpravo budou sbíhat; může to být asymptotická série . Lze jej však použít k získání přesnějších aproximací derivátu. Například zachováním prvních dvou členů řady vznikne aproximace druhého řádu na f '( x ) zmíněná na konci části Rozdíly vyššího řádu .
Analogické vzorce pro operátory zpětného a centrálního rozdílu jsou
Počet konečných rozdílů souvisí s umbrálním počtem kombinatoriky. Tato pozoruhodně systematická korespondence je dána identitou komutátorů umbrálních veličin k jejich analogům kontinua ( h → 0 limitů),
Velký počet formálních diferenciálních vztahů standardního počtu zahrnujícího funkce f ( x ) tak systematicky mapuje analogie umbrálních konečných rozdílů zahrnující f ( xT-1
h) .
Například umbrální analog monomiální x n je zobecněním výše uvedeného klesajícího faktoriálu ( Pochhammerův symbol k ),
aby
proto výše uvedený Newtonův interpolační vzorec (porovnáním koeficientů při expanzi libovolné funkce f ( x ) v takových symbolech) atd.
Například umbrální sinus je
Stejně jako v limitu kontinua, vlastní funkce Δ h/h je také exponenciální,
a proto jsou Fourierovy součty funkcí kontinua snadno mapovány na umbrální Fourierovy součty věrně , tj. zahrnující stejné Fourierovy koeficienty násobící tyto exponenciály umbrálního základu. Tato umbral exponenciální tak představuje exponenciální funkci generující z symbolů Pochhammer .
Funkce Dirac delta se tedy například mapuje na svého mozkového korespondenta, kardinální sinusovou funkci ,
a tak dále. Diferenční rovnice lze často řešit technikami velmi podobnými technikám pro řešení diferenciálních rovnic .
Inverzní operátor dopředného rozdílového operátoru, tedy pak umbrální integrál, je neurčitý součet nebo antidiferenční operátor.
Pravidla pro počet operátorů konečných rozdílů
Analogicky s pravidly pro hledání derivátu máme:
- Konstantní pravidlo : Pokud c je konstanta , pak
Všechna výše uvedená pravidla platí stejně dobře pro všechny rozdílové operátory, včetně ∇ jako Δ .
- nebo
Viz reference.
Zobecnění
- Zobecnit konečný rozdíl je obvykle definována jako
- Zobecněný rozdíl lze vidět jako polynomické prstence R [ T h ] . Vede k rozdílovým algebrám.
- Operátor rozdílu generalizuje inverzi Möbius přes částečně uspořádanou sadu .
- Jako konvoluční operátor: Prostřednictvím formalismu incidenčních algeber mohou být diferenční operátory a další Möbiova inverze reprezentovány konvolucí s funkcí na posetu, nazývanou Möbiova funkce μ ; pro diferenční operátor μ je posloupnost (1, −1, 0, 0, 0, ...) .
Vícenásobné konečné rozdíly
Konečné rozdíly lze uvažovat ve více než jedné proměnné. Jsou analogické parciálním derivacím v několika proměnných.
Některé částečné derivační aproximace jsou:
Alternativně pro aplikace, ve kterých je výpočet f nejnákladnějším krokem a musí být vypočítán první i druhý derivát, je účinnějším vzorcem pro poslední případ
protože jediné hodnoty pro výpočet, které již nejsou potřebné pro předchozí čtyři rovnice, jsou f ( x + h , y + k ) a f ( x - h , y - k ) .
Viz také
Reference
- ^ a b c Paul Wilmott; Sam Howison; Jeff Dewynne (1995). Matematika finančních derivátů: studentský úvod . Cambridge University Press. p. 137 . ISBN 978-0-521-49789-3.
- ^ a b c Peter Olver (2013). Úvod do parciálních diferenciálních rovnic . Springer Science & Business Media. p. 182. ISBN 978-3-319-02099-0.
- ^ a b c M Hanif Chaudhry (2007). Tok otevřeného kanálu . Springer. p. 369. ISBN 978-0-387-68648-6.
- ^ Jordán, op. cit., s. 1 a Milne-Thomson, s. xxi. Milne-Thomson, Louis Melville (2000): The Calculus of Finite Differences (Chelsea Pub Co, 2000) ISBN 978-0821821077
- ^ Fraser, Duncan C. (1. ledna 1909). „K grafickému vymezení interpolačního formuláře“ . Časopis Institutu pojistných matematiků . 43 (2): 235–241. doi : 10,1017/S002026810002494X . Citováno 17. dubna 2017 .
- ^ Newton, Isaac, (1687). Principia , kniha III, Lemma V, případ 1
- ^ Richtmeyer, D. a Morton, KW, (1967). Diferenční metody pro problémy počáteční hodnoty , 2. vydání, Wiley, New York.
- ^ Boole, George , (1872). Pojednání o počtu konečných rozdílů , 2. vydání, Macmillan a společnost. Online . Také [Dover edition 1960]
- ^ Jordan, Charles, (1939/1965). „Kalkul konečných rozdílů“, Chelsea Publishing. On-line: [1]
- ^ Zachos, C. (2008). „Umbrální deformace na diskrétním časoprostoru“. International Journal of Modern Physics A . 23 (13): 2005–2014. arXiv : 0710.2306 . Bibcode : 2008IJMPA..23.2005Z . doi : 10,1142/S0217751X08040548 . S2CID 16797959 .
- ^ Curtright, TL; Zachos, CK (2013). „Umbral Vade Mecum“ . Hranice ve fyzice . 1 : 15. arXiv : 1304.0429 . Bibcode : 2013FrP ..... 1 ... 15C . doi : 10,3389/fphy.2013.00015 . S2CID 14106142 .
- ^ Levy, H .; Lessman, F. (1992). Rovnice konečných rozdílů . Dover. ISBN 0-486-67260-3.
- ^ Ames, WF, (1977). Numerické metody pro parciální diferenciální rovnice , oddíl 1.6. Academic Press, New York. ISBN 0-12-056760-1 .
- ^ Hildebrand, FB , (1968). Rovnice a simulace konečných rozdílů , oddíl 2.2, Prentice-Hall, Englewood Cliffs, New Jersey.
- ^ Flajolet, Philippe; Sedgewick, Robert (1995). „Mellinovy transformace a asymptotika: konečné rozdíly a Riceovy integrály“ (PDF) . Teoretická počítačová věda . 144 (1–2): 101–124. doi : 10,1016/0304-3975 (94) 00281-M ..
- Richardson, CH (1954): Úvod do kalkulu konečných rozdílů (Van Nostrand (1954) online kopie
- Mickens, RE (1991): Diferenční rovnice: Teorie a aplikace (Chapman a Hall/CRC) ISBN 978-0442001360