《第10章仿真的输出分析课件》由会员分享,可在线阅读,更多相关《第10章仿真的输出分析课件(29页珍藏版)》请在金锄头文库上搜索。
1、第七章 仿真的输出分析,7-1 引言 7-2 性能测度及其估计 7-3 终态仿真的输出分析 7-4 稳态仿真的输出分析,引 言,什么是输出分析? 为什么要进行输出分析? 输出分析的两种状态,确定的输入激励一个确定的系统,得到的输出就是一个确定的输出。通过一次确定的仿真便可得出解。,随机的输入激励一个随机的系统,得到的输出是? 输出的表达形式如何? 需要经过多少次的仿真才能说明输出结果?,输出分析的目的在于预测一个系统的性能,或比较两个或多个不同系统设计的性能。估计系统的性能参数,以及性能参数估计的有效范围。用仿真统计得到的作为观察值的估计量。统计得到的方差S2 就是估计量的偏差范围。或确定出达
2、到给定精度所需的观察次数。,在离散事件仿真中,大多数仿真输出数据呈现出自相关的特征,即:前面的输出往往会影响到后面的输出数据。 如:库存系统中的初期库存、生产系统中的初始状态、排队系统中初始排队状态和初始服务状态等。,接下页,输出分析的输出状态,暂态(终 态) 稳态(非终态),终态仿真就是指在某个持续时间TE之内系统的仿真,这里E是停止仿真的一个指定的事件,这样被仿真系统在指定初始条件下于时刻0“打开”,并在停止时刻TE“闭合”。终态系统常被用来研究系统的固有特性,研究系统在初始条件作用下的响应。,非终态系统是指系统在持续循环运行时间内,前一时间结束的仿真结果影响到后一时间的仿真条件。非终态系
3、统是连续运行的系统,至少在很长一段时期内运行。稳态系统仿真常被用来研究系统对外界条件变化的响应能力。通常稳态系统的响应与系统的初始状态无关。,例 题,某一个通信系统由几个部件加上几个备用部件组成。其中一个分支环节由A、B、C、D四个部分组成,B和C呈并联方式连接。 在系统失效为止的 时间周期TE内考虑系统。 停止事件E定义为 E=A失效,或D失效, 或B与C同时失效 初始条件为各部件在时刻0都是新的(系统处于理想状态)。,例 题,研究内容:电器元件的平均寿命 研究方法:在相同的实验条件下,进行元件的寿命测量 即:在相同的实验环境下,从时刻0开始测量,一直进行到E事件变真。 结 论 :这样的仿真
4、我们称其为终态仿真 研究内容的变化:如果对于同样的系统,研究的是系统的特性,如通讯能力、通讯容量等,我们采用的是非终态仿真。 因此终态或非终态仿真是随研究要求的变化而改变。,性能测度估计的方法,假设系统性能可用参数(或)表示,系统仿真的目的是:通过仿真,希望得到表示系统性能(或)的值。如何得到或统计得到此值? 我们可以运用参数的估计方法:既要得到这个值点估计,又要得到这个值的精度范围区间估计。区间估计的范围(或长度)是点估计准确度的一个测度。 同样,仿真的数据也有两种: 离散性仿真:仿真输出数据具有离散形式Y1,Y2,Yn,用来估计 连续性仿真:仿真输出数据具有连续形式Y(t),0tT,用来估
5、计,点估计(离散随机变量),基于数据Y1,Y2,Yn的的点估计定义为 式中 是基于样本量n的样本均值。如果 的数学期望是,即 则 是的无偏估计。 称为离散仿真系统性能的平均测度。,点估计(连续随机变量),基于数据Y(t),0tT的的点估计定义为 式中 T是仿真的运行长度, 称为Y(t)在0,T上的时间平均值。 称为连续仿真系统性能的平均测度。,=,是的无偏估计。,区间估计,首先,确定在无偏估计下,估计点估计 (或 )的方差。 令 表示点估计 的真实方差 令 表示基于数据Y1,Y2,Yn的方差 的估计值。 假设 B称为在方差估计中的偏差系数。 如果 是近似无偏的(B1.0),那么取统计量,根据数
6、理统计的定理,可知,统计量,当B=1, 为点估计方差 的一个无偏估计。,区间估计,为了使达到近似的100(1-)%置信区间,必须满足 P(| t | t,)= f=n-1为 t检验的自由度数。n为样本数。 通过 t分布的标准统计表,可以查得自由度为f,满足置信区间的t,值 根据,问题:如何求得上式中点估计方差的近似无偏估计 ?,影响方差估计值 的因素,实验观测值的相关性将影响到方差估计值 ,及其区间估计的精度,区间估计,如果Y1,Y2,Yn是统计独立的观察值 由点估计定义式计算 ,然后计算样本方差 当Yi是独立的、相同分布时,那么样本方差S2就是总体方差2=var(Yi)(对所有i=1,2,n
7、,皆为常数)的无偏估计。由于 的方差为 ,那么2( )的无偏估计具有 f =n-1的自由度,只要 点估计是无偏的,那么置信区间 便是近似正确的,称为点估计 的标准偏差。标准偏差是点估计准确度的测度。,区间估计,如果Y1,Y2,Yn不是统计独立的观察值(相关的) 那么 是点估计 真实方差 的有偏估计。当Y1,Y2,Yn是从一次单独运行中得到的输出观察序列时,Y1,Y2,Yn是一个自相关序列(称之为时间序列)。 为了定量表示自相关对方差估计的影响,假设时间序列Y1,Y2,Yn是均值为的协方差平稳的时间序列(稳态过程),对协方差平稳时间序列Y1,Y2,Yn用 表示滞后k的协方差。,区间估计,依据协方
8、差的定义 是Y的自相关函数 又其对应的是平稳过程(或近似平稳过程),协方差函数与i无关 协方差函数是一个对称的函数 当k=0时,函数取最大值0 是序列的总体方差 定义 (-1 k 1) k=1,2,n 当对所有的k, k 0,序列是正自相关的。大部分仿真序列是正相关的。 当对所有的k, k 0,序列是负自相关的。,对具有式 所定义的采祥均值的平稳时间序列来讲, 的方差由下式给出 如果Yi是独立观察值,那么当k=1,2,3,时, =0, ,上式为 将其表示为如同独立观察值的形式 ,其期望值由下式给出,与独立观察值的估计区间比较,其取值范围的表示方式与B的值有关,B有可能大于1或小于1,如果相关系
9、数 大体上是正的,那么c1,这样n/c1,在这种情况下,用 来估计 则偏高。由上述不等式计算出的名义的100(1-)%置信区间将是太大,它的真实置信度将大于(1-)。换句话说,点估计的真实精度将比用它的方差估计 表示的精度来得高,即 这个误差不像第一种情况的误差那么严重,因为当呈现负相关时,估计值 的准确度将倾向于比用 方差估计表示的准确度更好。,独立重复运行法,一个终态仿真,它在仿真时间区间0,TE中运行,并由此得到观察值Y1,Y2,Yn。样本量n可以是固定数,也可以是随机变量。 终态仿真的目的是估计,设仿真共重复R次,每次运行都利用不同的随机数流和独立选择的初始条件(也包含所有含有相同初始
10、条件的情况)。令Yri是第r次重复运行的第i次观察i=1,2,n,以及r=1,2,R。,重复运行的数据记录,当固定r时,Yr1,Yr2,是自相关序列,但对不同的响应r和s,rs,Yri和Ysi是统计独立的。对每一次运行r,其样本均值 为 r=1,2,R,独立重复运行法,R个样本均值。 , , 是统计独立的,具有同一分布,并且是的无偏估计,于是可以应用经典的置信区间估计的方法。假设做了R次独立的重复运行,用 来计算整个的点估计 。 用 来估计 的方差 则 其100(1-)%的置信区间为 自由度 f=R-1,量 = 称为点估计 的标准偏差,它的大小反映了的点估计的 准确度。当R增加时,标准误差 倾
11、向于变得越来越小而趋于0。 即:独立重复运行次数的增加,可以减小标准误差的值,也就是使置信区间缩小,提高了性能测度的精度。,稳态仿真的作用,一个仿真模型的单次运行的目的在于估计系统的稳态或长期特征。 设该单次运行得到的观察值是Y1,Y2,一般情况下,它是一个自相关时间序列的采样值。所要估计的稳态(或长期)的均值性能测度由下式定义 上式表明:系统模型利用同一统计特征的不同随机数进行的仿真,都将产生样本均值收敛于的序列Yi,i=1,2,的值与初始条件无关。,稳态仿真中初始条件所引起的偏差,稳态仿真运行一般可以分成两段(目的是为了消除初始条件的影响) 第一段从时刻0到时刻T0为初始阶段; 第二段接着
12、从T0到停止时刻T0+ TE为数据收集阶段。 T0的选择是十分重要的,因为用I来表示系统在时刻T0的状态比用时刻0原来的初始条件I。更能反映出系统的稳态行为。(无动态偏差的状况下),暂 态,稳 态,系统在时间T0的状态I是随机变量,系统在此点已达到近似稳态,指在时刻T0的系统状态的概率分布充分接近稳态概率分布,从而使响应变量点估计的动态偏差可以忽略不计。,数据收集阶段的长度TE 应足够长以保证得到充分准确的系统稳态行为的估计。,稳态仿真重复运行方法,我们已经看到通过设定T0、TE可以将点估计中初始条件引起的偏差已被减少到可忽略的程度,另外,独立重复运行的方法同样可用来估计点估计值变化范围(适用
13、于无动态偏差),并构造置信区间。 如果,在点估计中有明显的动态偏差(如:系统存在交替出现的偏差),那么采用大量的重复运行来减少点估计值的变化范围,就会导致错误的置信区间。 大量的重复运行会使得置信区间围绕点发生 “偏移” ,使原本围绕着变短的置信区间“偏移”到围绕着“错误的点”(+b)变短。,在稳态仿真中样本量与准确度的关系,提高在100(1-)%置信区间内的系统估计性能测度的准确度在之内的方法主要有2种: 增加重复运行数R 增加运行长度TE,增加重复运行数R提高准确度,置信区间 不等式是均值的基于t分布的置信区间,其半长是,S是样本标准偏差,R是重复运行次数。,=,设给定一个准确度临界值,我
14、们希望用 来估计(具有准确度)这个事件的发生具有较高的概率,比如说至少为 。即:需要取足够大的样本量R来满足,增加重复运行数R提高准确度,仿真初期,假设初始仿真样本量R0(已经运行并得到了观察记录),即仿真者最初已做了R0次独立的重复运行。R0次重复运行将用来得到总体方差2的初始估计S02 。 为符合半长临界值必须选择更大的样本量R,使得RR0,且 R是满足RR0,以及 的正整数。,增加重复运行数R提高准确度(算法),由于t/2,R-1 RZ/2(这里Z/2为R的t/2,可查表得到),那么对R的初始估计量 其中:Z(x)=1-norm(x) 在确定出最终样本量R之后,再做R-R0次附加的观察,
15、并得到 的100(1- )%置信区间,和S是基于所有R次重复运行而计算出来的,若按上一种方法计算,需要附加重复运行次数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批。,在稳态仿真中区间估计的批平均值,设5个批平均值Y11,