To retrieve the phase from the noisy measured intensities in the diffraction planes, an iterative Wiener deconvolution based method is proposed. With the same iterative scheme as the iterative angular spectrum method (IAS), the propagation of the optical wave function between the input plane and the diffraction planes is calculated by Wiener deconvolution in this method. The angular spectrum convolution kernel used in the iterative angular spectrum method is incorporated into the Wiener filter. The simulation experiments show that the proposed method can reduce the impact of the noise on the retrieved phase and performed better than the pre-denoising method. Furthermore, the proposed method exhibits great advantage compared to IAS for retrieving the complicated phase distribution from two measured intensities.