仿真的输出分析

第六章仿真结果分析与模型校验引言什么是输出分析？为什么要进行输出分析？输出分析的两种状态系统（结构数量是确定的）（结构参数是随机的）输入（参数是随机的）输出？确定的输入激励一个确定的系统，得到的输出就是一个确定的输出。通过一次确定的仿真便可得出解。随机的输入激励一个随机的系统，得到的输出是……?输出的表达形式如何？需要经过多少次的仿真才能说明输出结果？输出分析的目的在于预测一个系统的性能，或比较两个或多个不同系统设计的性能。估计系统的性能参数，以及性能参数估计的有效范围。用仿真统计得到的‘作为观察值的估计量。统计得到的方差S2就是估计量的偏差范围。或确定出达到给定精度所需的观察次数。在离散事件仿真中，大多数仿真输出数据呈现出自相关的特征，即：前面的输出往往会影响到后面的输出数据。如：库存系统中的初期库存、生产系统中的初始状态、排队系统中初始排队状态和初始服务状态等。接下页输出分析的输出状态暂态（终态）稳态（非终态）终态仿真就是指在某个持续时间TE之内系统的仿真，这里E是停止仿真的一个指定的事件，这样被仿真系统在指定初始条件下于时刻0“打开”，并在停止时刻TE“闭合”。终态系统常被用来研究系统的固有特性，研究系统在初始条件作用下的响应。非终态系统是指系统在持续循环运行时间内，前一时间结束的仿真结果影响到后一时间的仿真条件。非终态系统是连续运行的系统，至少在很长一段时期内运行。稳态系统仿真常被用来研究系统对外界条件变化的响应能力。通常稳态系统的响应与系统的初始状态无关。例题某一个通信系统由几个部件加上几个备用部件组成。其中一个分支环节由A、B、C、D四个部分组成，B和C呈并联方式连接。在系统失效为止的时间周期TE内考虑系统。停止事件E定义为E={A失效，或D失效，或B与C同时失效}初始条件为各部件在时刻0都是新的（系统处于理想状态）。ABCD图7.1通讯系统例题•研究内容：电器元件的平均寿命•研究方法：在相同的实验条件下，进行元件的寿命测量即：在相同的实验环境下，从时刻0开始测量，一直进行到E事件变真。•结论：这样的仿真我们称其为终态仿真•研究内容的变化：如果对于同样的系统，研究的是系统的特性，如通讯能力、通讯容量等，我们采用的是非终态仿真。因此终态或非终态仿真是随研究要求的变化而改变。ABCD图7.1通讯系统性能测度估计的方法假设系统性能可用参数θ(或φ)表示，系统仿真的目的是：通过仿真，希望得到表示系统性能θ(或φ)的值。如何得到或统计得到此值？我们可以运用参数的估计方法：既要得到这个值——点估计，又要得到这个值的精度范围——区间估计。区间估计的范围（或长度）是点估计准确度的一个测度。同样，仿真的数据也有两种：•离散性仿真：仿真输出数据具有离散形式{Y1，Y2，…，Yn}，用来估计θ•连续性仿真：仿真输出数据具有连续形式{Y(t)，0≤t≤T}，用来估计φ§6-1终态(暂态)仿真的输出分析6.1.1重复运行法•一个终态仿真，它在仿真时间区间[0，TE]中运行，并由此得到观察值Y1，Y2，…，Yn。样本量n可以是固定数，也可以是随机变量。•终态仿真的目的是估计niiYnE11设仿真共重复R次，每次运行都利用不同的随机数流和独立选择的初始条件(也包含所有含有相同初始条件的情况)。令Yri是第r次重复运行的第i次观察i=1，2，…，n，以及r=1，2，…，R。6.1终态仿真结果分析§6-1终态(暂态)仿真的输出分析2008.11Y1n1i112y11y12y1iy1n1......Y2n2i12y21y22y2iy2n1......Yrnri12yr1yr2yriyrn1......重复运行的数据记录当固定r时，Yr1，Yr2，…是自相关序列，但对不同的响应r和s，r≠s，Yri和Ysi是统计独立的。对每一次运行r，其样本均值为r=1，2，…，RrˆrnirrirnY1ˆ§6-1终态(暂态)仿真的输出分析2008.11独立重复运行法1ˆ2ˆRˆRrrR1ˆ1ˆR个样本均值。，，…，是统计独立的，具有同一分布，并且是θ的无偏估计，于是可以应用经典的置信区间估计的方法。假设做了R次独立的重复运行，用来计算整个的点估计。用来估计的方差则其100(1-α)%的置信区间为自由度f=R-1RS22ˆˆˆRrrRRRS1222ˆˆ11ˆˆˆˆˆˆˆˆ,2,2fftt量=称为点估计的标准偏差，它的大小反映了θ的点估计的准确度。当R增加时，标准误差倾向于变得越来越小而趋于0。即：独立重复运行次数的增加，可以减小标准误差的值，也就是使置信区间缩小，提高了性能测度的精度。ˆˆˆˆ2ˆˆˆˆ§6-1终态(暂态)仿真的输出分析2008.116.1.2序贯程序法§6-2稳态仿真的输出分析2008.11稳态仿真的作用一个仿真模型的单次运行的目的在于估计系统的稳态或长期特征。设该单次运行得到的观察值是Y1，Y2，…，一般情况下，它是一个自相关时间序列的采样值。所要估计的稳态(或长期)的均值性能测度由下式定义上式表明：系统模型利用同一统计特征的不同随机数进行的仿真，都将产生样本均值收敛于θ的序列Yi，i=1，2，…，θ的值与初始条件无关。niinYn11lim§6-2稳态仿真的输出分析2008.11稳态仿真中初始条件所引起的偏差•稳态仿真运行一般可以分成两段(目的是为了消除初始条件的影响)–第一段从时刻0到时刻T0为初始阶段；–第二段接着从T0到停止时刻T0+TE为数据收集阶段。T0的选择是十分重要的，因为用I来表示系统在时刻T0的状态比用时刻0原来的初始条件I。更能反映出系统的稳态行为。(无动态偏差的状况下)暂态稳态I0IT0T0+TE指定初始条件稳态初始条件长度为T0的初始段长度为TE的数据收集段系统在时间T0的状态I是随机变量，系统在此点已达到近似稳态，指在时刻T0的系统状态的概率分布充分接近稳态概率分布，从而使响应变量点估计的动态偏差可以忽略不计。数据收集阶段的长度TE应足够长以保证得到充分准确的系统稳态行为的估计。§6-2稳态仿真的输出分析2008.11稳态仿真重复运行方法我们已经看到通过设定T0、TE可以将点估计中初始条件引起的偏差已被减少到可忽略的程度，另外，独立重复运行的方法同样可用来估计点估计值变化范围（适用于无动态偏差），并构造置信区间。如果，在点估计中有明显的动态偏差（如：系统存在交替出现的偏差），那么采用大量的重复运行来减少点估计值的变化范围，就会导致错误的置信区间。大量的重复运行会使得置信区间围绕点发生“偏移”，使原本围绕着θ变短的置信区间“偏移”到围绕着“错误的点”(θ+b)变短。§6-2稳态仿真的输出分析2008.11在稳态仿真中样本量与准确度的关系•提高在100(1-)%置信区间内的系统估计性能测度θ的准确度在ε之内的方法主要有2种：–增加重复运行数R–增加运行长度TE§6-2稳态仿真的输出分析2008.11增加重复运行数R提高准确度置信区间不等式是均值θ的基于t分布的置信区间，其半长是ˆˆˆˆˆˆ,2,2ffttˆˆ1,2RtS是样本标准偏差，R是重复运行次数。ˆˆRS=设给定一个准确度临界值ε，我们希望用来估计θ(具有准确度ε)这个事件的发生具有较高的概率，比如说至少为。即：需要取足够大的样本量R来满足ˆ11ˆP§6-2稳态仿真的输出分析2008.11增加重复运行数R提高准确度仿真初期，假设初始仿真样本量R0（已经运行并得到了观察记录），即仿真者最初已做了R0次独立的重复运行。R0次重复运行将用来得到总体方差2的初始估计S02。为符合半长临界值必须选择更大的样本量R，使得RR0，且R是满足RR0，以及的正整数。RStR01,2201,2StRR§6-2稳态仿真的输出分析2008.11若按上一种方法计算，需要附加重复运行次数R-R0，那么，我们可按同一比例R/R0把运行长度T0+TE增加到新的运行长度(R/R0)(T0+TE)。于是，要把附加的数据从时刻0删除到时刻(R/R0)T0，这样有更多的数据可用来计算点估计值。增加仿真运行的时间长度T0+TE初始阶段数据采集阶段初始阶段数据采集阶段0T00(R/R0)T0(R/R0)(T0+TE)T0+TE增加每次重复运行的总运行长度并删除总运行长度中固定比例T0/(T0+TE)部分的好处是：在点估计中，任何残留的动态偏差将由于附加删除了运行的初始阶段的数据而进一步减少。而该方法可能具有的缺点是：为了继续进行全部R次响应的仿真（从时刻T0+TE运行到(R/R0)(T0+TE)），必须记下模型在时刻T0+TE的状态，以便让模型重新开始并运行所需的附加时间。§6-2稳态仿真的输出分析2008.11在稳态仿真中区间估计的批平均值重复运行法的缺点是，必须在每一次重复运行中删除数据，在某种意义上，把数据删除掉是一种浪费或至少损失了信息。批平均值法是采用一次长运行，在其中数据仅需删除一次，输出数据被放到少(比如说5到10批)而大(数据量大)的批次中，并且分析员用这些似乎独立的少数批平均值进行统计分析。Law(1966)发现，对于固定的总样本量，最好利用很少量的批数但尽可能长的运行长度，他建议用5批。I0IT0T0+TE指定初始条件稳态初始条件长度为T0的初始段长度为TE的数据收集段§6-2稳态仿真的输出分析2008.11在稳态仿真中区间估计的批平均值设5个批平均值Y11，Y12，Y13，Y14，Y15，假设它们是统计独立的(虽然批平均值并不独立)，在这种假设下，可用通常的方法进行计算。•点估计是•样本方差是•标准偏差是•95%的置信区间是51151jjYY45251212YYSjj5.SYesYestYYestY..4,025.04,025.0§6-2稳态仿真的输出分析2008.116.2.2稳态序贯法6.2.3再生法