离散事件系统仿真基础优秀课件.ppt

C 语言中与随机数有关的函数 ? ? ? randomize(void); 初始化 x=random(int M); 产生 0 ~ M 之间的随机 数 x=rand(void); 产生 0 ~ 2 15 -1 之间的随 机数 Matlab 中与随机数有关的函数 ? ? ? ? Rand(n) ： n 个 [0,1] 之间均匀分布随机数 Rand(m,n) ： m*n 个 [0,1] 之间均匀分布 随机数 randn(n):n 个 N(0,1) 标准正态分布随机数 randn(m,n) ： m*n 个 N(0,1) 标准正态分布 随机数 例 射击命中率 在我方某前沿防守地域，敌人以一个炮排 （含两门火炮）为单位对我方进行干扰和破 坏．为躲避我方打击，敌方对其阵地进行了伪 装并经常变换射击地点．经过长期观察发现， 我方指挥所对敌方目标的指示有 50 ％是准确的， 而我方火力单位，在指示正确时，有 1/3 的射 击效果能毁伤敌人一门火炮，有 1/6 的射击效 果能全部消灭敌人． 现在希望能用某种方式把我方将要对敌人 实施的 20 次打击结果显现出来，确定有效射击 的比率及毁伤敌方火炮的平均值 说明 ? ? 这是一个概率问题，只涉及到简单的随机数的 产生和处理，可以通过理论计算得到相应的概 率和期望值。但这样只能给出作战行动的最终 静态结果，而显示不出作战行动的动态过程 为了能显示我方 20 次射击的过程，考虑采用计 算机模拟 1. 问题分析 需要模拟出以下两件事： ? (1) 观察所对目标的指示正确与否 ? ? 模拟试验有两种结果，每一种结果出现的概率都是 1/2 因此，可以用 投掷一枚硬币 的方式予以确定，当硬 币出现正面时为指示正确，反之为不正确 ? (2) 当指示正确时，我方火力的射击结果 模拟试验有三种结果： ? 毁伤一门火炮的可能性为 1/3( 即 2/6) ， ? 毁伤两门的可能性为 1/6 ， ? 没能毁伤敌火炮的可能性为 1/2( 即 3/6) ． 这时 可用投掷骰子的方法来确定： ? 出现１、２、３点：则认为没击中敌人； ? 出现４、５点： 则认为击毁敌人一门火炮； ? 出现６点： 则认为击毁敌人两门火炮 2. 符号假设 i ：要模拟的打击次数； k1 ：没击中敌人火炮的射击总数； k2 ：击中敌人一门火炮的射击总数； k3 ：击中敌人两门火炮的射击总数． E ：有效射击比率； E1 ： 20 次射击平均每次毁伤敌人的火炮数 3. 程序流图 1 ， 2 ， 3 初始化 :i=0,k 1 =0,k 2 =0,k 3 =0 i=i+1 Y 硬币正面 ? N 骰子点数 ? 4 ， 5 6 k 1 =k 1 +1 k 2 =k 2 +1 k 3 =k 3 +1 Y k 1 =k 1 +1 i ＜ 20? N E=(k 2 +k 3 )/20 E 1 =0*k 1 /20+1*k 2 /20+2*k 3 /20 停止 4. 模拟结果 试验 序号 １ ２ ３ ４ ５ ６ ７ ８ ９ １０ 投硬币 结 果 正 正 反 正 正 反 正 正 反 反 指示 正确 ∨ ∨ ∨ ∨ ∨ ∨ 指 示 不正确 ∨ ∨ ∨ ∨ 掷骰子 结 果 ４ ４ １ ２ ３ ６ 消灭敌人火炮数 ０ ∨ ∨ ∨ ∨ ∨ ∨ ∨ １ ∨ ∨ ２ ∨ 试验 序号 １１ １２ １３ １４ １５ １６ １７ １８ １９ ２０ 投硬币 结 果 正 反 正 反 正 正 正 正 反 正 指示 正确 ∨ ∨ ∨ ∨ ∨ ∨ ∨ 指 示 掷骰子 不正确 结 果 ∨ ∨ ∨ ２ ３ ６ ４ ２ ４ ６ 消灭敌人火炮数 ０ １ ２ ∨ ∨ ∨ ∨ ∨ ∨ ∨ ∨ ∨ ∨ ? 从以上模拟结果可以计算出： 有效射击 E =7/20=0.35 13 4 3 平均值 E 1 ? 0 ? ? 1