图1-5.m

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

%发射分集,选择性合并
clear;
for SNR=[1 sqrt(10) 10 sqrt(1000) 100 sqrt(100000) 1000]    %信噪比
    C=zeros(1,8);
    i=1;
    for Nr=[1 3 5 7 9 24 48 64]   %发射天线的数目
        ronglian=zeros(1,200);
        for n=1:200                         %重复次数
            H=(randn(Nr,1)+j*randn(Nr,1))/sqrt(2);     %瑞利信道，零均值复高斯随机变量，方差为0.5
            G=(abs(H)).^2;                             %选择性合并
            K=max(G);                                  %选择性合并
            ronglian(1,n)=log2(1+SNR*K);               %求容量
        end
        C(1,i)=sum(ronglian)/200;                      %求算术平均
        i=i+1;
    end
    Nr=[1 3 5 7 9 24 48 64];
    plot(Nr,C,'-*');
    grid on;
    hold on;
    axis([0,70,0,14]);
end
xlabel('接收天线Nr');
ylabel('容量/(bit/s/Hz)');
text(40,2.7,'SNR=0dB');
text(40,4.4,'SNR=5dB');
text(40,6.5,'SNR=10dB');
text(40,7.4,'SNR=15dB');
text(40,9.0,'SNR=20dB');
text(40,10.7,'SNR=25dB');
text(40,12.5,'SNR=30dB');