设计方案比较课件.ppt

1、第十一章 设计方案的比较与评价 11.1 引言 11.2 两个系统设计方案的比较 11.3 k个系统设计方案之间的比较 物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-12211.1 引言如何运用仿真对系统进行评价 系统参数的随机性对系统输出结果的影响 系统特征参数的确定与选择(系数、运行规则等)系统特征参数的评价 物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-12311.2 两个系统设计方案的比较生产系统的设计方案 系统的设备形式、数量、参数；系统作业的工艺流程；系统

2、的控制方法（手动、半自动、全自动）系统的布置形式；（包括生产物料的流动方式）系统的运行策略（库存保证、订单生产、供应链生产）；生产调度策略；生产系统的人力规划；物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-12411.2 两个系统设计方案的比较系统的比较系统的比较是基于系统的同一参数（设计参数、运行规则等同一定义下的系统特征）。这一（或这些）参数在系统的重复运行中可以得到的输出数据（可观测的）。对于两个系统的设计方案进行比较，可用i(i=1，2)来表示系统i的性能（系统均值性能）。如果是稳态仿真，保证i的估计是近似无偏的。仿真实验

3、的目标是要获得均值性能之间的差别，即仿真实验的目标是要获得均值性能之间的差别，即12的点估计及其区间估计。的点估计及其区间估计。物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-12511.2 两个系统设计方案的比较系统仿真的关键参数稳态仿真的关键参数有下列几个：J 仿真模型的稳态运行时间TEJ 模型的重复运行次数Ri 系统 i 的第 r 次重复运行产生均值性能测度i的一个估计Yri。假设估计值Yri是(至少近似是)无偏的，那么 111,2,1,RrYEr222,2,1,RrYEr 物流工程与管理系物流工程与管理系Logistics

4、Engineering&Management2023-5-12611.2 两个系统设计方案的比较系统性能的比较计算两个性能测度之间的差别12的置信区间，可用来回答以下两个问题：均值差别有多大，以及均值差别的估计有多准确?两个系统之间有显著的差别吗?12 物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-12711.2 两个系统设计方案的比较系统性能比较的三种可能 如果12的置信区间绝大部分在零的左侧，那么120或等价地10或等价地12的假设便有强的证据。如果12的置信区间包含零点，那么，根据现有的数据还没有强的统计证据表明一个系统设计

5、方案优于另一个。121212 物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-12811.2 两个系统设计方案的比较系统性能参数比较的置信区间形式 对12的置信区间有下列形式 21,221.YYestYYf2121.YYesYY ，iY是系统i在所有重复运行上的样本均值性能测度 iRrriiiYRY11f 是相应于方差估计的自由度，是在自由度为 f 的 t 分布中 点处的值s.e.()表示指定的点估计的标准偏差。ft,2%21100 物流工程与管理系物流工程与管理系Logistics Engineering&Management20

6、23-5-129例题：正态分布的标准偏差和自由度计算 具有相等方差的独立采样具有相等方差的独立采样 具有不相等方差的独立采样具有不相等方差的独立采样 相关采样相关采样 物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-1210具有相等方差的独立采样具有相等方差的独立采样独立采样是指用不同的且独立的随机数流来仿真两个系统。这意味着Yr1，r=1，2，R1与Yr2，r=1，2，R2是统计独立的。于是样本均值 的方差由下式给出：iYiiiriiRRYY2varvari=1，2利用独立采样的性质，与 是统计独立的，于是 1Y2Y2221212

7、121varvarvarRRYYYY 物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-1211具有相等方差的独立采样具有相等方差的独立采样如果两次独立采样的方差相等，均值性能差别的点估计是样本方差的无偏估计 依据方差相等条件 ，则 的联合估计由下式给出：它具有f=R1+R22个自由度。那么12的置信区间表达式的标准偏差为i=1，222212121YY2121221111iiRrriiRririiiYRYRYYRSii222212211212222112RRSRSRSp212111.RRSYYesp 物流工程与管理系物流工程与管理系L

8、ogistics Engineering&Management2023-5-1212具有不相等方差的独立采样具有不相等方差的独立采样 独立采样是指用不同的且独立的随机数流来仿真两个系统。这意味着Yr1，r=1，2，R1与Yr2，r=1，2，R2是统计独立的。于是样本均值 的方差由下式给出：iYiiiriiRRYY2varvari=1，2利用独立采样的性质，与 是统计独立的，于是 1Y2Y2221212121varvarvarRRYYYY 物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-1213具有不相等方差的独立采样具有不相等方差的

9、独立采样如果 ，那么12的 近似置信区间w为：均值性能差别的点估计是样本方差的无偏估计点估计的标准偏差自由度 f 的近似计算式%11002121YY2121221111iiRrriiRririiiYRYRYYRSii22212121.RSRSYYes1122222121212222121RRSRRSRSRSf2221 物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-1214相相 关关 采采 样样 相关采样指的是对每一次重复运行，利用相同的随机数来仿真两个系统。因此系统的仿真次数R1和R2相等，为了表达的方便，假设R1=R2=R。对每

10、个第r次重复运行，两个估计Yr1和Yr2不再是独立的，而是相关的。由于对任意两次不同的重复运行利用的是独立的随机数流，那么Yr1和Yr2之间仍然是独立的rs。利用相关采样的目的是让Yr1与Yr2产生正相关，并从而达到使均值差 的点估计的方差减小的目的。21YY 物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-1215相相 关关 采采 样样 方差的一般表达形式12为Yr1与Yr2之间的相关系数，与r无关。令相关采样的方差为Vcor 令独立采样的方差(设R1=R2=R)称之为Vind RRRYYYYYY211222212121212,c

11、ov2varvarvarVcorVindR21122=-如果相关采样是正相关，那么12将是正的，于是 indcorVV较小的方差意味着基于相关采样的估计更为准确。物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-1216相关采样数据的置信区间计算相关采样数据的置信区间计算 令Dr=Yr1Yr2，于是Dr(r=1，2，R)是独立的、具有相同分布的随机样本，其样本均值 样本方差 它具有自由度 f=R1。对12的置信区间估计中的标准偏差为 RrrDRD11式中 21YYD2121221111DRDRDDRSRrrRrrD RSYYesDes

12、D21.当工作正常时(即120)，在给定样本量下，相关采样产生的置信区间要比独立采样所产生的来得短。物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-121711.3 k个系统设计方案之间的比较多系统设计方案的比较方法多系统设计方案的比较方法较多，主要可以分为两种v 固定样本量法 v 时序采祥(或多阶段采样)法 1.预先确定出仿真的样本量(包括运行长度TE以及重复运行次数R)2.通过假设检验和/或置信区间作出论断固定样本量法的优点是在进行仿真实验前，花费计算机机时是已知的，适用于机时有限或作些初步研究。固定样本量法的主要缺点是不可能有

13、强有力的结论，例如，置信区间对实际应用来说可能太宽或假设检验可能导致不拒绝零假设。需要收集越来越多的数据一直到估计值达到预先给定的准确度，或者一直到从几个可供选用的假设中选出一个为止。物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-121811.3 k个系统设计方案之间的比较多系统设计方案的比较方法假设要计算总共C个置信区间，其中第i个置信区间具有的置信系数1-i。令第i个置信区间是一个命题，称之为Si，对给定的一组数据它可以为真或假，为真的概率为1-i。那么Bonferroni不等式是ECjj11CjjE1称为总误差概率。等价于

14、P(一个或多个命题Si为假，i=1，2，C)E于是 给出了结论为假的概率的上界。CjjE111,2,11jCjiCiSP为真，所有命题 物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-1219例题：总误差概率的分解当进行一个作C次比较的实验时，首先选择总误差概率，比如说E=0.05或0.11。单个的j可以选为相等()或不相等。由于j的值比较小，则第j个置信区间将比较宽。例如总置信水平要求1-E=95%，当要作11个比较时，那么对所关心的差数(或差别)去构造11个1-J=99.5%的置信区间。当进行大量比较时，Bonferroni法的

15、主要缺点是每一单个区间宽度增加。课本上的例题表明：Bonferroni法适用于少量设计方案进行比较，其上限不要超过11个方案为宜。CEj 物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-1220例题：总误差概率的分解43.143.1807.2025.00025.0zz99.1776.2598.54,025.04,0025.0tt对一组给定的数据和一个大的样本量，设样本量C=11，如果每个样本的误差区间j相同，若1-j=99.5%，j/2=0.0025，这样的样本置信区间宽度将是总置信区间宽度（1-E=95%）的的1.43倍，即：对小

16、样本量来说，比如当样本量C为5时，99.5%的置信区宽度将是95%的置信区间宽度的1.99倍 物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-122111.3 k个系统设计方案之间的比较Bonferroni不等式可达到的三个目标 单个置信区间单个置信区间 Bonferroni法给出的置信区间是最低可能的总置信水平。与现有系统进行比较与现有系统进行比较 Bonferroni法将所有其它系统方案的置信区间与现有系统的置信区间进行比较 所有可能的比较所有可能的比较 对所有的设计方案进行相互比较，即对任意两个系统的设计方案ij，构造ij的置

17、信区间。对于k个设计方案，则要计算的置信区间数是C=k(k1)/2。总置信系数的下限是 jiijE11 物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-122211.3 k个系统设计方案之间的比较试验统计设计的目的 试验的统计设计是设计并评价试验的一组原则。在统计学中，系统的输入变量，如决策变量、结构假设以及随机变量的参数，都称为因子。因子的每一可能的值叫因子的水平。全部因子在给定水平上的一个组合叫一个“处理”。若仿真在相同的处理下运行，但利用的是独立随机数流时，就认为是作了一次独立重复运行试验。统计设计的目的就是确定各种因子对响应

18、变量的影响。物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-122311.3 k个系统设计方案之间的比较因子的分类 定性因子定性因子是用一系列定性的策略描述的规则、逻辑等。例如：排队规则，如采用FIFO还是采用优先权；定量因子定量因子可以用数值来表示，如并行服务台数，到达速率及定货策略等。其它的因子分类：受控因子和非受控因子某些因子在管理控制之下并能随意变化，这些因子统称为决策变量或策略变量，如并行服务台数和定货策略。其它的因子，如随机到达速率或随机需求速率，都不能由管理人员来控制。然而，在仿真模型中，像需求速率这样的非策略变量也可

19、以由分析员来控制。物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-122411.3 k个系统设计方案之间的比较试验统计设计的目的就系统本身来讲，即使某些输入变量没有在策略制定者控制之下，但这些随机变量的特定值却可由分析员在一定程度上由指定所用的随机数种子和随机数流来控制。随机波动源被慎重地引进到模型中是为了准确地表示系统行为。另一方面,仿真实验没有受到具体试验的外部波动的影响，如量测误差的影响。物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-122511.3 k个系统设计

20、方案之间的比较单因子完全随机化试验设计 首先考察只有一个因子可以影响响应变量Y的情况。这是完全随机化试验设计中最为简单的问题。对一个排队系统来讲，单个因子可以是排队规则，它可能有三个水平，如先到先服务，或有优先级的服务，或轮流服务。排队规则是一个定性的策略因子的例子，当仅有一个因子，该因子具有k个水平时，该试验称为单因子试验。物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-122611.3 k个系统设计方案之间的比较单因子完全随机化试验设计的统计模型具有k个处理水平的单因子完全随机试验设计的分析所采用的统计模型是 Yrj=+j+rj

21、 r=1，2，.，Rj；j=1，2，.，k式中，Yrj是因子在第j个水平时，响应变量的第r个观察值，称之为总的平均影响，j是由于因子的第j个水平所引起的影响，rj是在水平j之下第r个观察值的“随机误差”，而Rj是在水平j时的观察次数。物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-122711.3 k个系统设计方案之间的比较单因子完全随机化试验设计统计模型的分析01kjj假设随机偏差项rj具有零均值、协方差为2的正态独立分布。参数及j被假设是固定的且满足 。当因子的水平可由分析人员选定时，这种模型称为固定影响模型。如果因子的水平不能

22、被选定，而是从某一总体中随机选择出来的，j假设是正态分布，那么得到的是随机影响模型。这里仅讨论固定影响模型。物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-1228单因子固定影响完全随机试验的初步分析由统计假设检验组成：H0：j=0，j=1，2，k即因子的水平对响应没有影响。单向方差分析可用于上述统计检验。该检验本身是由计算F统计量并把它的值与一适当的临界值进行比较组成。如果假设H0没有被拒绝，那么分析人员可得出结论：对因子的所有水平的平均响应是对因子的所有水平的平均响应是，即因子对响应变量没有明显即因子对响应变量没有明显的影响。如

23、果假设的影响。如果假设H0被拒绝，那么分析人员有理由相信因子的水被拒绝，那么分析人员有理由相信因子的水平对平均响应有某些影响。平对平均响应有某些影响。物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-122911.3 k个系统设计方案之间的比较F统计量的单向方差分析检验方法 该检验基本上是把观察到的Yrj的变动分成两个成分，其中一个成分是由因子的水平所引起的，而另一成分是由于被仿真的过程所固有的变动所引起的。首先，把观察到的响应Yry作成数据观测表，并计算总数Tj和T，以及第j个水平的样本均值 和整个样本均值或总均值 。jY Y 物流

24、工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-123011.3 k个系统设计方案之间的比较数据观测表第r次重复运行单 因 子 的 水 平 j 12jK1Y11Y12Y1jY1k2Y21Y22Y2 jY2kRj总数T1T2TjTk均值11RY22RYjRjYkRkY1Y2YjYkY第r次重复运行单 因 子 的 水 平 j 12jK1Y11Y12Y1jY1k2Y21Y22Y2 jY2kRj总数T1T2TjTk均值11RY22RYjRjYkRkY1Y2YjYkY第r次重复运行单 因 子 的 水 平 j 12jK1Y11Y12Y1jY1k2Y

25、21Y22Y2 jY2kRj总数T1T2TjTk均值11RY22RYjRjYkRkY1Y2YjYkY第r次重复运行单 因 子 的 水 平 j 12jK1Y11Y12Y1jY1k2Y21Y22Y2 jY2kRj总数T1T2TjTk均值11RY22RYjRjYkRkY1Y2YjYkY仿真实验的观测值记录第j个仿真实验处理时所有响应之和：jRrrjjYT1T 仿真实验所有响应之总和 kjRrrjkjjjYTT111kjjRR1重复运行的总次数Y 整个样本的均值 kjRrrjkjjjjRYRYRRTY111 物流工程与管理系物流工程与管理系Logistics Engineering&Managemen

26、t2023-5-123111.3 k个系统设计方案之间的比较方差估计的计算方法为了估计计算方差，我们先观察下列算式该式反映了仿真响应变量围绕整个样本均值的变化 jrjjrjYYYYYY这种变化由两部分组成 由于某个仿真处理的均值对总体均值之差 由于每个响应对该水平的采样均值响应之差 对上式两端平方，再对所有r及j求和，可以得到：物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-1232方差估计的计算方法 kjRrjrjkjjjkjjjjjjkjjjkjRrjrjkjRrjjrjjkjjjkjRrjrjkjRrjrjjkjjjkjRrj

27、rjkjRrjjrjkjRrjkjRrjrjkjRrjjrjjjrjkjRrjjrjkjRrrjjjjjjjjjjjjjYYYYRYRYRYYYYRYYYRYYYYYRYYYYYYYYRYYYYYYYYYYYYYYYYYYYYYYYY112121121121112112111211211112112112211211222222系统总的平方和可以简单表示为 SSTOTAL =SSTREAT +SSE SSTREAT是由处理引起的平方和 SSE 是误差平方和。注意，这里误差是指在水平j时单个响应Yrj与在水平j时采样平均响应 的偏差。jY 物流工程与管理系物流工程与管理系Logistics En

28、gineering&Management2023-5-123311.3 k个系统设计方案之间的比较方差的无偏估计如果rj是系统仿真的协方差，那么其均方差MSE=SSE/(Rk)是响应变量Y的方差2的无偏估计，即 EMSE=2。如果统计假设检验组成：H0：j=0，j=1，2，k成立，则MSTREAT=SSTREAT/(k1)是响应变量Y的方差2的无偏估计。在任何情况下，MSTREAT与MSE是统计独立的。SSTREAT/2 是自由度为(k1)的2分布。SSE/2是自由度为(Rk)的2分布。用以检验统计量的检验式为：kRMSkMSMSMSFETREATETREAT1这个检验统计量是具有自由度为(k

29、1)与(Rk)的F分布。物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-123411.3 k个系统设计方案之间的比较单向方差分析检验的判别 设检验的显著性指标110(l)%，这里l是具有自由度为k1以及Rk的F分布在临界值为F1-处的概率。注意：k 单因子的水平总数；R各因子水平的运行总次数。如果F F1-，则拒绝零假设H0 物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-123511.3 k个系统设计方案之间的比较两因子的析因设计模型 二因子模型，其统计模型为Yijr=

30、+Qi+Nj+QNij+ijr i=1，2，q；j=1，2，n；r=1，2，R式中Yijr是响应变量Y在第一因子为水平i，第二因子为水平j时的第r个响应的观察值。Yijr是一个全随机设计，即在一个处理里的所有重复运行以及所有处理上的所有重复运行都是统计独立地做出的。随机偏差项ijr，假设是服从零均值、协方差为2的独立、正态分布。是系统总的平均影响 Qi是因子Q在水平i下的影响 Nj是因子N在水平j下的影响 QNij是因子Q在水平i、因子N在水平j的相互起作用的影响 ijr是实际响应的随机波动 物流工程与管理系物流工程与管理系Logistics Engineering&Management202

31、3-5-1236两因子的析因设计模型的方差分析 假设因子Q有q个水平而因子N有n个水平，且每个处理水平有m次重复运行，总共有R=qmn次重复运行。qinjmrijrTOTALYYSS1112 qiiQYYmnSS12njjNYYqmSS12 qinjjiijQNYYYYmSS112由因子Q引起的响应的变化 由因子N引起的响应的变化 由因子Q、N共同引起的响应的变化 总体的响应变化 同样可以证明：SSTOTAL =SSQ +SSN +SSQN +SSE qinjmrijqinjmrijiijqinjmrjjiijqinjmriijijrqinjmrjijijrqinjmrjiijijijrqin

32、jmrqinjmrijqinjmrjiijqinjmrijijrqinjmrijjiijijijrqinjmrijrTOTALYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYSS111111111111111111111111221112111211121112222222 物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-1237两因子的析因设计模型的方差分析SSTOTAL =SSQ +SSN +SSQN +SSE上式右侧的各个平方和项是独立的分布，因此当每个平方和项用2去除，

33、得：SSQ/2符合自由度为(q-1)的2分布；SSN/2符合自由度为(n-1)的2分布；SSQN/2符合自由度为(q-1)(n-1)的2分布；SSE/2符合自由度为qn(m-1)的2分布。物流工程与管理系物流工程与管理系Logistics Engineering&Management2023-5-1238两因子的析因设计模型的方差分析变化的来源自由度SSMS F因子Qq1SSQ因子NN1SSNQN相互作用(q1)(n1)SSQN误差qn(m1)SSE总数qmn1 SSTOTAL1qSSMSQQ1nSSMSNN11nqSSMSQNQN1mqnSSMSEEEQMSMSENMSMSEQNMSMS3个假设可被检验：H01：Qi=0，对所有iH02：Nj=0，对所有jH03：(QN)ij=0，对所有i，j