%% *******************************************************************
% GPSA算法仿真(等步长,已知移相)
% 作者:James_Ray_Murphy
% 参考文献:高精度干涉测量随机移相技术研究_苏志德
%****************************%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;clc;close all;
delta = [0 pi/2 pi pi*3/2];
N = 512;
xmax = 1;
ymax = 1;
x = linspace(-xmax,xmax,N);
y = linspace(-ymax,ymax,N);
[X,Y] = meshgrid(x,y);
A = 145*exp(-1.8*(X.^2+Y.^2));
B = 100*exp(-0.2*(X.^2+Y.^2));
phi = 5*pi*(X.^2+Y.^2);
figure;mesh(phi);title('初始相位图');
I0 = A + B.*cos(phi+delta(1));
I1 = A + B.*cos(phi+delta(2));
I2 = A + B.*cos(phi+delta(3));
I3 = A + B.*cos(phi+delta(4));
figure;subplot(221);imshow(I0,[]);title('四步非等步长相位图de
% GPSA算法仿真(等步长,已知移相)
% 作者:James_Ray_Murphy
% 参考文献:高精度干涉测量随机移相技术研究_苏志德
%****************************%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;clc;close all;
delta = [0 pi/2 pi pi*3/2];
N = 512;
xmax = 1;
ymax = 1;
x = linspace(-xmax,xmax,N);
y = linspace(-ymax,ymax,N);
[X,Y] = meshgrid(x,y);
A = 145*exp(-1.8*(X.^2+Y.^2));
B = 100*exp(-0.2*(X.^2+Y.^2));
phi = 5*pi*(X.^2+Y.^2);
figure;mesh(phi);title('初始相位图');
I0 = A + B.*cos(phi+delta(1));
I1 = A + B.*cos(phi+delta(2));
I2 = A + B.*cos(phi+delta(3));
I3 = A + B.*cos(phi+delta(4));
figure;subplot(221);imshow(I0,[]);title('四步非等步长相位图de