clear
clc
%%?MMK∞和單隊(duì)列多通道—多服務(wù)臺(tái)混聯(lián)式（即在每一個(gè)通道上，每個(gè)區(qū)域的服務(wù)臺(tái)保持串聯(lián)，而同一個(gè)區(qū)域內(nèi)有多個(gè)服務(wù)臺(tái)并聯(lián)）
N=100;%乘客總數(shù)
M=4;%有4個(gè)串聯(lián)的區(qū)域
u=0.45;%大約有45%的乘客會(huì)選擇PC通道（Pre-Check）
lambda0=4.65;%乘客的到達(dá)率（每分鐘平均到達(dá)4.65人）
%%
%乘客的平均到達(dá)時(shí)間
l0=1/lambda0;
%常規(guī)區(qū)域的平均服務(wù)時(shí)間，時(shí)間單位均為分鐘
lambda1=11.21/60;%第一個(gè)區(qū)域有兩個(gè)服務(wù)臺(tái)但是這里為了簡(jiǎn)化模型，假設(shè)兩個(gè)服務(wù)臺(tái)的服務(wù)時(shí)間是完全一樣的
lambda2=13.2/60;
lambda3=27/60;
lambda4=13.2/60;

%PC通道的平均服務(wù)時(shí)間，由于乘客不用移除衣物和一些個(gè)人物品，所以將lambda2和lambda4的時(shí)間減少25%
lambda1_PC=11.21/60;%第一個(gè)區(qū)域有兩個(gè)服務(wù)臺(tái)但是這里為了簡(jiǎn)化模型，假設(shè)兩個(gè)服務(wù)臺(tái)的服務(wù)時(shí)間是完全一樣的
lambda2_PC=13.2/60*0.75;
lambda3_PC=27/60;
lambda4_PC=13.2/60*0.75;
%%
%按負(fù)指數(shù)分布產(chǎn)生各乘客到達(dá)的時(shí)間間隔
arr_interval=?exprnd（l01N）;
%各乘客的到達(dá)時(shí)刻等于時(shí)間間隔的累積和
arr_moment=?cumsum（arr_interval）;
%產(chǎn)生45個(gè)隨機(jī)數(shù)，代表選擇PC的乘客的編號(hào)
r=randperm（10045）;
r=sort（r）;
%按負(fù)指數(shù)分布產(chǎn)生每個(gè)乘客在常規(guī)通道各區(qū)域的服務(wù)時(shí)間serve（ij）表示第i個(gè)乘客在第j個(gè)區(qū)域的服務(wù)時(shí)間//此處第一個(gè)區(qū)域有兩個(gè)服務(wù)臺(tái)但是這里為了簡(jiǎn)化模型，假設(shè)兩個(gè)服務(wù)臺(tái)的服務(wù)時(shí)間是完全一樣的
serve（:1）=exprnd（lambda11N）‘;
serve（:2）=exprnd（lambda21N）‘;
serve（:3）=exprnd（lambda31N）‘;
serve（:4）=exprnd（lambda41N）‘;
%按負(fù)指數(shù)分布產(chǎn)生每個(gè)乘客在PC通道各區(qū)域的服務(wù)時(shí)間serve（ij）表示第i個(gè)乘客在第j個(gè)區(qū)域的服務(wù)時(shí)間//此處第一個(gè)區(qū)域有兩個(gè)服務(wù)臺(tái)但是這里為了簡(jiǎn)化模型，假設(shè)兩個(gè)服務(wù)臺(tái)的服務(wù)時(shí)間是完全一樣的
PC_serve（:1）=exprnd（lambda1_PC1N）‘;
PC_serve（:2）=exprnd（lambda2_PC1N）‘;
PC_serve（:3）=exprnd（lambda3_PC1N）‘;
PC_serve（:4）=exprnd（lambda4_PC1N）‘;
%%
%第i個(gè)乘客在常規(guī)通道第j個(gè)區(qū)域花費(fèi)的時(shí)間為A1（ij）（i=1:100j=1:4花費(fèi)的時(shí)間等于服務(wù)的時(shí)間加上等待的時(shí)間）
%第i個(gè)乘客離開(kāi)常規(guī)通道第j個(gè)區(qū)域的時(shí)刻為B1（ij）

%在SLM的基礎(chǔ)上做點(diǎn)修改簡(jiǎn)化下代碼//在這一節(jié)編程的時(shí)候就直接用“C”矩陣了

%第i個(gè)乘客在PC通道第j個(gè)區(qū)域花費(fèi)的時(shí)間為A2（ij）（i=1:100j=1:4花費(fèi)的時(shí)間等于服務(wù)的時(shí)間加上等待的時(shí)間）
%第i個(gè)乘客離開(kāi)PC通道第j個(gè)區(qū)域的時(shí)刻為B2（ij）
%在假設(shè)不計(jì)區(qū)域之間花費(fèi)的時(shí)間的條件下有這樣一個(gè)關(guān)系:到達(dá)第j個(gè)區(qū)域的時(shí)刻等于離開(kāi)第j-1個(gè)區(qū)域的時(shí)刻
%a1，a2分別用來(lái)記錄乘客離開(kāi)常規(guī)通道和PC通道內(nèi)第一區(qū)域的兩個(gè)服務(wù)臺(tái)的時(shí)刻為了簡(jiǎn)化，假設(shè)乘客遇到兩個(gè)服務(wù)臺(tái)時(shí)（即使兩個(gè)都是空的）會(huì)選擇乘客先行離開(kāi)的那個(gè)服務(wù)臺(tái)
a1=zeros（12）;
a2=zeros（12）;
A1=zeros（NM）;
A2=zeros（NM）;
C1=zeros（N5）;
C2=zeros（N5）;
C1（:1）=arr_moment‘;
C2（:1）=arr_moment‘;
if?r（1）~=1
?????A1（1:）=serve（1:）;
?????C1（12）=A1（11）+C1（11）;
?????a1（1）=C1（12）;
????for?j=2:4
???????C1（1j+1）=C1（1j）+A1（1j）;
????end
else?
????A2（1:）=PC_serve（1:）;
????C2（12）=A2（11）+C2（11）;
????a2（1）=C2（12）;
????for?j=2:4
???????C2（1j+1）=C2（1j）+A2（1j）;
????end
end


for?i=2:N
????f