[求助]相位恢复算法 [复制链接]

有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 Ya%-/u  
v-(Ry<fT9  
2d centered Fourier transform，居中的傅里叶变换 [,s{/32s  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% /| nZ)?  
% Citation for this code and algorithm: MLu@|Xgh  
% Tatiana Latychevskaia and Hans-Werner Fink aP#nK  
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", q_5hKipd\b  
% Appl. Optics 54, 2424 - 2434 (2015) 7W)*IJ  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Ia>07av  
% The code is written by Tatiana Latychevskaia, 2002 kOu C@~,  
% The version of Matlab for this code is R2010b {S+ $C  
*,hg+?lZ  
function [out] = FT2Dc(u0); x]e&G!|  
V^v?;f?  
[Nx, Ny] = size(1i,1j); _TGv"c@V  
OA_WjTwDs  
f1 = zeros(Nx,Ny); TzPG(f  
NCid`a$  
for ii = 1:Nx
dX720/R  
    for jj = 1:Ny "@P)  
         f1(ii,jj) = exp(1i*pi*(ii + jj)); &d!Q%  
    end |a>W9Ym  
end )<<}8Fs  
~
d9R:t1  
FT = fft2(f1.*in); M,uQ8SZA[  
!=-{$& {  
out = f1.*FT; | Ns-l (l  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

求大神指教