Unfortunately it's for 2D signals (image analysis), but I believe his conclusion would also apply to 1D signal.
J.F. Kirby, "Which wavelet best reproduces the Fourier power spectrum?", Computers & Geosciences 31 (2005) 846–864
Basically, his conclusion is to go with the Fan wavelet, which is a 2D rotated version of the Morlet wavelet. In 1D, I'd suggest the complex Morlet. It's the mix of real and complex part that allows for a good similarity to a Fourier power spectrum.
In better precision, here what it should look like, converted to 1D from Kirby (2005):
$$ \Psi = exp\Big(-\frac{ik_0x}{\lambda} - \frac{x^2}{2\lambda^2}\Big),$$
where $\lambda$ is the scale you're looking at, and $k_0=5.336$ is a constant selected to give the best "scale sampling" vs "frequency sampling". I didn't include the normalization constant because in every computational situation, it's better to just divide the final wavelet by its maximum value, and subtract its average. It gives pretty much the same result with less headache.
Basically, the complex Morlet wavelet is a Fourier transform "wave" ($exp(-i k_0 x/\lambda)$) bounded by a Gaussian kernel ($exp(-x^2/2)$). I suspect you might get a good power spectrum using only the real part (using $cos(x) \cdot exp(-x^2/2)$), but you would loose phase information.
Try comparing the spectrum obtained from a Fourier transform, from a complex Morlet and from a real Morlet. Watch out for bad/non-standard normalization found in many FFT algorithms.