%Background and Modulation terms

a = (255/2);

b = (255/2);

%phase map

P1 = 2*pi*10*sqrt(((X.^2+Y.^2)));

P2 = 2*pi*60*exp(((X.^2+Y.^2)));

P2 = P2+5.*peaks(768);

% P = P-mean2(P);

%fringe patterns

A_P1 = a+b.*(cos(P1)); A_P2 = a+b.*(cos(P1+2*pi/3)); A_P3 = a+b.*(cos(P1+4*pi/3));

B_P1 = a+b.*(cos(P2)); B_P2 = a+b.*(cos(P2+2*pi/3)); B_P3 = a+b.*(cos(P2+4*pi/3));

figure;

subplot(231); imshow(uint8(A_P1)); subplot(232); imshow(uint8(A_P2)); subplot(233); imshow(uint8(A_P3));

subplot(234); imshow(uint8(B_P1)); subplot(235); imshow(uint8(B_P2)); subplot(236); imshow(uint8(B_P3));