sdma_sigle

来自「无线多址技术SDMA的MATLAB仿真及其相关知识讲解！！！！」· 代码 · 共 49 行

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 A=1:8
       lambda1=20; 
       lambda2=A/lambda;
       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(1/lambda2);
        if  comeint>sevint;
          ichnum=ichnum+1;
          for i=1: ichnum-qchnum
              n=n+1;
              temp=rand;
              comeint =exprnd(lambda1);
              sevint=exprnd(1/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;
SNRindB=1:8;
semilogy(SNRindB,p1,'o');