busy_rate.m

这是一个基于M/M/S类型的排队论仿真程序。M/M/S表示到达时间为负指数分布

%阻塞率的理论值计算
function[busy_r]=busy_rate(mean_arr,mean_serv,server_num)
C=server_num;
A=mean_serv/mean_arr;  
%到达速率与服务速率之比等于单位服务时间与单位到达时间之比
sum=0;
for k=0:C-1
    sum=sum+A^k/jiecheng(k);
end
%前提mean_serv/mean_arr<C
busy_r=(A^C)/((A^C)+jiecheng(C)*(1-A/C)*sum);