busy_rate.m

M/M/1 模型,主要用来求解单对单服务台的模型

%阻塞率的理论值计算

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/prod(1:k);

end

%前提mean_serv/mean_arr<C

busy_r=(A^C)/((A^C)+prod(1:C)*(1-A/C)*sum);