fig26.m

Space-Time Codes and MIMO Systems

% fig26.m, for the figure 2.6, pp.29
% capacity for ergodic channel by Monte-Carlo method,
% the channel objection the rayleigh distribution,
% Date :  Dec.,2007
       clear; 
       j = sqrt(-1);
      mR = 2;
      mT = 2;
     SNR = 10;
 testnum = 2^(15);
     snr = 10^(SNR/10);
 thresld = 0.1;
  expect = 0.35;
    Rate = 15;
 Rateval = 0:14;   
 Vateval = 0:0.01:14;   
                 for ll = 1 : testnum,
                     matrix_A = rand(mR,mT);
                     H_matrix = pdf('rayl',expect,matrix_A);
                     M_Channl = H_matrix*H_matrix';
                     [Vect,D] = eig(M_Channl);
                     F_matrix = eye(size(D))+snr*D./mT;
                     G_mt = sum(log2(diag(F_matrix)));
                     out_p(ll) = G_mt;  
                 end
                 [A_vect B_vect]= hist(out_p,Rate);
                 prob_ccdf = cumsum(A_vect./testnum);
                 prob_ccdg = interp1(Rateval,prob_ccdf,Vateval);   
                 figure(1);
                 plot(Vateval,prob_ccdg);
                 hold on;
                 grid,
                 title('the channel is Rayleigh distribution')
                 xlabel('Rate(bit/s/Hz)');
                 ylabel('CDF');
                 hold off;