系统工程概论_系统仿真.ppt

第一章 系统工程概论 System Engineering 蒙特卡洛法仿真库存问题： 离散事件系统仿真 —— 库存系统仿真 —— 例 解：第二步：计算商品缺货的的概率和平均缺货的个数； (1) 缺货的概率： 订货点低，可以减少库存量，但是有可能发生缺货损失。 当订货点（OP）在20个时，即商品库存量在20个时就开始订货，则不可能发生缺货，缺货的概率为0。 当订货点（OP）在19个时，订货期中的需要量（DDLT）在20个的情况下，就会发生缺货，缺货的概率为从表5中可知为P（DDLT＞19）＝P（DDLT＝20）＝0.0004 蒙特卡洛法仿真库存问题： 离散事件系统仿真 —— 库存系统仿真 —— 例 解：第二步：计算商品缺货的的概率和平均缺货的个数； 当订货点（OP）在18个时，订货期中的需要量（DDLT）在20个及19个的情况下，就会发生缺货，缺货的概率为从表5中可知为：P（DDLT＞18）＝P（DDLT＝20）＋P（DDLT＝19）＝0.0004+0.0006＝0.0010 当订货点（OP）在17个时，订货期中的需要量（DDLT）在20个、19个和18个的情况下，就会发生缺货，缺货的概率为从表6中可知为：P（DDLT＞17）＝P（DDLT＝20）＋P（DDLT＝19）＋P（DDLT＝18）＝0.0004+0.0006+0.0012＝0.0022 蒙特卡洛法仿真库存问题： 离散事件系统仿真 —— 库存系统仿真 —— 例 解：第二步：计算商品缺货的的概率和平均缺货的个数； 依次类推，可以求出订货点（OP）从20个直到0个为止，所有发生缺货的概率，如表6-9所示。 表6-9 订货点OP和缺货的概率 0.2306 10 0.0000 20 0.3244 9 0.0004 19 0.4354 8 0.0010 18 0.5522 7 0.0022 17 0.6476 6 0.0046 16 0.7386 5 0.0094 15 0.8144 4 0.0230 14 0.8722 3 0.0452 13 0.9496 2 0.0832 12 0.9740 1 0.1438 11 0.9902 0 缺货的概率P（DDLT＞OP） 订货点OP 缺货的概率P（DDLT＞OP） 订货点OP 蒙特卡洛法仿真库存问题： 离散事件系统仿真 —— 库存系统仿真 —— 例 解：第二步：计算商品缺货的的概率和平均缺货的个数； (2) 平均缺货的个数： 当订货点OP＝20个时，不可能发生缺货的现象，缺货的个数为0。即：E（DDLT＞20）＝0 当订货点OP＝19个时，而订货期的需要量DDLT在20个的情况下，就会发生缺货，个数为：E（DDLT＞19）＝（20－19）×P（DDLT＝20）＝1×0.0004＝0.0004 离散事件系统仿真 —— 随机变量的产生 解：则 ，其中为雅可比行列式，即 从而得到 其中f(ρ)，f(θ)，分别为随机变量ρ，θ的密度函数 反变换法 —— 例6-3：正态分布随机变量的产生 离散事件系统仿真 —— 随机变量的产生 解：其中f(ρ)，f(θ)，分别为随机变量ρ，θ的密度函数 其分布函数为 反变换法 —— 例6-3：正态分布随机变量的产生 ρ,θ的分布函数具有封闭形式。可用反变换法，即独立产生两个[0,1]均匀分布的随机变量u1,u2，分别对f(ρ)和f(θ)进行反变换，得 根据x1, x2与ρ，θ的变换关系，可得 仿真时钟的推进：是离散事件系统仿真不可缺少的组成部分，是仿真的时间控制部件。离散事件系统仿真中仿真时钟的推进方法是系统仿真的基础。事件法和时间间隔法两种： 事件法：又称面向事件的仿真时钟或事件调度法。它按照下一最早发生事件的发生时间来推进仿真时钟，仿真以不等距的时间间隔。具体而言，该方法是在处理完当前事件所引起的系统变化后，从未来将发生的各类事件中挑选最早发生的任何一类事件，将仿真钟推进到该事件发生时刻，再作以上重复处理直到仿真运行满足某终止条件为止。 离散事件系统仿真 —— 离散系统的仿真策略 仿真时钟的推进： 时间间隔法：又称面向时间间隔的仿真时钟或固定增量推进法。采用该法之前需确定某一时间单位T作为仿真钟推进的固定时间增量。仿真从开始即按时间单位T等距推进，每次推进都需要扫描所有活动，检查该时间区间内