¿Cuál es la mejor manera de simular el corto de la tasa de $r(t)$ en un simple factor de Casco de proceso en Blanco? Supongamos que tengo
$$ dr(t) = (\theta(t)-\alpha r(t))dt+\sigma dW_t $$
donde $\theta(t)$ es calibrado para intercambiar curva, las constantes de $\alpha$ y $\sigma$ están calibrados para tapas cerrado utilizando el formulario de solución para el bono cupón cero opciones. La mejor manera que puedo pensar es una de Euler discretisation, que es:
$$ r(t+\Delta t) = r(t) + \theta(t)\Delta t - \alpha r(t) \Delta t + \sigma \sqrt {\Delta t} Z $$ donde $Z \sim N(0,1)$. En este caso, necesito $t$ para ir de 0 a 10 años, idealmente en incrementos de 0.25. Pero con Euler, que había necesidad de utilizar pequeñas $\Delta t$, así que tal vez 0,025 o menos? Una vez que tenga la cadena de $r(t)$, que fácilmente se puede calcular $P(t,T)$ bonos de cupón cero.
Agradezco cualquier otra idea o si alguien puede me apunte en la dirección correcta. Soy bastante nuevo en las tasas de modelado!