Eulerin menetelmä

Wikipedia
Loikkaa: valikkoon, hakuun
Tutkittavan sinisen kuvaajan yhtälö on tuntematon. Punainen kuvaaja on Eulerin menetelmän avulla tehty approksimaatio sinisestä kuvaajasta.

Numeerisessa analyysissä Eulerin menetelmällä tarkoitetaan Leonhard Eulerin mukaan nimettyä menetelmää differentiaaliyhtälöiden alkuarvo-ongelmien ratkaisuun. Se on yksinkertaisin numeerisista integrointimenetelmistä, mutta siitä huolimatta menetelmä on tärkeä. Sen kautta on johdettu monet paremmista menetelmistä ja sen ymmärtäminen on pohjana näille.

Johtaminen[muokkaa | muokkaa wikitekstiä]

Eulerin menetelmä johdetaan derivaatan määritelmästä

 f'(x) = \lim_{h\rightarrow 0}\frac{f(x+h) - f(x)}{h}.

Nyt ei kuitenkaan anneta muuttujan h lähestyä nollaa, vaan pidetään se pienenä. Näin saadaan yhtälöstä

 f'(x) = \frac{f(x + h) - f(x)}{h},

jossa h tarkoittaa menetelmässä käytettyä askelta. Olkoot alkuarvo-ongelmamme muotoa

 y'(t) = F(t, y(t)),\quad y(t_0) = y_0.

Nyt siis funktiolla F tarkoitetaan muuttujista t ja y(t) koostuvaa lauseketta, joka tulee ylläolevan yhtälön oikeaksi puoleksi. Määritetään nyt, että Eulerin menetelmällä askel ajanhetkestä tn ajanhetkeen tn + 1 voidaan laskea kaavalla

 y_{n+1} = y_n + h\cdot F(t_n, y_n).[1]


Virhe[muokkaa | muokkaa wikitekstiä]

Eulerin menetelmän virhe voidaan laskea selvittämällä Taylorin sarja muuttujalle y. Tulkitsemalla ylläolevaa tulosta voidaan menetelmä merkitä muodossa

 y(t_0 + h ) = y(t_0) + hy'(t_0).

Tässä on kaksi ensimmäistä termiä Taylorin sarjasta

 y(t_0 + h ) = y(t_0) + hy'(t_0) + \frac{1}{2}h^2y''(t_0) + \mathcal{O}(h^3)

ja yhtälöiden välinen virhe saadaan siis suuruudeltaan vastaamaan termejä

 \frac{1}{2}h^2y''(t_0) + \mathcal{O}(h^3).

Nähdään, että käytettäessä pientä askelen h arvoa on menetelmän virheessä dominoiva termi suhteessa verrannollinen h2:n. Ratkaistessa menetelmällä alkuarvotehtävää on tarvittavien askelien määrä kääntäen verrannollinen askeleeseen h. Näin on pääteltävissä, että kokonaisvirhe Eulerin menetelmässä on verrannollinen askeleeseen h (virhe askeleessa kertaa askelien määrä). Eulerin menetelmää kutsutaankin ensimmäisen kertaluokan menetelmäksi.

Tarkkuus on parannettavissa ottamalla huomioon Taylorin sarjasta lisää termejä. Tämä kuitenkin saattaa olla hankalaa, sillä niitä varten pitäisi ratkaista muuttujan y(t) derivaattoja. Tämän takia onkin kehitetty useaa ajankohtaa alkuarvoinaan käyttäviä multistep-menetelmiä tai yhden ajanhetken derivaattaa useammassa pisteessä approksimoivia menetelmiä kuten esim. Runge-Kutta-menetelmät.

Lähteet[muokkaa | muokkaa wikitekstiä]

  1. Adams, Robert A.: ”17.3.: "Existence, Uniqueness and Numerical Methods”, Calculus: A Complete Course, s. 911. Pearson: Adisson Wesley, 6. painos.

Katso myös[muokkaa | muokkaa wikitekstiä]