图1-7.m

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

clear;
for SNR=[1 sqrt(10) 10 sqrt(1000) 100 sqrt(100000) 1000]
    C=zeros(1,11);
    i=1;
    for n=[1 3 5 7 10 12 16 24 32 48 65]
        C(1,i)=n*(log2(SNR)-1);
        i=i+1;
    end 
    n=[1 3 5 7 10 12 16 24 32 48 65];
    plot(n,C,'-*');
    axis([0,70,0,600]);
    grid on;
    hold on;
end
xlabel('天线数量n');
ylabel('容量/(bit/s/Hz)');