sdma_blocking.asv

新一代无线通信核心技术OFDM的matlab仿真源程序！！！！！

lambda2=3;  
text_num=10000;%%%%the text num。
ichnum=0;%%%the num of seviceman in the main channel
qchnum=0;%%the num of  seviceman in  the complex channel
channelnum=5;  %%%the num  of  each channel
block_num=0;%%%用the num which is blocking
p1=0;    %%%to catch the probability of the blocking instance
ps=0.99;  %%%the probability of successful spatial allocation
%%%%%%%%%%%%%%FD  Blocking probability %%%%%%%%%%%%%%%%%%%%%
   for lambda1=1:8
       ichnum=0;
       qchnum=0;
       block_num=0;
       n=0;
    for i=1:text_num; 
       if ichnum<=channelnum;     
      %%%%%%%%%%%%%  
        if(ichnum==channelnum && qchnum==channelnum)
           block_num=block_num+1;
            ichnum=0;
           qchnum=0;

       end;
      %%%%%%%%%%%%%% 
        comeint =exprnd(lambda1);
        sevint=exprnd(lambda2);
        if  comeint>sevint;
          ichnum=ichnum+1;
          for i=1: ichnum-qchnum
              n=n+1;
              temp=rand;
               comeint = exprnd(lambda1);
               sevint=exprnd(lambda2);
              if(temp<ps && comeint>sevint)
                  qchnum=qchnum+1;
              end
              end
        end   
    else   block_num=block_num+1;
          ichnum=0;
          qchnum=0;
end;
end;
p1(lambda1)=block_num/(n+text_num);
end;
%%%%%%%%%%%%%%%%DL Blocking probability%%%%%%%%%%%%%%%
   for lambda1=1:8
       ichnum=0;
       qchnum=0;
       block_num=0;
       n=0;
    for i=1:text_num; 
        if ichnum<=channelnum;     
      %%%%%%%%%%%%%  
        if(ichnum==channelnum )
           block_num=block_num+1;
            ichnum=0;
           qchnum=0;

       end;
      %%%%%%%%%%%%%% 
        comeint =exprnd(lambda1);
        sevint=exprnd(lambda2);
        if  comeint>sevint;
          ichnum=ichnum+1;
              end
    
    else   block_num=block_num+1;
          ichnum=0;
          qchnum=0;
end;
end;
p2(lambda1)=block_num/(n+text_num);
end;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55555
SNRindB=1:8;
semilogy(SNRindB1,p1,'o');
hold
semilogy(SNRindB1,p2,'*');
xlabel('traffic');
ylabel('blocking probability');
title('blocking probability of FD ,DL  and the Theoretical');
legend( 'FD',' Rayleigh  'DL' ,' theoretical ');