% Experimental Setup
% clear all;

init;

guideExpSetupTrans; % set transmitters positions
guideExpSetupRece;  % set receivers positions
guideSetGrid; % set a grid

seti.model = 'helmholtz2D';
seti.k = 250; % wave number

seti = setIncField(seti); % Evaluate incident fields on ROI
seti = setMeasKer(seti); % set up measurement kernels (such that measurement = k^2 * kernel * solution * voxelVolume)

figure(1); imagesc(real(reshape(seti.incField(:,1),[seti.nROI seti.nROI]))); axis xy; colorbar; % incident field in ROI of 1st transmitter
% set(gca,'FontSize',20); axis square; print(1,'-depsc','guideExpSetup1.eps');

figure(2); imagesc(real(reshape(seti.measKer(5,:),[seti.nROI seti.nROI]))); axis xy; colorbar; % measurement kernel on ROI of 5th receiver
% set(gca,'FontSize',20); axis square; print(2,'-depsc','guideExpSetup2.eps');