Neste estudo, estimamos os futuros do VIX usando um modelo gaussiano, seguindo o artigo de Avellaneda e Papanicolau (2018). Utilizamos o modelo de dois fatores com ambas as variáveis seguindo um processo de Ornstein-Uhlenbeck que considera a estacionariedade, conforme observado a partir de dados empíricos. Adotamos duas estratégias diferentes: na primeira, seguimos rigorosamente a referência acima e incluímos o VIX à vista e todos os futuros do VIX como variáveis de observação (ou CMF - futuros de vencimento constante); na segunda abordagem, estimamos o modelo usando apenas os futuros do VIX. Nos dois casos, foi utilizada a metodologia do filtro de Kalman. Na segunda abordagem, encontramos erros menores na amostra. A estimação do modelo sem o VIX à vista envolveu um menor tempo computacional. A contribuição da dissertação com a estimação alternativa permite a simulação de instrumentos atrelados a contratos mais longos com menores erros. Uma vez estimados os hiperparâmetros, realizamos a simulação de alguns produtos, como ETFs e ETNs. Da perspectiva prática, essa modelagem é útil para entender o comportamento dos contratos futuros do VIX, pois todos os instrumentos derivados do VIX podem ser usados por agentes para estabelecer estratégias para proteger suas posições no mercado.