fig29.m

Space-Time Codes and MIMO Systems

% fig29.m, for the figure 2.9, pp.32
% capacity for ergodic channel by Monte-Carlo method,
% Date :  Dec.,2007
       clear
       j = sqrt(-1);
      Mr = [1 2 3 4];
      Mt = [1 2 3 4];
     LRT = length(Mr);
     SNR = 0 : 20;
    Lsnr = length(SNR);
 testnum = 2^(8);
 probnum = 10^4;
 out_num = 10^3;
 
         for ln = 1 : LRT;
             mR = Mr(ln);
             mT = Mt(ln);
             for lm = 1 : Lsnr;
                 snr = 10^(SNR(lm)/10);
                 for np = 1 : probnum,
                     for ll = 1 : testnum,
                            H_matrix = (randn(mT,mR) + j*randn(mT,mR))/sqrt(2);
                            M_Channl = H_matrix*H_matrix';
                            [Vect,D] = eig(M_Channl);
                            F_matrix = eye(size(D))+snr*D./mT;
                            G_mt(ll) = sum(log2(diag(F_matrix)));
                     end
                     F_mt(np) = sum(G_mt)/testnum;
                 end
                 H_mt = sort(F_mt);
                 D_mt(ln,lm) = H_mt(out_num); 
             end
         end
         
         bvar = C_mt(4,:);
         save capacity29.mat bvar -ascii; 
         clf,
         figure(1);
         plot(SNR,D_mt(1,:),'k-*');
         hold on;
         plot(SNR,D_mt(2,:),'k-p');
         plot(SNR,D_mt(3,:),'k-h');
         plot(SNR,D_mt(4,:),'k-v');
         xlabel('SNR(dB)');
         ylabel('10% outage capacity(bit/s/Hz)');
         legend('M=1','M=2','M=3','M=4');
         hold off