tu1-22.m

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

clear;
alfa_r=[1,5,10,20,40]./180.*pi;
d_bi_lamda=0:0.1:5;
SNR_dB=20;
n_R=4;
SNR=10.^(SNR_dB./10);
for s=1:length(alfa_r)
    for ss=1:length(d_bi_lamda)
        SitaR=sitaRfun(d_bi_lamda(ss),n_R, alfa_r(s));
        for sim=1:10000
            H=(randn(n_R,n_R)+j*randn(n_R,n_R))./sqrt(2);
            C(sim)=log2(det(eye(n_R,n_R)+SNR./n_R.*SitaR'*H*H'));%1-110
        end
        Cmean(s,ss)=real(mean(C));
        %[s, ss] %display progress
    end
end
plot(d_bi_lamda, Cmean');grid on;axis([0 5 8 24]);
xlabel('d/λ');
ylabel('平均容量/(bit/s/Hz)');
text(3.1,11,'αr=1°');
text(2.5,15.5,'αr=5°');
text(2.3,19,'αr=10°');
text(1.5,21,'αr=20°');
text(1.1,22.5,'αr=40°');