70 votos

¿Cuáles son algunas aproximaciones útiles a la fórmula de Black-Scholes?

Definamos la fórmula de Black-Scholes como la función $f(S, X, T, r, v)$ .

Tengo curiosidad por saber qué funciones son computacionalmente más sencillas que el Black-Scholes y que dan resultados que se aproximan a $f$ para un conjunto determinado de entradas $S, X, T, r, v$ .

Entiendo que "computacionalmente más simple" no está bien definido. Pero me refiero a más simple en términos de número de términos utilizados en la función. O incluso más específicamente, el número de pasos computacionales distintos que hay que completar para llegar al resultado de Black-Scholes.

Obviamente, Black-Scholes es computacionalmente sencillo, pero estoy dispuesto a cambiar algo de precisión por una función aún más sencilla que dé resultados que se aproximen a B&S.

¿Existe alguna aproximación más sencilla?

61voto

Tim Boland Puntos 4063

Esto es sólo para ampliar un poco en Respuesta de vonjd .

La fórmula aproximada mencionada por vonjd se debe a Brenner y Subrahmanyam ("A simple solution to compute the Implied Standard Deviation", Diario de los Analistas Financieros (1988), pp. 80-83). No dispongo de un enlace gratuito al artículo, así que me limitaré a hacer una derivación rápida y sucia.

Para la opción de compra at-the-money, tenemos $S=Ke^{-r(T-t)}$ . Introduciendo esto en la fórmula estándar de Black-Scholes $$C(S,t)=N(d_1)S-N(d_2)Ke^{-r(T-t)},$$ conseguimos que $$C(S,t)=\left[N\left(\frac{1}{2}\sigma\sqrt{T-t}\right)-N\left(-\frac{1}{2}\sigma\sqrt{T-t}\right)\right]S.\qquad\qquad(1)$$ Ahora, la fórmula de Taylor implica para pequeños $x$ que $$N(x)=N(0)+N'(0)x+N''(0)\frac{x^2}{2}+O(x^3).\qquad\qquad\qquad\qquad(2)$$ Combinando (1) y (2), obtendremos con algunas cancelaciones obvias que $$C(S,t)=S\left(N'(0)\sigma\sqrt{T-t}+O(\sigma^3\sqrt{(T-t)^3})\right).$$ Pero $$N'(0)=\frac{1}{\sqrt{2\pi}}=0.39894228...$$ por lo que finalmente tenemos, para los pequeños $\sigma\sqrt{T-t}$ que $$C(S,t)\approx 0.4S\sigma\sqrt{T-t}.$$ La fórmula modificada $$C(S,t)\approx 0.4Se^{-r(T-t)}\sigma\sqrt{T-t}$$

da una aproximación ligeramente mejor.

5 votos

Buena respuesta. Me pregunto si hay alguna aproximación para las opciones que son no ¿en el dinero? Mi enfoque simple sería asumir una delta del 50% para la opción de compra ATM e implicar un precio para la opción no ATM como Opción = Precio ATM + 0,5*(Strike - Forward). ¿Alguien tiene algo mejor?

0 votos

@Robert Estoy de acuerdo con la idea pero parece que hay un error de señal. Debería ser Option = ATM Price + 0.5*(Forward - Strike) ? Entonces, la opción out-of-the-money ( Strike>Forward ) debería ser menos valiosa que la opción del cajero automático.

46voto

penti Puntos 93

Esta es la mejor aproximación que he visto:

Si odias los ordenadores y los lenguajes informáticos no te rindas, todavía hay ¡esperanza! ¿Qué tal si tomas el Black-Scholes en tu cabeza? Si la opción es sobre at-the-money-forward y es un corto plazo hasta el vencimiento, entonces puede utilizar la siguiente aproximación:

call = put = Precio de la acción * 0,4 * volatilidad * Sqrt( Tiempo )

Fuente: http://www.espenhaug.com/black_scholes.html

0 votos

¿Por qué exigimos poco tiempo de maduración? ¿La aproximación falla cuando el tiempo hasta el vencimiento es largo? Además, ¿por qué la opción de compra y la de venta tienen el mismo valor?

0 votos

@Idonknow La aproximación se basa en la aproximación de la serie de Taylor de primer orden de una función de d1. Cuando d1 es "grande", los términos de segundo orden y superiores de la serie de Taylor se vuelven no despreciables, lo que hace que la aproximación de primer orden sea "mala". Dado que d1 para una opción a plazo es 0,5*sigma*sqrt(tiempo), un tiempo grande hace que d1 sea grande, lo que hace que la aproximación sea mala. Una volatilidad grande tendrá el mismo efecto.

30voto

Erlend Halvorsen Puntos 843

La fórmula Black-Scholes "normal-vol" conduce rápidamente a una aproximación similar a la descrita por olaker.

Haga clic en aquí para un artículo que contiene una derivación formal de los precios de compra y venta basada en un modelo normal (es decir, un movimiento browniano en lugar de un movimiento browniano geométrico).

La fórmula del precio de compra es:

$$\text{Call} = (F-K)N(d_1) + \frac{\sigma \sqrt{T-t}}{\sqrt{2 \pi}} e^{-\frac{1}{2} d_1^2},$$

donde

$$d_1 = \frac{F-K}{\sigma\sqrt{T-t}}.$$

Por cierto, trabajo en renta fija, por lo que siempre tiendo a escribir una versión que se adapte a las swaptions.

Opciones que son ATM

Puede ver que para $F=K$ esto se convierte en $\text{ATMCall} = \frac{\sigma \sqrt{T-t}}{\sqrt{2 \pi}} \approx 0.4 \, \sigma \, \sqrt{T-t}$ .

Opciones que no son ATM

Recientemente he descubierto una generalización de esta fórmula que funciona muy bien también para las huelgas que no están en el dinero. Ver mi blog para una discusión más larga, aquí están los puntos principales.

La descomposición estándar de una opción es:

$$\text{Option value} = \text{Intrinsic value} + \text{Time value}.$$

En la fórmula normal de Black-Scholes anterior, si se investiga el término $(F-K)N(d_1)$ en una hoja de cálculo, verás que para pequeños niveles de volatilidad y madurez (prueba, por ejemplo, $\sigma=0.0025$ , Madurez=1) se aproxima bastante a $\max(0,F-K)$ - que es el valor intrínseco de la llamada.

En consecuencia, la fórmula normal de la BS es casi :

$$\text{Call Price} = \text{Intrinsic Value} + \text{ATMPrice} \times e^{-\frac{1}{2} d_1^2},$$

donde

$$\text{ATMPrice} = \frac{\sigma \sqrt{T-t}}{\sqrt{2 \pi}}.$$

Sin embargo, si se compara esta aproximación con la verdadera fórmula de BS en una hoja de cálculo, se verá que alrededor de la huelga (especialmente para valores mayores de $\sigma$ ) da demasiado valor a la llamada: básicamente el término $\text{ATMPrice} \times e^{-\frac{1}{2}d_1^2}$ es demasiado grande cuando $d_1$ es distinto de cero y pequeño. Esto nos dice que la diferencia entre $(F-K)N(d_1)$ y $\max(0,F-K)$ se vuelve importante cerca de la huelga. Como digo, mira en una hoja de cálculo.

No obstante, esta fórmula sencilla pero errónea para el precio de la opción de compra nos indica la dirección correcta: muestra que el valor temporal de la opción debe escribirse en función del precio de la opción ATM .

Esta es mi solución.

Yo lo llamo La descomposición de Hardy :

$$\text{Call Price} = \text{Intrinsic Value} + \text{ATMPrice}\times\text{HardyFactor}$$ donde

$$\text{HardyFactor} = e^{-\frac{1}{2} d_1^2} + \frac{d_1}{0.4}N(d_1) – \max(\frac{d_1}{0.4},0).$$

Hasta aquí se trata de una reorganización de la fórmula original de Black-Scholes "normal-vol". El resultado clave es que el $\text{HardyFactor}$ se aproxima bien con una expresión sencilla:

$$\text{HardyFactor} \approx (1- 0.41 |d_1|) e^{-|d_1|}$$

Por lo tanto, se puede utilizar lo siguiente como una aproximación bastante buena para los precios de las llamadas:

$$\text{Call Price} = \text{Intrinsic Value} + \text{ATMPrice}\times (1- 0.41 |d_1|) e^{-|d_1|}.$$

Un resultado similar es válido para las opciones de venta.

Puede utilizar este Descomposición Hardy para calcular los precios de las opciones en su cabeza - sólo necesita recordar algunos valores:

Remember these few values

1 votos

No entiendo por qué, ATM, no hay el valor spot subyacente en su fórmula ?

1 votos

Porque estoy usando un normal vol, no un lognormal En la respuesta de vonjd su término de "volatilidad" es un vol lognormal, por lo que (S x sigma) es el correspondiente normal vol.

25voto

gary Puntos 4856

Además de lo que vonjd ya publicado le recomendaría que mirara el artículo de E.G. Haug - El genio de las opciones. Wilmott.com . Se pueden encontrar algunas aproximaciones de BS no sólo para las call y put europeas vainilla sino incluso para algunas exóticas. Por ejemplo:

  • opción de elección: call = put = $0.4F_{0} e^{-\mu T}\sigma(\sqrt{T}-\sqrt{t})$
  • opción asiática: call = put = $0.23F_{0} e^{-\mu T}\sigma(\sqrt{T}+2\sqrt{t})$
  • llamada de retorno a la huelga flotante = $0.8F_{0} e^{-\mu T}\sigma\sqrt{T} - 0.25{\sigma^2}T$
  • put con efecto de retroceso de la huelga flotante = $0.8F_{0} e^{-\mu T}\sigma\sqrt{T} + 0.25{\sigma^2}T$
  • spread a plazo: call = put = $0.4F_{1} e^{-\mu T}\sigma\sqrt{T}$

2 votos

El enlace wilmott no está disponible ahora. ¿Alguien conoce un nuevo enlace?

10voto

Muhammed Refaat Puntos 97

La distribución logística se aproxima a la función de distribución normal utilizada en el Black-Scholes. Los inconvenientes de la función de distribución acumulativa normal son que no puede calcularse exactamente a través de funciones elementales, no puede invertirse algebraicamente (es decir, la biyección inversa no puede resolverse algebraicamente) y es computacionalmente cara. Una alternativa es la distribución logística, que puede calcularse con funciones elementales, puede invertirse algebraicamente y es computacionalmente inexistente. Además, como resultado de la propiedad de inversión, los valores límite óptimos de las opciones americanas pueden estimarse fácilmente para la distribución logística.

La solución general de Black-Scholes viene dada por:

$$V_t = \Phi[d_1]S_t - \Phi[d_2]K e^{-r (T-t)} $$

donde $\Phi[X]$ es la FDA de la distribución normal.

Se puede demostrar que una aproximación numérica rápida (es decir, computacionalmente barata) es la siguiente:

$$V_t \approx \left(\frac{1}{2} + \frac{1}{2} \operatorname{Tanh} \left( \frac{\pi \, d_1}{2 \sqrt{3}} \right) \right) S_t - \left(\frac{1}{2} + \frac{1}{2} \operatorname{Tanh} \left( \frac{\pi \, d_2}{2 \sqrt{3}} \right) \right) K e^{-r(T-t)} = \left(1- \frac{1}{1+e^{\pi d_1 / \sqrt{3}}} \right)S_t - \left(1- \frac{1}{1+e^{\pi (d_1-\sigma \sqrt{T-t}) / \sqrt{3}}} \right) K e^{-r(T-t)}$$


La lógica

Una distribución logística $F$ -- que puede expresarse como una tangente hiperbólica reescalada, puede aproximarse a la función de distribución normal $\Phi$ . Asimismo, su función inversa -- la función "logit" $F^{-1}$ -- se puede reescalar para aproximar la CDF normal inversa -- la función "probit" $\Phi^{-1}$ .

En comparación, la distribución logística tiene colas más gruesas (lo que puede ser deseable). Mientras que la FCD y la FCD inversa de la distribución normal ("probit") no pueden expresarse mediante funciones elementales, las expresiones de forma cerrada para la FCD de la distribución logística y su inversa se derivan fácilmente y se comportan como funciones algebraicas elementales.

La distribución logística surge de la ecuación diferencial $\frac{d}{dx}f(x) = f(x)(1-f(x))$ . Intuitivamente, esta función se utiliza normalmente para modelar un proceso de crecimiento en el que la tasa se comporta como una curva de campana. En física, surge como la "Distribución límite de un movimiento aleatorio amortiguado de velocidad finita descrito por un proceso telegráfico en el que los tiempos aleatorios entre cambios de velocidad consecutivos tienen distribuciones exponenciales independientes con parámetros linealmente crecientes".

En comparación, la distribución normal surge de la siguiente ecuación diferencial: $ \frac{d \,f(x)}{dx}=f(x)\frac{(\mu-x)}{\sigma^2}$ ). La distribución normal se utiliza habitualmente para modelar procesos de difusión. Por ejemplo, un proceso de Wiener es un proceso estocástico que tiene incrementos independientes distribuidos normalmente con media $\mu$ y la varianza $\sigma^2$ . En el límite, se trata de un movimiento browniano.

Curiosamente, la distribución logística surge en un proceso físico análogo al movimiento browniano.

Obsérvese que la FCD de la distribución logística $F$ puede expresarse mediante la función tangente hiperbólica:

$F(x;\mu ,s)={\frac {1}{1+e^{{-{\frac {x-\mu }{s}}}}}}={\frac 12}+{\frac 12}\;\operatorname {Tanh}\!\left({\frac {x-\mu }{2s}}\right)$

Dado que la varianza de la distribución es ${\tfrac {s^{2}\pi ^{2}}{3}}$ la distribución logística puede escalarse para aproximarse a la distribución normal multiplicando su varianza $\frac{3}{\pi ^2}$ . La aproximación resultante tendrá los mismos primeros y segundos momentos que la distribución normal, pero tendrá una cola más gruesa (es decir, "platocurótica").

También, $\Phi$ está relacionada con la función de error (y su complemento) por: $\Phi (x)={\frac {1}{2}}+{\frac {1}{2}}\operatorname {erf} \left(x/{\sqrt {2}}\right)={\frac {1}{2}}\operatorname {erfc} \left(-x/{\sqrt {2}}\right)$


Así, para una distribución normal estándar con $\mu =0$ y $\sigma =1$ : $$\operatorname{erf}(\frac{x}{\sqrt{2}}) \approx \operatorname{Tanh}\left(\frac{x \, \pi}{2 \sqrt{3}} \right) \equiv \frac{e^{\frac{\pi\,x}{\sqrt{3}}}-1}{e^{\frac{\pi\,x}{\sqrt{3}}}+1} $$

$$\operatorname{erf}(x) \approx \operatorname{Tanh}\left(\frac{x \, \pi}{ \sqrt{6}} \right) \equiv \frac{e^{\pi\,x\frac{2}{\sqrt{3}}}-1}{e^{\pi\,x\frac{2}{\sqrt{3}}}+1} $$

$$\Phi \left( x \right) \approx \frac{1}{2} + \frac{1}{2} \operatorname{Tanh} \left( \frac{\pi \, x}{2 \sqrt{3}} \right) \equiv 1-\frac{1}{1+e^{\pi x \over\sqrt{3}}} $$

Y fácilmente, así: $$x \mapsto \Phi^{-1}\left(p\right) \approx -\frac{2\sqrt{3}\operatorname{ArcTanh}\left( 1-2p \right)}{\pi}$$


De nuevo, la principal ventaja de aproximar la normal con la distribución logística es que la FCD y la FCD inversa pueden expresarse fácilmente utilizando funciones elementales. He comprobado que esta propiedad es muy útil para estimar los valores de las opciones americanas o en el caso general de que no se conozcan las condiciones límite terminales/óptimas.

Sin embargo, hay que tener en cuenta que las distribuciones se comportan de forma diferente en "las colas", lo que podría dar lugar a residuos asintóticos si se espera que sean normales.

Comparison of the logit function with a scaled probit (i.e. the inverse CDF of the normal distribution), comparing {\displaystyle \operatorname {logit} (x)} \operatorname {logit} (x) vs. {\displaystyle \Phi ^{-1}(x)/{\sqrt {\frac {\pi }{8}}}} \Phi ^{-1}(x)/{\sqrt {\frac {\pi }{8}}}, which makes the slopes the same at the y-origin.

4 votos

-1. Parece que simplemente se aproxima basándose en el perfil de las dos funciones. No hay justificación numérica ni teórica. De hecho, como puedes ver en tu gráfico, la aproximación no es buena. De hecho, para grandes positivos $x$ , $1-\Phi(x)\sim e^{-\frac{x^2}2}$ en lugar de $e^{-ax}$ como su $1-\tanh(x)$ lo tendría. Tienes la asintótica mal. Hay aproximaciones rápidas y teóricamente justificadas. Busca la relación de Mills.

0 votos

He aquí una justificación teórica para modelar la dinámica de los precios de las acciones como si siguieran procesos telegráficos: link.springer.com/article/10.1023%2FA%3A1009437108439

1 votos

Un documento interesante. Sin embargo, está detrás de un muro de pago. Sin embargo, ¿cómo aborda el documento mi objeción, que es que su respuesta o el modelo vinculado no se basan en el modelo Black-Scholes, que es la preocupación de la pregunta del OP, en lugar de algún otro modelo? Se puede plantear un número infinito de otros modelos. Es irrelevante para una pregunta relativa a un modelo en particular, que en este caso es Black-Scholes.

Finanhelp.com

FinanHelp es una comunidad para personas con conocimientos de economía y finanzas, o quiere aprender. Puedes hacer tus propias preguntas o resolver las de los demás.

Powered by:

X