M_M_1.m

clearclc;
Total=100000;              %顾客的总数；
Lamda=10;               %顾客到达率；
Mu=4;                 %服务率；
Arrive_time=zeros（1Total）; %顾客到来时间变量加快速度；
Leave_time=zeros（1Total）;
Arrive_num=zeros（1Total）; %顾客到达人数变量
Leave_num=zeros（1Total）;
Arrive_interval=Createpoiss（Lamda1Total）; %顾客到达时间间隔，服从指数分布
Serve_interval=Createpoiss（Mu1Total）; %服务台服务时间，服从指数分布

%各顾客到来时间；
Arrive_time（1）=Arrive_interval（1）;%初始化
Arrive_num（1）=1;
for i=2:Total
    Arrive_time（i）=Arrive_time（i-1）+Arrive_interval（i）;%下一个顾客到达时间等于前一个顾客到达的时间加上到达间隔时间
    Arrive_num（i）=i;
end
%Arrive_time

%各顾客离开时间；
 Leave_time（1）=Arrive_time（1）+Serve_interval（1）;
 Leave_num（1）=1;
for i=2:Total
    if Leave_time（i-1）        Leave_time（i）=Arrive_time（i）+Serve_interval（i）;%则下一个顾客的离开时间等于他到达的时间加上服务间隔时间
    else         
         Leave_time（i）=Leave_time（i-1）+Serve_interval（i）;%如果前一个顾客在下一个顾客到来之时或之后才离开
    end                                                  %则下一个顾客的离开时间等于前一个顾客离开时间加上服务间隔时间
    Leave_num（i）=i; %离开顾客数增加一名
end
%Leave_time

Wait_time=Leave_time-Arrive_time;      %各个顾客在系统中的等待时间
Queue_time=Wait_time-Serve_interval;  %各个顾客的排队时间
Wait_time_avg=mean（Wait_time）;        %平均系统等待时间
Queue_time_avg=mean（Queue_time）;      %平均排队时间