实验课(十一).pptVIP

实验课（十一） 随机模拟 1. 随机变量的模拟 unifrnd(0,1,1,6) %产生1×6的[0,1]均匀随机数 unifrne(0,10,1,6) %产生1×6的[0,10]均匀随机数 unidrnd(10,1,8) %产生1×8的≤10的离散均匀随机数 normrnd(1,1/3,1,10) %产生1×10的N(1,1/3)正态随机数 normrnd(0,1,1,10) %产生1×10的N(0,1)正态随机数 exprnd(2,1,10) %产生1×10的λ=1/2的指数随机数 2. 随机现象的模拟 问题1. 顾客到达收款台的规律是: 每分钟内40%的时间没人来，30%的时间有1个人来，30%的时间有两个人来。 模拟30 分钟内顾客到达收款台的情况。 r=rand(1,30); n1=(r=0.4);n2=(r=0.7); n=n1+n2 顾客到达时间（n中非零元的位置） m1=find(n) 顾客到达情况（n中非零元的数值） m2=n(find(n)) 市场服务 超市有两个出口的收款台， 两项服务：收款、装袋。 两名职工在收款台处工作。 有两种安排方案： 开一个出口，一人收款、一人装袋； 开两个出口，每个人既收款又装袋。 问商店经理应选择哪一种收款台的服务方案 假设： 1. 顾客的到达是随机的。 2. 收款装袋的时间是相同的。 3. 第一种方案中，收款与装袋同时进行 建模分析开设两个收款台的方案，根据模拟结果评价两个方案。 (1). 假设对不同的顾客一位服务员收款并装袋所需的服务时间服从正态分布 N(2，2/3)， 两位服务员分别收款装袋所需的服务时间服从正态分布 N(1，1/3). (2). 假设顾客到达的平均时间间隔是0.9分钟,顾客在收款台前排成一队。? 模型I（时间步长法） 参量、变量： n(t) 到达人数，Δt 时间步长， τ 服务时间， L(t) 队列长, T1(t) 排队时间, 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 模型II 事件表法： 参量，变量： t(i): 第i位顾客到达时间， t2(i):第i位顾客受到的服务时间(随机变量)， T(i): 第i位顾客离去时间. 将第i位顾客到达作为第i件事发生; t(i+1)- t(i)= t1(i) (随机变量) 平衡关系： 当 t(i+1)?T(i) 时, T(i+1)=t(i+1)+t2(i+1); 否则, T(i+1)=T(i)+t2(i+1) 模拟20位顾客到收款台前的排队情况。 两个款台： 由于一位服务员既要收款又要装袋，所以设一个收款台的服务时间服从正态分布 N(2，2/3)。 平衡关系： 当 t(i+1)?max(T(i),T(i-1))，即前两位顾客都已离开收款台，或t(i+1)max(min(T(i),T(i-1)),max(T(1:i-2))) 即前两位顾客中至少有一位已经离开，两位之前的顾客都已离开时，第i+1位顾客到达收款处时至少有一个收款台是空闲着， 所以他马上接受服务 T(i+1)=t(i+1)+t2(i+1)，没有增加顾客等待的总时间; 否则, 他要排队等待前两位中还逗留在收款台的那位顾客离去或在两位前的某位因收款拖延较长时间的顾客离去，记他跟随的那位顾客的离去时间为k, 则 T(i+1)=k+t2(i+1)，顾客等待的总时间增加 k-t(i+1)。 作业 1. 假设如果今天是雨天，则明天有75%的可能是雨天，同样地如果今天是晴天，则明天有75%的可能是晴天。在假期的第一天是晴天（X1 = 0）。使用蒙特卡洛模拟估计一周内出现至少三个连续雨天的概率。 2.市场服务问题的综合分析，建模分析开设两个收款台的方案，根据模拟结果评价两个方案。 (1). 假设对不同的顾客一位服务员收款并装袋所需的服务时间服从正态分布 N(2，2/3)， 两位服务员分别收款装袋所需的服务时间服从正态分布 N(1，1/3). (2). 假设顾客到达的平均时间间隔是0.9分钟,顾客在收款台前排成一队。有两种安排方案：开一个出口，一人收款、一人装袋；开两个出口，每个人既收款又装袋。问商店经理应选择哪一种收款台的服务方案。 * * * * * * * * * * * *