基于matlab的网络排队仿真,Matlab的mmn型排队仿真源程序

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

if  lea_time_a(a)

%第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));

end

end

%连接两个状态矩阵

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;

for i=1:peo_num

t1=t1+state(3,i)-state(1,i);

t2=t2+lea_time2(i)-state2(1,i);

end

meant=(t1+t2)/peo_num;