随机模拟实验报告.docx

随机模拟实验报告目的：了解随即现象，分析随即过程找出规律，解决实际问题实验原理：1. 随机变量及其分布随机事件：在一定条件下有可能发生的事件概率：随机事件发生的可能性的度量P(A), 0 W P(A) W 1.随机变量：在一定的范围内随机取值的变量，“X=ak”(k=l,2,...,n),或 X e[a,b]随机发生.随机变量的分布： “PCX=ak'') = pk, £久=1 则称ak : ala2 … an玖 pl为随机变量X的分布列,p2 … pn简称X的分布若已知P(”x w[q,/M”)= f卩⑴力则称P(x)为随机变量X的分布密度，简称X的分布实验过程：假设：市场服务超市有两个出口的收款台，两项服务：收款、装袋两名职丁•在出 口处工作有两种安排方案：1. 开一个出口，一人收款、—•人装袋；2. 开两个出口，每个人既收款乂装袋问商店经理应选择哪一种收款台的服务方案随机变量x={0, I, 2}表示每分钟到达超市收款台的人数，有分布列xk 0 I 2pk 0.4 0.3 0.3利用matlab软件的rand产生的随即数模拟十分钟内顾客到达收款台的状况模型假设：1.顾客的到达是随机的。

2.收款装袋的时间是相同的第一种方案中，收款与装袋同时进行，用的时间相同代码如下：elfline=zeros(l,60);%等待的顾客人数 wait=zeros(l ,60);%等待时间的累加 serve=zeros(l,60); %服务时间的累加 num=zeros(l,60); %LI到达顾客人数累加 st=l;%服务时间r=rand( 1,60);% 随机数for i=l:60;if r(i)<0.4n(i)=0;elseif 0.4<=r(i)&r(i)<0.7n(i)=l；else n(i)=2;endend讦 n(l)==2line(l)=l;serve(l)=st;endfor t=2:60if line(t-l )==0&(n(t)==0ln(t)== 1) line(t)=0;serve(t)=serve(t-1 )+n(t); wait(t)=wait(t-1 )+n(t); num(t)=num(t-1 )+n(t);elseline(t)=line(t-1 )+n(t)-st; wait(t)=wait(t-1 )+line(t); serve(t)=serve(t-1)+1; num(t)=num(t-l )+n(t);endend m=l:60; hold on plotCm^/b1); ploKrnJine/r1);hold off server_average=serve(60)/60 totla_num=num(60) wait_average=wait(60)/60 运行结果如下： server_average =0.8500totla_num =53 wait_average =1.5500结果评价：试验多次总有等待的高峰期。

即是总有等待时候代码如2•第二种方案中，收款与装袋分开进行，用的时间增加一倍 elfline2=zeros( 1,60);%等待的顾客人数line 1 =zeros(l ,60);%等待时I'可的累加 serve=zeros(l,61); %服务时间的累加 num=zeros(l,61); %L1到达顾客人数累加 st=l/2;%服务时间r=rand( 1,60);% 随机数for i=l:60if r(i)<0.4n(i)=0;elseif 0.4<=r (i)&r(i)<0.7n(i)=l；else n(i)=2;endendifn(l)==2line l(l)=st;line2( l)=st;elseif n(l)==1.0linel(l)=st;endfor t=2:60if n(t)==l讦 line 1 (t-l)>line2(t-l)Iine2(t)=line2(t-1 )+st;elseline 1 (t)=line l(t-I)+st;endelseif n(t)==2Iine2(t)=line2(t-1 )+st; linel(t)=linel(t-l)+st; elseIine2(t)=line2(t-1 )-st; linel(t)=linel(t-l)-st; endiflinel(t)<=0 linel(t)=O; endif line2(t)<=0 line2(t)=0; endif linel>0&line2>0wait(t)=wait(t-1 )+n(t); endserve(t)=serve(t-1 )+n(t); num(t)=num(t-1 )+n(t);endm=l:60;hold onplot(m,n/br); plot(mjinel/rr); plot(mjine2/y,);plot(m,line 1 +line2,g);hold off server_average=serve(60)/60 totla_num=nuni(60) 运行结果：3 • I2.5 - !server_average =1.0833totla_num =65结果评价：多次运行结果表明，第二种方案也常有等待髙峰时期，和第一•种方案在处理速 度上相差不人，再考虑工人的熟练程度和其它方面的情况，笔者认为还是第一种方案效率最 咼。