来自「模拟随机过程求解超市实际问题。很具有代表性。matlab」· 代码 · 共 39 行

clear;
clc;
L=zeros(1,31);% L 等待的顾客人数，
T1=zeros(1,31);%T1总合等待时间，
T2=zeros(1,31); %T2总合服务时间,
L1=zeros(1,31); %L1总合到达顾客人数。
t=1;
tau=1;
x=0:30;
r=rand(1,30);%随机数
%随机模拟
for i=1:30; 
    t=t+1;
    if 0<=r(i) & r(i)<0.4
        n=0;
       elseif 0.4<=r(i) & r(i)<0.7
        n=1;
       else n=2;
    end;
%排队分析
if L(t-1)==0 & n==0
        L(t)=L(t-1);T1(t)=T1(t-1);%模型
        T2(t)=T2(t-1);L1(t)=L1(t-1);
        else
        L(t)=L(t-1)+n-1;T1(t)=T1(t-1)+L(t); 
        T2(t)=T2(t-1)+tau; L1(t)=L1(t-1)+n; 
    end;
end;
r=[0 r];
a=[x',r',L',L1',T1',T2']

eL=T2(end)/tau  %已被服务的人数
L2=(find(L1>eL))
L3=sum(L(L2))%未被服务的顾客等待时间总和
g1=(T1(end)-L3)/eL %平均等待时间
g2=g1+tau  %平均逗留时间
g3=eL/30  %平均每分钟服务的顾客人数