模拟40436.pptVIP

排队服务系统的模拟 1、模拟目的 2、系统假设与输入 3、系统的状态 4、初始条件和终止条件 5、模拟过程的运行 6、系统的性能指标 7、与分析模型的比较 8、其它排队系统 假设由概率分布产生了如下的数据： 3、系统的状态 队长：排队等候的顾客数目。记作L(t). 服务员的状态用S(t)表示，当服务员工作时，令S(t)=1，空闲时令S(t)=0. 事件：引起系统状态L(t)和S(t)改变的行为。 　　　 “顾客到达”和“顾客离开” 设第k位顾客到达和离开时刻分别为ak和dk，则有 　　 4、初始条件和终止条件 　　不妨假设t=0是第1位顾客到达的时刻，此 时服务员处于空闲状态。模拟的终止时刻是T （给定值），或模拟直到第n个（给定值）顾 客进入服务的时间。 5、模拟过程 　　进行模拟时，可以事先产生一批数据ik,sk,再计算ak,dk，然后让时钟t按ak,dk从小到大的顺序跳动；也可在时钟t跳到某一事件发生时，才产生一个所需要的i或s。这里采用后一种方法。 　　设当前时钟为t，队长L(t)记作WL，服务员状态S(t)记作SS，t以后下一个“顾客到达”事件的发生时刻记作AT, t以后下一个“顾客离开”事件的发生时刻记作DT. 6、系统的性能指标 平均队长 平均待时间 服务利用率 7、与分析模型的比较 　　　对于一般的排队服务系统，很难建立分析模型并求出结果，即使是单服务员系统，也只有在顾客到达间隔和服务时间服从某种特殊的概率分布时，才可以在稳态情况下由分析模型得到上述三个指标的解析形式的结果。 8、其它的排队服务系统 有两个（或多个）服务员 串联服务系统，每个顾客接受多次服务 从排队顾客中选服务时间最短的最先服务 队的长度有限制 带有优化目标的排队服务系统 排队问题 Email: cquhrb@ ATDT WL=0,SS=0,AT=0,DT=999 置t=AT SS=0 置t=DT WL0 WL=WL+1 令SS＝1 产生s 置DT=t+s 产生i 置AT＝t+i t=T 令SS＝0 WL＝WL－1 产生s 置DT＝t+s 停止 置DT＝999 是 否 是 否 是 否 否 是 分析模型基本假设： 1、到达间隔服从参数λ的指数分布； 2、服务时间服从参数为μ的指数分布； 3、排队规则是“先到先服务”，顾客数目和排队长度都无限制。 排队论中将这样的单服务员系统记作M/M/1。 　　 服务强度： 稳态（ρ1） 库存太多造成浪费或资金积压，库存太少不能满足需求也造成损失，需要进行决策：何时进货，进多少，使得平均费用最少，而收益最大。 库存系统 方案甲：按前一天的销售量作为当天的生产量； 方案乙：按前二天的平均销售量作为当天的生产量； 假定市场对该产品的每天需要量是一个随机变量，但从以往的统计分析得知它服从正态分布N(135,22.42）. 例1 某企业生产易变质的产品。（如蛋糕）当天生产的产品必须当天售出，否则就会变质。该产品单位成本为2.5元，单位产品售价为5元。企业为避免存货过多而造成损失，拟从以下两种货存方案中选出一个较优的方案.(如何决定当天的生产数量？) 模拟的基本思路： 需要获得市场对该产品的需要量的样本 值；(产生抽样值) ② 寻找货存量与需要量之间的关系； ③ 按照两种不同方案（货存量≥{≤}需要量）计算出经 T 天后企业的利润值(累计)； ④ 比较大小，从中选出一个较优的方案； 方案甲：货存量Y1 分析两种状态下 X＞Y1和X≤Y1 的利润值L1 置初始状态 获取产品的需求量X X= normrnd（135，22.5） N = 1 方案乙：货存量Y2 分析两种状态下 X＞Y2和X≤Y2 的利润值L2 累计利润值TL1 累计利润值TL2 N=N+1 判断：N≤T ( T天) 比较计算：max{TL1，TL2} 输出最大利润值及方案 MATLAB编制的程序： kucun.m。 运行程序：kutest.m 初始状态：T=100;S1=136;S21=126;S22=148; [TL1,TL2]=kucun(T,S1,S21,S22) 观察计算结果 如果多数情况有：TL1TL2, 则认为方案乙较优. 思考： 1）一般模拟结果波动性较大，如何减少这种波动？ 2）修改上述程序。 可靠性问题 一设备上有三个相同的轴承，每个轴承正常工作寿命为随机变量，概率分布如下： 0.05 0.05 0.06 0.02 0.12 0.09 0.13 0.25 0.13 0.10 概率 1900 1800 1700 1600 1500 1400 1300 1200 1100 1000 寿命/h 有轴承损坏→设备停止工作→检