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）;??????%平均排隊時間