9 votos

Último Paso en el cálculo de la opción de precios por debajo de la Heston Modelo de Volatilidad Estocástica

Vamos a: $$ \alpha = -\frac{u^2}{2}-\frac{ui}{2}+dci\\ \beta = \lambda\rho \eta, i, u - j \rho \eta\\ \gamma = \frac{\eta ^2}{2}\\ $$

donde $j \in \{0,1\}$ y $i^2=-1$, $g=\frac{r_-}{r_+}$ y $r_{\pm}=\frac{\beta\pm \sqrt{\beta^2-4\alpha\gamma}}{2\gamma}=\frac{\beta \pm d }{\eta^2}$. A continuación, vamos a: $$ D(u,\tau)=r_-\frac{1-e^{-d\tau}}{1-ge^{-d\tau}}\\ C(u,\tau)=\lambda \left[\tau r_--\frac{2}{\eta^2}log(\frac{1-ge^{-d\tau}}{1-g}) \right] $$ Donde $\lambda$ es una constante.

Para la función: $$ \hat{P}_j(u,v,\tau)=\frac{1}{ui}exp\left[C(u,\tau)\bar{v}+D(u,\tau)v \right] $$ $\bar{v}\in \mathscr R$.

Demostrar que: $$ P_j(x,v,\tau)=\int_{-\infty}^{\infty}\frac{e^{iux}}{2\pi} \hat{P_j}(u,v,\tau)du\\=\frac{1}{2}+\frac{1}{\pi}\int_0^\infty Re\left[\frac{exp[C_j(u,\tau)\bar{v}+D_j(u,\tau)v+iux]}{ui} \right]du $$

Este resultado viene en un texto que estoy leyendo en los modelos de volatilidad estocástica, se indica pero no se ha comprobado, las matemáticas involucradas es probablemente la forma por encima de mi nivel, pero sería genial ver a los pasos involucrados en probar este. Supongo que no equivale simplemente a mostrar la parte imaginaria de esta integral es igual a la mitad.

No me refiero a la pregunta perezoso, he trabajado a través de todas las otras partes de este relativamente larga de la prueba, esta técnica suele utilizarse mucho en la volatilidad estocástica y la liquidez de los modelos que estoy interesado en estudiar más, y esa es la razón por la que he publicado.

4voto

SDiv Puntos 788

Gatheral deja de lado algo importante y eso hace toda la diferencia para la derivación, pero no el resultado final. Específicamente, en la pg. 17 dice lo siguiente, $$\hat{P}(u,v,0)=\int_{-\infty}^\infty\Theta(x)e^{-iux}dx=\frac{1}{iu},$$ where $\Theta(x)$ is the Heaviside unit step function (as defined in his Eq. 2.7, take a look at his definition and convince yourself that it is basically the same as Wikipedias). But this is not true - there is a Dirac delta missing. See rule 313 column 3 on WP - Tables of important Fourier transforms. The equation should be: $$\hat{P}(u,v,0)=\int_{-\infty}^\infty\Theta(x)e^{-iux}dx=\frac{1}{iu}+\pi\delta(u),$$ where $\delta(x)$ is the Dirac delta distribution. Next, we have $$\hat{P}_j(u,v,t)=\exp(C[u,\tau]\bar{v}+D[u,\tau]v)\hat{P}_j(u,v,0)\\ =(\frac{1}{iu}+\pi\delta[u])\exp(C[u,\tau]\bar{v}+D[u,\tau]v).$$ Then $$\int_{-\infty}^\infty \frac{e^{iux}}{2\pi}\hat{P}_j(u)\text{d}u = \int_{-\infty}^\infty \frac{e^{iux}}{2\pi}(\frac{1}{iu}+\pi\delta[u])\exp(C[u,\tau]\bar{v}+D[u,\tau]v)\text{d}u.$$ Lets split this into two parts. Lets consider the integral over the delta function first. Noting that $C(0,\tau)=0$ and $D(0,\tau)=0$ we have $$ \int_{-\infty}^\infty \frac{e^{iux}}{2\pi}(\pi\delta[u])\exp(C[u,\tau]\bar{v}+D[u,\tau]v)\text{d}u=\frac{1}{2}.$$ Now consider the part with $\frac{1}{iu}$. Let $$\hat{P}'_j:=\frac{1}{iu}\exp(C[u,\tau]\bar{v}+D[u,\tau]v)$$

La cosa importante a observar es, asumiendo $j,\rho,\eta,\lambda$ $\tau$ son reales, que $r_\pm(u)$ es conjugado-incluso en $u$. Es decir, $r_\pm(-u)=r_\pm^\ast(u)$ donde $\ast$ denota el complejo conjugado. A continuación, g y d y, por tanto, $D$ $C$ y, por tanto, $\hat{P}'_j$ son conjugadas. Hago la suposición adicional de que $\bar{v}$ indica la media del parámetro real $v$ y no el conjugado complejo. Entonces, $$\int_{-\infty}^\infty \frac{e^{iux}}{2\pi}\hat{P}'_j(u)\text{d}u = \int_{-\infty}^0 \frac{e^{iux}}{2\pi}\hat{P}'_j(u)\text{d}u + \int_{0}^\infty \frac{e^{iux}}{2\pi}\hat{P}'_j(u)\text{d}u \\ = \int_{0}^\infty \frac{e^{-iux}}{2\pi}\hat{P}'_j(-u)\text{d}u + \int_{0}^\infty \frac{e^{iux}}{2\pi}\hat{P}'_j(u)\text{d}u \\ = \frac{1}{2\pi}\int_{0}^\infty \left [ e^{-iux}\hat{P'}_j^\ast(u) + e^{iux}\hat{P}'_j(u) \right ]\text{d}u \\ = \frac{1}{\pi}\int_{0}^\infty \text{Re}\left[e^{iux}\hat{P}'_j(u) \right ]\text{d}u. $$ On the last line I have used the identity $2\texto{Re}(z)=z+z^\ast$.

Poniendo todo junto da $$P_j(x,v,\tau)=\frac{1}{2}+\frac{1}{\pi}\int_{0}^\infty \text{Re}\left[e^{iux}\hat{P}'_j(u) \right ]\text{d}u,$$, que es la solución en el libro.

Actualización: Si no se queda satisfecho y se muestran escépticos de que Gatheral omitido el delta, considere esto: En la pagina 19 Gatheral da las condiciones $C(u,0)=0$$D(u,0)=0$. Observe de nuevo que de acuerdo a su definición, $\hat{P}(u,v,0)=1/(iu)$. Luego, utilizando la inversa de la transformación como se define en la ecuación. 2.8, esto da $$\Theta(x)=\lim_{\tau \rightarrow 0}P_j(x,v,\tau)=\lim_{\tau \rightarrow 0}\int_{-\infty}^\infty \frac{1}{2\pi iu}\exp(C[u,\tau]\bar{v}+D[u,\tau]v)e^{iux}du\\=\int_{-\infty}^\infty \frac{1}{2\pi iu}e^{iux}du \\ = \frac{1}{2}\text{sgn}(x) \neq \Theta(x)$$ where on the last line we used the inverse transform version of rule 309 on the WP fourier transform page. Notice that this answer is off by exactly $1/2$! That is, $1/2+ (1/2)\texto{sgn}(x)=\Theta(x)$. En otras palabras, Gatheral omitido el delta y su derivación como impreso no es auto-consistente. Cabe señalar, sin embargo, que la omisión no afecta su resultado final (que sigue siendo válido).

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X