图1-4.m

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

%%%%%绘制的是《空时编码技术》上第一章的图1-4,基于蒙特卡罗的思想，书上的求期望，
%%%%%在程序里面是通过重复求n多次，然后求算术平均的方法实现的(接收分集,最大比合并)
clear;
for SNR=[1 sqrt(10) 10 sqrt(1000) 100 sqrt(100000) 1000]    %信噪比
    C=zeros(1,11);
    i=1;
    for Nr=[1 3 6 8 10 12 16 24 32 48 65]   %发射天线的数目
        ronglian=zeros(1,200);
        for n=1:200                         %重复次数
            H=(randn(Nr,1)+j*randn(Nr,1))/sqrt(2);     %瑞利信道，零均值复高斯随机变量，方差为0.5
            G=H'*H;                                    %最大比合并
            ronglian(1,n)=log2(1+SNR*G);               %求容量
        end
        C(1,i)=sum(ronglian)/200;                      %求算术平均
        i=i+1;
    end
    Nr=[1 3 6 8 10 12 16 24 32 48 65];
    plot(Nr,C,'-*');
    grid on;
    hold on;
    axis([0,70,0,17]);
end
xlabel('接收天线Nr');
ylabel('容量/(bit/s/Hz)');
text(40,6.3,'SNR=0dB');
text(40,8.0,'SNR=5dB');
text(40,9.3,'SNR=10dB');
text(40,11.2,'SNR=15dB');
text(40,12.8,'SNR=20dB');
text(40,14.5,'SNR=25dB');
text(40,16.3,'SNR=30dB');