tu1-12.m

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

SNR_dB=15;
SNR=10.^(SNR_dB./10);
n=1;
simtimes=10000;
C_E(simtimes)=0;%init
for s=1:simtimes
    H=(randn(n,n)+i.*randn(n,n))./sqrt(2);
    C_E(s)=1./n.*log2(det(eye(n,n)+SNR./n.*(H'*H)));
end
j=0;k=0;l=0;m=0;a=0;b=0;c=0;d=0;e=0;f=0;
for s=1:10000
    if C_E(s)>=0
        j=j+1;
    end
    if C_E(s)>=0.2
        c=c+1;
    end
    if C_E(s)>=0.5
        k=k+1;
    end
    if C_E(s)>=0.8
        d=d+1;
    end
    if C_E(s)>=1.0
        l=l+1;
    end
    if C_E(s)>=1.3
        e=e+1;
    end
    if C_E(s)>=1.5
        m=m+1;
    end
    if C_E(s)>=1.7
        f=f+1;
    end
    if C_E(s)>=2.0
        a=a+1;
    end
    if C_E(s)>=2.2
        b=b+1;
    end
end
p1=j/10000;p2=c/10000;p3=k/10000;p4=d/10000;p5=l/10000;p6=e/10000;p7=m/10000;p8=f/10000;p9=a/10000;p10=b/10000;
p=[p1 p2 p3 p4 p5 p6 p7 p8 p9 p10];
r=[0 0.2 0.5 0.8 1.0 1.3 1.5 1.7 2.0 2.2];
plot(r,p);
axis([0,4,0.90,1]);
hold on;
grid on;