sdma_single.asv

无线多址技术SDMA的MATLAB仿真及其相关知识讲解！！！！

text_num=10000;%%%%the text num。
channelnum1=5;  %%%the num  of   channel
channelnum2=10;
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
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  for A=0.5:7.5;
     lambda1=40; 
       lambda2=A/ lambda1;
       ichnum=1;
       qchnum=0;
       block_num=0;
       n=0;
    for i=2:text_num; 
     %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  
     
         if ichnum<=channelnum;     
        comeint(i)=exprnd(1/lambda2);
        sevint(i)=exprnd(lambda1);
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
        for j=i:(-1):(i-ichnum-qchnum);
         ltime(j)=sevint(j)-comeint(i); 
         sevint(j)=ltime(j);
     %%%%%%%%%%%%%%%%%%%%%%%%%%%%
   if  ltime(j)<=0;
           if ichnum>qchnum
           ichnum=ichnum-1;       
           else if ichnum==qchnum
               qchnum=qchnum-1;
                  end
           end
         end 
                if qchnum<0;
                  qchnum=0;             
                  if ichnum<0;
                  ichnum=0;
         end   
                end
        end
     
      %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
       if  all(ltime(j))>0;
              temp=rand; 
           if ichnum>qchnum && temp<ps                      
                  qchnum=qchnum+1;
           else    ichnum=ichnum+1;
            end
        end  
        end;
   
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%555555
        
        if(ichnum==channelnum && qchnum==channelnum)
           block_num=block_num+1;
            ichnum=1;
           qchnum=0;  
    
      else   if  ichnum>channelnum;  
            block_num=block_num+1;
          ichnum=1;
          qchnum=0;
      
          end
    end
           end 
p1(A+0.5)=block_num/text_num
  end
  SNRindB=1:8;
semilogy(SNRindB,p1,'g');
xlabel('traffic');
ylabel('blocking probability');
title('blocking probability of FD  and the Theoretical');