σμimplied vol · KKmax(S−K,0)udβρᵢⱼfat tailsPnL · drawdownyield · TdSₜ = μSₜ dt + σSₜ dWₜC = S·N(d₁) − K·e⁻ʳᵀ·N(d₂)σₚ² = wᵀΣwd₁ = [ln(S/K) + (r + σ²⁄2)T] ⁄ σ√T∂V/∂t + ½σ²S²·∂²V/∂S² + rS·∂V/∂S − rV = 0Sharpe = (Rₚ − R_f) ⁄ σₚVaR₉₉ = μ − 2.33σρ = Cov(X,Y) ⁄ σₓσᵧE[Rₚ] = Σ wᵢ·E[Rᵢ]φ(x) = e^(−x²/2) ⁄ √2πE[Rᵢ] = R_f + βᵢ(E[Rₘ] − R_f)dXₜ = θ(μ − Xₜ)dt + σ dWₜσ²ₜ = ω + α·ε²ₜ₋₁ + β·σ²ₜ₋₁C − P = S − K·e⁻ʳᵀΔ = ∂V/∂S Γ = ∂²V/∂S²p = (e^(rΔt) − d) ⁄ (u − d)rₜ = ln(Pₜ ⁄ Pₜ₋₁)E[Xₜ₊₁ | ℱₜ] = Xₜf* = (bp − q) ⁄ bz = (x − μ) ⁄ σP(N=k) = λᵏe^(−λ) ⁄ k!Θ = −∂V/∂Tν = S·φ(d₁)·√Tmin wᵀΣw s.t. wᵀμ = μ*∫Σ∂∞λ
uantCraft AcademyHojas de rutaRutas
Ctrl K
AO
+1,2 k estudiantes
Iniciar sesiónUnirse ↗
Hojas de rutaRutas
AO
+1,2 k estudiantes
Iniciar sesión
  • English
  • Français
  • Español
  • العربية
  • Português
  • Deutsch
  • 日本語
  • 中文
σμimplied vol · KKmax(S−K,0)udβρᵢⱼfat tailsPnL · drawdownyield · TdSₜ = μSₜ dt + σSₜ dWₜC = S·N(d₁) − K·e⁻ʳᵀ·N(d₂)σₚ² = wᵀΣwd₁ = [ln(S/K) + (r + σ²⁄2)T] ⁄ σ√T∂V/∂t + ½σ²S²·∂²V/∂S² + rS·∂V/∂S − rV = 0Sharpe = (Rₚ − R_f) ⁄ σₚVaR₉₉ = μ − 2.33σρ = Cov(X,Y) ⁄ σₓσᵧE[Rₚ] = Σ wᵢ·E[Rᵢ]φ(x) = e^(−x²/2) ⁄ √2πE[Rᵢ] = R_f + βᵢ(E[Rₘ] − R_f)dXₜ = θ(μ − Xₜ)dt + σ dWₜσ²ₜ = ω + α·ε²ₜ₋₁ + β·σ²ₜ₋₁C − P = S − K·e⁻ʳᵀΔ = ∂V/∂S Γ = ∂²V/∂S²p = (e^(rΔt) − d) ⁄ (u − d)rₜ = ln(Pₜ ⁄ Pₜ₋₁)E[Xₜ₊₁ | ℱₜ] = Xₜf* = (bp − q) ⁄ bz = (x − μ) ⁄ σP(N=k) = λᵏe^(−λ) ⁄ k!Θ = −∂V/∂Tν = S·φ(d₁)·√Tmin wᵀΣw s.t. wᵀμ = μ*∫Σ∂∞λ
uantCraft AcademyHojas de rutaRutas
Ctrl K
AO
+1,2 k estudiantes
Iniciar sesiónUnirse ↗
Hojas de rutaRutas
AO
+1,2 k estudiantes
Iniciar sesión
  • English
  • Français
  • Español
  • العربية
  • Português
  • Deutsch
  • 日本語
  • 中文
Ruta de Calculus›Módulo 3 · Aproximación y dinámica›Tema 3.2·~50 minIntermedio

Ecuaciones diferenciales: Euler y RK4

Posición en la rutaMódulo 3 de 4 · Tema 2 de 2
3.2 Ecuaciones diferenciales: Euler y RK4
—Puntuación del lab— / 100XP total—Racha—

Seguir una tasa de cambio

Una ecuación diferencial dice con qué rapidez cambia una cantidad: dydt=f(t,y)\frac{dy}{dt}=f(t,y)dtdy​=f(t,y). El dinero que crece a una tasa continua rrr cumple dydt=ry\frac{dy}{dt}=rydtdy​=ry, con solución y0erty_0e^{rt}y0​ert. Una tasa corta atraída hacia un nivel de largo plazo θ\thetaθ cumple drdt=κ(θ−r)\frac{dr}{dt}=\kappa(\theta-r)dtdr​=κ(θ−r). Pocas ecuaciones tienen una solución en fórmula, así que avanzamos en el tiempo numéricamente.

yn+1=yn+h f(tn,yn)y_{n+1}=y_n+h\,f(t_n,y_n)yn+1​=yn​+hf(tn​,yn​)(1)

El método de Euler sigue la tangente durante un paso corto hhh: valor nuevo = valor anterior + paso × pendiente. Es simple pero solo de primer orden: reducir hhh 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+h6(k1+2k2+2k3+k4)y_{n+1}=y_n+\frac{h}{6}(k_1+2k_2+2k_3+k_4)yn+1​=yn​+6h​(k1​+2k2​+2k3​+k4​), y su error disminuye como h4h^4h4: reducir hhh a la mitad lo divide por 16.

Crecimiento a tasa 1, de dos manerasResolvamos y′=yy'=yy′=y, y(0)=1y(0)=1y(0)=1 hasta t=1t=1t=1 con h=0,5h=0{,}5h=0,5. Euler: y1=1+0,5⋅1=1,5y_1=1+0{,}5\cdot1=1{,}5y1​=1+0,5⋅1=1,5, y2=1,5+0,5⋅1,5=2,25y_2=1{,}5+0{,}5\cdot1{,}5=2{,}25y2​=1,5+0,5⋅1,5=2,25, frente a e=2,7183e=2{,}7183e=2,7183, un error de 0,4683. RK4 para un paso hasta t=0,5t=0{,}5t=0,5: k1=1k_1=1k1​=1, k2=1,25k_2=1{,}25k2​=1,25, k3=1,3125k_3=1{,}3125k3​=1,3125, k4=1,65625k_4=1{,}65625k4​=1,65625, así que y=1+0,56(1+2,5+2,625+1,65625)=1,6484y=1+\frac{0{,}5}{6}(1+2{,}5+2{,}625+1{,}65625)=1{,}6484y=1+60,5​(1+2,5+2,625+1,65625)=1,6484, frente al verdadero e0,5=1,6487e^{0{,}5}=1{,}6487e0,5=1,6487.
  • Euler: yn+1=yn+h f(tn,yn)y_{n+1}=y_n+h\,f(t_n,y_n)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
0/ 7 respondidos
Ejercicios
1Respuesta numérica15 ptsSin responder

Considere y′=yy'=yy′=y con y(0)=1y(0)=1y(0)=1 y paso h=0,5h=0{,}5h=0,5. Dé dos pasos de Euler y escriba y2y_2y2​ (la estimación de y(1)y(1)y(1)).

y2=y_2 =y2​=aceptado: ± 0.001
2Respuesta numérica15 ptsSin responder

Considere y′=yy'=yy′=y con y(0)=1y(0)=1y(0)=1 y paso h=0,5h=0{,}5h=0,5. La solución exacta es y(1)=e=2,7183y(1)=e=2{,}7183y(1)=e=2,7183. ¿Cuánto se aparta la estimación de Euler y2=2,25y_2=2{,}25y2​=2,25?

error=\text{error} =error=aceptado: ± 0.001
3Respuesta numérica20 ptsSin responder

Considere y′=yy'=yy′=y con y(0)=1y(0)=1y(0)=1 y paso h=0,5h=0{,}5h=0,5. Dé un paso de RK4 desde t=0t=0t=0 y escriba y1y_1y1​, la estimación de y(0,5)y(0{,}5)y(0,5), con 4 decimales. Para f(t,y)=yf(t,y)=yf(t,y)=y: k1=ynk_1=y_nk1​=yn​, k2=yn+h2k1k_2=y_n+\frac{h}{2}k_1k2​=yn​+2h​k1​, k3=yn+h2k2k_3=y_n+\frac{h}{2}k_2k3​=yn​+2h​k2​, k4=yn+hk3k_4=y_n+hk_3k4​=yn​+hk3​.

y1=y_1 =y1​=aceptado: ± 0.0002
4Elección10 ptsSin responder

¿Cuáles son los órdenes de precisión del método de Euler y de RK4?

5Respuesta numérica15 ptsSin responder

El método de Euler tiene un error de 0,08 con paso h=0,1h=0{,}1h=0,1. ¿Qué error espera aproximadamente con h=0,05h=0{,}05h=0,05?

error=\text{error} =error=aceptado: ± 0.005
6Ordenar los pasos10 ptsSin responder

Ordene los pasos de una iteración de Euler.

Arrastra para reordenar, o enfoca un elemento, pulsa Espacio para levantarlo, las flechas para moverlo y Espacio para soltarlo.

  1. 1Actualizar y←y+h⋅pendientey\leftarrow y+h\cdot\text{pendiente}y←y+h⋅pendiente
  2. 2Partir de los valores iniciales t0t_0t0​ y y0y_0y0​
  3. 3Avanzar el tiempo t←t+ht\leftarrow t+ht←t+h y repetir
  4. 4Calcular la pendiente f(t,y)f(t,y)f(t,y) en el punto actual
7Explicar con palabras15 ptsSin responder

En dos o tres frases, explique por qué RK4 suele ser más preciso que Euler con el mismo paso y qué cuesta.

Modo matemáticasEntrada por teclado o LaTeXEjercicio 1 de 7 · Ctrl ⏎ para comprobar

Inicia sesión para ejecutar tu código.

Todas las clases se pueden leer gratis. Crea una cuenta o inicia sesión para ejecutar código, lanzar las pruebas y enviar este lab.

  • IDE en el navegador con terminal bash — sin instalar nada
  • Corrección instantánea, XP y tu lugar en la ruta
Crear una cuentaIniciar sesión
O
Continuar con GitHub