Articulo de referencia

Movimiento browniano geométrico

Para la simulación que genera las realizaciones, ver a continuación. Un movimiento browniano geométrico (GBM) (también conocido como movimiento browniano exponencial ) es un pro...

Para la simulación que genera las realizaciones, ver a continuación.

Un movimiento browniano geométrico (GBM) (también conocido como movimiento browniano exponencial ) es un proceso estocástico de tiempo continuo en el que el logaritmo de la cantidad que varía aleatoriamente sigue un movimiento browniano (también llamado proceso de Wiener ) con deriva . [1] Es un ejemplo importante de procesos estocásticos que satisfacen una ecuación diferencial estocástica (SDE); en particular, se utiliza en finanzas matemáticas para modelar los precios de las acciones en el modelo de Black-Scholes .

Definición técnica: la SDE

Se dice que un proceso estocástico S t sigue un GBM si satisface la siguiente ecuación diferencial estocástica (SDE):

d S a = micras S a d a + σ S a d Yo a {\displaystyle dS_{t}=\mu S_{t}\,dt+\sigma S_{t}\,dW_{t}}

donde es un proceso de Wiener o movimiento browniano , y ('la deriva porcentual') y ('la volatilidad porcentual') son constantes. Yo a Estilo de visualización Wt micras {\estilo de visualización \mu} σ {\estilo de visualización \sigma}

El primer parámetro se utiliza para modelar tendencias deterministas, mientras que el último parámetro modela eventos impredecibles que ocurren durante el movimiento.

Solución de la SDE

Para un valor inicial arbitrario S 0 la SDE anterior tiene la solución analítica (según la interpretación de Itô ):

S a = S 0 exp ( ( micras σ 2 2 ) a + σ Yo a ) . {\displaystyle S_{t}=S_{0}\exp \left(\left(\mu -{\frac {\sigma ^{2}}{2}}\right)t+\sigma W_{t}\right).}

La derivación requiere el uso del cálculo de Itô . La aplicación de la fórmula de Itô conduce a

d ( En S a ) = ( En S a ) " d S a + 1 2 ( En S a ) " d S a d S a = d S a S a 1 2 1 S a 2 d S a d S a {\displaystyle d(\ln S_{t})=(\ln S_{t})'dS_{t}+{\frac {1}{2}}(\ln S_{t})''\,dS_{t}\,dS_{t}={\frac {dS_{t}}{S_{t}}}-{\frac {1}{2}}\,{\frac {1}{S_{t}^{2}}}\,dS_{t}\,dS_{t}}

donde es la variación cuadrática de la SDE. d S a d S a {\displaystyle dS_{t}\,dS_{t}}

d S a d S a = σ 2 S a 2 d Yo a 2 + 2 σ S a 2 micras d Yo a d a + micras 2 S a 2 d a 2 {\displaystyle dS_{t}\,dS_{t}\,=\,\sigma ^{2}\,S_{t}^{2}\,dW_{t}^{2}+2\sigma S_{t}^{2}\mu \,dW_{t}\,dt+\mu ^{2}S_{t}^{2}\,dt^{2}}

Cuando , converge a 0 más rápido que , ya que . Por lo tanto, el infinitesimal anterior se puede simplificar mediante d a 0 {\displaystyle dt\to 0} d a {\estilo de visualización dt} d Yo a Estilo de visualización: dW_{t} d Yo a 2 = Oh ( d a ) {\displaystyle dW_{t}^{2}=O(dt)}

d S a d S a = σ 2 S a 2 d a {\displaystyle dS_{t}\,dS_{t}\,=\,\sigma ^{2}\,S_{t}^{2}\,dt}

Introduciendo el valor de en la ecuación anterior y simplificando obtenemos d S a Estilo de visualización dS_{t}

En S a S 0 = ( micras σ 2 2 ) a + σ Yo a . {\displaystyle \ln {\frac {S_{t}}{S_{0}}}=\left(\mu -{\frac {\sigma ^{2}}{2}}\,\right)t+\sigma W_{t}\,.}

Tomando la exponencial y multiplicando ambos lados por obtenemos la solución indicada anteriormente. S 0 Estilo de visualización S_{0}

Movimiento browniano aritmético

El proceso para satisfacer la SDE incógnita a = En S a S 0 {\displaystyle X_{t}=\ln {\frac {S_{t}}{S_{0}}}}

d incógnita a = ( micras σ 2 2 ) d a + σ d Yo a , {\displaystyle dX_{t}=\left(\mu -{\frac {\sigma ^{2}}{2}}\,\right)dt+\sigma dW_{t}\,,}

o más generalmente el proceso de resolución de la SDE

d incógnita a = metro d a + en d Yo a , {\displaystyle dX_{t}=m\,dt+v\,dW_{t}\,,}

donde y son constantes reales y para una condición inicial , se denomina Movimiento Browniano Aritmético (MBA). Este fue el modelo postulado por Louis Bachelier en 1900 para los precios de las acciones, en el primer intento publicado de modelar el movimiento browniano, conocido hoy como modelo de Bachelier . Como se mostró anteriormente, la SDE del MBA se puede obtener a través del logaritmo de un MBG mediante la fórmula de Itô. De manera similar, un MBA se puede obtener mediante la exponenciación de un MBA mediante la fórmula de Itô. metro {\estilo de visualización m} en > 0 {\displaystyle v>0} incógnita 0 Estilo de visualización X_{0}

Propiedades del GBM

La solución anterior (para cualquier valor de t) es una variable aleatoria distribuida log-normalmente con valor esperado y varianza dados por [2] S a Estilo de visualización S_{t}

mi ( S a ) = S 0 mi micras a , {\displaystyle \operatorname {E} (S_{t})=S_{0}e^{\mu t},}
Variedad ( S a ) = S 0 2 mi 2 micras a ( mi σ 2 a 1 ) . {\displaystyle \operatorname {Var} (S_{t})=S_{0}^{2}e^{2\mu t}\left(e^{\sigma ^{2}t}-1\right).}

Se pueden derivar utilizando el hecho de que es una martingala , y que O a = exp ( σ Yo a 1 2 σ 2 a ) {\displaystyle Z_{t}=\exp \left(\sigma W_{t}-{\frac {1}{2}}\sigma ^{2}t\right)}

mi [ exp ( 2 σ Yo a σ 2 a ) F s ] = mi σ 2 ( a s ) exp ( 2 σ Yo s σ 2 s ) , 0 s < a . {\displaystyle \operatorname {E} \left[\exp \left(2\sigma W_{t}-\sigma ^{2}t\right)\mid {\mathcal {F}}_{s}\right]=e^{\sigma ^{2}(ts)}\exp \left(2\sigma W_{s}-\sigma ^{2}s\right),\quad \forall 0\leq s<t.}

La función de densidad de probabilidad de es: S a Estilo de visualización S_{t}

F S a ( s ; micras , σ , a ) = 1 2 π 1 s σ a exp ( ( En s En S 0 ( micras 1 2 σ 2 ) a ) 2 2 σ 2 a ) . {\displaystyle f_{S_{t}}(s;\mu ,\sigma ,t)={\frac {1}{\sqrt {2\pi }}}\,{\frac {1}{s\sigma {\sqrt {t}}}}\,\exp \left(-{\frac {\left(\ln s-\ln S_{0}-\left(\mu -{\frac {1}{2}}\sigma ^{2}\right)t\right)^{2}}{2\sigma ^{2}t}}\right).}

Al derivar otras propiedades de GBM, se puede hacer uso de la SDE de la cual GBM es la solución, o se puede usar la solución explícita dada anteriormente. Por ejemplo, considere el proceso estocástico log( S t ). Este es un proceso interesante, porque en el modelo de Black-Scholes está relacionado con el retorno logarítmico del precio de las acciones. Usando el lema de Itô con f ( S ) = log( S ) se obtiene

d log ( S ) = f ( S ) d S + 1 2 f ( S ) S 2 σ 2 d t = 1 S ( σ S d W t + μ S d t ) 1 2 σ 2 d t = σ d W t + ( μ σ 2 / 2 ) d t . {\displaystyle {\begin{alignedat}{2}d\log(S)&=f'(S)\,dS+{\frac {1}{2}}f''(S)S^{2}\sigma ^{2}\,dt\\[6pt]&={\frac {1}{S}}\left(\sigma S\,dW_{t}+\mu S\,dt\right)-{\frac {1}{2}}\sigma ^{2}\,dt\\[6pt]&=\sigma \,dW_{t}+(\mu -\sigma ^{2}/2)\,dt.\end{alignedat}}}

De lo cual se deduce que . E log ( S t ) = log ( S 0 ) + ( μ σ 2 / 2 ) t {\displaystyle \operatorname {E} \log(S_{t})=\log(S_{0})+(\mu -\sigma ^{2}/2)t}

Este resultado también se puede obtener aplicando el logaritmo a la solución explícita de GBM:

log ( S t ) = log ( S 0 exp ( ( μ σ 2 2 ) t + σ W t ) ) = log ( S 0 ) + ( μ σ 2 2 ) t + σ W t . {\displaystyle {\begin{alignedat}{2}\log(S_{t})&=\log \left(S_{0}\exp \left(\left(\mu -{\frac {\sigma ^{2}}{2}}\right)t+\sigma W_{t}\right)\right)\\[6pt]&=\log(S_{0})+\left(\mu -{\frac {\sigma ^{2}}{2}}\right)t+\sigma W_{t}.\end{alignedat}}}

Tomando la expectativa obtenemos el mismo resultado que el anterior: . E log ( S t ) = log ( S 0 ) + ( μ σ 2 / 2 ) t {\displaystyle \operatorname {E} \log(S_{t})=\log(S_{0})+(\mu -\sigma ^{2}/2)t}

Simulación de trayectorias de muestra

# Código Python para la trama

importar  numpy  como  np
importar  matplotlib.pyplot  como  plt

mu  =  1
n  =  50
dt  =  0,1
x0  =  100
np . aleatorio . semilla ( 1 )

sigma  =  np . rango ( 0.8 ,  2 ,  0.2 )

x  =  np . exp (
    ( mu  -  sigma  **  2  /  2 )  *  dt
    +  sigma  *  np . aleatorio . normal ( 0 ,  np . sqrt ( dt ),  tamaño = ( len ( sigma ),  n )) . T
)
x  =  np . vstack ([ np . ones ( len ( sigma )),  x ])
x  =  x0  *  x . cumprod ( eje = 0 )

plt . parcela ( x )
plt . leyenda ( np . ronda ( sigma ,  2 ))
plt .xlabel ( "$t$ " )
plt . ylabel ( "$x$" )
plt . título (
    "Realizaciones del movimiento browniano geométrico con diferentes varianzas \n $\mu=1$"
)
plt . mostrar ()

Versión multivariada

El GBM se puede extender al caso en que existen múltiples trayectorias de precios correlacionadas. [3]

Cada trayectoria de precios sigue el proceso subyacente

d S t i = μ i S t i d t + σ i S t i d W t i , {\displaystyle dS_{t}^{i}=\mu _{i}S_{t}^{i}\,dt+\sigma _{i}S_{t}^{i}\,dW_{t}^{i},}

donde los procesos de Wiener están correlacionados de tal manera que donde . E ( d W t i d W t j ) = ρ i , j d t {\displaystyle \operatorname {E} (dW_{t}^{i}\,dW_{t}^{j})=\rho _{i,j}\,dt} ρ i , i = 1 {\displaystyle \rho _{i,i}=1}

Para el caso multivariado, esto implica que

Cov ( S t i , S t j ) = S 0 i S 0 j e ( μ i + μ j ) t ( e ρ i , j σ i σ j t 1 ) . {\displaystyle \operatorname {Cov} (S_{t}^{i},S_{t}^{j})=S_{0}^{i}S_{0}^{j}e^{(\mu _{i}+\mu _{j})t}\left(e^{\rho _{i,j}\sigma _{i}\sigma _{j}t}-1\right).}

Una formulación multivariada que mantiene independientes los movimientos brownianos impulsores es W t i {\displaystyle W_{t}^{i}}

d S t i = μ i S t i d t + j = 1 d σ i , j S t i d W t j , {\displaystyle dS_{t}^{i}=\mu _{i}S_{t}^{i}\,dt+\sum _{j=1}^{d}\sigma _{i,j}S_{t}^{i}\,dW_{t}^{j},}

donde la correlación entre y ahora se expresa a través de los términos. S t i {\displaystyle S_{t}^{i}} S t j {\displaystyle S_{t}^{j}} σ i , j = ρ i , j σ i σ j {\displaystyle \sigma _{i,j}=\rho _{i,j}\,\sigma _{i}\,\sigma _{j}}

Uso en finanzas

El movimiento browniano geométrico se utiliza para modelar los precios de las acciones en el modelo Black-Scholes y es el modelo más utilizado para el comportamiento de los precios de las acciones. [4]

Algunos de los argumentos para utilizar GBM para modelar los precios de las acciones son:

  • Los rendimientos esperados de GBM son independientes del valor del proceso (precio de la acción), lo que coincide con lo que esperaríamos en la realidad. [4]
  • Un proceso GBM sólo asume valores positivos, como los precios reales de las acciones.
  • Un proceso GBM muestra el mismo tipo de “rugosidad” en sus trayectorias que vemos en los precios de las acciones reales.
  • Los cálculos con procesos GBM son relativamente fáciles.

Sin embargo, el GBM no es un modelo completamente realista; en particular, se aleja de la realidad en los siguientes puntos:

  • En los precios de las acciones reales, la volatilidad cambia con el tiempo (posiblemente de forma estocástica ), pero en GBM se supone que la volatilidad es constante.
  • En la vida real, los precios de las acciones a menudo muestran saltos causados ​​por eventos o noticias impredecibles, pero en GBM, la trayectoria es continua (sin discontinuidad).

Además de modelar los precios de las acciones, el movimiento browniano geométrico también ha encontrado aplicaciones en el seguimiento de estrategias comerciales. [5]

Extensiones

En un intento de hacer que el GBM sea más realista como modelo para los precios de las acciones, también en relación con el problema de la sonrisa de volatilidad , se puede abandonar el supuesto de que la volatilidad ( ) es constante. Si suponemos que la volatilidad es una función determinista del precio de la acción y el tiempo, esto se llama un modelo de volatilidad local . Una extensión directa del GBM de Black Scholes es una SDE de volatilidad local cuya distribución es una mezcla de distribuciones de GBM, la dinámica de mezcla lognormal, que resulta en una combinación convexa de precios de Black Scholes para opciones. [3] [6] [7] [8] Si en cambio suponemos que la volatilidad tiene una aleatoriedad propia, a menudo descrita por una ecuación diferente impulsada por un movimiento browniano diferente, el modelo se llama modelo de volatilidad estocástica , véase por ejemplo el modelo de Heston . [9] σ {\displaystyle \sigma }

Véase también

Referencias

  1. ^ Ross, Sheldon M. (2014). "Variaciones sobre el movimiento browniano". Introducción a los modelos de probabilidad (11.ª ed.). Ámsterdam: Elsevier. pp. 612–14. ISBN 978-0-12-407948-9.
  2. ^ Øksendal, Bernt K. (2002), Ecuaciones diferenciales estocásticas: una introducción con aplicaciones , Springer, pág. 326, ISBN 3-540-63720-6
  3. ^ ab Musiela, M. y Rutkowski, M. (2004), Métodos de martingala en modelos financieros, segunda edición, Springer Verlag, Berlín.
  4. ^ ab Hull, John (2009). "12.3". Opciones, futuros y otros derivados (7.ª ed.).
  5. ^ Rej, A.; Seager, P.; Bouchaud, J.-P. (enero de 2018). "Estás en una fase de reducción. ¿Cuándo deberías empezar a preocuparte?". Wilmott . 2018 (93): 56–59. arXiv : 1707.01457 . doi :10.1002/wilm.10646. S2CID  157827746.
  6. ^ Fengler, MR (2005), Modelado semiparamétrico de la volatilidad implícita, Springer Verlag, Berlín. DOI: https://doi.org/10.1007/3-540-30591-2
  7. ^ Brigo, Damiano ; Mercurio, Fabio (2002). "Dinámica de mezclas lognormales y calibración para sonrisas de volatilidad del mercado". Revista Internacional de Finanzas Teóricas y Aplicadas . 5 (4): 427–446. doi :10.1142/S0219024902001511.
  8. ^ Brigo, D, Mercurio, F, Sartorelli, G. (2003). Dinámica de precios de activos alternativos y volatilidad, QUANT FINANC, 2003, Vol: 3, Páginas: 173 - 183, ISSN  1469-7688
  9. ^ Heston, Steven L. (1993). "Una solución de forma cerrada para opciones con volatilidad estocástica con aplicaciones a opciones sobre bonos y divisas". Review of Financial Studies . 6 (2): 327–343. doi :10.1093/rfs/6.2.327. JSTOR  2962057. S2CID  16091300.
  • Modelos de movimiento browniano geométrico para el movimiento de acciones excepto en eventos raros.
  • Simulación en Excel de un movimiento browniano geométrico para simular precios de acciones
  • "Aplicación web interactiva: Procesos estocásticos utilizados en finanzas cuantitativas". Archivado desde el original el 20 de septiembre de 2015. Consultado el 3 de julio de 2015 .
  • Sitio web sobre cálculo no newtoniano
  • Monitoreo de la estrategia comercial: modelado del PnL como un movimiento browniano geométrico
Retrieved from "https://en.wikipedia.org/w/index.php?title=Geometric_Brownian_motion&oldid=1210860085"