管理系统模拟5仿真结果与系统方案分析.pptVIP

5 仿真结果与系统方案分析5.1简介; 设y11，y12…，y1m是随机变量Y1，Y2，…，Ym单次仿真运行的结果，观测长度为m，进行仿真的时候所用的随机数为 u 11，u 12…，u 1m；u 21，u 22…，u 2m；u n1，u n2…，u nm；（在第j次仿真运行时用的第i个随机数记为uji）。 假定进行了n次独立的重复运行，即每次仿真运行的随机数不同、初始条件相同，每次仿真运行开始时计数器重置。得到以下结果： ;在同一行上的数值来自一次重复运行，不是独立同分布（IID）。在同一列上的数值，y11，y21…，yn1，是变量Y1的观测值满足独立同分布。 系统仿真结果分析就是用多次独立仿真运行的观测值yij（i=1，…，m；j=1，…，n），估计随机变量Y1，Y2，…，Ym的参数。 ; 例5.1：某银行有5位出纳，到达银行的顾客排成一个队列，每位出纳员一次为一位顾客服务。银行上午9点开门，下午5点关门，但继续为在下午5点时已经在银行内的顾客服务完毕。要求确定顾客在银行办理业务需要等待的时间。 银行仿真模型的10次独立重复运行结果如表 ;5.2 仿真结果的瞬态与稳态特征; 不同时刻的随机变量Yi1，Yi2，Yi3，Yi4的瞬时分布的概率密度函数如图 ;系统存在稳态并不表示在某次仿真运行中系统进入稳态后，不同时刻的随机变量取相同的数值，而是进入稳态后不同时刻的随机变量服从相同的分布。 这些随机变量也可能是不独立的。 稳态分布F（y）不依赖于初始条件I，但是瞬时分布收敛于稳态分布的速率会依赖于初始条件I。; 例5.2：单服务台的排队系统 设Di为第i个顾客的排队等待时间。初始队列长度对顾客排队等待时间的影响如图5.2所示.;5.3 系统仿真的类型;5.3.1终止型仿真;例5.3：某飞机制造商接到了生产100架飞机的订单，要求在18个月内交货。公司用仿真方法来确定满足交货期要求的、成本最小的生产方案。 例5.4：某公司只销售一种产品，要确定在120个月内需要维持多少库存。给定初始库存水平，系统仿真的目标为：确定每个月的采购量使得平均每个月的库存维护成本最低。 例5.5：某制造公司每天运行16个小时（分2个班次），当天未完成的工作留在第二天继续进行。用仿真方法确定每个班次的平均产量。;5.3.2 非终止型仿真;5.4 区间估计与置信区间;; 例：正态分布具有未知的均值μ，10个观测结果为1.20，1.50，1.68，1.89，0.95，1.49，1.58，1.55，0.50，1.09。构造μ的具有90%置信度的置信区间。;5.5 终止型仿真的结果分析;5.5.1 固定样本数量法; 例5.8：对于例中的银行，我们希望知道在一天当中顾客的平均排队时间是多少。 由观测结果计算样本均值和方差，构造90％置信度的置信区间，;;;用固定样本数量法进行仿真运行，仿真运行的次数分别为5、10、20，在表5．3中给出了输出分析结果。;5.5.2 序贯法;序贯法进行仿真试验的步骤： （1）预定重复运行的次数 n≥3，建议 n=5；; 例5.10 对于例中的排队系统，要求统计出的顾客平均等待时间的绝对误差小于0.60分钟，用序贯法计算仿真运行次数。 先执行 5次仿真运行，得到样本的方差为 S2（5）=1.0758，置信区间的半长为 1.29，不满足绝对误差要求。估算达到绝对误差要求所需的仿真运行次数， ;5.6 稳态仿真的结果分析;5.6.1重复／删除法;; 采用重复／删除法获得输出参数的点估计和置信区间的方法与固定样本数量法相似。假设仿真运行的总长度为m，预热时段长度为l，独立仿真运行的次数为n。系统输出变量的点估计为， ;例 第1次： Y11 Y12 Y1l Y1 l+1 Y1 l+2 … Y1m 第2次： Y21 Y22 Y2l Y2 l+1 Y2 l+2 … Y2m … 第n次： Yn1 Yn2 Ynl Yn l+1 Yn +2 … Ynm;5.6.2 批均值法;例;1.53 1.66 1.24 2.34 2.00 1.69 2.69 2.86 1.70 2.60;;;假定一公司生产某种产品的速度为 100件／天，生产每批产品的设置成本为100元，每件产品每天的存储成本为0.1元，该产品每天的需求量为25件。若最大生产批量为400件，年初该产品的初始库存量为90件，试用模拟法分析： （1）100天内累积生产量与累积需求量的变动情况。 （2）100天内库存量的变动情况。 （3）生产批量为多少时可使年总成本最低（一年按300天计算）。 若产品每天的需求量为15 ~ 40件的均匀分布，其他不变，则又如何？ 要求： （1）说明利用 Excel