随机过程模拟.pptVIP

回忆：;二、随机型动态系统仿真;;目标：选择什么方案？

依据是什么？

主体：顾客和商家

顾客的满意度——排队时间尽可能的短，效劳时间短。

商家希望付出的人力少，得到更好的顾客满意度。

排队问题——有随机因素参加——随机模拟;研究开一个收款台的情况

假设：

1.顾客的到达是随机的。

2.收款装袋的时间是相同的。

3.收款与装袋同时进行

时间步长法：

把整个仿真过程分为许多相等的时间间隔

每个间隔为一个时间单位——时间步长。

在每个时间步长内模拟系统的动态。

仿真时钟：用以控制时间步进的过程

〔每一次步进一个步长〕;参量、变量：

n(t)到达顾客人数，L(t)队列长,τ效劳时间

T1总合排队时间,T2总合效劳时间，Δt时间步长

模型：平衡关系

当L(t)=0且n(t)=0时

L(t+Δt)=L(t);

T1(t+Δt)=T1(t);T2(t+Δt)=T2(t);

否那么

L(t+Δt)=L(t)+n(t)-1;

T1(t+Δt)=T1(t)+L(t+Δt)Δt;

T2(t+Δt)=T2(t)+τ

取时间步长Δt=1，收款或装袋的时间τ=1。

仿真30分钟内收款台处的排队情况;

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);%随机数;%随机模拟——顾客到来情况

fori=1:30;

t=t+1;

if0=r(i)r(i)0.4

n=0;

elseif0.4=r(i)r(i)0.7

n=1;

elsen=2;

end;

;%排队分析

ifL(t-1)==0n==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)Δt;

T2(t)=T2(t-1)+tau;L1(t)=L1(t-1)+n;

end;

end;;r=[0r];

a=[x,r,L,L1,T1,T2]

eL=T2(end)/tau%已被效劳的人数

L2=(find(L1eL))

L3=sum(L(L2))%未被效劳的顾客等待时间总和

g1=(T1(end)-L3)/eL%平均等待时间

g2=g1+tau%平均逗留时间

g3=eL/30%平均每分钟效劳的顾客人数;a=[x‘,r’,L‘,L1’,T1‘,T2’]

000000

1.00000.01530000

2.00000.74681.00002.00001.00001.0000

3.00000.44511.00003.00002.00002.0000

4.00000.93182.00005.00004.00003.0000

5.00000.46602.00006.00006.00004.0000

6.00000.41862.00007.00008.00005.0000

7.00000.84623.00009.000011.00006.0000

8.00000.52523.000010.000014.00007.0000

9.00000.20262.000010.000016.00008.0000

10.00000.67212.000011.000018.00009.0000

11.00000.83813.0000