图1-6.m

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

clear all;
for SNR=[1 sqrt(10) 10 sqrt(1000) 100 sqrt(100000) 1000]    %信噪比
    C=zeros(1,8);
    i=1;
    for NT=1:8   %发射天线的数目
        ronglian=zeros(1,200);
        for n=1:200                         %重复次数
            H=(randn(NT,1)+j*randn(NT,1))/sqrt(2);     %瑞利信道，零均值复高斯随机变量，方差为0.5
            G=H'*H;                                    %最大比合并
            ronglian(1,n)=log2(1+SNR*G/NT);               %求容量
        end
        C(1,i)=sum(ronglian)/200;   %求算术平均
        C1(1,i)=log2(1+SNR);
        i=i+1;
    end
    NT=1:8;
    plot(NT,C,'-*');
    hold on;
    plot(NT,C1);
    grid on;
    hold on;
    axis([1,8,0,11]);
end
xlabel('发射天线NT');
ylabel('容量/(bit/s/Hz)');
text(5,1.3,'SNR=0dB');
text(5,2.3,'SNR=5dB');
text(5,3.7,'SNR=10dB');
text(5,5.3,'SNR=15dB');
text(5,6.9,'SNR=20dB');
text(5,8.6,'SNR=25dB');
text(5,10.3,'SNR=30dB');