MMN的排队源程序matlab.doc

旅客登机前需要进行安检。首先接受安检人员的检查，有两名安检人员，旅客根据安检人员前的等待人数排队；然后通过一道安检门完成检查。如下图所示（时间单位分钟）。请编制仿真程序，估计顾客通过安检所需时间。添加函数newMM2andMM1，参数的解释也在源程序里说明了。则按题意输入meant=newMM2andMM1(1,1.5,2,0.2,1,300) ，得到结果为时间为meant =2.2418；各个人得到的结果是不太一样的，是正确的，因为程序模拟的随机分布的数，是符合实际的。源程序如下：function meant=newMM2andMM1(mean_arr,mean_serv,mean_serv1,mean_serv21,mean_serv22,peo_num) %mean_arr,到达的时间参数%mean_serv,服务台a的服务时间参数%mean_serv1,服务台b的服务时间参数%mean_serv21,mean_serv22,二级服务台的均匀分布的参数%peo_num,总服务人数nt=exprnd(mean_arr,1,peo_num);%各顾客到达时间间隔服从指数分布state_a=zeros(3,peo_num);%用一个三行矩阵表示a台每个顾客的状态 %三行依次为：到达时间间隔，服务时间，等待时间 state_b=zeros(3,peo_num); %用一个三行矩阵表示b台每个顾客的状态 %三行依次为：到达时间间隔，服务时间，等待时间 state_a(2,:)=exprnd(mean_serv,1,peo_num);%生成a台各顾客服务时间的矩阵state_b(2,:)=exprnd(mean_serv1,1,peo_num); %生成b台各顾客服务时间的矩阵state_a(3,1)=0; state_b(3,1)=0; a=1;%a台服务的人数b=1;%b台服务的人数arr_time=cumsum(nt); %到达时间由时间间隔变成连续时间state_b(1,1)=arr_time(1);state_a(1,1)=arr_time(2); %state(1,:)=arr_time;lea_time_a(1)=sum(state_a(:,1); %先计算前1名顾客的离开时间 lea_time_b(1)=sum(state_b(:,1); %先计算第2名顾客的离开时间 for i=3:peo_num iflea_time_a(a)lea_time_b(b)%第i个顾客到达，服务台满，等待时间为 %当时服务台最早离开的顾客的离开时间减去第i个顾客的到达时间 a=a+1; state_a(1,a)=arr_time(i); if state_a(1,a)=state_a(3,a-1)+state_a(2,a-1) state_a(3,a)=state_a(3,a-1)+state_a(2,a-1)-state_a(1,a-1); else state_a(3,a)=0; end lea_time_a(a)=sum(state_a(:,a); else b=b+1; state_b(1,b)=arr_time(i); if state_b(1,b)=state_b(3,b-1)+state_b(2,b-1) state_b(3,b)=state_b(3,b-1)+state_b(2,b-1)-state_b(1,b-1); else state_b(3,b)=0; end lea_time_b(b)=sum(state_b(:,b); endend%连接两个状态矩阵state=state_a(:,1:a),state_b(:,1:b);state(3,:)=lea_time_a(1:a),lea_time_b(1:b);%连接两个离开时间%g,m=min(lea_time_a);%h,n=min(lea_time_b);lea_time=lea_time_a,lea_time_b;%按离开时间的先后顺序排队guodu1=lea_time;guodu2=zeros(1,peo_num);for i=1:peo_num guodu2(i),j=min(guodu1); guodu1(j)=max(guodu1);end state2=zeros(3,peo_num); %用一个三行矩阵表示二级服务台每个顾客的状态 %三行依次为：到达时间间隔，服务时间，等待时间 state2(2,:)=unifrnd(mean_serv21,mean_serv22,1,peo_num);%产生二级服务台的服务时间分布 state2(1,:)=guodu2; for i=2:peo_num if state2(1,i)=state2(3,i-1)+state2(2,i-1)%需要等待，更新等待时间 state2(3,i)=state2(3,i-1)+state2(2,i-1)-state2(1,i); else state2(3,i)=0; end; end; %arr_time2=cumsum(state2(1,:); %state2(1,:)=arr_time2; lea_time2=sum(state2);%计算平均时间t1=0;t2=0