Hw = (randn(n_r, n_t, nsample) + sqrt(-1) * randn(n_r, n_t, nsample)) / sqrt(2);
for i = 1:nsample
H = R_r_sqrt * Hw(:,:,i) * R_t_sqrt;
lambda = eig(H * H');
[mu pl] = WFL(lambda, SNR);
cwf(i) = sum(log2(1 + pl .* lambda));
cep(i) = sum(log2(1 + SNR/n_t * lambda));
end
Ce = mean(
1