计算机试验设计均匀布点两算法.pdfVIP

计算机试验设计均匀布点两算法 钱 进 中南财经大学，武汉，430064 摘要 为推广均匀设计并用于计算机试验设计．给均匀设计布点法则一设计了加速算法，绐布 点法则二和三研究了累数宿集．可计算出素数的特征数备快速戢用。 关键字均匀设计，计算机试验设计．布点法则，加速算法，素数宿集 1 引言 计分会第二届学术交流会(青岛)上，张建方以[3]中观点多次在会上质疑而回答者说服力不 够强，致使有与会者认为：E3j的缘起多半源于[5]中计算结果所提出问题之存在，[4]能超脱 [33“限于多因素最优化的范围”．给[23开辟了新的应用领域，但在目前大气候下应用的可能性 不大。实质上[33用的正是[5]用过的计算机验证方法，类似方法被[7]发挥尽致，变成了生产 力，1996年李开坤用类似方法有力地辩驳了[3]。[2]中介绍了国外在计算机试验设计方向做 了许多工作，可惜没有简明有效的采样方法而失于相当复杂。鉴于上述原因，以及[4]改成为 [12]的进展．特把[2]中三个布点法则的计算机自动实现的关键算法归纳于下。 2求均匀点列的加速算法 考虑s个因素，q个水平的试验，用均匀设计法布点。[6]给出了q在3l以内的均匀设计 参量表。为了求得高维的参量，[5]引入了加速算法。 [引理][G]中证得均勺废公式： (1) D(q，i)=二∑Ⅱ(1一．二Ln(2．Sin(Ⅱ．二=二)))一1+o(qt(Lnq)s [引理二]有同余逆关系的P重序元有相同的均匀度。(前人类似结果，[5]简单证明。) 纳。) 定它在不大于N的多重序元集中的序号M。即 、l M= 气。 10 ＼1／ 态对应于第一个P重序元(1，…，P—l，P)。分析历经的状态数，用高阶等整级数可推导出第 1／ 式。 [算法二](组合调度算法D)(2 枚举证明，并在[83中提出“杨辉表”快速计算二项式系数。 用(i。j)表示第i行第j列的表项，“杨辉表”有如下特点： 『l+i一21 ／n＼ 4、(i，j)=l。．I，5、(1．})一(n—i+1．i+1)。 ¨一l／ ＼1／ [算法三]均匀设计参量加速算法(UDV) 从80以内的某数中取若干个组合数的各种情况。 ／n——1＼ u2．(查组合总数)由杨辉表的特点4，查表项(q--s+l，s)得到M—ljl； 、S—l／ 把二级双精度级数展开计算集中在此，以后不再重复计算。 U4．(求缩剩余系)以下针对[6]中布点一。欧拉函数9(0)合于s≤9(0)q； U5．(建J挚表)存放1．2…，M。并选一个足够大的初始fx。 U6．(顺取一序数)由序表中取当前尚未划去的最小数M。， U7．(产生组态i)据[算法D]求得M。对应的S重序元i=(a，，a。．m)； U8．(计算均匀度)据?式计算均匀度f。--D(q，i)， U9．(弃大留小)当Lf。时，记下：；=i，f；一f。； U10．(求同余逆E)由数论知一般有S一1个不同组态6一(b。，b。…，b；)， U11．(求序号)据(2)式用[算法c]求对应序号Mb，将它从序表中划去。 U12．(判断一)若还有同余逆则转U10．，否则继续u13．； U1 3．(判断二)若序表中还有数转U6．。否则继续U14．； U1 了“均匀设计”算法原理、数据结构、机上实现、FORTRANⅣ源程序及框图。课题背景是原七