阻塞率matlab仿真

采用MATLAB编程，这个程序是理论和仿真呼叫阻塞率的对比

                   %本程序时间单位是小时
clear;                
k=0;               %设定初始阻塞次数为0                    
s=3;               %信道数为3
a=1/100;           %服务时间   
num=1000;          %呼叫次数
c=zeros（1s）;      %信道初始化，全部为空闲
y=zeros（1s）;      %y为某一信道的剩余时间
g=rand（1num）;
pr=zeros（6001）;
for q=1:num        %随机按指数分布产生服务时间
    ser（q）=-a*log（1-g（q））;
end

for r=1:600        %r为信号到达率
    k=0;
    x=rand（1num）;
    for n=1:num    %产生信号时间间隔
       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;         %如果信道全部占用则阻塞次数加1       
        end        
        for t=1:s          %分析信号到达间隔期间的信道情况
            y（t）=y（t）-arrive（i）;    
            if y（t）>0
                c（t）=1;
            else
                c（t）=0;
            end
        end
    end
                   
    pr（r）=k/num;           %输出数据
    pr（r）
  
end
save simulate.txt pr -ASCII 



  
    
    


                        
        
            
                
                
    

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1418  2006-09-03 15:56  yaojin.m

----------- ---------  ---------- -----  ----

                 1636                    2