FrogapenaEDAaurreratua

Dedukzioa: 4. ordenako Runge-Kutta

Simpson-en erregela PVIaren forma integralari aplikatzeak eta malda ezezagunak kateatutako barne-ebaluazioekin hurbiltzeak RK4 klasikoa sortzen du eta bere 1, 2, 2, 1 pisuak azaltzen ditu.

Simpson forma integralaren gainean

RK4-ren lau etapek Simpson-en malda zehatzak kateatutako barne-ebaluazioekin ordezkatzen dituzte.

Lehenik ekuazio diferentziala $[t_k,t_{k+1}]$ azpitartean integratzen dugu:
$\int_{t_k}^{t_{k+1}} y'(\tau)\,d\tau=\int_{t_k}^{t_{k+1}} f\bigl(\tau,y(\tau)\bigr)\,d\tau$
Ezkerreko aldeari Kalkuluaren Oinarrizko Teorema aplikatuz, PVIaren forma integral zehatza geratzen da:
$y(t_{k+1})=y(t_k)+\int_{t_k}^{t_{k+1}} f\bigl(\tau,y(\tau)\bigr)\,d\tau$
Integrala Simpson-en erregelarekin hurbiltzen dugu, muturrak eta erdiko puntua $t_{k+\frac12}=t_k+\frac{h}{2}$ erabiltzen dituena:
$\int_{t_k}^{t_{k+1}} f\,d\tau\approx\frac{h}{6}\Bigl(f(t_k,y_k)+4f\bigl(t_{k+\frac12},y_{k+\frac12}\bigr)+f(t_{k+1},y_{k+1})\Bigr)$
Arazoa: ez dakigu ez $y_{k+\frac12}$ ez $y_{k+1}$ . Runge-Kuttaren irtenbidea kateatutako barne-maldekin estimatzea da, bakoitza aurrekoarekin aurrera eginez eraikia:
$\begin{aligned} k_1&=f(t_k,\,y_k)\\ k_2&=f\bigl(t_k+\tfrac{h}{2},\,y_k+\tfrac{h}{2}k_1\bigr)\\ k_3&=f\bigl(t_k+\tfrac{h}{2},\,y_k+\tfrac{h}{2}k_2\bigr)\\ k_4&=f(t_k+h,\,y_k+h\,k_3) \end{aligned}$
Simpson-en erdiko malda erdiko puntuko bi estimazioen batez bestekoaz hurbiltzen da, eta amaierakoa $k_4$ -rekin:
$f\bigl(t_{k+\frac12},y_{k+\frac12}\bigr)\approx\frac{k_2+k_3}{2},\qquad f(t_{k+1},y_{k+1})\approx k_4$
Simpson-en ordezkatuz, erdiko puntuaren $4$ pisua honela banatzen da: $4\cdot\frac{k_2+k_3}{2}=2k_2+2k_3$ , eta RK4 klasikoa agertzen da bere $1,2,2,1$ pisuekin:
$y_{k+1}=y_k+\frac{h}{6}\bigl(k_1+2k_2+2k_3+k_4\bigr)$
Heun-enaren antzeko Taylor-en analisi batek (baina laugarren ordenaraino) baieztatzen du barne-hurbilketa horiek $\mathcal{O}(h^5)$ errore lokala mantentzen dutela, eta beraz metodoa 4. ordenakoa dela.