Este trabalho tem como objetivo desenvolver uma metodologia de análise probabilística do comportamento de um sistema massa-mola-amortecedor submetido a carregamentos aleatórios, modelados como processos estocásticos estacionários, abordando o problema tanto de forma numérica quanto analítica. Essa abordagem se justifica pela relevância de vibrações aleatórias em diversos problemas de engenharia. Com o auxílio do MATLAB, foram implementados algoritmos baseados na metodologia de Monte Carlo, envolvendo a geração de valores pseudoaleatórios de variáveis aleatórias, a construção de modelos estatísticos e a análise da convergência entre estatísticas teóricas e amostrais. As ferramentas de análise utilizadas incluem gráficos envelope, histogramas univariados e bivariados, distâncias de Wasserstein, densidades espectrais da resposta e investigações sobre a influência do amortecimento do sistema na convergência da resposta para um processo estocástico estacionário. Por fim, foram comparados os resultados das abordagens analítica e numérica. A solução analítica foi obtida por meio da derivação da equação da dinâmica, enquanto as aproximações numéricas foram obtidas com o método de Runge-Kutta. A comparação avaliou vantagens como o custo computacional e a aplicabilidade da metodologia para sistemas mais complexos.