正交实验教材第三章.DOCVIP

第三章 均匀设计表的构造和运用 本章介绍均匀设计表的构造和使用表的来源，其中均匀性度量──偏差将起关键作用，我们将介绍偏差的定义，并给出正交设计与均匀设计各自偏差的比较，从中可以了解为什么均匀设计可以比正交设计节省试验次数，本章还介绍拟水平在均匀设计中的使用和有关表的构造，熟悉本章内容对于正确理解和使用均匀设计有很大帮助。 3.1 均匀设计表的构造 定义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 试验将是所有最高水平相组合，在有些试验中，例如在化工试验中，所有最高水平组合在一起可能使反应过分剧烈，甚至爆炸。反之，若每个因素的水平都是由高到低排列，则 表中最后一号试验将是所有低水平的组合，有时也会出现反常现象，甚至化学反应不能进行。表则没有类似现象，比较容易安排试验。 iii）若n 为偶数， 表比表有更多的列。如上面讨论过的表只有2列，而表可以有6列。 iv）若n 为奇数，则表列数通常少于表。 v）表比表有更好的均匀性，应优先采用 表，其细节将在下节讨论。 vi）若将或的元素组成一个矩阵的秩最多分别为 及。 本书附录I，列出了2≤s≤7，5≤n≤31，及n 37的表或表，供使用时选择，为了节省篇幅，凡使用表中没有推荐的列我们就没有列出。 3.2 均匀性准则和使用表的产生 在第一章1.6 节我们曾指出均匀设计在