Este documento se enfoca en el problema inverso de identificar la función fuente dependiente del espacio y el valor inicial de la ecuación de onda de difusión no homogénea fraccional en el tiempo a partir de datos ruidosos medidos en el tiempo final en un caso multidimensional. Se presenta un método de regularización de molificación basado en un núcleo exponencial bilateral para resolver la mal-posedness del problema por primera vez. Se obtienen estimaciones de error con una estrategia a priori y una regla de elección a posteriori para encontrar el parámetro de regularización. Experimentos numéricos de interés muestran que nuestro método propuesto es efectivo y robusto con respecto al ruido de perturbación en los datos.