tu1-13.m

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

SNR_dB=0:5:30;
SNR=10.^(SNR_dB./10);
n=8;
for k=1:length(SNR)
    simtimes=10000;
    C_E(simtimes)=0;%init
    for s=1:simtimes
        H=(randn(n,n)+j.*randn(n,n))./sqrt(2);
        C_E(s)=1./n.*log2(det(eye(n,n)+SNR(k)./n.*(H'*H)));
    end
    dx=0.01;
    [c,xout] = hist(C_E, 0:dx:10);
    pdf=c./simtimes./dx; %pdf fun
    %plot(xout,pdf);
    %sum(pdf).*dx %verify pdf fun
    plot(xout, 1-cumsum(pdf).*dx);hold on;
    axis([0 9 0.9 1]);grid on;
end
xlabel('容量n/(bit/s/Hz)');
ylabel('大于横坐标的概率 [容量n>横坐标的概率]');
text(0.8,0.975,'0dB');
text(1.5,0.962,'5dB');
text(2.4,0.975,'10dB');
text(3.8,0.962,'15dB');
text(5.1,0.975,'20dB');
text(6.7,0.962,'25dB');
text(8.1,0.975,'30dB');