Seguir una tasa de cambio
Una ecuación diferencial dice con qué rapidez cambia una cantidad: dtdy=f(t,y). El dinero que crece a una tasa continua r cumple dtdy=ry, con solución y0ert. Una tasa corta atraída hacia un nivel de largo plazo θ cumple dtdr=κ(θ−r). Pocas ecuaciones tienen una solución en fórmula, así que avanzamos en el tiempo numéricamente.
yn+1=yn+hf(tn,yn)(1) El método de Euler sigue la tangente durante un paso corto h: valor nuevo = valor anterior + paso × pendiente. Es simple pero solo de primer orden: reducir h a la mitad reduce el error a la mitad. Runge–Kutta de orden 4 (RK4) muestrea la pendiente cuatro veces por paso y las combina, yn+1=yn+6h(k1+2k2+2k3+k4), y su error disminuye como h4: reducir h a la mitad lo divide por 16.
Crecimiento a tasa 1, de dos manerasResolvamos
y′=y,
y(0)=1 hasta
t=1 con
h=0,5. Euler:
y1=1+0,5⋅1=1,5,
y2=1,5+0,5⋅1,5=2,25, frente a
e=2,7183, un error de 0,4683. RK4 para un paso hasta
t=0,5:
k1=1,
k2=1,25,
k3=1,3125,
k4=1,65625, así que
y=1+60,5(1+2,5+2,625+1,65625)=1,6484, frente al verdadero
e0,5=1,6487.
- Euler: yn+1=yn+hf(tn,yn), precisión de primer orden.
- RK4: cuatro pendientes por paso combinadas con pesos 1, 2, 2, 1; precisión de cuarto orden.
- Error de truncamiento: lo que se pierde al avanzar por pasos en lugar de resolver exactamente.
Un paso menor no es gratisReducir el paso a la mitad duplica el número de pasos. Más allá de cierto punto, el error de redondeo crece mientras el de truncamiento ya no importa, y algunas ecuaciones rígidas divergen si el paso no es lo bastante pequeño. Compare dos tamaños de paso para estimar el error antes de fiarse de un resultado.
Lab: qué se califica
- Dos pasos de Euler15 pts
- El error de Euler15 pts
- Un paso de RK420 pts
- Órdenes de precisión10 pts
- Reducir el paso a la mitad15 pts
- El bucle de Euler10 pts
- Por qué RK4 supera a Euler15 pts