排队论matlab代码

function out=MMSkteam（skmu1mu2T）
%多服务台
%s——服务台个数
%k——最大顾客等待数
%T——时间终止点
%mu1——到达时间间隔服从指数分布
%mu2——服务时间服从指数分布
%事件表：
%   arrive_time——顾客到达事件
%   leave_time——顾客离开事件
%mintime——事件表中的最近事件
%current_time——当前时间
%L——队长
%tt——时间序列
%LL——队长序列
%c——顾客到达时间序列
%b——服务开始时间序列
%e——顾客离开时间序列
%a_count——到达顾客数
%b_count——服务顾客数
%e_count——损失顾客数

%初始化

arrive_time=exprnd（mu1）;
leave_time=[];
current_time=0;
L=0;
LL=[L];
tt=[current_time];
c=[];
b=[];
e=[];
a_count=0;
b_count=0;
e_count=0;
%循环
while min（[arrive_timeleave_time]）    current_time=min（[arrive_timeleave_time]）;
    tt=[ttcurrent_time];    %记录时间序列
    if current_time==arrive_time          %顾客到达子过程
        arrive_time=arrive_time+exprnd（mu1）;  % 刷新顾客到达事件
        a_count=a_count+1; %累加到达顾客数
        if  L