Ratkaiseminen : Numeerisen ratkaisemisen menetelmät

Parannettu Eulerin menetelmä

Integroimalla differentiaaliyhtälö y' = f(x, y(x)) puolittain välin [xk, xk+1] yli saadaan

y(xk+1) - y(xk) =  integral  xk+1

 xkf(x, y(x)) dx.

Kun funktion y arvot korvataan niiden approksimaatioilla ja integraalin arvioimiseen käytetään puolisuunnikassääntöä, päädytään kaavaan

yk+1 = yk + 1
2h[f (xk,yk) + f(xk+1,yk+1)] .

Tämä ei suoranaisesti anna mahdollisuutta differentiaaliyhtälön numeeriseen ratkaisemiseen, koska yk+1 esiintyy myös kaavan oikealla puolella.

Kaavaa voidaan kuitenkin soveltaa siten, että arvolle yk+1 lasketaan ensin karkeampi approksimaatio tavallisella Eulerin menetelmällä ja tätä käytetään eo. kaavan oikealla puolella. Kaavaa käytetään siis karkean arvon tarkentamiseen.

Laskeminen tapahtuu kahdessa vaiheessa:

~yk+1 = yk + hf(xk, yk),
yk+1 = yk + 12h[f(xk,yk) + f(xk+1,y~k+1)] .
Edellistä kaavaa kutsutaan ennustajaksi ja jälkimmäistä korjaajaksi. Menetelmän nimenä on parannettu Eulerin menetelmä tai myös Heunin menetelmä. Voidaan osoittaa, että se on jossain määrin tarkempi kuin tavallinen Eulerin menetelmä. Kumpikaan ei kuitenkaan ole käyttökelpoinen todellisissa ongelmissa. Parannettu Eulerin menetelmä voidaan johtaa myös Taylorin kaavan avulla. Alkuarvoprobleeman

y' = f(x, y),   y(x0) = y0

ratkaisun toisen asteen Taylorin kehitelmä pisteessä xk on

y(xk + h)= y(xk) + hy'(xk) + 12y''(xk)h2 + O(h3)
= y(xk) + [hf (x  ,y(x )) + 1h2(f (x ,y(x  )) + f (x ,y(x ))f (x ,y(x )))]
      k    k     2    x  k     k     y  k     k     k    k + O(h3).
On mahdollista (vaikka ei aivan yksinkertaista) osoittaa, että hakasulkulauseketta voidaan tarkkuutta menettämättä approksimoida lausekkeella

1
2h[f (x ,y(x )) + f(x  + h,y(x  ) + hf (x ,y(x )))]
    k    k        k         k        k    k .

Tämä johtaa samaan kaksivaiheiseen kaavaan kuin edellä.


Teoria: numeerisen ratkaisemisen perusidea
Ratkaiseminen: tavallinen Eulerin menetelmä
Esimerkki: alkuarvoprobleeman numeerinen ratkaiseminen
Esimerkki: Airyn yhtälön numeerinen ratkaiseminen

SKK 15.5.2001