Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- %%
- 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);
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement