0803-第三节 均匀设计表的构造和运用.doc

第三节 均匀设计表的构造和运用 本节介绍均匀设计表的构造和使用表的来源，其中均匀性度量──偏差将起关键作用，我们将介绍偏差的定义，并给出正交设计与均匀设计各自偏差的比较，从中可以了解为什么均匀设计可以比正交设计节省试验次数，本节还介绍拟水平在均匀设计中的使用和有关表的构造，熟悉本节内容对于正确理解和使用均匀设计有很大帮助。 均匀设计表的构造 定义1 每一个均匀设计表是一个方阵，设方阵有n行m列，每一行是{1，2，...，n}1，2，…，n的重新排列)，{1，2，…，n}的一个子集，但不一定是真子集。 显然，第一章表4-6列举的U(64)，U(74)和U(7)都符合上述定义。 符合定义1的均匀设计表数量太多，本节仅介绍用好格子点法(good lattice point)构造的均匀设计表，其方法如下： 1) 给定试验数n，寻找比n小的整数h，且使n和h的最大公约数为1。符合这些条件的正整数组成一个向量h＝(h，…，h )。 2) 均匀设计表的第j列下法生成 [mod n] (3.1) 这里[mod n] 表示同余运算，若jh超过n，则用它减去n的一个适当倍数，使差落在[1，n] 之中。U可以递推来生成 (3.2) 例如，当n＝9 时，符合条件1)的h有1，2，4，5，7，8；而h=3 或h=6 时不符合条件1)，因为最大公约数(3，9)=3 ，(6，9)=3，1.所以最多只可能有6列，又如当时，用公式 (3.2) 来生成该列时其结果依次如下: 其结果列于表16的第三列。 表16 1 2 3 4 5 6 1 1 2 4 5 7 8 2 2 4 8 1 5 7 3 3 6 3 6 3 6 4 4 8 7 2 1 5 5 5 1 2 7 8 4 6 6 3 6 3 6 3 7 7 5 1 8 4 2 8 8 7 5 4 2 1 9 9 9 9 9 9 9 用上述步骤生成的均匀设计表记作 ，向量h称为该表的生成向量，有时为了强调h 的作用，可将 记成. 给定n ，相应的h 可以象上例那样方便地求得，从而m 也就确定.所以m 是n 的一个函数，这个函数曾由大数学家欧拉研究过，称为欧拉函数，记为E(n) .这个函数告诉我们均匀设计表最多可能有多少列.下面的结果来自数论： i)当n为素数时 ，E(n-1)=n-1所谓素数就是一个正整数，它与其所有比它小的正整数的最大公约数均为1.如2，3，4，5，11，13， ii)当n 为素数幂时，即n 可表成n= ，这里p为素数l，l为正整数，这时 (3.3) 例如n=9 可表为 ，于是 即至多可以有6列。 iii)若n 不属于上述两种情形，这时n一定可以表为不同素数的方幂积，即 (3.4) 这里 为不同的素数， 为正整数，这时 … (3.5) 例如n=12 可表为n= ，于是 即 最多只可能有4列。 上述三种情形中，以素数情形为最好，我们最多可以获得n-1列，而非素数情形，在上述表的结构中永远不可能有n-1 列，例如n=6= ，此时，这说明，当n=6 时，用上述办法生成的均匀设计表只有2列，即最多只能安排两个因素，这是太少了，为此，王元，方开泰(1981)建议，可将 表的最后一行去掉来构造 ，为了区别于由(3.2) 生成的均匀设计表，我们记它为 ，在U 的右上角加一个“*”号，表列于表17，对照表16我们看到U 表和表之间的关系和各自特点： i)所有的表是由 表中划去最后一行而获得； ii)表的最后一行全部由水平n 组成，表的最后一行则不然。若每个因素的水平都是由低到高排列，表中最后一号 表17 No. 1 2 3 4 5 6 1 1 2 3 4 5 6 2 2 4 6 1 3 5 3 3 6 2 5 1 4 4 4 1 5 2 6 3 5 5 3 1 6 4 2 6 6 5 4 3 2 1 试验将是所有最高水平相组合，在有些试验中，例如在化工试验中，所有最高水平组合在一起可能使反应过分剧烈，甚至爆炸。反之，若每个因素的水平都是由高到低排列，则 表中最后一号试验将是所有低水平的组合，有时也会出现反常现象，甚至化学反应不能进行。表则没有类似现象，比较容易安排试验。