计算机仿真-精.pptVIP

一、计算机仿真的基本概念 仿真常用术语: 仿真研究步骤 例3 随机变量 x = {0，1，2}表示每分钟到达超市收款台的人数，有分布列 xk 0 1 2 pk 0.4 0.3 0.3 模拟十分钟内顾客到达收款台的状况. 三、时间步长法 时间步长法 应用举例-池水含盐量问题 池水含盐量仿真流程 池水含盐量仿真结果 四、事件步长法 事件步长法 例7 收款台前的排队过程的仿真 假设： 1、顾客的到达收款台是随机的，平均时间间隔 为0.5分钟，即间隔时间服从?=2的指数分布。 2、对不同的顾客收款和装袋的时间服从正态分布 N(1，1/9)。 现要求模拟20位顾客到收款台前的排队情况， 主要关心问题：每个顾客的平均等待时间,队长,服务 员的工作效率. 假设 t(i): 第i位顾客到达时刻， t2(i):第i位顾客受到的服务时间(随机变量)， T(i): 第i位顾客离去时刻. 将第i位顾客到达作为第i件事发生; t(i+1)- t(i)= r(i) (随机变量) 平衡关系： 当 t(i+1)?T(i) 时, T(i+1)=t(i+1)+t2(i+1); 否则, T(i+1)=T(i)+t2(i+1) 五、Monte Carlo方法 应用举例-报童的策略 初始状态 时间步长加1 在当前步长内， 考察分析，计算和 记录系统的活动 仿真时间到否？ 结束 输出结果 是 否 例6:池水含盐量问题 某水池有2000m3水，其中含盐2kg，以每分钟6m3的速率向水池内注入含盐率为0.5kg/m3的盐水，同时又以每分钟4m3的速率从水池流出搅拌均匀的盐水．试用计算机仿真该水池内盐水的变化过程，并每隔10min计算水池中水的体积、含盐量和含盐率．欲使池中盐水的含盐率达到0.2kg/m3，需经过多少时间？ 池水含盐量随时间变化的示意图 分析:系统中，实体是水，属性是水的体积、含盐量和含盐率，活动是水的注入与流出，由于注入和流出活动的作用，池中水的体积与含盐量、含盐率均随时间变化，初始时刻含盐率为0.001kg/m3,以后每分钟注入含盐率为0.5 kg/m3的水6m3，流出混合均匀的盐水4m3,当池中水的含盐率达到0.2kg/m3时,结束仿真过程. 为了能定量地考察系统在任一时刻的属性，引入下列记号： 注水速度： WI＝6 m3/min 排水速度： WO＝4 m3/min 注入水的含盐率：SI＝0.5 kg/m3 最终含盐率：SF＝0.2 kg/m3 T时刻水的体积：VTm3 T时刻水的含盐量：ST kg T时刻水的含盐率：SR=ST/VT kg/m3 对于这样一个连续系统仿真时，必须在一系列离散时间点t0 t1 t2… tn上来进行考察，这些离散时间点之间的间隔T= ti- ti-1(i=1,2,…,n)就是时间步长．若取步长为1min时，就要每隔1min考察一次系统的状态特性，并作出相应的记录与分析．在本题中每隔1min池水的动态变化过程是： 每分钟水的体积增加 6-4＝2 m3； 每分钟向池内注入盐 6×0.5＝ 3 kg； 每分钟向池外流出盐 4×SR kg； 每分钟池内增加盐 3-4×SR kg． 取T＝0作为系统仿真的初始时刻，取T=0时系统的属性为系统的初始状态．取池中盐水的含盐率达到0.2kg/m3的时刻作为仿真的终止时刻． 初始化 按时间步长前进1min 计算池水体积VT,含盐量ST,含盐率SR 是否到达10min？ 结束 输出结果 是 否 含盐率达到SF? 打印次数加1 输出时间T,水体积VT 含盐量ST,含盐率SR 下一个10min开始，计时单元清零 是 否 0.19643 463.5802 2360 180 0.20026 474.6208 2370 185 0.18857 441.2655 2340 170 0.18044 418.6306 2320 160 0.17203 395.6642 2300 150 0.16331 372.3546 2280 140 0.15429 348.6897 2260 130 0.14494 324.6566 2240 120 0.13524 300.242 2220 110 0.12519 275.4319 2200 100 0.11478 250.2118 2180 90 0.10397 224.5663 2160 80 0.09275 1