%总仿真时间
Total_time = 10;
%队列最大长度
N = 20;
%到达率与服务率
lambda = 10;
mu = 6;
%平均到达时间与平均服务时间
arr_mean = 1/lambda;
ser_mean = 1/mu;
%可能到达的最大顾客数(round:四舍五入求整数)
arr_num = round(Total_time*lambda*2);
%顾客事件表初始化
events = [];
%按负指数分布产生各顾客达到时间间隔
events(1,:) = exprnd(arr_mean,1,arr_num);
%各顾客的到达时刻等于时间间隔的累积和
events(1,:) = cumsum(events(1,:));
%按负指数分布产生各顾客服务时间
events(2,:) = exprnd(ser_mean,1,arr_num);
%计算仿真顾客个数,即到达时刻在仿真时间内的顾客数len_sim = sum(events(1,:)<= Total_time);
%*****************************************
% 计算第1 个顾客的信息
%*****************************************
%第1 个顾客进入系统后直接接受服务,无需等待
events(3,1) = 0;
%其离开时刻等于其到达时刻与服务时间之和
events(4,1) = events(1,1)+events(2,1);
%其肯定被系统接纳,此时系统内共有1 个顾客,故标志位%置1
events(5,1) = 1;
%其进入系统后,系统内已有成员序号为1
member = [1];
%*****************************************
% 计算第i 个顾客的信息
%*****************************************
for i = 2:arr_num
%如果第i 个顾客的到达时间超过了仿真时间,则跳出循环if events(1,i)>Total_time
break;
%如果第i 个顾客的到达时间未超过仿真时间,则计算在其%到达时刻系统中已有的顾客个数
else number = sum(events(4,member) > events(1,i));