sumcfun.m

我自己在学习空时编码一书时编写的第一章的仿真程序

function sumC=sumCfun(n_T, n_R, SNR_dB, threshold)
SNR=10.^(SNR_dB./10);
simtimes=20000;
C_E(simtimes)=0;%init
for s=1:simtimes
    H=(randn(n_R,n_T)+j.*randn(n_R,n_T))./sqrt(2);
    C_E(s)=log2(det(eye(n_T,n_T)+SNR./n_T.*(H'*H)));
end
dx=0.01;
[c,xout] = hist(C_E, 0:dx:50);
pdf=c./simtimes./dx; %pdf fun
cdf=1-cumsum(pdf).*dx;
for s=1:length(xout)
    if cdf(s)<threshold
        sumC=xout(s);
        break;
    end
end