系统工程(第5章).pptxVIP

系统仿真;一、系统仿真概述;仿真的实质和作用;仿真的实质和作用;仿真的应用;系统仿真方法;连续系统仿真;离散事件系统仿真;仿真的一般步骤;仿真语言;蒙特卡罗法简介;MC 的起源和发展 ;;思路：;MC的应用举例 ;假如有一台计算机产生[0，1]区间内均匀分布的随机数，让产生的这个随机数表示正方形边长的系数。用两个随机数表示一个点在正方形上的具体位置。然后判断该点是否落在不规则图形内。通过记录落在不规则图形内的点数和试验的总次数，这两个数值的比值就是他们的面积之比。使用计算机的好处就是能够在很短的时间内产生大量的随机数，而且还能记录统计结果。由此可见，蒙特卡罗方法依靠计算机产生的随机数来模拟实际系统。同时，计算机也负责相关的计数工作，以便获得最终的统计结果。;在计算积分上，MC的实用场合是计算重积分 其中 是 维空间的点，当 较大时，用MC方法比一般的数值方法有优点，主要是它的误差与维数 无关。;可以得到 的一个估计 具体试验步骤为 ;注1 随机投点法的思想简单明了，且每次投点结果服从二项分布， ，其中 注2 可证 是 的无偏估计。;;蒙特卡罗法应用——PERT网络计划仿真问题;;蒙特卡罗法应用;随机数的生成;随机数的生成;均匀随机数的产生;均匀随机数的产生;均匀随机数的产生;随机数;例：;;? ?用混合同余法产生的随机整数是模M的余数，必定小于M，显然的取值范围为M，最多只能有M个不同的取值。因此，产生的随机数有周期性，重复周期TM。 为了获得足够长周期的随机数序列，应尽可能取大的M值，并合理选择其他参数，使得到的是全周期的。一般认为，要得到最大周期，参数选择应满足如下条件： c0,且c与M互质； 乘子a满足(a-1)是4的倍数。 ;;1) 逆变换法 以指数分布为例 分布函数： 逆变换步骤： 计算分布函数。 在x取值范围内，置F(x)=R，即 用R表示x，即 x=F-1(R)=(-1/λ)ln(1-R) 产生均匀分布随机数R1,R2,…，计算所需随机变量， xi=(-1/λ)ln(Ri);2) 产生泊松分布的接受－拒绝技术 若相继两个事件出现的时间间隔服从指数分布，则单位时间内事件出现的次数服从泊松分布，其概率密度函数： 若到达时间间隔为xi，则有 化简为;产生泊松分布的步骤： 1.置n=0，S=1。 2.产生一个随机数Rn+1，S=S*Rn+1。 3.若Se-λ，则接受n，否则，拒绝现在的n，令n=n+1，转步骤2. ;3) 正态分布的直接变换 从两个独立随机变量R1,R2来产生两个独立的标准正态分布随机变量Z1,Z2的方法： 通过变换Y=μ+σZ，可以得到均值为μ，方差为σ 的正态分布随机变量。;4) 查表法产生离散分布随机变量 某随机变量取值只能为0、1、2，其概率分布如下： P(0)=p(X=0)=0.5 P(1)=p(X=1)=0.3 P(2)=p(X=2)=0.2 则累积分布为：;5）三角分布;α;输入数据分析;输出数据分析;离散系统动态仿真;输入数据分析;离散事件系统仿真——排队系统;排队系统有关概念;排队规则 损失制 等待制 先进先出（FIFO） 后进先出（LIFO） 随机服务 按优先级服务 最短处理时间先服务 混合制 限制队长 限制等待时间 限制逗留时间;排队系统分析;排队系统的性能研究;平均队长Lq和平均顾客数L： 其中，qi为[Ti-1，Ti]内排队人数。 Ti为第i个任意一类事件发生的时间。 忙期/闲期：忙期是指服务台全部处于非空闲状态的时间段，否则称为非忙期。闲期指服务台全部处于空闲状态的时间段。;排队系统仿真的建模;仿真钟的推进;面向事件的仿真时钟，又称事件调度法，是按下一最早发生事件的发生时间来推进仿真时钟的方法。采用这种方法，每当某一事件发生，即将仿真钟推进到该事件发生时刻并立即计算其直接后续事件的发生时间，在处理完当前事件所引起的系统变化之后，从未来将发生的各类事件中挑选最早发生的任何一类事件，将仿真钟推进到该事件发生时刻，再作以上重复处理。这个过程中，仿真以不等距的时间间隔向前推进，直到仿真运行满足某终止条件为止。;面向时间间隔的仿真时钟，也称固定增量推进法。采用这种方法，首先要确定某一时间单位T作为仿真钟推进的固定时间增量，即仿真钟从开始即按时间单位T等距推进(跳跃)，每次推进都需要扫描所有的活动，检查该时间区间内是否有事件发生。若无，则仿真钟继续等距推进；若有，则记录此时间区间，从而得到有关事件的时间参考数。若有若干事件同时发生，除了记录该事件的时间参数