MMN的排队源程序matlab

旅客登机前需要进行安检。首先接受安检人员的检查，有两名安检人员，旅客根据安检人员前的等待人数排队；然后通过一道安检门完成检查。如下图所示（时间单位分钟）。请编制仿真程序，估计顾客通过安检所需时间。添加函数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名顾客的离开时间fori=3:peo_numiflea_time_a(a)lea_time_b(b)%第i个顾客到达，服务台满，等待时间为%当时服务台最早离开的顾客的离开时间减去第i个顾客的到达时间a=a+1;state_a(1,a)=arr_time(i);ifstate_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);elsestate_a(3,a)=0;endlea_time_a(a)=sum(state_a(:,a));elseb=b+1;state_b(1,b)=arr_time(i);ifstate_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);elsestate_b(3,b)=0;endlea_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);fori=1:peo_num[guodu2(i),j]=min(guodu1);guodu1(j)=max(guodu1);endstate2=zeros(3,peo_num);%用一个三行矩阵表示二级服务台每个顾客的状态%三行依次为：到达时间间隔，服务时间，等待时间state2(2,:)=unifrnd(mean_serv21,mean_serv22,1,peo_num);%产生二级服务台的服务时间分布state2(1,:)=guodu2;fori=2:peo_numifstate2(1,i)=state2(3,i-1)+state2(2,i-1)%需要等待，更新等待时间state2(3,i)=state2(3,i-1)+state2(2,i-1)-state2(1,i);elsestate2(3,i)=0;end;end;%arr_time2=cumsum(state2(1,:));%state2(1,:)=arr_time2;lea_time2=sum(state2);%计算平均时间t1=0;t2=0;fori=1:peo_numt1=t1+state(3,i)-state(1,i);t2=t2+lea_time2(i)-state2(1,i);endmeant=(t1+t2)/peo_num;在此感谢一篇微博的主人分享