%% CF plot from WF
clear
clc

load('Hy_CF3(x1,x2,y1,y2)_10cm.mat')

L=60;  %sets the number of sample points 
D=0.295/2;    %sets the physical distance range

x = linspace(-D,D,60);
y = linspace(-D,D,60);


figure
pcolor(x,y,squeeze(abs(Hy_CF3(30,:,30,:)))'); %CF2(x1,x2,y1,y2)
axis tight
shading interp
title('Correlation function at 10 cm - WF ')
set(gca,'YDir','normal')
set(gca,'FontSize',13) % Creates an axes and sets its FontSize to 18
%caxis([0 3.15e-3])
colormap('jet'); % set the colorscheme
colorbar
xlabel('$$x_2 [m]$$','Interpreter','latex');
ylabel('$$y_2 [m]$$','Interpreter','latex');

%% CF plot from WF
clear
clc

load('Hy_MoM(x1,y1,x2,y2)_10cm.mat')

L=60;  %sets the number of sample points 
D=0.295/2;    %sets the physical distance range

x = linspace(-D,D,60);
y = linspace(-D,D,60);


figure
pcolor(x,y,squeeze(abs(Hy_MoM(30,30,:,:)))'); %prop_4D(x1,y1,x2,y2)
axis tight
shading interp
title('Correlation function at 10 cm - MoM ')
set(gca,'YDir','normal')
set(gca,'FontSize',13) % Creates an axes and sets its FontSize to 18
colormap('jet'); % set the colorscheme
colorbar
xlabel('$$x_2 [m]$$','Interpreter','latex');
ylabel('$$y_2 [m]$$','Interpreter','latex');

%% CF plot from measured data at 10 cm
clear
clc

load('Hy_R3(x1,y1,x2,y2)_10cm.mat')

L=60;  %sets the number of sample points 
D=0.295/2;    %sets the physical distance range

x = linspace(-D,D,60);
y = linspace(-D,D,60);

figure
pcolor(x,y,squeeze(abs(Hy_R3(30,30,:,:)))'); % R3(x1,y1,x2,y2)
axis tight
shading interp
title('Correlation function at 10 cm - Measured ')
set(gca,'YDir','normal')
set(gca,'FontSize',13) % Creates an axes and sets its FontSize to 18
colormap('jet'); % set the colorscheme
colorbar
xlabel('$$x_2 [m]$$','Interpreter','latex');
ylabel('$$y_2 [m]$$','Interpreter','latex');

%% CF plot from measured data at 10 cm

clear
clc

load('Hy_R1(x1,y1,x2,y2)_1cm.mat')

L=60;  %sets the number of sample points 
D=0.295/2;    %sets the physical distance range

x = linspace(-D,D,60);
y = linspace(-D,D,60);

figure
pcolor(x,y,squeeze(abs(Hy_R1(30,30,:,:)))'); % R3(x1,y1,x2,y2)
axis tight
shading interp
title('Correlation function at 1 cm - Measured (Source)')
set(gca,'YDir','normal')
set(gca,'FontSize',13) % Creates an axes and sets its FontSize to 18
colormap('jet'); % set the colorscheme
colorbar
xlabel('$$x_2 [m]$$','Interpreter','latex');
ylabel('$$y_2 [m]$$','Interpreter','latex');