function R = calculate_covariance_matrix(r)

% r = SV * interference_data; % Interference and noise only
% r_zero_mean = r - diag(mean(r.'))*ones(size(r));
% R = (1/size(r_zero_mean,2))*(r_zero_mean*r_zero_mean');% + 5e10 * eye(size(r,1));    %R_smi_int_only = cov(r');

R = (1/size(r,2))*(r*r');% + 5e10 * eye(size(r,1));    %R_smi_int_only = cov(r');