This paper concerns the random source problems for the time-harmonic acoustic and elastic wave equations in two and three dimensions. The goal is to determine the compactly supported external force from the radiated wave field measured in a domain away from the source region. The source is assumed to be a microlocally isotropic generalized Gaussian random function such that its covariance operator is a classical pseudo-differential operator. Given such a distributional source, the direct problem is shown to have a unique solution by using an integral equation approach and the Sobolev embedding...