V numerické matematice je numerické řešení obyčejných diferenciálních rovnic postup, kterým můžeme získat přibližné řešení obyčejných diferenciálních rovnic. Používá se v případech, kdy by bylo nalezení přesného (analytického) řešení náročné nebo v případech, kdy analytické řešení nelze najít.
Diferenciální rovnice a její počáteční podmínky bývají často uváděny v tomto tvaru:
![](https://wikimedia.org/api/rest_v1/media/math/render/svg/8499f43d20ed2abf557d5a6989eb40b868b1823d)
![](https://wikimedia.org/api/rest_v1/media/math/render/svg/3ba6cc4b1d5692b5d384b9e1f93705ae00f05cbc)
Metody řešení
Funkce f(t,y) (někdy se nazývá stavová rovnice) může být obecně velmi komplikovaná, proto je nutné řešit rovnici numericky. V takovém případě probíhá řešení v diskrétních časových krocích
:
![](https://wikimedia.org/api/rest_v1/media/math/render/svg/05d137771a82437e1c5e729a27032cc152ec21a9)
je funkce (někdy též směrová funkce), která se snaží aproximovat
tak, aby
bylo co nejpřesnější.
Eulerova metoda
Existuje více metod, jak v daném čase získat co nejlepší aproximaci derivace, nejjednodušší je Eulerova metoda:
![](https://wikimedia.org/api/rest_v1/media/math/render/svg/23c6d2d34c52747e656a0b2f8ba20860680fcbc0)
Rungeovy–Kuttovy metody
Obecně lze Rungeovy–Kuttovy metody zapsat následovně:
![](https://wikimedia.org/api/rest_v1/media/math/render/svg/1a9f5cb628eda139d979c70b0eb03fa1dbd7433b)
![](https://wikimedia.org/api/rest_v1/media/math/render/svg/20930be5b1560202a1e8bd7446a16289aea8e1fa)
Koeficienty u těchto metod jsou vypočteny tak, aby metoda řádu
odpovídala Taylorovu polynomu funkce
stejného řádu. (Eulerova metoda je vlastně metodou prvního řádu.)
Často se používá čtyřbodová metoda Runge-Kutta (RK4), která je čtvrtého řádu.
![](https://wikimedia.org/api/rest_v1/media/math/render/svg/9c29c02f9db39c5af84fa1f1fcfa64aa6069b2ad)
![](https://wikimedia.org/api/rest_v1/media/math/render/svg/c05fe7e25aad6c7a55aafaf97aa86a67f1422d0e)
![](https://wikimedia.org/api/rest_v1/media/math/render/svg/d8d07ecfc6e2c0d19352badb655a1e3f78895fcd)
![](https://wikimedia.org/api/rest_v1/media/math/render/svg/2b553256b19a3d88fca0c12c76315ed303c1c178)
![](https://wikimedia.org/api/rest_v1/media/math/render/svg/9c6c4ccafc1de85e4638275f5160abaea4889b2a)
- (Korespondence různých způsobů zápisu:
;
;
;
. Korespondence s obecným vzorcem:
;
;
;
;
.)
Vícekrokové metody
U vícekrokových metod je hodnota
vypočtena z předchozích hodnot
(respektive
,
) proložených interpolačním polynomem. Řád metody zde odpovídá řádu interpolačního polynomu. (Eulerova metoda je v podstatě jednokrokovou metodou.)
Obecnou vícekrokovou metodu lze zapsat následovně:
![](https://wikimedia.org/api/rest_v1/media/math/render/svg/55653b0b7ba285e319b5f628a5c7817111cf4812)
Explicitní metody
Pokud je
, lze hodnotu
určit z
předchozích hodnot
(respektive z
předchozích hodnot
) a jedná se o metodu explicitní.
Příklad 1, explicitní metoda Adams-Bashford druhého řádu:
![](https://wikimedia.org/api/rest_v1/media/math/render/svg/f1036e1a22e655e785cba04fb9e5112947aefc5f)
- (Korespondence s obecným vzorcem:
;
;
;
;
;
.)
Příklad 2, explicitní metoda Adams-Bashford čtvrtého řádu:
![](https://wikimedia.org/api/rest_v1/media/math/render/svg/baf28e252beb46c544ba82813c9ecf6957059b08)
Implicitní metody
Pokud je
různé od nuly, je pro výpočet
nutná znalost
a jedná se o metodu implicitní.
Příklad, implicitní metoda Adams-Moulton čtvrtého řádu:
![](https://wikimedia.org/api/rest_v1/media/math/render/svg/cb95e2d16a00dbd3222414d7b2af94c163e62b16)
Metody prediktor-korektor
Metody prediktor-korektor jsou sloučením explicitních a implicitních metod. Nejprve je použita explicitní metoda pro odhad nového
. V tomto bodě je vypočtena derivace
, která je následovně použita v implicitní metodě pro výpočet přesnější aproximace
.
Související články