GS法计算全息重现

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）;%二维离散傅里叶变换