clear?all
close?all
clc

phi0=peaks（200）*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