计算机仿真100714.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 1) 建立概率模型。 当粒子在某内点B时，统计与其相连的临点数k,并将临点依次用1到k编号，然后从1到k中等概率地任选一个数，则粒子移到对应此编号的点。若此点为边界点，则游动结束，记录边界吸收的能量，然后开始下一轮游动。 例2 二维随机游动问题 { A i i ＝1,…,10}组成，粒子在任何一个内点以等概率向其临近点移动，若到达边界点Ai则被吸收。边界获得能量f i, 求由Q点出发的粒子被边界吸收的能量。 如图所示的网格，由内点 {Q, P i, i =1,2,3,4,5 } 及边界点 2)抽样计算 在第m轮游动中，当粒子游动到某内点B时，生成均匀分布随机数 r，若与B其相连的临点数为k，则将区间(0,1) 分成k个小区间，若r落在第i个小区间，则粒子游动到B的第i个临点，如此继续，直至粒子被某边界点A i吸收，则本轮游动结束，记录边界吸收的能量Em = Ei，然后开始下一轮游动。 则由Q点出发的粒子被边界吸收的能量的近似值为 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 2) 模型建立 输入参数 订货点 P ；订货量 Q 注意 P Q 输入常量 从发出订货单到收到货物见隔d=3天 每辆自行车每天的保管费为c1=0.75元 每辆自行车每天的缺货费c2=1.8元 每次的订货费c3=75元 原始库存s=115辆 仿真总时间 T＝150天 变量 到货日期D，需求量R,日期t，预定到货量q 输出参数 总费用 C ①数据整理与分析 ②系统状态 系统状态变量：需求量，库存量，总费用 ③有关费用的计算 1) 发出订货单需要计算订货费； 2)当需求量小于库存量时需计算保管费； 3)当需求量大于库存量时需计算缺货费； 初 始 化: c1=0.75,c2=1.8,c3=75,s=115, T=150,d=3,D=0,q=0, C=0, t=1 s←s-R, C ← C + c1×s 是 仿真是否完毕 t = T ? 否 R ≤ s ? 否 结束 输出结果 是 产生今天的随机需求量R 否 s ← s＋Q, q=0 t = D ？ 是 C ← C＋(R-s) ×c2 D = t + 3, q=Q C ← C＋c3 s + q ≤P ? 是 t ← t + 1 否 初 ④ 仿 真 程 序 流 程 输入参变量 P,Q 上述仿真过程做了如下假设 ①仿真开始前没有订货； ②在上次订货到达前，如果缺货再订货； 显然，假设2实际上是一种订货管理策略，可修改。 3) 程序运行、分析 4) 设计格式，输出仿真结果 比较这五种方案仿真得到的总费用，可以看出，方案4最好，即仓库管理员应取最低订货点为175辆，每次订250辆自行车的方案，这时在150天的总费用为26094元。 方案编号 方案1 方案2 方案3 方案4 方案5 总费用 38679.75 31268.25 29699.25 26094.00 27773.25 五种方案各运行100次的平均结果 思考 ①程序中变量q的作用？ ②如何修改假设2？写出修改后的程序。 ？ 3 事件表法 3.1 事件表 事件表好像一本记事本，干完一件事后就把它从记事本中勾销，而把新的要完成的工作再登记到记事本中。 事件表法的主要思想是：将系统的仿真过程看成一个事件序列，根据事件出现的时序，用事件表来调度事件执行的顺序，以此使得系统的仿真过程有条不紊地进行下去。这种方法要求对系统的各种事件进行详细的描述，因此，当事件之间没有太多的相互作用和事件数目不是太多时，应用事件表法比较有效。 初 始 化 3.2 仿真流程 仿真时钟步进 事件A类型？ 调用处理第i类事件的子程序 i类 产生条件事件？ 事件表 否 更新状态数据 产生新事件，删除已处理的事件 否 仿真是否完毕 仿真结束 输出结果 是 是 调用处理相关条件 事件的子程序 扫描事件表 处理最近事件A 更新 3.3 事件表法特点 1) 仿真时钟每次前进的步长取决于事件之间的间隔； 2) 在一个步长内，系统的状态不变； 3) 只有在有事件发生时，才对系统进行一次全面的考察。 3.4 实例分析 例1 排队系统 （考虑一个收款台的排队系统） 某杂货店只有一个收款台，顾客到达收款台的间隔服从均值为4.5的指数分布，每个顾客的服务时间服从均值为3.2、标准差为0.6的正态分布。时间的单位为分，服务时间不取负值。 试对100个顾客去收款台缴款排队过程进行仿真，并估计系统的以下特征： 1）顾客的平