系统工程-第四章探讨.pptVIP

第四章 系统仿真 第四章 系统仿真 4.1 基本概念 4.2 蒙特卡罗仿真 4.3 均匀随机数的产生方法 4.4 随机变量的产生方法 4.5 离散事件随机系统的仿真方法 4.6 随机分布的确定 4.7 系统动力学仿真 4.8 仿真语言 4.9 可视化仿真技术 4.10 虚拟现实技术 4.11 仿真技术在矿业中的应用 作业 4.1 基本概念 在下图中，已知采区来煤量q1是一个随机变量，根据统计结果，它服从正态分布，已知其均值μ=1.45t/ min，其标准差σ=0.095t/min。大巷是胶带输送机，其运输量是常数：q2=1.33t/min。求8h后，采区最大煤仓容量。 投掷硬币随机模拟实验 建立一个计算机能够运算的模拟模型,问题的关键 是用随机数去模拟投掷硬币的实验。 所谓随机数就是一组随机出现的数列，通常用的是 [0，1]区间的随机数。这时要求[0，1]区间的随机数必 须具备下列性质：一是它的均匀性，即这些随机数落在 [0，1]区间内任一位置的概率是相同的；另一个是它的 随机性，即这些数在 [0，1]区间内任一位置出现与否 是随机的、独立的。如以R表示随机数，则出现0≤R≤0.5 和 0.5≤R≤1的概率相等，都为0.5.这样如以随机数R≤0.5模拟出现正面，以R≥0.5模拟出现反面，则模拟 模型就建立起来了。 设 I 是投掷次数 N 为正面出现次数 R 为随机数 P是正面出现频率 4.2 蒙特卡罗仿真 设f(x)是区间（a,b）内的连续函数，且已知在该区 间内f(x)=0,要求计算 在（a，b）和（0，c） 区间中分别产生两个均 匀分布的随机数xi和yi， 由此构成的坐标点 Pi(xi，yi)必定在矩形 abde内。 若yi≤f(xi)，则接受该点并 作计数统计；否则舍弃该点 不作统计. 设共产生N个随机仿真点，如果共接受个m点， （即恰在曲线上或在曲线以下的点），当N充分大时 有 曲线下面积 4.3 均匀随机数的产生方法 4.3.1 均匀随机数 [0，1]区间上的均匀随机数，它的密度函数f(x) 和分布F(x)分别为 4.3.2 均匀随机数的产生方法 用数学公式产生的均匀随机数，称为伪随机数,常 用的公式方法是线性同余法 . 若有a、b两个整数被某个整数M除后，它们的余数 相同，则称a、b两数对模M同余，记作 a=b(ModM) 求余数的公式为 例如, 产生随机数的线性同余法是一个递推公式： 即 通过下式得到落在[0，1]区间内的随机数γ： 两个较好的计算公式如下 4.3.3 均匀随机数检验方法 设N个随机数从小到大排列起来并分成互不重叠 的等区间的k段.从理论上讲，这些随机数应该均匀分 布在每一段上，在每一段上理论频数为 如果统计出随机数实际落在每段的频数为ni，则令 x2值越大表明随机数列越不均匀。 给定判定标准称α为显著度,根据自由度r=k-1，从x2表中查得 值，如果计算所得的x2 ，就认为该随机数列在显著度α的情况下，其均匀性符合要求。 例如，假设由公式求得了100个随机数，将其分为 5段（k=5）。找出该随机数列的最大值为31252，最 小值为1488，并求出每段区间数值为5952.8，便可找 出各段的范围。已知各段的理论频数mi=100/50=20， 统计出各段的实际频数ni，最后求出x2=2。查x2表， 根据自由度r=5-1=4，显著度α=0.05查分布临界值表， 得 。由于 ,则认为在显著度0.05情 况下，该随机数列的均匀性符合要求。 随机数的独立性 前后距离为j个数的相关系数为： 当N充分大时（N-j50），若取零假设H0：ρ=0，并 取统计量 (近似服从N（0，1）标准正 态分布)，同时选定显著度α=0.05，则根据概率论， 当|u|1.96时为差异显著，拒绝H0假设；当|u|≤1.96 时接受假设，即检验的随机数符合独立性要求。 4.3.4 计算机函数 计算机高级语言有产生随机数的标准函数，可直接调用。 4.4 随机变量的产生方法 4.4.1 逆变换法 已知某随机变量的分布函数F（x） ，并知其定义域为 [0，1]区间，这时可用[0，1]区间的均匀分布随机