Embedded data are used to retrieve phases quicker with high accuracy in phase-modulated holographic data storage (HDS). We propose a method to design an embedded data distribution using iterations to enhance the intensity of the high-frequency signal in the Fourier spectrum. The proposed method increases the antinoise performance and signal-to-noise ratio (SNR) of the Fourier spectrum distribution, realizing a more efficient phase retrieval. Experiments indicate that the bit error rate (BER) of this method can be reduced by a factor of one after 10 iterations.