最新计算机仿真.pptVIP

计 算 机 仿 真 1.系统仿真概论 2. 时间步长法 3. 事件表法 4. Monte Carlo 方法 1. 系统仿真概论 1.1 系统 1) 系统 一些具有特定的功能、相互之间以一定的规律联系的对象所组成的总体。 2) 系统边界 每个系统都存在于一定的环境中, 系统与环境之间的分界称为边界. ?? 边界确定了系统的范围，边界以外对系统的作用称为系统的输入，系统对边界以外环境的作用称为系统的输出. 随机型活动通常用概率分布加以描述，如电路系统中的噪声、电话系统中通话时间的长短等，属于随机系统。 2 时间步长法 2.1 仿真时钟 在进行系统仿真时，可以把整个仿真过程分为许多相等的时间间隔，时间间隔的长度可根据实际问题分别取作秒、分、小时、天、周、月、年等。程序中按此间隔前进的虚拟时钟就是仿真的时钟，此时间间隔称为仿真时钟的步长。 选取系统的一个初始状态作为仿真时钟的零点，仿真时钟每步进一次，就对系统的所有实体、属性和活动进行一次前面的扫描考察，按照预定的计划和目标进行分析、计算和记录系统状态的变化，这个过程一直进行到仿真时钟结束为止。此即为时间步长仿真法。 3 事件表法 3.1 事件表 事件表好像一本记事本，干完一件事后就把它从记事本中勾销，而把新的要完成的工作再登记到记事本中。 事件表法的主要思想是：将系统的仿真过程看成一个事件序列，根据事件出现的时序，用事件表来调度事件执行的顺序，以此使得系统的仿真过程有条不紊地进行下去。这种方法要求对系统的各种事件进行详细的描述，因此，当事件之间没有太多的相互作用和事件数目不是太多时，应用事件表法比较有效。 4 Monte Carlo 方法 4.1 Monte Carlo 方法 Monte Carlo 是世界著名赌城－摩洛哥的蒙特卡罗，在二次世界大战中，美国军方将一项绝密研究计划的代号命名为Monte Carlo ，其目的是研究铀裂变过程中链式反应的能量计算问题。 Monte Carlo 方法是随机模型的计算机仿真方法，在用传统方法难以解决的问题中，有很大一部分可以用概率模型描述，由于这类问题含有不确定的随机因素，难以用定量分析法得到解析结果，在这种情况下， Monte Carlo 方法是非常有效。 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 上述仿真过程做了如下假设 ①仿真开始前没有订货； ②在上次订货到达前，如果缺货