第九章 模拟模型

1、 输入区:输入区:试验区:试验区: 变量1_标题直接输入值的固定参数1试验次数 =结果单元 变量2_标题直接输入值的固定参数21试验结果1 变量3_标题用控件可以调整的参数12试验结果2 变量4_标题用控件可以调整的参数23试验结果3 变量5_标题随机数分布参数14试验结果4 变量6_标题随机数分布参数25试验结果5 变量7_标题随机数分布参数36试验结果6 7试验结果7 生成区:生成区:8试验结果8 随机数1_标题=f1-1(rand(),分布参数1,分布参数2)9试验结果9 随机数2_标题=f2-1(rand(),分布参数1,分布参数3)10试验结果10 11试验结果11 输出区:输出区:

2、12试验结果12 结果_标题=f(随机数1,随机数2,变量1,变量2,变量3,变量4)13试验结果13 14试验结果14 统计区:统计区:15试验结果15 样本平均值=average(试验结果1,试验结果30)16试验结果16 样本标准差=stdev(试验结果1,试验结果30)17试验结果17 95%置信下限 =norminv(0.025,样本平均值,样本标准差)18试验结果18 95%置信上限 =norminv(0.975,样本平均值,样本标准差)19试验结果19 偏斜系数=skew(试验结果1,试验结果30)20试验结果20 峰态系数=kurt(试验结果1,试验结果30)21试验结果21

3、22试验结果22 图形区:图形区:23试验结果23 接收区间频数分布24试验结果24 分段统计点1=frequency(统计点1:统计点8,结果1:结果30)25试验结果25 分段统计点2=frequency(统计点1:统计点8,结果1:结果30)26试验结果26 分段统计点3=frequency(统计点1:统计点8,结果1:结果30)27试验结果27 分段统计点4=frequency(统计点1:统计点8,结果1:结果30)28试验结果28 分段统计点5=frequency(统计点1:统计点8,结果1:结果30)29试验结果29 分段统计点6=frequency(统计点1:统计点8,结果1:结

4、果30)30试验结果30 蒙特卡洛模拟模型的一般框架蒙特卡洛模拟模型的一般框架 步骤 步骤 步骤 步骤 步骤 步骤 步骤 单价2345678 概率5%10%20%30%20%10%5% 0.0%0.0% 1.0%1.0% 2.0%2.0% 3.0%3.0% 4.0%4.0% 5.0%5.0% 6.0%6.0% 7.0%7.0% 8.0%8.0% 9.0%9.0% 10.0%10.0% 概概 率率 值值 -16.5-16.5-11.4-11.4-6.3-6.3-1.2-1.23.93.99.09.014.114.119.319.324.424.429.529.534.634.6 项目净现值(百万

5、元)项目净现值(百万元) 投资项目净现值的概率分布图投资项目净现值的概率分布图 投投资资项项目目净净现现值值的的累累计计概概率率分分布布 0 0. .0 0% % 1 10 0. .0 0% % 2 20 0. .0 0% % 3 30 0. .0 0% % 4 40 0. .0 0% % 5 50 0. .0 0% % 6 60 0. .0 0% % 7 70 0. .0 0% % 8 80 0. .0 0% % 9 90 0. .0 0% % 1 10 00 0. .0 0% % - -2 20 0. .0 0- -1 15 5. .0 0- -1 10 0. .0 0- -5 5. .0

6、 00 0. .0 05 5. .0 01 10 0. .0 01 15 5. .0 02 20 0. .0 02 25 5. .0 03 30 0. .0 0 项项目目净净现现值值( (百百万万元元) ) 累累 计计 概概 率率 值值 投投资资项项目目大大于于某某净净现现值值的的概概率率 4 43 3. .4 4% % 0 0. .0 00 0% % 1 10 0. .0 00 0% % 2 20 0. .0 00 0% % 3 30 0. .0 00 0% % 4 40 0. .0 00 0% % 5 50 0. .0 00 0% % 6 60 0. .0 00 0% % 7 70 0.

7、.0 00 0% % 8 80 0. .0 00 0% % 9 90 0. .0 00 0% % 1 10 00 0. .0 00 0% % - -2 20 0. .0 0- -1 15 5. .0 0- -1 10 0. .0 0- -5 5. .0 00 0. .0 05 5. .0 01 10 0. .0 01 15 5. .0 02 20 0. .0 02 25 5. .0 03 30 0. .0 0 项项目目净净现现值值( (百百万万元元) ) 概概 率率 值值 13.2 净现值大于5.6 概率值等于43.4% 输入区:输入区:输出区:输出区:试验区:试验区: 变量1_标题标题标题=

8、f2() 关键参数1关键参数2 变量2_标题=f1()=f2()1试验结果1试验结果1 变量3_标题工作区:工作区:2试验结果2试验结果2 变量4_标题序列3试验结果3试验结果3 变量5_标题4试验结果4试验结果4 变量6_标题关键参数5试验结果5试验结果5 变量7_标题6试验结果6试验结果6 7试验结果7试验结果7 生成区:生成区:8试验结果8试验结果8 随机数1_标题=f1-1()9试验结果9试验结果9 随机数2_标题=f2-1()10试验结果10 试验结果10 11试验结果11 试验结果11 统计区:统计区:12试验结果12 试验结果12 样本1平均值13试验结果13 试验结果13 样本

9、1标准差14试验结果14 试验结果14 95%置信下限15试验结果15 试验结果15 95%置信上限16试验结果16 试验结果16 偏斜系数17试验结果17 试验结果17 峰态系数18试验结果18 试验结果18 19试验结果19 试验结果19 图形区:图形区:20试验结果20 试验结果20 接收区间频数分布21试验结果21 试验结果21 22试验结果22 试验结果22 23试验结果23 试验结果23 24试验结果24 试验结果24 25试验结果25 试验结果25 26试验结果26 试验结果26 27试验结果27 试验结果27 28试验结果28 试验结果28 29试验结果29 试验结果29 30

10、试验结果30 试验结果30 活动列统计列 系统模拟模型的一般框架系统模拟模型的一般框架 步骤 步骤 步骤 步骤 步骤 步骤 步骤 2 )/ 1 ( 1 2 3 4 5 6 7 8 9 10 bcd 离离散散型型随随机机数数的的生生成成 累积分布f概率分布p随机变量x 0.000.1050 0.100.2560 0.350.3570 0.700.2080 0.900.1090 掷骰子0.16 随机数60 2 2 1 )( a s n q q q r e 提前期 下订单收到货 时间 数量 t1t2 库存状况 库存水平 库存状况=库存水平 随机需求下总成本随订购量的变化随机需求下总成本随订购量的变化

11、 22002200 24002400 26002600 28002800 30003000 32003200 34003400 36003600 38003800 40004000 125125150150175175200200225225250250275275300300325325350350375375400400425425450450475475500500 订购量订购量 总成本总成本 1000次模拟平均值 eoq理论公式结果 随机需求下平均缺货次数随订购量的变化随机需求下平均缺货次数随订购量的变化 0.000.00 0.200.20 0.400.40 0.600.60 0.80

12、0.80 1.001.00 1.201.20 1.401.40 1.601.60 1.801.80 2.002.00 125125150150175175200200225225250250275275300300325325350350375375400400425425450450475475500500 订购量订购量 平均缺货次数平均缺货次数 400 450 500 550 600 650 700 750 800 400 450 500 550 600 650 700 750 800 总成本总成本 订购数量q订购数量q 再订购点r再订购点r 再订购点订购量对总成本的模拟结果再订购点订购量

13、对总成本的模拟结果 9000 -10000 8000 -9000 7000 -8000 6000 -7000 5000 -6000 单单服服务务台台等等待待线线系系统统的的时时空空图图 0 0 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 1 10 0 1 11 1 1 12 2 1 13 3 0 02 20 04 40 06 60 08 80 01 10 00 0 1 12 20 0 1 14 40 0 1 16 60 0 1 18 80 0 2 20 00 0 2 22 20 0 2 24 40 0 2 26 60 0 2 28 80 0 3 30 00 0 3 32 20 0 3 34 40 0 3 36 60 0 3 38 80 0 4 40 00 0 时时间间 系系统统人人数数 多多服服务务台台等等待待线线系系统统的的时时空空图图 0 0 1 1 2 2 3 3 4 4 5 5 6