tu1-12-1.m

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

clear all;
n_R=1;n_T=1;
SNR_dB=15;
SNR=10.^(SNR_dB./10);
P=1;
sigma2=P./SNR;
for simutimes=1:10000
    H=(randn(n_R, n_T)+i*randn(n_R, n_T))./sqrt(2);
    lamda=eig(H*H');
    Pri=lamda.*P./n_T;%eq1-20
    C(simutimes)=sum(log2(1+Pri./sigma2))./n_T;%eq1-21
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(s)>=0
        j=j+1;
    end
    if C(s)>=0.2
        c=c+1;
    end
    if C(s)>=0.5
        k=k+1;
    end
    if C(s)>=0.8
        d=d+1;
    end
    if C(s)>=1.0
        l=l+1;
    end
    if C(s)>=1.3
        e=e+1;
    end
    if C(s)>=1.5
        m=m+1;
    end
    if C(s)>=1.7
        f=f+1;
    end
    if C(s)>=2.0
        a=a+1;
    end
    if C(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;
p0=[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,p0);
axis([0,4,0.90,1]);
hold on;
grid on;