%%
matrixR = zeros(1000,1000);
matrixR = polyval([-62.09 -22.3 33.51],r); %% evaluate parabolic fit of log(|S|^2) for a grid of radii
matrixR = sqrt(exp(matrixR)); %% exponentiate and square root to get |S|

matrixCR = zeros(1000,1000);
matrixCR = sqrt(exp(F))./matrixR; %% matrix F is: F = log(abs(fftshift(fft2(m15gn))).^2);

matrixPhi = exp(F)./(exp(F) + exp(22.74)); %% sqrt(22.74) is the white noise component of log(|S|^2)

%%
matrixUf = matrixCR.*matrixPhi;

matrixU = ifftshift(ifft2(matrixUf));

matrixUabs = abs(matrixU);
figure();
imagesc(matrixUabs);