%本程序時(shí)間單位是小時(shí)
clear;                
k=0;               %設(shè)定初始阻塞次數(shù)為0                    
s=3;               %信道數(shù)為3
a=1/100;           %服務(wù)時(shí)間   
num=1000;          %呼叫次數(shù)
c=zeros(1,s);      %信道初始化，全部為空閑
y=zeros(1,s);      %y為某一信道的剩余時(shí)間
g=rand(1,num);
pr=zeros(600,1);
for q=1:num        %隨機(jī)按指數(shù)分布產(chǎn)生服務(wù)時(shí)間
    ser(q)=-a*log(1-g(q));
end

for r=1:600        %r為信號(hào)到達(dá)率
    k=0;
    x=rand(1,num);
    for n=1:num    %產(chǎn)生信號(hào)時(shí)間間隔
       arrive(n)=-(1/r)*log(1-x(n));
    end
    for i=1:num           
        for j=1:s                %依次分析信道占用情況
            if c(j)==0
                c(j)=1;                               
                y(j)=ser(i);
                break                
            end
        end        
        if c(3)==1
            k=k+1;         %如果信道全部占用則阻塞次數(shù)加1       
        end        
        for t=1:s          %分析信號(hào)到達(dá)間隔期間的信道情況
            y(t)=y(t)-arrive(i);    
            if y(t)>0
                c(t)=1;
            else
                c(t)=0;
            end
        end
    end
                   
    pr(r)=k/num;           %輸出數(shù)據(jù)
    pr(r)
  
end
save simulate.txt pr -ASCII