3 votos

¿Cómo derivas esta ecuación tipo Carr-Madan?

¿Cómo se deriva la ecuación (3) a continuación? La ecuación está etiquetada como ecuación (11) en este documento: http://janroman.dhis.org/finance/IR/Heston%E2%80%93Hull%E2%80%93White%20Model%20Part%20I.pdf

Hay partes de este documento que no entiendo. Sospecho que hay algunos pequeños errores, que hacen que las cosas sean más difíciles de desentrañar. Así es como estoy leyendo la ecuación (3).

  1. $C(t)$ es el precio de una opción de compra en el tiempo $t$ con precio de ejercicio $K$ y tiempo de vencimiento $T$.
  2. $\varphi(z)$ es la función característica de $x=\log(S(T))$, donde $S$ es el precio del activo en el tiempo $T$. Es decir,

$$\varphi(z) = \mathbb{E}[e^{xt}|F(t)]\tag{1}$$

donde $F$ es la filtración. El autor no lo indica, pero también creo que este valor esperado se calcula bajo la medida neutral al riesgo con respecto a $F(t)$. Si $q_T$ es la función de densidad de $x$ en el tiempo $T$ bajo la medida neutral al riesgo, entonces

$$\varphi(z) = \int_{-\infty}^{\infty}e^{ixz}q_T(x)\mathrm{d}x\text{.}\tag{2}$$

  1. $f(x)$ es el pago de la opción de compra cuando el precio del activo en el tiempo $T$ es $e^{x}$. $e^{-cx}f(x)$ es una transformación de $f$ necesaria para hacer la transformada de Fourier. $\hat{f}$ es la transformada de Fourier de $e^{-cx}f(x)$.

Finalmente, nuestra ecuación es la integral compleja de línea

$$C(t) = \frac{e^{-r(T-t)}}{2\pi}\int_{c-\infty}^{c+\infty} \varphi(-z) \hat{f}(z)\mathrm{d}z\tag{3}\text{.}$$

Esta ecuación parece ser similar a la ecuación (5) del documento de Carr & Madan aquí. Esa ecuación es simplemente la transformada de Fourier inversa

$$C(t) = \frac{e^{-ck}}{2\pi}\int_{-\infty}^{\infty}e^{ivk}\psi_T(v)\mathrm{d}v\tag{4}$$

donde $k=\log(K)$ y $\psi_T$ es la transformada de Fourier de $e^{ck}C(t)$. Una gran diferencia aquí con mi ecuación (3) es el intercambio de $\varphi(-z)$ por $e^{ivk}$.

¿Cómo derivo mi ecuación (3)? ¿Cómo está relacionada con la ecuación (4)? ¿Qué es $f$ y $\hat{f}$?

5voto

drN Puntos 571

La ecuación (11) en el artículo de Kammeyer y Kienitz es una fórmula de valoración de opciones muy conocida y popular. Se remonta al trabajo de Lewis (2001), ver Teorema 3.2 en el artículo de Lewis.


Fórmula Original de Lewis (2001)

La fórmula en Lewis para el valor de un derivado de estilo europeo es $$V(S_0) = \frac{e^{-rT}}{2\pi} \int_{\color{red}{i}\nu-\infty}^{\color{red}{i}\nu+\infty}\varphi_T(-z)\hat{w}(z)\text{d}z,$$ donde

  • $\nu$ es un número real. Define el camino a lo largo del cual integramos en el plano complejo: $\{z\in\mathbb{C}:\text{Im}(z)=\nu\}$. Aporto más información sobre esto a continuación.
  • $\varphi_T$ es la función característica (generalizada) de $\ln(S_T)$, el precio de la acción logarítmico terminal en el cual depende nuestro pago de la opción
  • $w$ es la función de pago (como función de $\ln(S_T)$). Para una opción de compra estándar, $w(s)=\max\{e^s-K,0\}$. La función $\hat{w}$ es la transformada de Fourier (generalizada) de la función $w$.

Nota: tanto $\varphi_T$ como $\hat{w}$ se evalúan en puntos en el plano complejo (¡no necesariamente en la línea real!)


Tu Fórmula

Kammeyer y Kienitz afirman que el valor en el tiempo $t$ de una opción de compra es $$C(t)=\frac{e^{-r(T-t)}}{2\pi} \int_{c-\infty}^{c+\infty} \varphi(-z)\hat{f}(z)\text{d}z.$$ Primero, dos puntos importantes

  • Hay un pequeño error tipográfico en la fórmula. Los límites de la integral deben ser $\color{red}{i}c-\infty$ y $\color{red}{i}c+\infty$.
  • Tu fórmula de valoración de opciones es para el precio de la opción en el tiempo $t$. Por lo tanto, todo es condicional en $\mathcal{F}_t$, la filtración generada hasta el tiempo $t$, ver esta respuesta. Lewis (2001) simplemente establece $t=0$.

El resto es idéntico a la fórmula original de Lewis. $\varphi$ es la función característica de $\ln(S_T)$, condicional en $\mathcal{F}_t$ y $f$ es la función de pago y $\hat{f}$ es su transformada de Fourier (generalizada).


Transformada de Fourier de la función de pago

Sea $f(x)=\max\{e^x-e^k,0\}$ el pago de una opción de compra estándar con precio de ejercicio $K=e^k$. ¡Esta función no está en $L^1$ y no tiene una transformada de Fourier tradicional! Sin embargo, tiene una transformada de Fourier generalizada. Normalmente, si $f:\mathbb{R}\to\mathbb{R}$, entonces definimos la transformada de Fourier (en finanzas) como $\hat{f}:\mathbb{R}\to\mathbb{C}, u\mapsto \int_\mathbb{R} e^{iux}f(x)\text{d}x$. Para que esta integral exista, $f$ necesita decaer rápidamente o tener un soporte compacto. La función de pago no cumple con esto.

La transformada de Fourier generalizada de $f$ es $\hat{f}:\mathcal{S}_f\subset\mathbb{C}\to\mathbb{C}, u\mapsto \int_\mathbb{R} e^{iux}f(x)\text{d}x$. ¡Por lo tanto, está definida para un subconjunto de los números complejos! Resulta que este $\mathcal{S}_f$ es una franja horizontal en el plano complejo. Podemos calcular la transformada para la función de pago de la siguiente manera \begin{align} \hat{f}(u) &= \int_{-\infty}^\infty e^{iux} \left(e^x-e^k\right)^+ \text{d}x \\ &= \int_k^\infty \left(e^{x(iu+1)} - Ke^{iux}\right) \text{d}x \\ &= \left[ \frac{e^{x(iu+1)}}{iu+1} - K\frac{e^{iux}}{iu}\right]_{x=k}^{x=\infty} \\ &= \frac{e^{ik(ui)}}{u(u-i)}. \end{align} Este último paso solo es válido si el término realmente se desvanece cuando $x\to\infty$. Esto solo ocurre si $\text{Im}(z)>1$. Por lo tanto, la franja para el pago de la opción de compra es $\mathcal{S}_f=\{z\in\mathbb{C}:\text{Im}(z)>1\}$. De manera similar, para una opción de venta, tenemos $\mathcal{S}_f=\{z\in\mathbb{C}:\text{Im}(z)<0\}$. Estas son las franjas de integración en las cuales las funciones de pago tienen una transformada de Fourier válida. Tenga en cuenta que ambas franjas excluyen la línea real (es decir, no hay una transformada de Fourier estándar).


Prueba de la Fórmula de Valoración de Opciones de Lewis

Comenzando con la valoración estándar de riesgo neutral, \begin{align} V &= e^{-rT}\mathbb{E}^\mathbb{Q}[w(\ln(S_T)] \\ &=e^{-rT}\mathbb{E}^\mathbb{Q}\left[\frac{1}{2\pi}\int_{i\nu-\infty}^{i\nu+\infty}e^{-iz\ln(S_T)}\hat{w}(z)\text{d}z\right] \\ &=\frac{e^{-rT}}{2\pi}\int_{i\nu-\infty}^{i\nu+\infty}\mathbb{E}^\mathbb{Q}\left[e^{i(-z)\ln(S_T)}\right]\hat{w}(z)\text{d}z \\ &=\frac{e^{-rT}}{2\pi}\int_{i\nu-\infty}^{i\nu+\infty}\varphi_T(-z)\hat{w}(z)\text{d}z \\ \end{align} Aquí, solo estamos utilizando la definición de transformadas de Fourier (inversas) generalizadas y el teorema de Fubini. También es posible una prueba utilizando el teorema de Plancherel (o el teorema de Parseval). Para que Fubini aplique y las integrales estén bien definidas, necesitamos integrar a lo largo de un camino en el complejo donde todos los términos estén bien definidos, de ahí la condición de $\nu\in\mathcal{S}_V=\mathcal{S}_w\cap\mathcal{S}_f^*$.

3voto

ir7 Puntos 435

Solo una nota para agregar a la respuesta anterior. El parámetro de amortiguación $c$, número real, se convierte en la parte imaginaria de un número complejo debido a esta simple observación:

$$f_{c}(x) := {\rm e}^{-c x}f(x)$$

$$ \hat{f_c}(x) = \int {\rm e}^{ixy}f_c(y) dy = \int {\rm e}^{i(x+ic)y}f(y) dy = \hat{f}(x+ic) $$

(el sombrero se coloca en dos funciones diferentes, $f$ y $f_c$). Entonces (genéricamente):

$$ E[f(X)] = \int {\rm e}^{cy} f_c(y) q_X(y) dy = \int {\rm e}^{cy} \left(1/2\pi \int {\rm e}^{-ixy} \hat{f_c}(x) dx \right) q_X(y) dy $$

$$ \stackrel{Fubini}{=} 1/2\pi\int \left( \int {\rm e}^{-i(x+ic)y} q_X(y) dy \right) \hat{f_c}(x) dx $$

$$ \stackrel{observación}{=} 1/2\pi \int \phi_X(-(x+ic))\hat{f}(x+ic) dx $$

$$ = 1/2\pi \int_{-\infty+ic}^{\infty +ic} \phi_X(-z)\hat{f}(z) dz $$

(Con el cálculo de $\hat{f}$ para el payoff de una llamada en la respuesta anterior y la relación (6) entre $\phi_T$ y $\psi_T$ en el paper de Carr-Madan, deberíamos obtener la conciliación.)

0 votos

¡Buena adición, +1! ¿No usaría la transformada inversa de Fourier $\frac{1}{2\pi}$ en lugar de $\frac{1}{\pi}$?

0 votos

@Kevin Dejé las integrales sin límites :) (dije 'genéricamente'). La raíz de tal dilema podría ser: si tanto $h$ como $\hat{h}$ son integrables y $h$ es una función real, entonces $x \rightarrow e^{-ixy}\hat{h}(x)$ es par, por lo que los límites pueden ser $(-\infty$, $\infty)$ y se usa $2\pi$, o $(0$, $\infty)$ y se usa $\pi$. Lo editaré a $2\pi$.

0 votos

Totalmente de acuerdo, cada vez que expresamos un precio de opción de valor real como $\frac{1}{2\pi}\int_\mathbb{R} e^{-iux}\phi(u)\text{d}u$, entonces la simetría normalmente entra en acción y obtenemos $\frac{1}{\pi}\int_0^\infty \text{Re}\left(e^{-iux}\phi(u)\right)\text{d}u$. Solo preguntaba porque la última línea en la respuesta sugiere que integramos a lo largo de toda la línea real (desplazada), no solo en su parte positiva. ¡Muchas gracias por la aclaración! :)

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