本软件下载地址未录入，如想要软件请评论，我们会补上！最小二乘解决方案使用代码进行包装, 包括 lsunwrap、unwrapphase、wrapphase
PP=fft2(lou);%对方程进行傅里叶变换
for ii=1:2*m-1 ii
for jj=1:2*n-1
k1=2*cos((ii)*pi/(m));
k2=2*cos((jj)*pi/(n));
PH(ii,jj)=PP(ii,jj)/(k1+k2-4);%在DCT域内两者关系
end
end
phi3=ifft2(PH);%反傅里叶变换得到真实相位分布
