系统工程蒙特卡罗

系统工程用蒙特卡罗法对下面报童问题进行仿真。对新闻的重要度分布概率进行随机数编码。随机数采用两位数字（0099） ，如下表1 2 3 4重要程度 概率 累积概 率 随机编 码重要 0.35 0.35 0034一般 0.45 0.80 3579平淡 0.20 1.00 8999对报纸的销量分布概率进行随机数编码。随机数采用两位数字（0099） ，如下表累积分布报纸需求量 重大新 闻 随机编 码 一般新 闻 随机编码 平淡新 闻 随机编 码40.00 0.03 0002 0.10 0009 0.44 004350.00 0.08 0307 0.28 1027 0.66 446560.00 0.23 0822 0.68 2867 0.82 668170.00 0.43 2342 0.88 6887 0.94 829380.00 0.78 4377 0.96 8895 1.00 949990.00 0.93 7892 1.00 9699 1.00 100.00 1.00 9399 1.00 1.00 用 excel 随机产生 100 个数，重要度 报纸需求模拟次数 随机数 程度 随机数 报纸需求量1 5 重要 84 902 41 一般 84 703 73 一般 33 604 53 一般 22 505 29 重要 1 406 9 重要 71 807 67 一般 7 408 69 一般 73 709 56 一般 45 6010 93 平淡 86 7011 61 一般 25 5012 3 重要 3 5013 34 重要 47 8014 55 一般 47 6015 15 重要 17 6016 51 一般 63 6017 98 平淡 10 4018 10 重要 14 6019 52 一般 0 4020 77 一般 91 8021 25 重要 15 6022 71 一般 70 7023 88 平淡 20 4024 57 一般 2 4025 58 一般 94 8026 66 一般 71 7027 75 一般 18 5028 28 重要 13 6029 33 重要 61 8030 19 重要 55 8031 6 重要 68 8032 73 一般 72 7033 97 平淡 97 8034 23 重要 68 8035 61 一般 76 7036 68 一般 81 7037 15 重要 84 9038 66 一般 66 6039 58 一般 85 7040 98 平淡 90 7041 62 一般 74 7042 24 重要 24 7043 27 重要 48 8044 62 一般 87 7045 26 重要 27 7046 4 重要 7 5047 58 一般 42 6048 13 重要 93 10049 89 平淡 8 4050 24 重要 53 8051 20 重要 84 9052 35 一般 70 7053 48 一般 18 5054 91 平淡 84 7055 80 平淡 32 4056 41 一般 13 5057 40 一般 5 4058 15 重要 0 4059 88 平淡 85 7060 48 一般 23 5061 33 重要 18 6062 88 平淡 9 4063 24 重要 92 9064 20 重要 42 7065 70 一般 58 6066 60 一般 92 8067 53 一般 1 4068 75 一般 57 6069 54 一般 94 8070 64 一般 91 8071 20 重要 3 5072 94 平淡 6 4073 93 平淡 63 5074 41 一般 36 6075 2 重要 92 9076 4 重要 86 9077 84 平淡 57 5078 42 一般 35 6079 12 重要 67 8080 7 重要 54 8081 33 重要 36 7082 22 重要 51 8083 96 平淡 14 4084 54 一般 63 6085 32 重要 8 6086 16 重要 51 8087 91 平淡 93 7088 23 重要 45 8089 17 重要 97 10090 55 一般 6 4091 52 一般 24 5092 49 一般 12 5093 18 重要 46 8094 80 平淡 96 8095 70 一般 24 5096 70 一般 65 6097 51 一般 77 7098 20 重要 27 7099 5 重要 37 70100 27 重要 51 80根据上述的结果得到报纸的需求量，如下表报纸需求量的个数 个数 概率值40 15 0.1550 14 0.1460 18 0.1870 23 0.2380 22 0.2290 6 0.06100 2 0.02总数 100第二步 计算报纸过多的概率和平均需求量。当需求点 A=40，则不可能发生报纸多，过多的概率为 0当需求点 A=50，报纸需求量为 40 的情况下，就会发生多货，多货的概率从上表中可知为P（DDLT50）=P（DDLT=40）=0.15当需求点 A=60，报纸需求量为 40.50 的情况下，就会发生多货，多货的概率从上表中可知为P（DDLT60）=P（DDLT=40）+P （DDLT=50）=0.29当需求点 A=70，报纸需求量为 40.50.60 的情况下，就会发生多货， 、多货的概率从上表中可知为P（DDLT70）=P（DDLT=40）+P （DDLT=50）+P（DDLT=60 ）=0.47当需求点 A=80，报纸需求量为 40.50.60.70 的情况下，就会发生多货，多货的概率从上表中可知为P（DDLT60）=P（DDLT=40）+P （DDLT=50）+P（DDLT=60 ）+P（DDLT=70）=0.70当需求点 A=90，报纸需求量为 40.50.60.70.80 的情况下，就会发生多货，多货的概率从上表中可知为P（DDLT0 ） =P（DDLT=40）+P （DDLT=50）+P（DDLT=60 ）+P（DDLT=70）+P（DDLT=80）=0.92当需求点 A=100，报纸需求量为 90,80,70,60,50.40 的情况下，就会发生缺货，缺货的概率从上表中可知为P（DDLT100）=P（DDLT=40）+P （DDLT=90）+P（DDLT=80 ）+P（DDLT=70）+P（DDLT=60）+P （DDLT=50）=0.98所有发生多货的概率如表40 0.50 0.1560 0.2970 0.4780 0.7090 0.92100 0.98平均多货的个数。报纸需求量的个数 个数 概率值40 15 0.1550 14 0.1460 18 0.1870 23 0.2380 22 0.2290 6 0.06100 2 0.02总数 100当需求点 A=40，不可能发生多货，多货个数为 0，即E（DDLT40）=0当需求点 A=50，报纸需求量为 40 的情况下，就会发生多货个数为 E（DDLT50）=P （DDLT=40）*（50-40）=1.5当需求点 A=60，报纸需求量为 40.50 的情况下，就会发生多货个数为 E（DDLT60）=P （DDLT=40）*（60-40）+P（DDLT=50）*（60-50）=4.4当需求点 A=70，报纸需求量为 40.50.60 的情况下，就会发生多货个数为 E（DDLT70）=P （DDLT=40）*（70-40）+P（DDLT=50）*（70-50）+P（DDLT=60 ）*（ 70-60）=9.1当需求点 A=80，报纸需求量为 40.50.60.70 的情况下，就会发生多货个数为 E（DDLT80）=P （DDLT=40）*（80-40）+P（DDLT=50）*（80-50）+P（DDLT=60 ）*（ 80-60）+P（DDLT=70）*（80-70）=16.1当需求点 A=90，报纸需求量为 80,70,60.50.40 的情况下，就会发生多货个数为 E（DDLT90）=P （DDLT=40）*（90-40）+P（DDLT=50）*（90-50）+P（DDLT=60 ）*（ 90-60）+P（DDLT=70）*（90-70）+P （DDLT=80）*（90-80 ）=25.3当需求点 A=100，报纸需求量为 90，80,70,60，50.40 的情况下，就会发生多货个数为 E（DDLT100）=P （DDLT=40）*（100-40）+P（DDLT=50）*（100-50）+P（DDLT=60 ）*（100-60 ）+P （DDLT=70）*（100-70）+P（DDLT=80）*（100-80 ）+P（DDLT=90）*（100-90）=35.140 050 1.560 4.470 9.180 16.190 25.3100 35.1