He implementado el modelo de mercado del Libor en Matlab. Cuando genero una serie de trayectorias, observo que algunas de ellas explotan. ¿Alguien tiene una idea de lo que puede causar esto?
Ya he intentado resolver el problema disminuyendo el paso de tiempo (hasta dt=0,001) para reducir el error y también simulando con el esquema log-Euler en lugar del Euler "normal". En ambos casos no se resolvió el problema, ya que algunas de las trayectorias de los tipos Libor siguen siendo divergentes.
Específicos:
Simulo los tipos Libor a plazo bajo la medida spot, cuya dinámica viene dada por: $$dL_n\left(t\right)=\sigma_n\left(t\right)L_n\left(t\right)\sum_{j=q\left(t\right)}^n \frac{\tau_j \rho_{j,n} \sigma_j\left(t\right)L_j\left(t\right)}{1+\tau_j L_j\left(t\right)}dt + \sigma_n\left(t\right)L_n\left(t\right)dW\left(t\right)$$ donde $$L_n\left(t\right):=L\left(t;T_n,T_{n+1}\right),$$ $$\tau_n = T_{n+1}-T_n,$$ $$\sigma_n\left(t\right) = k_n \left[\left(a+b\left(T_n-t\right)\right)e^{-c\left(T_n-t\right)}+d\right],$$ función de índice $q\left(t\right)$ se define por $$T_{q\left(t\right)-1}\leq t < T_{q\left(t\right)},$$ $W$ es un movimiento browniano bajo la medida del punto.