clear;
close?all;
clc;

phi0=peaks（512）*3;
figure
imagesc（phi0）
phi=angle（exp（j*phi0））;
[mn]=size（phi）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
phi1=zeros（2*m-12*n-1）;
phi1（1:m1:n）=phi;
phi1（m+1:2*m-11:n）=phi（m:-1:2:）;
phi1（1:mn+1:2*n-1）=phi（:n:-1:2）;
phi1（m+1:2*m-1n+1:2*n-1）=phi（m:-1:2n:-1:2）;
figure
imagesc（phi1）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
phidx=zeros（2*m-12*n-1）;
phidy=zeros（2*m-12*n-1）;
phidx（1:2*m-2:）=angle（exp（j*（phi1（2:2*m-1:）-phi1（1:2*m-2:））））;
phidy（:1:2*n-2）=angle（exp（j*（phi1（:2:2*n-1）-phi1（:1:2*n-2））））;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
lou=zeros（2*m-12*n-1）;
loudx=zeros（2*m-12*n-1）;
loudy=zeros（2*m-12*n-1）;
loudx（2:2*m-1:）=phidx（2:2*m-1:）-phidx（1:2*m-2:）;
loudy（:2:2*n-1）=phidy（:2:2*n-1）-phidy（:1:2*n-2）;
lou=loudx+loudy;
%

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----
?????目錄???????????0??2017-12-10?17:54??LSunwrap\
?????文件????????1171??2017-12-12?17:01??LSunwrap\LSunwrap.m
?????文件??????309084??2009-08-07?09:36??LSunwrap\unwrapphase.fig
?????文件??????310476??2009-08-07?09:36??LSunwrap\wrapphase.fig