Integración acumulativa con respecto al precio del bono del modelo Vasicek y su precio a plazo

Question

(Mi pregunta)

Por favor, muéstrame cómo calcular la siguiente expectativa con su proceso de cálculo. Además, $B_t$ es S.B.M.

$$E\left[ \exp \left( - \int^T_t \int^u_0 \sigma e^{-b(u-s)} d B_s du \right) \right]$$

(Gracias por su ayuda de antemano).

---

(Enlace cruzado)

He publicado la misma pregunta en https://math.stackexchange.com/questions/3343590/cumulative-integration-with-regard-to-vasicek-models-bond-price-and-its-forward

---

(Pregunta original)

Resolver $P(t, T)$ con el siguiente modelo

$$dr_t=-br_t dt + \sigma dB_t$$

---

(Mi consideración)

• Puño,

$$r_u=e^{-bu} r_0 + \int^u_0 \sigma e^{-b(u-s)} dB_s$$

• Segundo,

\begin{eqnarray} P(t, T) &=& E \left[ \exp \left( - \int^T_t r_u du \right) \middle| \mathcal{F}_t \right] \\ &=& E\left[ \exp \left( - \int^T_t \left(e^{-bu} r_0 + \int^u_0 \sigma e^{-b(u-s)} dB_s \right) du \right) \middle| \mathcal{F}_t \right] \\ &=& E\left[ \exp \left( - \int^T_t e^{-bu} r_0 du - \int^T_t \int^u_0 \sigma e^{-b(u-s)} dB_s du \right) \middle| \mathcal{F}_t \right] \\ &=& \frac{r_0}{b} (e^{-bT}-e^{-bt}) E\left[ \exp \left(- \int^T_t \int^u_0 \sigma e^{-b(u-s)} dB_s du \right) \middle| \mathcal{F}_t \right] \end{eqnarray}

• En tercer lugar, supongo que utilizar la siguiente fórmula, pero no puedo tener ninguna idea para reemplazar el orden de integración.

$$E\left[ \exp \left( \int^T_t f(s) dB_s \right) \middle| \mathcal{F}_t \right] = \exp \left( \frac{1}{2} \int^T_t f(s)^2 ds \right) $$

Gracias por su ayuda de antemano.

Answer 1

(Mi respuesta)

el precio del bono Vasicek y su precio a plazo

• Recuerda el resultado del ejercicio 5.2.(1) o del ejercicio 4.5.(10). \begin{eqnarray} P(t, T) &=& E \left[ \exp \left( - \int^T_t r_u du \right) \middle| \mathcal{F}_t \right] \\ &=& E \left[ \exp \left( - \int^T_t \left( e^{-bu} r_0 + \sigma \int^u_0 e^{-b(u-s)} dB_s \right) du\right) \middle| \mathcal{F}_t \right] \\ &=& E \left[ \exp \left( - r_0 \int^T_t e^{-bu} du - \int^T_t\int^u_0 \sigma e^{-b(u-s)} dB_s du \right) \middle| \mathcal{F}_t \right] \\ &=& E \left[ \exp \left( - r_0 \left[ -\frac{1}{b} e^{-bu} \right]^T_t - \int^T_t\int^u_0 \sigma e^{-b(u-s)} dB_s du \right) \middle| \mathcal{F}_t \right] \\ %&=& E \left[ \exp \left( \frac{r_0}{b} (e^{-bT} -e^{-bt} ) - \int^T_t\int^u_0 \sigma e^{-b(u-s)} dB_s du \right) \middle| \mathcal{F}_t \right] \\ &=& \exp \left( \frac{r_0}{b} (e^{-bT} -e^{-bt} ) \right) \nonumber \\ && \qquad E \left[ \exp \left( - \int^T_t\int^u_0 \sigma e^{-b(u-s)} dB_s du \right) \middle| \mathcal{F}_t \right] \end{eqnarray}

• Aquí, se calcula la parte exponencial dentro de la expectativa utilizando la siguiente fórmula (Fórmula de integración acumulativa sustituyendo $du$ y $dB_s$ . (Desarrollado el 5 de septiembre de 2019 por Kojima)) \begin{eqnarray} \int^t_0 \int^u_0 dB_s \ du &=& \int^t_0 \int^u_s du \ dB_s \\ \int^T_t \int^u_0 dB_s \ du &=& \int^T_0 \int^u_s du \ dB_s - \int^t_0 \int^u_s du \ dB_s \end{eqnarray}

• Por lo tanto, \begin{eqnarray} && - \int^T_t\int^u_0 \sigma e^{-b(u-s)} dB_s \ du \nonumber \\ && \qquad = - \int^T_0 \int^u_0 \sigma e^{-b(u-s)} dB_s \ du + \int^t_0 \int^u_0 \sigma e^{-b(u-s)} dB_s \ du \\ && \qquad = - \int^T_0 \int^T_s \sigma e^{-b(u-s)} du \ dB_s + \int^t_0 \int^t_s \sigma e^{-b(u-s)} du \ dB_s \\ && \qquad = - \int^T_0 \left[ -\frac{ \sigma }{b} e^{-b(u-s)} \right]^T_s dB_s + \int^t_0 \left[ -\frac{ \sigma }{b} e^{-b(u-s)} \right]^t_s dB_s \\ && \qquad = \int^T_0 \frac{ \sigma }{b} \left( e^{-b(T-s)} -1 \right) dB_s - \int^t_0 \frac{ \sigma }{b} \left( e^{-b(t-s)} -1 \right) dB_s \\ && \qquad = \int^T_t \frac{ \sigma }{b} \left( e^{-b(T-s)} -1 \right) dB_s \nonumber \\ && \qquad \qquad + \int^t_0 \frac{ \sigma }{b} \left( e^{-b(T-s)} -1 \right) dB_s - \int^t_0 \frac{ \sigma }{b} \left( e^{-b(t-s)} -1 \right) dB_s \end{eqnarray}

• Sustituye los resultados anteriores en la parte exponencial. \begin{eqnarray} P(t, T) &=& \exp \left( \frac{r_0}{b} (e^{-bT} -e^{-bt} ) \right) \nonumber \\ && \qquad E \left[ \exp \left( - \int^T_t\int^u_0 \sigma e^{-b(u-s)} dB_s du \right) \middle| \mathcal{F}_t \right] \\ &=& \exp \left( \frac{r_0}{b} (e^{-bT} -e^{-bt} ) \right) \nonumber \\ && \ \cdot \exp \left( \int^t_0 \frac{ \sigma }{b} \left( e^{-b(T-s)} -1 \right) dB_s - \int^t_0 \frac{ \sigma }{b} \left( e^{-b(t-s)} -1 \right) dB_s \right) \\ && \ E \left[ \exp \left( \int^T_t \frac{ \sigma }{b} \left( e^{-b(T-s)} -1 \right) dB_s \right) \middle| \mathcal{F}_t \right] \end{eqnarray}

• En este caso, utilice la siguiente fórmula. \begin{eqnarray} E\left[ \exp \left( \int^T_t f(s) dB_s \right) \middle| \mathcal{F}_t \right] = \exp \left( \frac{1}{2} \int^T_t f(s)^2 ds \right) \end{eqnarray}

• Por lo tanto, \begin{eqnarray} P(t, T) &=& \exp \left( \frac{r_0}{b} (e^{-bT} -e^{-bt} ) \right) \nonumber \\ && \ \cdot \exp \left( \int^t_0 \frac{ \sigma }{b} \left( e^{-b(T-s)} -1 \right) dB_s - \int^t_0 \frac{ \sigma }{b} \left( e^{-b(t-s)} -1 \right) dB_s \right) \\ && \ E \left[ \exp \left( \int^T_t \frac{ \sigma }{b} \left( e^{-b(T-s)} -1 \right) dB_s \right) \middle| \mathcal{F}_t \right] \\ &=& \exp \left( \frac{r_0}{b} (e^{-bT} -e^{-bt} ) \right) \nonumber \\ && \ \cdot \exp \left( \int^t_0 \frac{ \sigma }{b} \left( e^{-b(T-s)} -1 \right) dB_s - \int^t_0 \frac{ \sigma }{b} \left( e^{-b(t-s)} -1 \right) dB_s \right) \\ && \ \cdot \exp \left( \frac{1}{2} \int^T_t \frac{ \sigma^2 }{b^2} \left( e^{-b(T-s)} -1 \right)^2 ds \right) \end{eqnarray}

• Se calcula la parte exponencial dejando que $\mathscr{P}(t, T, r)$ . \begin{eqnarray} \mathscr{P}(t, T, r) &=& \frac{r_0}{b} (e^{-bT} -e^{-bt} ) \nonumber \\ && \qquad + \int^t_0 \frac{ \sigma }{b} \left( e^{-b(T-s)} -1 \right) dB_s - \int^t_0 \frac{ \sigma }{b} \left( e^{-b(t-s)} -1 \right) dB_s \nonumber \\ && \qquad + \frac{1}{2} \int^T_t \frac{ \sigma^2 }{b^2} \left( e^{-b(T-s)} -1 \right)^2 ds \\ &=& \frac{r_0}{b} (e^{-bT} -e^{-bt} ) \nonumber \\ && \qquad + \int^t_0 \frac{ \sigma }{b} e^{-b(T-s)} dB_s - \int^t_0 \frac{ \sigma }{b} e^{-b(t-s)} dB_s \nonumber \\ && \qquad + \frac{ \sigma^2 }{2b^2} \int^T_t \left( e^{-2b(T-s)} -2e^{-b(T-s)} + 1 \right) ds \\ &=& \frac{r_0}{b} (e^{-bT} -e^{-bt} ) \nonumber \\ && \qquad + \frac{1}{b} e^{-bT} e^{bt}\int^t_0 \sigma e^{-b(t-s)} dB_s - \frac{1}{b} \int^t_0 \sigma e^{-b(t-s)} dB_s \nonumber \\ && \qquad + \frac{ \sigma^2 }{2b^2} \left[ \frac{1}{2b} e^{-2b(T-s)} - \frac{2}{b} e^{-b(T-s)} + s \right]^T_t \end{eqnarray}

• Aquí, utiliza el resultado. \begin{eqnarray} \sigma \int^t_0 e^{-b(t-s)} dB_s &=& r_t - e^{-bt} r_0 \\ \mathscr{P}(t, T, r) &=& \frac{r_0}{b} (e^{-bT} -e^{-bt} ) \nonumber \\ && \qquad + \frac{1}{b} e^{-bT} e^{bt}\int^t_0 \sigma e^{-b(t-s)} dB_s - \frac{1}{b} \int^t_0 \sigma e^{-b(t-s)} dB_s \nonumber \\ && \qquad + \frac{ \sigma^2 }{2b^2} \left[ \frac{1}{2b} e^{-2b(T-s)} - \frac{2}{b} e^{-b(T-s)} + s \right]^T_t \\ &=& \frac{r_0}{b} (e^{-bT} -e^{-bt} ) \nonumber \\ && \qquad + \frac{1}{b} e^{-b(T-t)} (r_t - e^{-bt} r_0 )- \frac{1}{b} (r_t - e^{-bt} r_0 ) \nonumber \\ && \qquad + \frac{ \sigma^2 }{2b^2} \frac{1}{2b} ( 1- e^{-2b(T-t)} ) \nonumber \\ && \qquad - \frac{ \sigma^2 }{2b^2} \frac{2}{b} ( 1- e^{-b(T-t)} ) + \frac{ \sigma^2 }{2b^2} (T-t) \end{eqnarray}

• Sustituye los resultados anteriores en la parte exponencial. De este modo, se llega al siguiente resultado para $P(t, T)$ . \begin{eqnarray} P(t, T) &=& \exp \left( \frac{r_0}{b} (e^{-bT} -e^{-bt} ) \right) \nonumber \\ && \qquad \cdot \exp \left( \frac{1}{b} e^{-b(T-t)} (r_t - e^{-bt} r_0 )- \frac{1}{b} (r_t - e^{-bt} r_0 ) \right) \nonumber \\ && \qquad \cdot \exp \left(\frac{ \sigma^2 }{4b^3} ( 1- e^{-2b(T-t)} ) \nonumber \right) \\ && \qquad \cdot \exp \left( -\frac{ \sigma^2 }{b^3} ( 1- e^{-b(T-t)} ) + \frac{ \sigma^2 }{2b^2} (T-t) \right) \end{eqnarray}

• Recordemos la definición del tipo de interés a plazo $f(t, T, S)$ . \begin{eqnarray} f(t, T, S) &\equiv& - \frac{ \log P(t, S) - \log P(t, T)}{S-T} \\ f(t, T_1, T_2) &\equiv& - \frac{ \log P(t, T_2) - \log P(t, T_1)}{T_2-T_1} \end{eqnarray}

• Aquí se llega a $-\log P(t, T_2)$ . \begin{eqnarray} -\log P(t, T_2) &=& - \frac{r_0}{b} e^{-bT_2} + \frac{r_0}{b} e^{-bt} -\frac{1}{b} e^{-b(T_2-t)} r_t + \frac{1}{b} e^{-b(T_2-t)} e^{-bt} r_0 \nonumber \\ && \qquad + \frac{1}{b} r_t - \frac{1}{b} e^{-bt} r_0 - \frac{\sigma^2}{4b^3} + \frac{\sigma^2}{4b^3} e^{-2b(T_2 -t)} \nonumber \\ && \qquad + \frac{\sigma^2}{b^3} -\frac{\sigma^2}{b^3}e^{-b(T_2 -t)} - \frac{\sigma^2}{2b^2} T_2 + \frac{\sigma^2}{2b^2} t \\ &=& - \frac{r_0}{b} e^{-bT_2} + \frac{r_0}{b} e^{-bt} -\frac{1}{b} e^{-b(T_2-t)} r_t + \frac{1}{b} e^{-bT_2} r_0 \nonumber \\ && \qquad + \frac{1}{b} r_t - \frac{1}{b} e^{-bt} r_0 - \frac{\sigma^2}{4b^3} + \frac{\sigma^2}{4b^3} e^{-2b(T_2 -t)} \nonumber \\ && \qquad + \frac{\sigma^2}{b^3} -\frac{\sigma^2}{b^3}e^{-b(T_2 -t)} - \frac{\sigma^2}{2b^2} T_2 + \frac{\sigma^2}{2b^2} t \end{eqnarray}

• Además, se llega a $\log P(t, T_1)$ . \begin{eqnarray} \log P(t, T_1) &=& \frac{r_0}{b} e^{-bT_1} - \frac{r_0}{b} e^{-bt} + \frac{1}{b} e^{-b(T_1-t)} r_t - \frac{1}{b} e^{-bT_1} r_0 \nonumber \\ && \qquad - \frac{1}{b} r_t + \frac{1}{b} e^{-bt} r_0 + \frac{\sigma^2}{4b^3} - \frac{\sigma^2}{4b^3} e^{-2b(T_1 -t)} \nonumber \\ && \qquad - \frac{\sigma^2}{b^3} + \frac{\sigma^2}{b^3}e^{-b(T_1 -t)} + \frac{\sigma^2}{2b^2} T_1 - \frac{\sigma^2}{2b^2} t \end{eqnarray}

• Por lo tanto, se llega a la siguiente ecuación. \begin{eqnarray} -\log P(t, T_2) + \log P(t, T_1) &=& - \frac{r_t}{b} ( e^{-b(T_2-t)} - e^{-b(T_1-t)} ) \nonumber \\ && \qquad + \frac{\sigma^2}{4b^3} ( e^{-2b(T_2-t)} - e^{-2b(T_1-t)} ) \nonumber \\ && \qquad - \frac{\sigma^2}{b^3} ( e^{-b(T_2-t)} - e^{-b(T_1-t)} ) \nonumber \\ && \qquad - \frac{\sigma^2}{2b^2} (T_2 -T_1) \\ &=& - \frac{\sigma^2}{2b^2} (T_2 -T_1) \nonumber \\ && \qquad - \left( \frac{r_t}{b} + \frac{\sigma^2}{b^3} \right) ( e^{-b(T_2-t)} - e^{-b(T_1-t)} ) \nonumber \\ && \qquad + \frac{\sigma^2}{4b^3} ( e^{-2b(T_2-t)} - e^{-2b(T_1-t)} ) \end{eqnarray}

• Por lo tanto, sustituya el resultado anterior en la definición del tipo de interés a plazo $f(t, T_1, T_2)$ . \begin{eqnarray} f(t, T_1, T_2) &\equiv& - \frac{ \log P(t, T_2) - \log P(t, T_1)}{T_2-T_1} \\ &=& - \frac{\sigma^2}{2b^2} - \frac{1}{T_2 -T_1 } \left( \frac{r_t}{b} + \frac{\sigma^2}{b^3} \right) ( e^{-b(T_2-t)} - e^{-b(T_1-t)} ) \nonumber \\ && \qquad + \frac{1}{T_2 -T_1 } \frac{\sigma^2}{4b^3} ( e^{-2b(T_2-t)} - e^{-2b(T_1-t)} ) \end{eqnarray}

$\square$

Answer 2

(Fórmula de integración acumulativa que sustituye a $du$ y $dB_s$ )

• ¡He desarrollado fórmulas para resolver esto por mí mismo!

\begin{eqnarray} \int^t_0 \int^u_0 dB_s \ du &=& \int^t_0 \int^u_s du \ dB_s \\ \int^T_t \int^u_0 dB_s \ du &=& \int^T_0 \int^u_s du \ dB_s - \int^t_0 \int^u_s du \ dB_s \end{eqnarray}

• Por lo tanto, podemos utilizar la siguiente fórmula como supongo.

$$E\left[ \exp \left( \int^T_t f(s) dB_s \right) \middle| \mathcal{F}_t \right] = \exp \left( \frac{1}{2} \int^T_t f(s)^2 ds \right) $$

• He resuelto el cálculo del precio de los bonos Vasicek utilizando las fórmulas.