Generalized propagation formulas for the elements of the cross-spectral density matrix of stochastic electromagnetic beams on inverse propagation through an axially symmetrical or nonsymmetrical optical system are derived with the help of Fourier transform and inverse Fourier transform. As an example, we apply the formula to the inverse source problem of stochastic electromagnetic Gaussian Schell-model beams.