实现MMN排队算法的源程序

算法提出问题旅客登机前需要进行安检首先接受安检人员的检查，有两名安检人员，旅客根据安检人员前的等待人数排队；然后通过一道安检门完成检查估计顾客通过安检所需时间添加函数newMM2andMMl,参数的解释也在源程序里说明了则按题意输入meant=newMM2andMMl1,

1.5,2,

0.2,1,300,得到结果为时间为meant=

2.2418；各个人得到的结果是不太一样的，是正确的，因为程序模拟的随机分布的数,是符合实际的源程序如下function[meant]=newMM2andMMlmean_arr,mean_serv,mean_servl,mean_serv21,mean_serv22,peo_num%mean_arr,到达的时间参数%mean_serv,服务台a的服务时间参数%mean_servl,服务台b的服务时间参数%mean_serv21,mean_serv22,二级服务台的均匀分布的参数%peo_nu叫总服务人数nt=exprndmean_arr,1,peo_num;%各顾客到达时间间隔服从指数分布state_a=zeros3,peo num;%用一个三行矩阵表示a台每个顾客的状态%三行依次为到达时间间隔，服务时间，等待时间state_b=zeros3,peo_num;%用一个三行矩阵表示b台每个顾客的状态%三行依次为到达时间间隔，服务时间，等待时间state a2,:=exprndmean serv,1,peo num;%生成a台各顾客服务时间的矩阵state_b2,:=exprndmean_servl,1,peo_num;%生成b台各顾客服务时间的矩阵state a3,1=0;state_b3,1=0;a=l1%a台服务的人数b=l;%b台服务的人数arr_t ime=cumsum nt;%到达时间由时间间隔变成连续时间state_b1,l=arr_time1;state_al,l=arr_time2;%statel,:=arr_time;lea_time_a1=sum state_a:,1;%先计算前1名顾客的离开时间lea_time_b1=sumstate_b:,1;%先计算第2名顾客的离开时间for i=3:peo_numif lea_time_aalea_time_bb%第1个顾客到达，服务台满，等待时间为%当时服务台最早离开的顾客的离开时间减去第i个顾客的到达时间a=a+1;state_al,a=arr_time i;if state_al,a=state_a3,a-1+state_a2,a-1state_a3,a=state_a3,a-l+state_a2,a-l-state_al,a-l;elsestate_a3,a=0;endlea_time_aa=sumstate_a:,a;elseb=b+l;state_b1,b=arr_timei;if state_b1,b=state_b3,b-1+state_b2,bTstate b3,b=state_b3,b-l+state_b2,b-1-state_b1,b-l;elsestate_b3,b=0;end lea_time_bb=sumstate_b:,b;endend%连接两个状态矩阵state=[state_a:,1:a,state_b:,l:b];state3,:=[lea time_al:a,lea timeb1:b];%连接两个离开时间%[g,m]=minlea_time_a;%[h,n]=minlea_time_b;1ea time=[1ea time_a,lea timeb];%按离开时间的先后顺序排队guodul=lea_time;guodu2=zeros1,peo num;for i=l:peo_num[guodu2i,j]=minguodul;guodulj=maxguodul;end state2=zeros3,peo num;%用一个三行矩阵表示二级服务台每个顾客的状态%三行依次为到达时间间隔，服务时间，等待时间state22,:=unifrndmean_serv21,mean_serv22,1,peo num;%产生二级服务台的服务时间分布state2l,:=guodu2;for i=2:peo_numif state21,i〈二state23,i-1+state22,i-1%需要等待，更新等待时间state23,i=state23,i-1+state22,i-l-state21,ielsestate23,i=0;%arr_time2=cumsunistate21,:;%state2l,:=arr_time2;Iea_time2=sumstate2;%计算平均时间tl=0;t2=0;for i=l:peo numtl=tl+state3,i-state1,i;t2=t2+lea_time2i-state21,i;endmeant=tl+t2/peo_num;在此感谢一篇微博的主人分享。