capacity_22.m

计算信道容量的matlab程序

%close all;
%clc;
hold on
Nt=2;
Nr=2;
[alph11,alph12,alph21,alph22]=correlation_cofficiency();
a11=alph11(10);
a12=alph11(20);
a21=alph12(10);
a22=alph12(20);
a31=alph21(10);
a32=alph21(20);
a41=alph22(10);
a42=alph22(20);

    h1=sqrt(1/2)*(a11+j*a12); 
    h2=sqrt(1/2)*(a21+j*a22);
    h3=sqrt(1/2)*(a31+j*a32);
    h4=sqrt(1/2)*(a41+j*a42);
    H=[h1,h2;h3,h4];
    snr_dB=0:1:20;
    snr=10.^(snr_dB/10);
    for i=1:length(snr)
        C2(i)=log2(det(eye(Nr)+snr(i)*H*H'));
    end

plot(snr_dB,C2,'');
xlabel('SNR(db)');
ylabel('信道容量(bits/s/Hz)');