第五章均匀的设计 《试验的设计与建模》课件.pptVIP

均匀设计 给定试验的诸元素 (n, s,X) 及均匀性测度 D，若一个设计 P? ∈PX 具有最小的偏差值，即 D(P?) = min D(P), 则称P? 为在 (n, s,X,D) 下的均匀设计，或简称为均匀设计。 求解均匀设计的方法分类： (A) 理论求解； (B) 优化数值求解； (C) 求近似解。 均匀设计的理论解 单因素试验 中心化偏差 可卷偏差 不同的偏差，其相应的均匀设计可以不同也可以相同 均匀设计的近似解 U 型设计 若 n× s 矩阵U = (uij) 中第j 列的取值为1, 2, · · · , qj 且这 qj 个数出现的次数相同，该设计称为U 型设计，记为 U(n; q1, · · · , qs)。 缩小设计空间 利用U 型设计，寻求均匀设计的设计空间可大大缩小如下： ?(n; ns) 为设计 U(n; ns) 的导出矩阵，即 U = (uij) = U(n; ns) ?(n; ns) = (xij) 其中 xij = (uij - 0.5) / n 给定n 和s，若设计的第一列选取为自然顺序1, 2, · · · , n，则总共有 (n!)s?1 个U 型设计 U(n; ns)，即使 n 和 s 不算太大，(n!)s?1 也是难以用穷举法来寻找均匀设计。 好格子点法(glpm) 是有效的 quasi-Monte Carlo 方法 (N. M. Korobov (1959), E. Hlawka (1962), H. Niederreiter (1978), L. K. Hua, Y. Wang (1981), J. E. N. Shaw (1988) , K. T. Fang and Y. Wang (1994). 5.5 好格子点法 Glpm 构造 Un(ns) 的步骤 寻找备选的正整数集 Hn={h: hn, gcd(n, h) = 1.}，若 Hn 中元素个数不小于s，转步骤2，否则glpm 无法构造需要的均匀设计，退出程序 在 Hn 中选取 h1, h2,?, hs, 产生一个 n?s 的矩阵 U=(uij)， 其中 uij=ihj (mod n) , 其中模运算要求 1?uij ? n. 记 U 为 U(n,h), 其中 h=(h1,?,hs) 称为 U 的生成向量. 记 Un,s 为所有 U(n,h) 的设计. 在 Un,s 中选择使得偏差值最小的设计即为所求的(近似) 均匀设计Un(ns). 1 2 3 例 5.5 对于 n = 15，s = 2 g.c.d {1,15} = 1 g.c.d {2,15} = 1 g.c.d {3,15} = 3 在中心化偏差意义下，当生成向量 h* = (1,11) 时，U(15, h*) 的偏差最小，故得均匀设计 U15(152). h* = (1,11) CL22 = 0.001600 好格子点法的注意事项 生成向量中元素与 n 互质的的目的是使得向量 hj =(hj , 2hj , · · · , nhj) (mod n) 构成 {1, 2, · · · , n} 的一个置换； 给定 n，其唯一的素数分解为 ， 与其互质的正整数的个数为欧拉函数 ?(n) 欧拉函数 ?(n) 的个数多少对最后选择出来的设计的均匀性有很大影响。 例 当 n = 6, H6 = {1,5}，其元素个数只有 2. 则好格子点法得到的设计 U6(62) 均匀性不好. 因此, 由 H6 更不可能得到 U6(6s), s 2 修正的好格子点法 寻找备选的正整数集 Hn+1 在 Hn+1 中选择 s 个不同元素 h=(h1, h2,?, hs), 产生一个 (n +1)? s 矩阵 U(n+1, h1,?, hs). 去掉第 n+1 行得到设计 U(n,h), 其全体记为 Un,s . 选择 使得 有最小偏差 则 U* 为(近似) 均匀设计. 1 2 3 U* U* ? Un,s 例 (续) 当 n = 6, 得 H7 = {1,2,3,4,5,6} 修正的好格子点法 U6*(62) 好格子点法 U6 (62) 所有 U6(6s), s ? 6 的近似均匀设计可用 H7 得到. 近似均匀设计 U6(62) 方幂好格子点法构造 Un(ns) 的步骤 寻找正整数集 An,s = {a : a n, gcd(a, n) = 1, 且a, a2 · · · , as 在