第3章静态仿真11.pptVIP

3.7 保险费模型 赔付率模型四种方案的仿真结果 方案 1 2 3 4 平均损失占保险费的比例rc 5.0 10.0 20.0 80.0 索赔的概率pc 0.16 0.08 0.04 0.01 仿真重复次数 500 500 500 500 均值 0.8111 0.8201 0.8099 0.7809 标准差 0.3057 0.4420 0.6068 1.2179 最小值 0.1496 0.0497 0.0000 0.0000 最大值 1.9745 2.5801 4.1687 10.5220 赔付率1.0 0.2440 0.2960 0.3440 0.2740 赔付率的标准误差 0.0268 0.0387 0.0532 0.1068 制表软件的扩展 加载项 PopTools Numerical Technologies tools @RISK Crystal Ball R 外部和内部工具 VBA 统计和数学工具 PopTools 免费软件网址 http://www.cse.csiro.au/poptools/index.htm 社会生态人口建模 产生随机数 产生随机变量 数据分析工具 其他功能 @RISK Commercial software Website Random number generation Random variate generation Data analysis and display of results Distribution fitting Optimization 不适宜电子制表软件仿真的问题 数据结构精细的问题 离散事件仿真。 树,列表, 多维数组。 要求非常精确结果问题 长时间仿真运行。 具有随机迭代次数的迭代计算问题。 应用中应注意的问题 多数电子制表软件的随机数发生器设计较差且无正式文档说明，不能用于实际问题分析。 很难实现复杂的数据分析，例如批均值 制表软件运算相对较慢，不适宜大规模问题的建模与仿真。 电子制表软件仿真参考文献： 1.Andrew F. Seila：SPREADSHEET SIMULATION. Proceedings of the 2006 Winter Simulation Conference. 2.Seila, A. F., V. Ceric, and P. Tadikamalla. 2003. Applied simulation modeling. * * * * * * * * * 3.4 蒙特卡罗方法 3.4.2 Monte－Carlo方法的收敛性与样本量 3.4 蒙特卡罗方法 3.4.3 估计π值的模型 园的面积 设圆的半径和正方形的边长都是1，则四分之一圆的面积为π/4，正方形的面积为1。 随机试验：在正方形中随机选择一点(x,y)，正方形中所有的点被选中的概率相同，则选择的点落在四分之一圆之内的概率p就等于四分之一圆面积与正方形面积之比，即p=π/4，估计出p值乘以4得到π的估计值。 0 1 1 3.4 蒙特卡罗方法 3.4.3 估计π值的模型 nq表示n次实验中点落在四分之一圆内的次数。p的估计值为： 利用仿真估计π值的算法总结如下： 1、初始化n（仿真次数），初始化nq和计数器count，令nq=0，count=0； 2、重复：生成变量X，Y；X和Y都服从均匀分布U(0,1)； 如果X2+Y21，则nq=nq+1；count=count+1；直到count=n； 3、计算π的估计值： 3.4 蒙特卡罗方法 3.4.4 结果分析 3.5评价决策—单周期库存模型 考虑单周期库存模型，购货量是随机需求，需求分布已知。例如，发出订单，订购一定数量的“饮料”供应足球比赛。需求数据有很大不确定性，根据历史数据，设需求D服从均值为5.0的指数分布。如果订单的订购量大于需求，过剩的每箱成本为Ce=$40；如果存货量小于需求，缺货成本每箱为Cs=$160。设D表示需求，为随机变量；x表示订货的箱数，是决策变量。 仿真这一模型，估计预期成本，确定订货量x，使成本最小。 3.5.1 单周期库存模型 =Mean_Demand*(-LN(RAND())) =IF(DemandSupply,Ce*(Supply-Demand),Cs*(Demand-Supply)) 库存模型 重复计算 G15-G18 G15+G18 ) 975 . 0 ( ) 2 / 1 ( NORMSINV Z = - a 产生多方案仿真结果步骤 在电子表格中选定方形区域，本例中I14：K22； 在I列（I15：I22）输入订货量x——决策变量值 （5,6,7,8,9,10,11,12）； 在“输入引用