iterative=300;???????????%設迭代次數為300次吧
imagename=‘test.jpg‘;????%你想要提取相位的圖像名稱
phaseimage=‘phase.png‘;??%要保存的相位圖像名稱


%空域輸入圖像的幅度（是已知的，也就是清晰的圖像，它的灰度就是幅值）和相位圖像（待恢復）
known_abs_spatial=imread（imagename）;???????????%作為輸入圖像的幅度，是已知的
known_abs_spatial?=rgb2gray（known_abs_spatial）;%注意要用單通道圖像做實驗，如果你讀取的是彩色圖像，那就吧這行取消注釋變成灰度圖像吧
known_abs_spatial=im2double（known_abs_spatial）;?%將圖像灰度映射到0～1

unknown_phase=known_abs_spatial;????????????????%Peppers圖像作為輸入圖像的相位，也即為待恢復的數據，
????????????????????????????????????????????????%要求它和known_abs_spatial大小一致，所以這里直接賦值就好了
unknown_phase=im2double（unknown_phase）;?????????%將圖像灰度映射到0～1
unknown_phase2=unknown_phase*2*pi;??????????????%相位范圍映射到0-2*pi
unknown_phase2（unknown_phase2>pi）=unknown_phase2（unknown_phase2>pi）-2*pi;%進一步映射至[-pi+pi]

[widthlength]=size（known_abs_spatial）;?????????%獲取Lena圖像的大小
input=known_abs_spatial.*exp（1i*unknown_phase2）;?%最終輸入圖像:幅度*e^（i*相位角度），它是復數圖像
known_abs_fourier=abs（fft2（input））;?????????????%先將input圖像進行傅立葉變換，然后取模，就是傅氏變換后的幅度
%以下開始迭代求相位
phase_estimate=pi*rand（widthlength）;???????????%這是生成了一副大小為（width*length）的圖像
??????????????????????????????????????