第4章 正交表的构造.ppt

优化试验设计及数据分析,第4章正交表的构造,4.1概述,1.基本术语不同类型的正交表的构造方法差异很大,甚至有些类型正交表的存在和构造至今还是一个未解决的数学问题。Lu(tq)型表是一类特殊的正交表,其中:t表示水平数。它限定为某数(即除1和它自己以外,不能被任何其他数整除的大于1的正整数。如2、3、5、7、11、13、等)或系数幂(如22、32、23、等)。u表示基本数列,可为任意正整数。q表示总列数。t、u为基本参数,当t、u给定后,则试验次数为tu次,列数为q=(tu-1)/(t-1),一般常用的正交表,如L4(23)、L16(215)、L9(34)L16(45)、L25(56)等属于此类型,其基本参数为(t=2、u=2)、(t=2、u=4)、(t=3、u=2)、(t=2、u=4)、(t=5、u=2)。2.正交表的同等变换正交表的行间置换、列间置换和同一列水平记号的置换,叫做正交表的三种初等变换。经过初等变换所能得到的一切表称为等价的(或同构的)。可以根据不同试验要求,把一个表变成与此等价的其他特殊型式的表。,9.22二水平正交表的构造,1.二水平运算法则构造此类正交表要用到有限域的理论(所谓有限域,大致说就是对有限个元素组成的集合定义了加法、乘法和除法的运算)。我们用0和1表示二水平记号,这个有限域只有两个元素,它们的加法和乘法定义为:a加法法则:0+0=00+1=11+1=0b乘法法则:00=001=011=1,用这种规则定义加法和乘法,是有限域理论所要求的。构造正交表时,将用到上面加法和乘法法则。以后凡是讲到加法和乘法都指这种有限域中的加法和乘法而言。2.正交表与交互作用列表的构造(1)L4(23)的构造L4(23)表是二水平中最小的一个表。它的两个基本参数是t=2、u=2,列数q=(4-1)/(2-1)=3。第一列是将4个试验分成两半,前一半是“0”水平,后一半是“1”水平,称为二分列。第二列是将第一列的两个“0”水平试验和两个“1”水平试验分别再分成一个“0”水平和1个“1”水平,称为四分列。二分列和四分列称为L4(23)的基本列第三列是将第一列与第二列的相应水平,按“加法规则”相加所得如表1所示。,表-1L4(23)表的构造1231000+0=12010+1=13101+0=14111+1=0列号abab,构造交互作用列表时,一般引进“列名”和“列名运算规则”来进行。用a、b分别一记L4(23)的两个基本列,称a为第1列的列号,b为第二列的列名。第3列是由第1列和第2列相加得到的,它的列名可用第一列列名与第二列列名相乘得到的,列名的运算规则是一种指数运算,指数的相加或相乘按加法表或乘法表给出的规则进行。二列交互作用列为其列名相乘。如第1,2两列的交互作用列为ab=ab,即第三列;第1、3列的交互作用列为aab=a1+1b=a0b=b列,即第二列。由此可见,当给出一组完备列的列名后,二列的交互作用列可由列名运算得到。(2)L8(27)的构造L8(27)的参数为t=2、u=3,它有三个基本列,分别置于第1、2、4列,如表2所示。第1列是:二分列,列名为a,这列8个试验被分成两半。第2列的列名为b,是一个四分列。第四列的列名为C,是8分列。其他4列通过列间运算才能得到。,表4L8(27)表的构造12345671000+0=000+0=00+0=00+0=02000+0=010+1=10+1=10+1=13010+1=100+0=01+0=11+0=14010+1=110+1=11+1=01+1=05101+0=101+0=10+0=01+0=16101+0=111+1=00+1=11+1=07111+1=001+0=11+0=10+0=08111+1=011+1=01+1=00+1=1,将列名列成表,它就是L8(27)的一组完备列名表如下:列号1234567列名ababcacbcabc,可以验证,L8(27)中任意二列的交互作用列是七列中的某一列,并可通过列名运算得到。如1、7两列的交互作用aabc=a2bc=bc列,即第6列。因此,可根据列名运算构造交互作用表(如表3所示)供直接查用。表-3L8(27)交互作用表1234567列号(1)3254761(2)167452(3)16543(4)1234(5)325(6)16(7)7,（2）Lu(2q)型正交表与交互作用列表的构造根据上面的方法,可以类似地构造任意基本列数为u的二水平正交表Lu(2q)和交互作用列表。a基本列的构造在Lu(2q)的正交表中,有u个基本列,分别置于第1列,第2列,第4列,第2u-1列上,基本列的列名分别用字母a、b、c来表示。b交互作用列表的构造通过上述a、b两个步骤,就可得到正交表的q=(2u-1)/(2-1)=2u-1个列,这q列即组成Lu(2q)的完备列。任意两列的交互作用列则是q列中的某一列,这一列可用这两列的列名相乘得到据此即可构造出交互作用列表。,下面再以L16(215)加以说明:它的基本列数u=4,即有4个基本列,分别置于1、2、4、8列,其列名分别为a、b、c、d。a为二分列；b为四分列；c为八分列；d为十六分列。第2至第4列之间的列(即第3列)为第2列加第1列,列名ab,而第4列与第8列间的列(即第5至第7列)是第4列依次与第1、2、3列相加而得。第9至15列(共7列)为第8列与前7列依次相加而得。L16(215)中的15列的列名表如下：列号123456789101112131415列名ababcacbcabcdadbdabdcdacdbcdabcd,这是一组完备列名表,任何两列的交互列均在这15列中,可用列名运算找出交互列,如第3列与第13列的交互列为abacd=a1+1bcd=bcd即14列。,二水平正交表的构造原理,可推广到三水平的情况,所不同的只是交互作用列和列名运算有些差别。1.三水平运算规则用0、1、2表示三个水平,其加法和乘法规则规定如下:+0120012112022012.正交表与交互作用列表的构造三水平表的最小一个表是L9(34),它的两个基本参数是t=3,u=2,从而得到其列数为q=(9-1)/(3-1)=4.第1列是三分列，记列名为a。第2列是称为九分列、列名记为b。第3列是由第1列与第2列按加法规则相加而得,其列名称为第1列的列名与第2列的列名相乘,即ab。第四列是将第一列的每个水平按乘法规则乘以“2”,然后与第2列相加得到的,其列名为a2b。,4.3三水平正交表的构造,01200010122021,L9(34)的四列a、b、ab、a2b就是一组标准化完备列名。可以验证,这四列的任两列的交互作用列就是另外两列。例如,第3列ab与第4列a2b的交互作用列为aba2b=a1+2b1+1=b2(b2)2=b22=b(第2列)(ab)2a2b=a2+2b2+1=ab0=a(第1列)反之,若列名最后一个字母的指数不是1,则均称为非标准化列名,如ab2。但是非标准化列名可以通过一定的列名运算就成与它等价的标准化列名。在三水平的情况下,只要将非标准化列名平方一下,即可达到目的。例如:ab2=(ab2)2=a2b22=a2b此外,如果一组列名不仅是一组标准化列名,而且又是一组完备列名,则称它为一组标准经完备列名。,通过对L9(34)表构造的讨论,可将L3(3q)型表和它的交互作用列表的构造法则归纳如下:(1)L3(3q)型正交表的构造给出基本参数u后,L3(3q)型表的总列数,q=(3u-1)/(3-1)=(3u-1)/2其中u列为基本列,分别置第1、2、5,(3u-1)/2+1列(共u列)。基本列的列号分别用字母a、b、c、d命名。第1列为三分列,第2列为九分列,第5列为二十七分列。第(3u-1)/2+1列为3u分列。在每个基本列后(除第一个基本列外)依次安排该基本列与该列前所有的交互作用(共两列),交互作用列的列名用乘法规则得到,这样所得到的q=(3u-1)/2个列名,就是一组标准化完备列名。,(2)交互作用列表的构造L3(3q)正交表的交互作用列表构造,按照列名运算法则可得到。下面看一下L27(313)。它是u=3的表。它的构造法和L9(34)相同,它有三个基本列,置于第1、2、5列,列名分别记为a、b、c。第1列a这三分列;第2列是九分列;第5列c是二十七分列。其余各列可按下列规则得出:前一列加后列以及前一列每水平乘2加后列,列名用乘积表示。L27(313)的列名表如下：列号12345678910111213列名ababa2bcaca2cbcb2cabca2b2ca2bcab2c读者可根据列名具体写出L27(313)