最小二乘法解包裹

最小二乘法解包裹的使用代码，包括LSunwrap.m，unwrapphase.m，wrapphase.m

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;
%

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           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