x=imread（‘lena.bmp‘）;
x=rgb2gray（x）;
x=double（x）;
mask=zeros（size（x））;
mask（200:300200:300）=1;
x1=x.*mask;
xm=uint8（x1）;
xf=fft2（x1）;
absfx=abs（xf）;
y=randn（size（x））;
iteration=400;
beta=0.97;
i=1:iteration;
for?iter=1:iteration
????y=myHIO（y?absfx?mask?beta）;
%??????????y=myER（y?absfx?mask）;
??????????x1=uint8（x1）;
??????????y=uint8（real（y））;
??????????figure（1）
??????????subplot（121）?imshow（x1）?axis?image?title（‘ER‘）;
??????????subplot（122）?imshow（y）?axis?image?title（iter）;
??????????y=double（y）;???????????
end
PSNR?=?Psnr（x1y）;
y=uint8（real（y））;
imwrite（y‘poisson?HIO?400.bmp‘‘bmp‘）;
figure（2）
subplot（121）imshow（xm[]）title（‘原始圖像‘）;
subplot（122）imshow（y[]）title（‘相位恢復后圖象‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2015-03-22?22:53??HIO?algorithm\
?????文件??????786486??2012-04-09?19:55??HIO?algorithm\lena.bmp
?????文件?????????750??2015-03-22?22:56??HIO?algorithm\main.m
?????文件?????????106??2014-04-21?17:34??HIO?algorithm\myHIO.m
?????文件??????????88??2014-04-21?17:21??HIO?algorithm\Pm1.m
?????文件??????????50??2014-06-10?16:16??HIO?algorithm\Ps1.m
?????文件?????????398??2014-05-05?09:50??HIO?algorithm\Psnr.m