tu1-17.m

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

SNR_dB=0:20;
%nTnR=[1 1; 2 1;2 2;4 2;4 4;8 8];
threshold=99/100; % that is 99%
for cases=1:6
    for s=1:length(SNR_dB)
        %n_T=nTnR(cases, 1);n_R=nTnR(cases, 2);
        sumC(s)=sumCfun(n_T, n_R, SNR_dB(s), threshold);
    end
    plot(SNR_dB, sumC);hold on; grid on;
    axis([0 20 0 40]);
    drawnow;
end
xlabel('信噪比/dB');
ylabel('容量/(bit/s/Hz)');
legend('T1R1','T2R1','T2R2','T4R2','T4R4','T8R8');