- Čo je to Eulerova metóda?
- Riešené cvičenia
- Cvičenie 1
- Riešenie
- Cvičenie 2
- Riešenie
- Cvičenie 3
- Riešenie
- Newtonovská dynamika a Eulerova metóda
- Cvičenie 4
- Riešenie
- Navrhované cvičenia pre domácnosť
- Cvičenie 1
- Cvičenie 2
- Referencie
Eulerova metóda je najjednoduchšie a jednoduché postupy používané na nájdenie numerická riešeniu približné na obyčajné diferenciálne rovnice v prvom poradí, za predpokladu, že počiatočný stav je známy.
Bežná diferenciálna rovnica (ODE) je rovnica, ktorá spája neznámu funkciu jednej nezávislej premennej s jej derivátmi.

Postupné aproximácie podľa Eulerovej metódy. Zdroj: Oleg Alexandrov
Ak je najväčšou deriváciou, ktorá sa objavuje v rovnici, stupeň jedna, potom je to obyčajná diferenciálna rovnica prvého stupňa.
Najbežnejším spôsobom, ako napísať rovnicu prvého stupňa, je:

x = x 0
y = y 0
Čo je to Eulerova metóda?
Eulerova metóda spočíva v hľadaní numerického riešenia diferenciálnej rovnice v intervale medzi X 0 a X f .
Najprv je interval diskretizovaný v n + 1 bodoch:
x 0 , x 1 , x 2 , x 3 …, x n
Získajú sa takto:
x i = x 0 + ih
Kde h je šírka alebo krok podintervalov:

Pri počiatočnom stave je potom možné derivát poznať aj na začiatku:
y '(x o ) = f (x o , y o )
Tento derivát predstavuje sklon dotyčnice k krivke funkcie y (x) presne v bode:
Ao = (x o , y o )
Potom sa urobí približná predpoveď hodnoty funkcie y (x) v nasledujúcom bode:
y (x 1 ) ≈ y 1
y 1 = y o + (x 1 - x o ) f (x o , y o ) = y o + hf (x o , y o )
Získal sa ďalší približný bod riešenia, ktorý by zodpovedal:
A 1 = (x 1 , y 1 )
Postup sa opakuje, aby sa získali nasledujúce body
2 , A 3 …, x n
Na obrázku znázornenom na začiatku predstavuje modrá krivka presné riešenie diferenciálnej rovnice a červená krivka predstavuje približné približné body získané Eulerovou procedúrou.
Riešené cvičenia
Cvičenie 1
I ) Nech je diferenciálna rovnica:

Pri počiatočných podmienkach x = a = 0; a a = 1
Pomocou Eulerovej metódy získame približné riešenie y na súradnici X = b = 0,5, pričom interval rozdelíme na n = 5 častí.
Riešenie

Numerické výsledky sú zhrnuté takto:

Z čoho vyplýva, že riešenie Y pre hodnotu 0,5 je 1,4851.
Poznámka: Na vykonanie výpočtov sa použil program Smath Studio, bezplatný program na bezplatné použitie.
Cvičenie 2
II ) Pokračovanie v diferenciálnej rovnici z cvičenia I), nájdite presné riešenie a porovnajte ho s výsledkom získaným Eulerovou metódou. Nájdite chybu alebo rozdiel medzi presným a približným výsledkom.
Riešenie

Presné riešenie nie je ťažké nájsť. Derivát funkcie sin (x) je známy ako funkcia cos (x). Preto riešením y (x) bude:
y (x) = sin x + C
Aby sa splnila počiatočná podmienka a (0) = 1, musí sa konštanta C rovnať 1. Presný výsledok sa potom porovná s približným:

Dospelo sa k záveru, že v vypočítanom intervale má aproximácia tri významné čísla presnosti.
Cvičenie 3
III ) Zoberme diferenciálnu rovnicu a jej počiatočné podmienky uvedené nižšie:
y, (x) = - y 2
Pri počiatočnej podmienke x 0 = 0; a 0 = 1
Pomocou Eulerovej metódy nájdite približné hodnoty riešenia y (x) v intervale x =. Použite krok h = 0,1.
Riešenie
Eulerova metóda je veľmi vhodná na použitie s tabuľkovým procesom. V tomto prípade použijeme tabuľku geogebra, bezplatný a open-source program.
Tabuľka na obrázku zobrazuje tri stĺpce (A, B, C), prvý je premenná x, druhý stĺpec predstavuje premennú y a tretí stĺpec je derivát y '.
Riadok 2 obsahuje počiatočné hodnoty X, Y, Y '.
Krok hodnoty 0,1 bol umiestnený do bunky absolútnej polohy ($ D $ 4).
Počiatočná hodnota y0 je v bunke B2 a y1 je v bunke B3. Na výpočet y 1 sa použije vzorec:
y 1 = y o + (x 1 - x o ) f (x o , y o ) = y o + hf (x o , y o )
Tento tabuľkový vzorec by mal byť číslo B3: = B2 + $ D $ 4 * C3.
Podobne by y2 bol v bunke B4 a jeho vzorec je uvedený na nasledujúcom obrázku:

Obrázok tiež ukazuje graf presného riešenia a body A, B, …, P približného riešenia pomocou Eulerovej metódy.
Newtonovská dynamika a Eulerova metóda
Klasickú dynamiku vyvinul Isaac Newton (1643 - 1727). Pôvodnou motiváciou Leonarda Eulera (1707 - 1783) k rozvoju tejto metódy bolo práve riešenie rovnice druhého Newtonovho zákona v rôznych fyzikálnych situáciách.
Newtonov druhý zákon sa zvyčajne vyjadruje ako diferenciálna rovnica druhého stupňa:

Kde x predstavuje polohu objektu v čase t. Tento predmet má hmotnosť ma je vystavený pôsobeniu sily F. Funkcia f sa týka sily a hmotnosti nasledovne:


Na použitie Eulerovej metódy sú potrebné počiatočné hodnoty času t, rýchlosti v a polohy x.
Nasledujúca tabuľka vysvetľuje, ako možno z počiatočných hodnôt t1, v1, x1 dosiahnuť aproximáciu rýchlosti v2 a polohy x2 v okamihu t2 = t1 + Δt, kde Δt predstavuje malé zvýšenie a zodpovedá kroku v metóde Euler.

Cvičenie 4
IV ) Jedným zo základných problémov v mechanike je to, že blok hmoty M je viazaný na pružinu (alebo pružinu) elastickej konštanty K.
Newtonov druhý zákon o tomto probléme by vyzeral takto:

V tomto príklade vezmeme pre jednoduchosť M = 1 a K = 1. Približné riešenia polohy xa rýchlosti v pomocou Eulerovej metódy nájdite v časovom intervale rozdelením intervalu na 12 častí.
Zoberte 0 ako počiatočný okamih, počiatočnú rýchlosť 0 a počiatočnú polohu 1.
Riešenie

Numerické výsledky sú uvedené v nasledujúcej tabuľke:

Zobrazujú sa aj grafy polohy a rýchlosti medzi časmi 0 a 1,44.


Navrhované cvičenia pre domácnosť
Cvičenie 1
Na určenie približného riešenia pomocou Eulerovej metódy pre diferenciálnu rovnicu použite tabuľku:
y '= - Exp (-y) s pôvodnými podmienkami x = 0, y = -1 v intervale x =
Začnite krokom 0,1. Vyneste výsledok.
Cvičenie 2
Pomocou tabuľky nájdite numerické riešenia nasledujúcej kvadratickej rovnice, kde y je funkciou nezávislej premennej t.
y '' = - 1 / y² s pôvodnými podmienkami t = 0; a (0) = 0,5; y '(0) = 0
Nájdite riešenie v intervale pomocou kroku 0,05.
Vyneste výsledok: y vs t; y 'vs t
Referencie
- Metóda Eurler Prevzatá z wikipedia.org
- Eulerov riešiteľ. Prevzaté z en.smath.com
