方开泰（香港浸会大学香港九龙塘）.PDF

均匀试验设计的理论﹑方法和应用 — 历史回顾 方开泰 （香港浸会大学 香港九龙塘） 摘要 本文回顾计算机仿真试验设计的主要两种方法：拉丁超立体抽样和均匀设计， 在过去二十五年的发展，特别是均匀设计的发展，包括均匀设计的优良性研究、 新的均匀性测度、均匀设计表的构造，以及均匀性在因子设计中的应用。 关键词：均匀设计，拉丁超立体抽样，因子设计，正交性，均匀性。 1 历史回顾 廿世纪七十年代，在系统工程、高科技发展的推动下，计算机仿真 (仿真) 试 验 (computer experiments) 的需求十分强烈，迫切要求高质量的试验设计。于是 计算机仿真试验设计 (Design of computer experiments) 在那时成为一个最有挑战 的课题。在北美洲，三位学者 (McKay, M.D., Beckman, R.J. and Conover, W.J. (1979))在 “Technometrics” 提出了“拉丁超立方体抽样” (Latin Hypercube Sampling) (简称 LHS) 的方法，并立即得到广泛的应用，一批学者对其理论和方法作了系统 地研究和发展，形成了一个独立的分枝。差不多在同一时间，在中国，我和王元 院士提出了“均匀设计” (Uniform Design) (简称UD) 。文章最初在1978 年发表在 中国科学院数学研究所的内部通讯，后来中、英文稿分别发表在 《应用数学学报》 和 《科学通报》。那时，中国正处于文化大革命刚结束，百废待兴的时代，学术上 与世界几乎隔绝。有趣的是，LHS 和UD 有异曲同工之处。表现于: (A) 两种方法均将试验点均匀地散布于输入参数空间，故在文献中广泛使用术语 “充满空间的设计” (space filling design) LHS 给出的试验点带有随机性，故称为抽 样；而UD 是通过均匀设计表来安排试验，不带有随机性。 (B) 两种方法的最初理论均来自“总均值模型” (Overall Mean Model) LHS 希望试 验点对输出变量的总均值提供一个无偏估值，且方差较小，而 UD 是希望试验点 能给出输出变量总均值离实际总均值的偏差最小。 (C) 两种设计均基于U-型设计。 (D) 两种设计能应用于多种多样的模型，且对模型的变化有稳健性。 经过了廿多年的发展，两种不同思路的方法是分道扬镳，还是相互补充，相 互融合呢？本文想作一些回顾和讨论。为此，我们需要介绍两种方法的思路、模 型、方法和应用。 11 2 总均值模型 设输入变量x , L, x 与输出变量有一个确定性的关系 1 s y f (x , L, x ), x (x , L, x )∈Cs . (2.1) 1 s 1 s s s [ ] s 这里假定试验区域为单位立方体C 0,1 ，变量y 在C 上的总均值为 ( ) ( ) E y f x , L, xs dx , L, dxs (2.2) ∫s 1