第七章仿真的输出分析.pptVIP

第六章 仿真结果分析与模型校验 引 言 什么是输出分析？ 为什么要进行输出分析？ 输出分析的两种状态 输出分析的输出状态 暂态（终 态） 稳态（非终态） 例 题 某一个通信系统由几个部件加上几个备用部件组成。其中一个分支环节由A、B、C、D四个部分组成，B和C呈并联方式连接。 在系统失效为止的 时间周期TE内考虑系统。 停止事件E定义为 E={A失效，或D失效， 或B与C同时失效} 初始条件为各部件在时刻0都是新的（系统处于理想状态）。 例 题 研究内容：电器元件的平均寿命 研究方法：在相同的实验条件下，进行元件的寿命测量 即：在相同的实验环境下，从时刻0开始测量，一直进行到E事件变真。 结 论 ：这样的仿真我们称其为终态仿真 研究内容的变化：如果对于同样的系统，研究的是系统的特性，如通讯能力、通讯容量等，我们采用的是非终态仿真。 因此终态或非终态仿真是随研究要求的变化而改变。 性能测度估计的方法 假设系统性能可用参数θ(或φ)表示，系统仿真的目的是：通过仿真，希望得到表示系统性能θ(或φ)的值。如何得到或统计得到此值？ 我们可以运用参数的估计方法：既要得到这个值——点估计，又要得到这个值的精度范围——区间估计。区间估计的范围（或长度）是点估计准确度的一个测度。 同样，仿真的数据也有两种： 离散性仿真：仿真输出数据具有离散形式{Y1，Y2，…，Yn}，用来估计θ 连续性仿真：仿真输出数据具有连续形式{Y(t)，0≤t≤T}，用来估计φ 6.1.1重复运行法 一个终态仿真，它在仿真时间区间[0，TE]中运行，并由此得到观察值Y1，Y2，…，Yn。样本量n可以是固定数，也可以是随机变量。 终态仿真的目的是估计 重复运行的数据记录 独立重复运行法 R个样本均值。 ， ，…， 是统计独立的，具有同一分布，并且是θ的无偏估计，于是可以应用经典的置信区间估计的方法。假设做了R次独立的重复运行，用 来计算整个的点估计 。 用 来估计 的方差 则 其100(1-α)%的置信区间为 自由度 f=R-1 稳态仿真的作用 一个仿真模型的单次运行的目的在于估计系统的稳态或长期特征。 设该单次运行得到的观察值是Y1，Y2，…，一般情况下，它是一个自相关时间序列的采样值。所要估计的稳态(或长期)的均值性能测度由下式定义 上式表明：系统模型利用同一统计特征的不同随机数进行的仿真，都将产生样本均值收敛于θ的序列Yi，i=1，2，…，θ的值与初始条件无关。 稳态仿真中初始条件所引起的偏差 稳态仿真运行一般可以分成两段(目的是为了消除初始条件的影响) 第一段从时刻0到时刻T0为初始阶段； 第二段接着从T0到停止时刻T0+ TE为数据收集阶段。 T0的选择是十分重要的，因为用I来表示系统在时刻T0的状态比用时刻0原来的初始条件I。更能反映出系统的稳态行为。(无动态偏差的状况下) 稳态仿真重复运行方法 我们已经看到通过设定T0、TE可以将点估计中初始条件引起的偏差已被减少到可忽略的程度，另外，独立重复运行的方法同样可用来估计点估计值变化范围（适用于无动态偏差），并构造置信区间。 如果，在点估计中有明显的动态偏差（如：系统存在交替出现的偏差），那么采用大量的重复运行来减少点估计值的变化范围，就会导致错误的置信区间。 大量的重复运行会使得置信区间围绕点发生 “偏移” ，使原本围绕着θ变短的置信区间“偏移”到围绕着“错误的点”(θ+b)变短。 在稳态仿真中样本量与准确度的关系 提高在100(1-?)%置信区间内的系统估计性能测度θ的准确度在ε之内的方法主要有2种： 增加重复运行数R 增加运行长度TE 增加重复运行数R提高准确度 置信区间 不等式是均值θ的基于t分布的置信区间，其半长是 增加重复运行数R提高准确度 仿真初期，假设初始仿真样本量R0（已经运行并得到了观察记录），即仿真者最初已做了R0次独立的重复运行。R0次重复运行将用来得到总体方差?2的初始估计S02 。 为符合半长临界值?必须选择更大的样本量R，使得R?R0，且 R是满足R?R0，以及 的正整数。 增加仿真运行的时间长度T0+TE 若按上一种方法计算，需要附加重复运行次数R-R0，那么，我们可按同一比例R/R0把运行长度T0+TE增加到新的运行长度(R/R0)(T0+TE)。于是，要把附加的数据从时刻0删除到时刻(R/R0)T0，这样有更多的数据可用来计算点估计值 。 在稳态仿真中区间估计的批平均值 重复运行法的缺点是，必须在每一次重复运行中删除