第6章仿浩真的输出分析z.pptVIP

* * * * * * * * * * * * * * * * * 若按上一种方法计算，需要附加重复运行次数R-R0，那么，我们可按同一比例R/R0把运行长度T0+TE增加到新的运行长度(R/R0)(T0+TE)。于是，要把附加的数据从时刻0删除到时刻(R/R0)T0，这样有更多的数据可用来计算点估计值 。 增加仿真运行的时间长度T0+TE 增加每次重复运行的总运行长度并删除总运行长度中固定比例T0/( T0+TE)部分的好处是：在点估计中，任何残留的动态偏差将由于附加删除了运行的初始阶段的数据而进一步减少。而该方法可能具有的缺点是：为了继续进行全部R次响应的仿真（从时刻T0+TE运行到(R/R0)( T0+TE)），必须记下模型在时刻T0+TE的状态，以便让模型重新开始并运行所需的附加时间。 在稳态仿真中区间估计的批平均值 重复运行法的缺点是，必须在每一次重复运行中删除数据，在某种意义上，把数据删除掉是一种浪费或至少损失了信息。 批平均值法是采用一次长运行，在其中数据仅需删除一次，输出数据被放到少(比如说5到10批)而大(数据量大)的批次中，并且分析员用这些似乎独立的少数批平均值进行统计分析。Law(1977)发现，对于固定的总样本量，最好利用很少量的批数但尽可能长的运行长度，他建议用5批。 I0 I T0 T0+ TE 指定初始条件 稳态初始条件 长度为T0的初始段 长度为TE的数据收集段 在稳态仿真中区间估计的批平均值 设5个批平均值Y11，Y12，Y13，Y14，Y15，假设它们是统计独立的(虽然批平均值并不独立)，在这种假设下，可用通常的方法进行计算。 点估计是 样本方差是 标准偏差是 95%的置信区间是 例 一机床加工零件的仿真系统，仿真的目的是分析机床的利用率和一个工作日内机床每加工一个零件的平均时间。在相同的初始条件下经过4次独立的仿真运行，得到结果如下表所示。 举例 试计算在α=0.05 时机床利用率v的置信区间和零件平均加工时间t的置信区间。 解：机床利用率v的样本均值为： 举例 机床利用率v的样本方差为： 举例 对于置信水平95%，查t-分布表，求双边临界值 。 下表是不同自由度下对应于不同显著水平t-分布的双边临界值。 举例 查t-分布表得 故 即利用率v的置信区间为： 同理，零件加工时间t的置信区间为： 举例 例 考察一个稳态仿真系统，对系统进行R次独立重复仿真运行，每次运行得到n个输出点，设Yij为第i次运行的第j 个观测值，如下表所示。 举例 每次重复运行都是对E[Y]估计的一个样本，对于每次运行剔除前d个观测值（初始阶段观察值），然后计算均值。对于第 r 次重复运行，其均值的计算公式为： 即每次运行都删去初始 d 个仿真数据，以消除初始条件的影响；同时每次仿真运行均用不同的随机数流，以保证输出数据的独立性，在t=0 时置相同的初始条件，则 举例 是独立同分布的随机样本，总的点估计为： 则 就是E[Y]的近似无偏估计。 相应的样本方差为： 举例 由此，可确定稳态仿真的置信区间为： 重复-删除法思路非常简洁，在实际常用，但存在一些缺点： （1）难以正确选择要删除的观察值数目； （2）数据的使用效率低，每次都剔除前d 个数据； （3）必须通过人为中断仿真运行来收集数据。 举例 * * * * * * * * * * * * * * * * * * * * *An Introduction to Database System * * 第六章 仿真的输出分析 §6-1 引言 §6-2 性能测度及其估计 §6-3 终态仿真的输出分析 §6-4 稳态仿真的输出分析 什么是输出分析？ 为什么要进行输出分析？ 输出分析的两种状态 系 统 （结构数量是确定的） （结构参数是随机的） 输 入 （参数是随机的） 输 出 ？ 确定的输入激励一个确定的系统，得到的输出就是一个确定的输出。通过一次确定的仿真便可得出解。 随机的输入激励一个随机的系统，得到的输出是……? 输出的表达形式如何？ 需要经过多少次的仿真才能说明输出结果？ 输出分析的目的在于预测一个系统的性能，或比较两个或多个不同系统的设计性能。估计系统的性能参数，以及性能参数估计的有效范围。用仿真统计得到的?作为观察值的估计量。统计得到的方差S2 就是估计量的偏差范围。或确定出达到给定精度所需的观察次数。 在离散事件仿真中，大多数仿真输出数据呈现出自相关的特征，即：前面的输出往往会影响到后面的输出数据。 如：库存系统中的初期库存、生产系统中的初始状态、排队系统中初始排队状态和初始服务状态等。 §6-1 引言 输出分析的输出状态 暂态（终 态） 稳态（非终态）