En este trabajo, proporcionamos un método numérico para discretizar osciladores estocásticos lineales con altas frecuencias constantes impulsados por una fuerza no lineal variable en el tiempo y una fuerza aleatoria. El método presentado se construye a partir de la fórmula de variación de constantes, en la cual aparecen integrales altamente oscilantes. Para proporcionar una discretización adecuada de este tipo de integrales, proponemos reglas de cuadratura basadas en expansiones asintóticas. Se introducen consideraciones teóricas y experimentos numéricos comparando el método con un enfoque estándar en modelos físicos.