Runge-Kutta meetodid

Allikas: Vikipeedia

Runge-Kutta meetodid on arvutusmatemaatikas algoritmide pere harilike diferentsiaalvõrrandite (ja harilike diferentsiaalvõrrandite süsteemide) ilmutatud või ilmutamata ligikaudse lahendi numbriliseks leidmiseks algtingimustega ülesande korral. Nad põhinevad iteratsioonil.

Meetodi algse kuju töötas välja saksa matemaatik Carl Runge 1895 ning seda üldistas Martin Wilhelm Kutta 1901.

Klassikaline 4. järku Runge-Kutta meetod[muuda | redigeeri lähteteksti]

Next.svg Pikemalt artiklis Klassikaline Runge-Kutta meetid

4. järku klassikaline Runge-Kutta meetod on nii laialt levinud, et seda nimetatakse sageli lihtsalt Runge-Kutta meetodiks.

Olgu meil Cauchy ülesanne \textbf{y}'=\textbf{f}(x,\textbf{y}), \textbf{y}(x_0)=\textbf{y}_0. Siis funktsiooni väärtus järgmises punktis arvutatakse järgmise valemi järgi:

 \textbf{y}_{n+1} = \textbf{y}_n + {h \over 6} (\textbf{k}_1 + 2\textbf{k}_2 + 2\textbf{k}_3 + \textbf{k}_4)

kus

 \textbf{k}_1 = \textbf{f} \left( x_n, \textbf{y}_n \right),
 \textbf{k}_2 = \textbf{f} \left( x_n + {h \over 2}, \textbf{y}_n + {h \over 2} \textbf{k}_1 \right),
 \textbf{k}_3 = \textbf{f} \left( x_n + {h \over 2}, \textbf{y}_n + {h \over 2} \textbf{k}_2 \right),
 \textbf{k}_4 = \textbf{f} \left( x_n + h, \textbf{y}_n + h\textbf{k}_3 \right),
 \textbf{ }h — võrgu sammu suurus \textbf{ }x järgi.

See meetod on 4. järku, st viga igal sammul on O(h^5) ja summaarne viga integreerimise lõppintervallil on O(h^4).

Otsesed Runge-Kutta meetodid[muuda | redigeeri lähteteksti]

Otseste Runge-Kutta meetodite pere on 4. järku Runge-Kutta meetodi üldistus. See on antud valemitega

 \textbf{y}_{n+1} = \textbf{y}_n + h\sum_{i=1}^s b_i \textbf{k}_i,

kus

 \textbf{k}_1 = f(x_n, \textbf{y}_n), \,
 \textbf{k}_2 = f(x_n+c_2h, \textbf{y}_n+a_{21}h\textbf{k}_1), \,
 \textbf{k}_3 = f(x_n+c_3h, \textbf{y}_n+a_{31}h\textbf{k}_1+a_{32}h\textbf{k}_2), \,
 \vdots
 \textbf{k}_s = f(x_n+c_sh, \textbf{y}_n+a_{s1}h\textbf{k}_1+a_{s2}h\textbf{k}_2+\cdots+a_{s,s-1}h\textbf{k}_{s-1}).

Konkreetse meetodi määravad arv s ning koefitsiendid b_i, a_{ij} ja c_i. Need koefitsiendid paigutatakse sageli tabelisse

0
 c_2  a_{21}
 c_3  a_{31}  a_{32}
 \vdots  \vdots  \ddots
 c_s  a_{s1}  a_{s2}  \cdots  a_{s,s-1}
 b_1  b_2  \cdots  b_{s-1}  b_s

Runge-Kutta meetodi koefitsiendid peavad rahuldama tingimusi \sum_{j=1}^{i-1} a_{ij} = c_i ( i=2, \ldots, s). Kui me tahame, et meetod oleks p-järku, siis tuleb tagada ka tingimus \bar\textbf{y}(h+x_0)-\textbf{y}(h+x_0)=O(h^{p+1}), kus \bar\textbf{y}(h+x_0) on Runge-Kutta meetodil saadud lähendus. Pärast mitmekordset diferentseerimist muutub see tingimus polünomiaalvõrrandite süsteemiks, mille lahendid on meetodi koefitsiendid.

Vaata ka[muuda | redigeeri lähteteksti]

Kirjandus[muuda | redigeeri lähteteksti]