yaojin.m

这个程序是理论和仿真呼叫阻塞率的对比!希望大家能多多交流学习!

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

for r=1:600        %r为信号到达率
    k=0;
    x=rand(1,num);
    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