物流仿真第3讲输入数据的分析课件

输入数据的分析2011-9-5第三讲

输入数据的分析

§1数据的收集§2分布的识别§3参数估计§4拟合度检验§5相关性分析

输入数据的分析2011-9-5引言输入数据是仿真实验的动力系统名称典型的输入数据排队系统顾客到达的间隔时间顾客被服务时间的分布库存系统需求顾客的分布顾客需求量的分布物料订货的提前期分布生产系统作业到达的间隔时间作业类型的概率每种作业每道工序服务时间的分布可靠性系统

生产无故障作业时间系统的仿真依靠这些原型系统的运行数据，缺乏这些数据的实验和实验值的提取，仿真也就毫无意义。

输入数据的分析2011-9-5引言收集原始数据

基本统计分布的辨识参数估计

拟合度检验

可信否？否是是输入数据分析的基础，需要分析的经验，对收集的方法、数据需要做预先的设计和估算。因此这是一个关键的、细致的工作。通过统计的数学手段（计数统计、频率分析、直方图制作等），得出统计分布的假设函数（如：正态分布、负指数分布、Erlang分布等）根据统计特征，计算确定系统的假设分布参数。运用统计分布的检验方法，对假设的分布函数进行可信度检验。通常采用的是2检验。输入数据模型确定的基本方法正确输入数据

-1

数据的收集2011-9-5数据收集过程中的注意事项

做好仿真计划，详细规划仿真所需要收集的数据在收集数据过程中要注意分析数据数据的均匀组合收集的数据要满足独立性的要求数据自相关性的检验根据问题的特征，进行仿真的前期研究。分析影响系统的关键因素。从相关事物的观察入手，尽量收集相关的数据。为此可以事先设计好调研表格，并注意不断完善和修改调研方式，使收集的数据更符合仿真对象的数据需要。数据的收集与仿真的试运行是密切相关的，应当是边收集数据、边进行仿真的试运行。然而系统仿真是一项专业性很强的工作，要正确认识“仿真”的含义，抓住仿真研究的关键，避免求全、求精。确信所收集的数据足以确定仿真中的输入分量，而对仿真无用或影响不显著的数据就没有必要去多加收集。针对仿真所收集的各个数据需要进行相关性检验。为了确定在两个变量之间是否存在相关。要建立两个变量的散布图。通过统计方法确定相关的显著性。尽量把均匀数据组合在一组里。校核在相继的时间周期里以及在相继日子内的一时间周期里的数据的均匀性。当校核均匀性时，初步的检验是看一下分布的均值是相同。考察一个似乎是独立的观察序列数据存在自相关的可能性。自相关可能存在于相继的时间周期或相继的顾客中。例如，第i个顾客的服务时间与(i+n)个顾客的服务时间相关。-2

分布的识别2011-9-5直方图对于离散系统的统计分析中，一般用频率统计的分析方法来计算分布函数。其图形描述用的就是直方图。直方图构筑方法取值区间划分水平区坐间标标轴注的计区算间确内定的每发一生数垂直标坐注标频轴数上绘上制的各发个生区频间数绘制直方图-2

分布的识别2011-9-5直方图分组区间数量的选取分组区间的组数依赖于观察次数以及数据的分散或散布的程度。一般分组区间组数近似等于样本量的平方根。即：如果区间太宽（m太小），则直方图太粗或呈短粗状，这样，它的形状不能良好地显示出来。如果区间太窄，则直方图显得凹凸不平不好平滑合适的区间选择（m值）是直方图制作，分布函数分析的基础。-3

参数估计2011-9-5参数估计的作用上一节通过对随机过程的样本值的直方图分析，我们已经得到了随机过程的分布假设，即假设随机过程的概率分布符合某一种标准随机分布。这是一种定性分析的结果。在给定了一种随机分布函数后，需要进一步获取这一分布函数的特征参数，这一标准分布函数的参数需通过参数估计来求得。因此，参数估计在这里是为了对随机分布函数参数求取的一个工具。-3

参数估计2011-9-5样本均值和样本方差

设某一个随机过程X，其n个抽样样本为x1，x2，…，xn，该样本的均值为该样本的方差为如果离散数据已按频数分组，则k是X中不相同数值的个数即分组数，f是X中数值Xj的观察频数-3

参数估计2011-9-5参数估计量分布参数建议使用的估计量泊松指数在（0，b）上的均匀分布正态，仿真中常用的一些分布参数建议值-4

拟合度检验2011-9-5可以证明：02近似服从具有自由度f=k-s-1的2分布。这里s表示由采样统计量所估计的假设分布的参数个数。假设检验作零假设：H0：观察值Xi是一组属于分组分布函数F的独立相同分布的随机变量。若2太大则拒绝H0，若拟合是好的，则期望值2很小。首先划分区间，定义k值计算各组的观察频数计算0查阅2表，得到该值反映的是假设分布的拟合程度。拟合度检验步骤拟合程度的判定-4

拟合度检验2011-9-5指定拟合度的检验我们可以根据拟合度检验的要求，设定一个拟合度的显著性指数，根据设定的显著性指数以及2分布的自由度数f=k-s-1，可以查2表得到，f2。如果则检验未通过，H0不成立。如果则检验通过，H0成立。在应用这个检验时，如果期望的频数太小，将对检验的有效性有所影响。一般情况下区间的个数k宜在30~40以下，并能使最小期望频数Ei≥5。如果Ei值太小，可以把它和相邻分组区间的期望频数相合并，对应的Oi值也应该合并起来，同时每当合并一个单元，k值应该减去1。-5

相关性分析2011-9-5相关性分析系统运行过程中，随机变量有多个，如存在多种因素的影响；系统参数的变化等。这些随机变量之间可能是独立的，也有可能是相互有牵连的，牵连程度的强弱有所不同。需要进行相关性分析。相关性分析的目的：更好地了解系统以及系统随机变量的关联性，更正确地把握问题的关键。相关性分析的方法：通常采用的是回归分析的统计方法-5

相关性分析2011-9-5单变量线性回归假设y的每一个观察值可用下式表示

y=0+1x+

式中是均值为0，方差为2的随机误差。假设存在n对观察值(xi，yi)，i=1，2，……，n，通常采用最小二乘法来估计上式中的yi。设yi=0+1xi

+ii=1，2，……，n，则i=yi-

0-1xi假设是不相关的随机变量。

-5

相关性分析2011-9-5随机变量偏差的平方和为(最小二乘法函数形式)为了使L（偏差）极小，可求出和，并置它们为0，从而可以得到0、1的线性代数方程，既有：-5

相关性分析2011-9-5单变量线性回归的显著性检验检验统计量的构造方法β1的均方误差：在xi处观测值yi与回归值yi之间的误差为均方误差值为也称为回归的剩余方差，它是误差方差的无偏估计量。^-5

相关性分析2011-9-5多变量线性回归X、Y变量可以是多种形式的变量，如，X、Y为非线性变量。-5

相关性分析2011-9-5多变量线性回归假设y={y1,y2,……,ym}T

的由m个变量构成的向量，每一个向量观察值可用下式表示

y=0+1x+

式中={

1,

2,……,

m}T是均值为0。

x={x1,x2,……,xn}Tn个影响观察值的控制变量。式中0={

1,

2,……,

m}T待求的相关系数（常数项）。1=[ij]为mn阶的系数矩阵。-5

相关性分析2011-9-5多变量线性回归为了计算的方便将上述表达形式改写为：y=x

B+

式中={

1,

2,……,

m}T

x={1,x1T,x2T,……,xnT}

n个影响观察值的控制变量。B=[0,

1]为n(m+1)阶待求的系数矩阵。用最小二乘法来估计上式中的y。设yi=B

xi

+