W=imread（‘D:\black.png‘）;%讀入原始圖像
O?=rgb2gray（W）;
figure;
subplot（121）;
imshow（O）%title（‘原始圖像‘）;%顯示原始圖像
Q=im2double（O）;
known_abs_fourier=abs（fft2（O））;
[MN]=size（known_abs_fourier）;
phase_estimate=rand（MN）;

P=Q.*exp（1i*2*pi*phase_estimate）;%數字化圖像乘以隨機相位因子平滑傅里葉變換譜
F=fftshift（fft2（P））;%二維離散傅里葉變換并移譜
C=real（F）;%二維離散傅里葉變換實部
D=imag（F）;%二維離散傅里葉變換