Processing math: 26%

数学物理学报, 2025, 45(2): 584-603

具有重复行的强度 2 的最优正交表的构造

庞善起,*, 路又维,, 王静,

河南师范大学数学与统计学院 河南新乡 453007

The Construction of Optimal Orthogonal Arrays with Repeated Rows and Strength 2

Pang Shanqi,*, Lu Youwei,, Wang Jing,

School of Mathematics and Statistics, Henan Normal University, Henan Xinxiang 453007

通讯作者: * 庞善起,E-mail:pangshanqi@263.net

收稿日期: 2024-06-18   修回日期: 2024-09-24  

基金资助: 国家自然科学基金(11971004)
国家自然科学基金(12201189)
国家自然科学基金(12471245)
中国博士后科学基金(2022M721055)
河南省自然科学基金(242300421164)

Received: 2024-06-18   Revised: 2024-09-24  

Fund supported: NSFC(11971004)
NSFC(12201189)
NSFC(12471245)
CPSF(2022M721055)
Natural Science Foundation of Henan(242300421164)

作者简介 About authors

路又维,E-mail:1277686343@qq.com;

王静,E-mail:wangjing19@126.com

摘要

具有重复行的正交表已经得到广泛应用, 它可以降低试验的复杂性和成本, 同时提高试验结果的可靠性. 具有重复行的最优正交表有更好的统计性质和组合结构, 然而目前对此类的最优正交表知之甚少. 该文主要研究各种水平变换和列变换构造具有重复行的最优正交表的方法. 首先介绍了一种独立列构造饱和的正交表的方法, 然后利用各种水平变换和列变换提出了具有重复行的最优正交表和m-最优正交表的构造方法, 最后作为这些方法的应用为使用人员提供了具有重复行的最优正交表的无穷类.

关键词: 正交表; 具有重复行的正交表; 具有重复行的最优正交表; 水平变换和列变换

Abstract

Orthogonal arrays with repeated rows have been widely used, which can reduce the complexity and cost of experiments and improve the reliability of experimental results. The optimal orthogonal arrays with repeated rows has better statistical properties and combinatorial structure, but little is known about this kind of optimal orthogonal arrays. In this paper, we mainly study the methods of constructing the optimal orthogonal arrays with repeated rows by various permutation of levels and permutation of columns. Firstly, a method of constructing saturated orthogonal arrays by using independent columns is introduced. Then, the construction methods of optimal orthogonal arrays with repeated rows and m-optimal orthogonal arrays are proposed by using various permutation of levels or columns. Finally, several infinite class of optimal orthogonal arrays with repeated rows are obtained and provided for users.

Keywords: orthogonal array; orthogonal array with repeated rows; optimal orthogonal array with repeated rows; permutation of levels and permutation of columns

PDF (636KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

庞善起, 路又维, 王静. 具有重复行的强度 2 的最优正交表的构造[J]. 数学物理学报, 2025, 45(2): 584-603

Pang Shanqi, Lu Youwei, Wang Jing. The Construction of Optimal Orthogonal Arrays with Repeated Rows and Strength 2[J]. Acta Mathematica Scientia, 2025, 45(2): 584-603

1 引言

正交设计是统计学中一个重要的部分, 主要应用于设计试验, 确定影响试验结果的因素和因素水平. 正交设计主要是用正交表来安排试验, 因此正交表是一种统计工具, 可用于优化试验设计, 减少试验次数, 节省时间和资源[1],[2],[3],[4],[5], 在因子分析实验中, 通常采用强度为 2 的正交表作为正交主效应设计, 它们对构建其他有用的正交表具有根本的重要性[6],[7],[8].

目前, 研究人员正致力于开发新的正交表构造方法, 提高正交设计的效率和可靠性, 拓展其在不同领域的应用, 解决实际应用问题, 如密码学[9],[10],[11],[12]、编码理论[13],[14],[15]、计算机试验[16],[17],[18],[19],[20],[21]以及量子信息理论[22],[23],[24],[25],[26],[27]. 此外, 研究人员正在努力开发相关软件和工具, 包括正交设计功能[28],[29], 使研究人员和工程师更容易使用正交表和正交试验设计[30],[31].

正交表存在一行重复出现 m 次, 其中 m2, 则称该正交表为具有重复行的正交表[32],[33]. 具有重复行的正交表是一种特殊类型的试验设计, 目前已经得到广泛研究, 旨在通过重复行来提高试验的稳定性和可靠性[34]. 一般情况下正交表的行是不相同的, 此时每组试验均无重复. 在试验设计中我们更希望使用具有重复行的正交表, 这样我们可以做纯误差估计[35]. 此外, 它也被应用于数学组合设计、物理学量子纠缠态[36]和量子纠错码以及计算机科学理论等许多方面. 例如在组合设计中, 类似于 Fisher 不等式的 Plackett-Burman 界, 给出强度为 2 的重复行正交表的期望不等式和相应证明, 并指出强度 t 的重复行正交表所需不等式的一般结果[37],[38], 而 Ray-Chaudhuri 和 Wilson[39] 的不等式又对 Fisher 不等式进行了推广. 之后, Wilson[40],[41] 又通过正交投影的方法, 证明了具有重复块的 t- 设计的存在性, 这与重复行的正交表是密切相关. 总的来说, 具有重复行的正交表因为其优良的组合性质和统计性质在不断扩大其应用领域.

具有重复行的最优正交表在某些情况下具有更好的性质[32],[33], 在 1980 年 Cheng[42] 开创了这一领域的研究, 证明了强度 2 正交表的普遍最优性. 之后研究人员就开始对其进行各个领域的研究: 在 2007 年 Yu 等[43]提出了基于 COAs 词型模式的最优性准则来选择最优单阵列; 在 2010 年 El-Zanati 等[44]提出了一些构造最优正交阵列的方法; 在 2014 年 Lin[45] 提出了具有重复行的最优正交表还可以应用于区组设计, 有效地减轻干扰变量的影响, 并使分析结果更加可靠; 在 2024 年 Du 等[46]提出了用 2-数对分布矩阵构造最优正交表的方法; 此外平衡不完全区组设计 (BIBD) 的研究与具有重复行的最优正交表有着密切联系[33],[47],[48],[49]. 虽然具有重复行的最优正交表在各个领域都有广泛的应用[50],[51],[52], 但是正交表的构造就可以称为 21 世纪的公开问题, 而具有重复行的正交表的构造更加困难, 具有重复行的最优正交表的构造则是难上加难. 目前我们对此类正交表的构造与研究少之又少.

本文还通过对正交表的汉明距离的研究[53],[54],[55],[56],[57], 将得到的具有重复行的最优正交表进行删除 x 列的操作, 并且通过该操作可以得到对应的 m-最优正交表[32].

本文主要研究通过水平变换和列变换构造具有重复行的最优正交表的方法. 首先介绍了一种由独立列构造饱和正交表的方法, 然后利用通过水平变换和列变换提出了具有重复行的最优正交表和 m-最优正交表的构造方法, 最后作为这些方法的应用为使用者提供了具有重复行的最优正交表的无穷类.

2 预备知识

首先我们介绍一些本文需要用到的符号、定义.

符号说明

1. AT 表示矩阵 A 的转置.

2. (s)=(0,1,2,,s1)T, 1s 表示 s×1 的全 1 向量.

3. 如果 A=(aij)m×n,B=(bij)u×v, 则 Kronecker 积 AB 表示 AB=(aijB)mu×nv, 其中 aijB 表示 u×v 矩阵且元素为 aijbrs(1ru,1sv).

4. 如果二整数 abm 除, 所得的余数相同, 就称 ab 对模 m 同余, 并且记为 ab(mod m)[58].

5. 对于二整数 ab, gcd(a,b)=1 表示 ab 的最大公约数等于 1.

6. F2 是一个二阶有限域, Fn2 是一个 n 维空间.

定义2.1[59] 一个第 j 列的元素为 0,1,2,,qj1n×m 矩阵 A=(a1,a2,,am), 称为一个强度 2 的正交表, 如果满足以下两个条件

1. 每一列中每个元素出现的次数相同.

2. 在任意两列 ai,aj(1i,jm) 中, 每一数对 (0,0),,(0,qj1),(1,0), ,(1,qj1), ,(qj1,qj1) 出现的次数相同.

我们用 Ln(q1qm) 表示一个正交表. 这里 n 称为试验次数, m 称为因素数, qj 称为第 j 个因素的水平数 (j=1,2,,m). 本文只研究所有的列具有相同水平的正交表, 即 q1=q2==qm=q, 记为 Ln(qm), 若 n=λq2λ 为该正交表的指数.

如果 (q1)m=n1, 称 Ln(qm) 是一个饱和正交表.

定义2.2 [2] 表中的列间置换、行间置换、同一列水平记号的置换称为正交表的三种初等变换, 经过初等变换所能得到的一切表称为同构的.

定义2.3 [32]k2,s2,λ1 为整数. 如果存在一个正交表 Lλs2(sk), 其中包含一行重复 m 次, 则

mλs2k(s1)+1.

如果等式成立, 则该正交表 Lλs2(sk) 被称为最优的.

如果

m=λs2k(s1)+1,

则该正交表被称为 m-最优正交表.

引理2.1 [32]k2,s2,λ1 为整数. 如果存在一个最优正交表 Lλs2(sk), 其中包含一个全 0 行重复 m 次, 则该最优正交表的其他的每一行恰好包含 k1s 个 0, 并且 k1(mod s).

定义2.4 [2]B1,,Bmm×k 矩阵 A 的所有行, 且第 i(i=1,,k) 列的元素来自阶为 qi(qi2) 的集合 Gi={0,1,,qi1}, 则 A 中两行 Bu=(aui,,auk)Bv=(avi,,avk) 的汉明距离 d(Bu,Bv) 定义为

d(Bu,Bv)=|{r:1rk,auravr}|.

如果 d 对于任意两行都是固定的, 则称矩阵 A 有固定的汉明距离.

引理2.2 [53] 饱和正交表 Ln(sk) 有固定的汉明距离ns.

下面我们引入一种运用 n 个独立列构造正交表 Lsn(ssn1+sn2++s+1)=Lsn(ssn1s1) 的方法, 其中 s2 为素数或素数幂, n2.

引理2.3n 个独立列为 a1=(s)1sn1,a2=1s(s)1sn2,a3=1s2(s)1sn3,,an=1sn1(s), 则该正交表的全部列均可由这两个独立列构造, 方法如下

1. 单独的独立列 a1 为第 1 列

2. x1a1+a2 的全部组合, 共 s

3. x1a1+x2a2+a3 的全部组合, 共 s2

n. x1a1+x2a2++xn1an1+an, 其中 0x1,x2,,xn1s1, 共 sn1 列.

该方法中的加法按有限域 GF(s) 上的加法进行, 综上, 共 (sn1+sn2++s+1)=sn1s1 列, 最终可得到期望的正交表

Lsn(ssn1s1)= (a1,a2,,an,a1+a2,,(s1)a1+a2,,(s1)a1++(s1)an1+an).

例2.1 s=2 时, 可以根据引理 2.3 构造 L4(23),L8(27),L16(215),,L2n(22n1); 当 s=3 时, 可以根据引理 2.3 构造 L9(34),L27(313),L81(340),,L3n(33n12); 当 s=4 时, 可以根据引理 2.3 构造 L16(45),L64(421),,L4n(44n13).

s=2, n=4 时,

a1=(2)123=[0000000011111111]T,a2=12(2)122=[0000111100001111]T, a3=122(2)12=[0011001100110011]T,a4=123(2)=[0101010101010101]T,

因此

L16(215)=(a1,a2,a1+a2,a3,a1+a3,a2+a3,a1+a2+a3,a4,a1+a4,a2+a4,a1+a2+a4,a3+a4,a1+a3+a4,a2+a3+a4,a1+a2+a3+a4)=[000000001111111100001111000011110000111111110000001100110011001100110011110011000011110000111100001111001100001101010101010101010101010110101010010110100101101001011010101001010110011001100110011001101001100101101001011010010110100110010110]T.

s=4,n=2 时,

a1=(4)14=[0000111122223333]T, a2=14(4)=[0123012301230123]T,

因此

L16(45)=(a1,a2,a1+a2,2a1+a2,3a1+a2)=[00001111222233330123012301230123012310322301321001232301321010320123321010322301]T.

3 具有重复行的 2 水平的最优正交表的构造

本节主要针对具有重复行的 2 水平的最优正交表的构造. 本节直接通过列变换构造期望的正交表, 并且这些列是需要进行选择的. 首先对引理 2.3 构造的 s=2 的正交表 L2n(22n1) 中选中的列进行列变换, 将变换后得到的正交表与原正交表组合在一起, 得到期望的具有重复行的最优正交表. 首先给出下面的定理.

定理3.1 存在只重复 2 个全 0 行的最优正交表 L2n+1(22n1), 其中 n4.

用引理 2.3 构造出正交表

B1=L2n(22n1)=(a1,a2,,an,a1+a2,,a1++an),

其中a1=(2)12n1,a2=12(2)12n2,a3=14(2)12n3,,an=12n1(2). 我们对 B1 进行列变换, 得到一个新的正交表 B2, 下证

L2n+1(22n1)=(B1B2)

为期望的正交表.

因为 B1n 个独立列构造, (a1,a2,,an)=Fn2 为全排列, 所以只有 (a1,a2,,an) 的最后一行为全 1 行, 而 (a1,a2,,an) 并上 (a1+a2), 得到的 (n+1) 列中不包含全 1 行. 为了得到 B2, 对 B1 中这不包含全 1 行的 (n+1) 列进行如下列变换 f=[1234(n1)n(n+1)], 并且 f 表示长度为 (n+1) 的循环变换, 此时除了全 0 和全 1 向量外, 其他任意一个 (n+1) 维向量变换前后必不相同. 所以 f 保证了上述 (n+1) 列中的非零 (n+1) 元组变换后一定不同, 也就是保证了 B2B1 的第 i 行在变换之后不同, 2i2n.

又因为 B1 的汉明距离为 2n1, 且 n4, 所以去掉 (n+1) 列之后的汉明距离 2n1n1 仍大于 1, 所以 B1 的另外 2nn2 列保持不变且每行都不相同, 因此变换之后的 B2 的第 j 行与 B1 的第 i 行不同, ij2i,j2n.

综上, 进行该列变换之后两个正交表 B1B2 除了全 0 行外没有重复行, 由定义 2.3 可知上述 L2n+1(22n1) 满足最优正交表的条件, 且只重复两个全 0 行.

注3.1 B1 中选择不包含全 1 行的 (n+1) 列时, 只需选择 (a1,a2,,an) 并上任意偶数个独立列生成的列即可, 如: (a1,a2,,an) 并上 (a1+a2+a3+a4); (a1,a2,,an) 并上 (a2+a3) 等.

注3.2 上述定理中的循环变换 f 可以有 n!

[1234(n1)n(n+1)], [1324(n1)n(n+1)], , [1(n+1)n(n1)32].

注3.3 n=3 时, 我们通过引理 2.3 构造出的正交表的汉明距离为 4, 不满足定理 3.1 的证明过程, 因此给出推论 3.1 说明 n=3 时的情况.

推论3.1 存在只重复 2 个全 0 行的最优正交表 L16(27).

用引理 2.3 构造出正交表

B1=L8(27)=[00000000001111011001101111001010101101101011001101101001],

我们对 B1 进行列变换, 得到一个新的正交表 B2, 使(B1B2)为期望的正交表. 首先根据定理 3.1 可知, 我们在 B1 中选择不包含全 1 行的 4 列 (a1,a2,a3,a1+a2) 进行如下列变换, 记为 f=[1234]. 并且 f 保证了 B2B1 的第 i 行在变换之后不同, 2i8.

又易证此时剩余三列

(a1+a3,a2+a3,a1+a2+a3)=[010110100110011001101001]T

为全排列, 因此 B2 的第 j 行与 B1 的第 i 行不同, ij,2i,j8.

综上, 进行该列变换之后两个正交表 B1B2 除了全 0 行外没有重复行, 由定义 2.3 可知上述

L16(27)=(B1B2)=[0000111100110011001100110101010101010101001111000011110000001111010110100101101001100110011001100110100101101001]T,

满足最优正交表的条件, 且只重复两个全 0 行.

注3.4 为了得到 B2, 对上述 B1 中不包含全 1 行的 4 列进行的列变换可以有如下 6 种

[1234], [1243], [1324], [1342], [1423], [1432].

f 取其中之一就能保证 B2B1 的第 i 行在变换之后不同.

注3.5 B1 中选择不包含全 1 行的 4 列时, 还可以选择 (a1,a2,a3,a1+a3)(a1,a2,a3,a2+a3), 此时剩余三列分别为

(a1+a2,a2+a3,a1+a2+a3)=[001111000110011001101001]T, (a1+a2,a1+a3,a1+a2+a3)=[001111000101101001101001]T,

并且均为全排列, 能保证 B2 的第 j 行与 B1 的第 i 行在变换之后不同.

定理3.2 存在只重复 m 个全 0 行的最优正交表 L2nm(22n1), 其中 n4, 且 n+1 为奇素数, m=3,4,,n+1.

B1,B2 均仍为定理 3.1 中的由 n 个独立列生成的正交表 L2n(22n1), 且已知(B1B2)为只重复两个全 0 行的 L2n+1(22n1). 下面我们设对 Bn 中不包含全 1 行的这 (n+1) 列进行的列变换为 [1234(n1)n(n+1)].

由定理 3.1 可知, 对 B2 的这 (n+1) 列进行的列变换 [1234(n1)n(n+1)] 可得到一个新的正交表, 设为 B3; 同样的, 对 B3 的这 (n+1) 列进行的列变换 [1234(n1)n(n+1)] 可得到一个新的正交表, 设为 B4; 同样的, 对 B4 的这 (n+1) 列进行的列变换 [1234(n1)n(n+1)] 可得到一个新的正交表, 设为 B5; ; 最后, 对 Bn 的这 (n+1) 列进行的列变换 [1234(n1)n(n+1)] 可得到一个新的正交表, 设为 Bn+1.

取任意两个正交表 Bi,Bj, 1i<jn+1, 通过上述构造过程我们可知 Bj 是由 Bi 的这 (n+1) 列进行 (ji) 次列变换 [1234(n1)n(n+1)] 得到的, 即 BiBj 的列变换可写为 [1234(n1)n(n+1)]ji, 又因为 (n+1) 为奇素数, 所以 [1234(n1)n(n+1)]ji 是一个距离为 (ji), 长度为 (n+1) 的循环变换, 且属于注 3.2 中 n! 种变换之一.

综上, B_1B_{n+1}(n+1) 个正交表 L_{2^{n}}(2^{2^{n}-1}) 两两之间都只重复一个全 0 行.

所以由定义 2.3 可知

L_{2^{n}m}(2^{2^{n}-1})=\left(\begin{array}{c} B_1\\ \vdots\\ B_m\\ \end{array}\right)

为期望的正交表, 其中 m=3,4,\cdots,n+1.

例3.1 最优正交表 L_{16m}(2^{15}) 的构造, 其中m=1,2,3,4,5.

首先, 用引理 2.3 构造正交表 B_1, 选定 B_1 中不包含全 1 行的 5 列为 12348 列, 进行列变化 [12345], 得到新的正交表 B_2

B_1=L_{16}(2^{15})=\left[\begin{array}{cccccccccccccccc} 0000000011111111\\ 0000111100001111\\ 0000111111110000\\ 0011001100110011\\ 0011001111001100\\ 0011110000111100\\ 0011110011000011\\ 0101010101010101\\ 0101010110101010\\ 0101101001011010\\ 0101101010100101\\ 0110011001100110\\ 0110011010011001\\ 0110100101101001\\ 0110100110010110\\ \end{array}\right]^{T}, \ \ B_2=\left[\begin{array}{cccccccccccccccc} 0000111100001111\\ 0000111111110000\\ 0011001100110011\\ 0101010101010101\\ 0011001111001100\\ 0011110000111100\\ 0011110011000011\\ 0000000011111111\\ 0101010110101010\\ 0101101001011010\\ 0101101010100101\\ 0110011001100110\\ 0110011010011001\\ 0110100101101001\\ 0110100110010110\\ \end{array}\right]^{T},

根据定理 3.1 可知 B_1,B_2 之间除了全 0 行外没有重复行, 组合在一起可以得到一个新的最优正交表

L_{32}(2^{15})=\left(\begin{array}{c} B_1\\ B_2\\ \end{array}\right)=\left[\begin{array}{cccccccccccccccc} 0000000011111111&0000111100001111\\ 0000111100001111&0000111111110000\\ 0000111111110000&0011001100110011\\ 0011001100110011&0101010101010101\\ 0011001111001100&0011001111001100\\ 0011110000111100&0011110000111100\\ 0011110011000011&0011110011000011\\ 0101010101010101&0000000011111111\\ 0101010110101010&0101010110101010\\ 0101101001011010&0101101001011010\\ 0101101010100101&0101101010100101\\ 0110011001100110&0110011001100110\\ 0110011010011001&0110011010011001\\ 0110100101101001&0110100101101001\\ 0110100110010110&0110100110010110\\ \end{array}\right]^{T},

其中只重复 2 个全 0 行.

根据定理 3.2 可以得到 B_3,B_4,B_5

B_3=\left[\begin{array}{cccccccccccccccc} 0000111111110000\\ 0011001100110011\\ 0101010101010101\\ 0000000011111111\\ 0011001111001100\\ 0011110000111100\\ 0011110011000011\\ 0000111100001111\\ 0101010110101010\\ 0101101001011010\\ 0101101010100101\\ 0110011001100110\\ 0110011010011001\\ 0110100101101001\\ 0110100110010110\\ \end{array}\right]^{T}\!\!\!, \ \ B_4=\left[\begin{array}{cccccccccccccccc} 0011001100110011\\ 0101010101010101\\ 0000000011111111\\ 0000111100001111\\ 0011001111001100\\ 0011110000111100\\ 0011110011000011\\ 0000111111110000\\ 0101010110101010\\ 0101101001011010\\ 0101101010100101\\ 0110011001100110\\ 0110011010011001\\ 0110100101101001\\ 0110100110010110\\ \end{array}\right]^{T}\!\!\!, \ \ B_5=\left[\begin{array}{cccccccccccccccc} 0101010101010101\\ 0000000011111111\\ 0000111100001111\\ 0000111111110000\\ 0011001111001100\\ 0011110000111100\\ 0011110011000011\\ 0011001100110011\\ 0101010110101010\\ 0101101001011010\\ 0101101010100101\\ 0110011001100110\\ 0110011010011001\\ 0110100101101001\\ 0110100110010110\\ \end{array}\right]^{T}\!\!\!,

因此 B_1,B_2,B_3,B_4,B_5 任意 m 个组合在一起可以得到一个新的最优正交表, 其中只重复 m 个全 0 行, 如

{\scriptsize L_{80}(2^{15})\!=\! \left[\begin{array}{ccccc} 0000000011111111&0000111100001111&0000111111110000&0011001100110011&0101010101010101\\ 0000111100001111&0000111111110000&0011001100110011&0101010101010101&0000000011111111\\ 0000111111110000&0011001100110011&0101010101010101&0000000011111111&0000111100001111\\ 0011001100110011&0101010101010101&0000000011111111&0000111100001111&0000111111110000\\ 0011001111001100&0011001111001100&0011001111001100&0011001111001100&0011001111001100\\ 0011110000111100&0011110000111100&0011110000111100&0011110000111100&0011110000111100\\ 0011110011000011&0011110011000011&0011110011000011&0011110011000011&0011110011000011\\ 0101010101010101&0000000011111111&0000111100001111&0000111111110000&0011001100110011\\ 0101010110101010&0101010110101010&0101010110101010&0101010110101010&0101010110101010\\ 0101101001011010&0101101001011010&0101101001011010&0101101001011010&0101101001011010\\ 0101101010100101&0101101010100101&0101101010100101&0101101010100101&0101101010100101\\ 0110011001100110&0110011001100110&0110011001100110&0110011001100110&0110011001100110\\ 0110011010011001&0110011010011001&0110011010011001&0110011010011001&0110011010011001\\ 0110100101101001&0110100101101001&0110100101101001&0110100101101001&0110100101101001\\ 0110100110010110&0110100110010110&0110100110010110&0110100110010110&0110100110010110\\ \end{array}\right]^{T}\!\!,}

其中只重复 5 个全 0 行.

4 具有重复行的 s 水平最优正交表的构造, 其中 {s\geq3}

上一节我们构造了具有重复行的 2 水平最优正交表, 本节我们将构造具有重复行的 s 水平最优正交表. 首先构造最优正交表 L_{s^{2}m}(s^{s+1}), 先对引理 2.3 中 n=2,s\geq3 时构造的正交表 L_{s^{2}}(s^{s+1}) 里面的独立列 a_{1} 进行水平变换, 再对整个正交表进行列变换, 然后将变换前后的正交表组合在一起, 构成一个只重复 m 个全 0 行的最优正交表.

其次构造 L_{2s^{n}}(s^{\frac{s^{n}-1}{s-1}}), 用引理 2.3 中 n\geq3,s\geq3 时构造的正交表 L_{s^{n}}(s^{\frac{s^{n}-1}{s-1}}) 通过列变换构造期望的只重复 2 个全 0 行的最优正交表, 与上述构造方法不同的是, 该方法不用进行水平变换, 而是只需要进行列变换, 当然这些列也是需要进行选择的. 并且本节所有的定理要求 s 为素数或素数幂, 首先给出下面的定理.

定理4.1 存在只重复 m 个全 0 行的最优正交表 L_{s^{2}m}(s^{s+1}), 其中 m=1,2,\cdots,s+1s\geq4, s 为素数或素数幂.

用引理 2.3 构造出正交表

B_1=L_{s^{2}}(s^{s+1})=(a_{1},a_{2},a_{1}+a_{2},\cdots,(s-1)a_{1}+a_{2}),

其中 a_1=(s)\otimes1_{s},a_2=1_{s}\otimes(s).

B_1 的第一列作水平变换: (1,2,\cdots,s-1) 变为 (2,3,\cdots,s-1,1) 得到一个新的正交表, 设为 B=(a,a_{2},a_{1}+a_{2},\cdots,(s-1)a_{1}+a_{2})=(a,b_{1},b_{2},\cdots,b_{s}). 再将 B 的第一列分别与后面 s 列位置互换得到 s 个新的表, 设为 B_2,B_3,\cdots,B_{s+1}. 因为符号变换与列变换都是正交表的同构变换, 不改变其正交性, 所以 B_2,B_3,\cdots,B_{s+1} 仍是正交表. 下证 L_{s^{2}m}(s^{s+1})=\left(\begin{array}{c} B_1\\ \vdots\\ B_m\\ \end{array}\right) 为期望的正交表.

首先证明 B_1B_2 的交集只有全 0 行.

B_1 分为两块, 第一块为前 s 行, 第二块为后 (s^{2}-s) 行, 因为 B_2B_1 交换前两列得到的, 所以 B_2 的前 s 行也是 B_1 交换前两列得到的. 明显 B_1 的前 s 行与 B_2 的前 s 行除全 0 行外都不相同, 因为 B_1s 行的 0 列与 B_2s 行的 0 列所在的位置不同.

因为 B_1 的指数为 1, 所以后 (s^{2}-s) 行的第三、四列的数对都不相同, 且在两种变换之后保持不变, 因此 B_1 的第 i 行与 B_2 的第 j 行不相同, 其中 i\neq js+1\leq i,j\leq s^{2}. 因此只要证明 B_1B_2(s^{2}-s) 行第一、二列的第 i 行的数对在两种变换之后一定不同, 则后 (s^{2}-s) 行的 (s+1) 元组都不相同. 设 B_1 中后 (s^{2}-s) 行第一、二列第 i 行数对为 (a_i',b_{1i}), 在 B 中变为 (a_i,b_{1i}), 且 a_i\neq a_i', 在 B_2 中变为 (b_{1i},a_i), 若 (a_i',b_{1i})=(b_{1i},a_i), 需要 a_i'=b_{1i}=a_i, 矛盾.

B_1B_i 的证明类似 B_1B_2, 其中 3\leq i\leq m.

再证 B_2B_3 的交集只有全 0 行.

B_2,B_3 分为两块, 第一块为前 s 行, 第二块为后 (s^{2}-s) 行, 明显 B_2B_3 的第一块除全 0 行外都不相同.

i 行取自第二块时, 设 B_2 的第 i 行前三列为 (b_i,a_i,c_i), B_3 的第 i 行前三列为 (c_i,b_i,a_i), 若 (b_i,a_i,c_i)=(c_i,b_i,a_i), 需要 a_i=b_i=c_i, 这种情况只会在第一块中出现, 矛盾, 因此 B_2 的第 i 行与 B_3 的第 i 行不同.

ij 行都取自第二块, 且 i\neq j 时, 设 B_2 的第 i 行前五列为 (b_i,a_i,c_i,d_i,e_i), B_3 的第 j 行前五列为 (c_j,b_j,a_j,d_j,e_j), 若 (b_i,a_i,c_i,d_i,e_i)=(c_j,b_j,a_j,d_j,e_j), 需要 (d_i,e_i)=(d_j,e_j), 而 \lambda=1, 矛盾.

B_iB_j 的证明类似 B_2B_3, 其中 i\neq j, 且 2\leq i,j\leq m.

综上, B_1,B_2,\cdots,B_{s+1} 除了全 0 行外没有别的重复行, 由定义 2.3 可知这 (s+1) 个正交表任意 m 个可以组合成期望的最优正交表

L_{s^{2}m}(s^{s+1})=\left(\begin{array}{c} B_1\\ \vdots\\ B_m\\ \end{array}\right),

其中 m=1,2,\cdots,s+1, 得证.

注4.1 具有重复行的最优正交表 "只重复 m 个全 0 行" 是指该正交表除了全 0 行外没有别的重复行. 当 s=3 时, 我们通过引理 2.3 构造出的正交表只有 4 列, 不满足定理 4.1 的证明过程, 因此给出推论 4.1 说明 s=3 时的情况.

推论4.1 存在只重复 m 个全 0 行的最优正交表 L_{9m}(3^{4}), 其中 m=1,2,3,4.

用引理 2.3 构造出正交表

B_1=L_9(3^4)=(a_{1},a_{2},a_{1}+a_{2},2a_{1}+a_{2})=(a',b,c,d)=\left[\begin{array}{ccccccccc} 0&0&0&1&1&1&2&2&2\\ 0&1&2&0&1&2&0&1&2\\ 0&1&2&1&2&0&2&0&1\\ 0&1&2&2&0&1&1&2&0\\ \end{array}\right]^{T},

B_1 的第一列作水平变换: (1,2) 变为 (2,1) 得到一个新的正交表, 设为

B=(a,b,c,d)=\left[\begin{array}{ccccccccc} 0&0&0&2&2&2&1&1&1\\ 0&1&2&0&1&2&0&1&2\\ 0&1&2&1&2&0&2&0&1\\ 0&1&2&2&0&1&1&2&0\\ \end{array}\right]^{T},

再将 B 的第一列分别与后面 3 列位置互换得到 3 个新的正交表, 设为 B_2,B_3,B_4

B_2=\left[\begin{array}{ccccccccc} 0&1&2&0&1&2&0&1&2\\ 0&0&0&2&2&2&1&1&1\\ 0&1&2&1&2&0&2&0&1\\ 0&1&2&2&0&1&1&2&0\\ \end{array}\right]^{T}, \ \ B_3=\left[\begin{array}{ccccccccc} 0&1&2&1&2&0&2&0&1\\ 0&1&2&0&1&2&0&1&2\\ 0&0&0&2&2&2&1&1&1\\ 0&1&2&2&0&1&1&2&0\\ \end{array}\right]^{T}, \ \ B_4=\left[\begin{array}{ccccccccc} 0&1&2&2&0&1&1&2&0\\ 0&1&2&0&1&2&0&1&2\\ 0&1&2&1&2&0&2&0&1\\ 0&0&0&2&2&2&1&1&1\\ \end{array}\right]^{T}.

下证 \left(\begin{array}{c} B_1\\ \vdots\\ B_m\\ \end{array}\right) 为期望的正交表, 其中 m=1,2,3,4.

B_1B_2, B_1B_3, B_1B_4 之间只重复全 0 行, 证明过程与定理 4.1 类似.

下证 B_2B_3 的交集只有全 0 行.

B_2,B_3 分为三块, 第一块为前三行, 第二块为中间三行, 第三块为后三行, 明显 B_2B_3 的前三行除 (0000) 外都不相同. 当 i 行取自第二块或者第三块时, 证明过程与定理 4.1 类似.

ij 行取自第二块或者第三块中的同一块, 且 i\neq j 时, 设 B_2 的第 i 行为 (b_i,a_i,c_i,d_i), B_3 的第 j 行为 (c_j,b_j,a_j,d_j), 若 (b_i,a_i,c_i,d_i)=(c_j,b_j,a_j,d_j), 需要 d_i=d_j, 而 ij 行取自同一块, 矛盾. 而当 B_2 的第 i 行与 B_3 的第 j 行取自第二块或者第三块中的不同块, 且 i\neq j 时, 设 B_2 的第 i 行为 (b_i,a_i,c_i,d_i), B_3 的第 j 行为 (c_j,b_j,a_j,d_j), 若 (b_i,a_i,c_i,d_i)=(c_j,b_j,a_j,d_j), 有 b_i=c_j,a_i=b_j,c_i=a_j,d_i=d_j, 因为 a_i\neq a_j, 设 b_i=c_j=\alpha,d_i=d_j=\beta,a_i=1,a_j=2, 所以 B_2i 行对应 B_1 中的 (a_j,b_i,c_i,d_i)=(2,\alpha,2,\beta), B_3j 行对应 B_1 中的 (a_i,b_j,c_j,d_j)=(1,1,\alpha,\beta), 而第二、三块中一行最多两个 1 或 2, 且最多只有一个 0, 矛盾, 因此 B_2 的第 i 行与 B_3 的第 j 行不同.

B_2B_4, B_3B_4 的证明类似 B_2B_3.

综上, B_1,B_2,B_3,B_4 除了全 0 行外没有别的重复行, 由定义 2.3 可知这四个正交表任意 m 个可以组合成期望的最优正交表 L_{9m}(3^{4}), 其中 m=1,2,3,4, 得证.

推论4.2 s\geq4 为素数或素数幂, 对每一个 m=1,2,\cdots,s+1 都存在 s-2 个最优正交表 L_{s^{2}m}(s^{s+1}), 其中每个都只重复 m 个全 0 行, 任意两个有 s^{2}+s(m-1) 个相同行.

用引理 2.3 构造出正交表

L_{s^{2}}(s^{s+1})=(a_{1},a_{2},a_{1}+a_{2},\cdots,(s-1)a_{1}+a_{2})=B_1, 其中 a_1=(s)\otimes1_{s},a_2=1_{s}\otimes(s).

B_1 的第一列中水平 (1,2,\cdots,s-1) 分别作水平变换为 (2,3,\cdots,s-1,1),(3,4,\cdots,s-1,1,2),\cdots,(s-1,1,2,\cdots,s-3,s-2) 得到 s-2 个新的正交表, 记为 C_1,C_2,\cdots,C_{s-2}. 将 C_1 的第一列分别与后面 s 列位置互换得到 s 个新的正交表, 设为 B_2,B_3,\cdots,B_{s+1}, 根据定理 4.1 可得

A_1=L_{s^{2}m}(s^{s+1})=\left(\begin{array}{c} B_1\\ B_2\\ \vdots\\ B_m\\ \end{array}\right)

为一个最优正交表.

同理, 由 C_2 可得第二个最优正交表

A_2=L_{s^{2}m}(s^{s+1})=\left(\begin{array}{c} B_1\\ B_2'\\ \vdots\\ B_m'\\ \end{array}\right).

C_3,C_4,\cdots,C_{s-2} 做同样的操作, 可以得到最优正交表 A_3,A_4,\cdots,A_{s-2}.

下证这 s-2 个最优正交表 A_1,A_2,\cdots,A_{s-2} 中任意两个有 s^{2}+s(m-1) 个相同行.

我们不妨只证 A_{1},A_{2} 满足要求, 明显 A_{1},A_{2} 的前 s^{2} 行是重复的.

再证 B_{2}B_{2}' 除前 s 行没有相同行. 将 B_2,B_{2}' 分为两块, 第一块为前 s 行, 第二块为后 (s^{2}-s) 行, 根据定理 4.1 可知 B_{2}B_{2}'s 行一定是相同的. 当 ij 行取自第二块时, 因为 B_{2}B_{2}' 的第二列相同位置的水平数一定不同, 所以 B_{2} 的第 i 行和 B_{2}' 的第 i 行一定不同. 又因为 B_{2}B_{2}' 都是饱和正交表, 根据引理 2.2 可知他们的汉明距离都为 s, 因此除第二列外 B_{2} 的第 i 行和 B_{2}' 的第 j 行一定不同. 因此, B_{2}B_{2}' 除前 s 行外没有相同行, 同理 B_{i}B_{i}' 的证明类似 B_{2}B_{2}', 其中 3\leq i\leq m.

最后证 B_{2}B_{3}' 的交集只有全 0 行. 将 B_2,B_3' 分为两块, 第一块为前 s 行, 第二块为后 (s^{2}-s) 行, 明显 B_2B_3' 的第一块除全 0 行外都不相同. 当 i 行取自第二块时, 设 B_2 的第 i 行前三列为 (b_i,a_{i1},c_i), B_3' 的第 i 行前三列为 (c_i,b_i,a_{i2}), 若 (b_i,a_{i1},c_i)=(c_i,b_i,a_{i2}), 需要 a_{i1}=a_{i2}, 而 A_{1},A_{2} 是对 B_1 的第一列中元素作了不同的水平变换, 所以矛盾, 因此 B_2 的第 i 行与 B_3' 的第 i 行不同. 当 ij 行都取自第二块, 且 i\neq j 时, 设 B_2 的第 i 行前五列为 (b_i,a_{i1},c_i,d_i,e_i), B_3' 的第 j 行前五列为 (c_j,b_j,a_{j1},d_j,e_j), 若 (b_i,a_{i1},c_i,d_i,e_i)=(c_j,b_j,a_{j1},d_j,e_j), 需要 (d_i,e_i)=(d_j,e_j), 而 B_{2}B_{3}' 的指数 \lambda=1, 矛盾. 因此, B_{2}B_{3}' 的交集只有全 0 行, 同理 B_iB_j' 的证明类似 B_2B_3', 其中 i\neq j, 且 2\leq i,j\leq m.

综上, A_{1},A_{2} 中有 s^{2}+s(m-1) 个相同行.

例4.1 最优正交表 L_{16m}(4^{5}) 的构造, 其中 m=1,2,3,4,5.

首先, 用引理 2.3 构造正交表 B_1, 将 B_1 的第一列作水平变换, (1,2,3) 变为 (2,3,1) 得到一个新的正交表, 设为 B

B_1=L_{16}(4^{5})=\left[\begin{array}{cccccccccccccccc} 0000111122223333\\ 0123012301230123\\ 0123103223013210\\ 0123230132101032\\ 0123321010322301\\ \end{array}\right]^{T}, \ \ B=\left[\begin{array}{cccccccccccccccc} 0000222233331111\\ 0123012301230123\\ 0123103223013210\\ 0123230132101032\\ 0123321010322301\\ \end{array}\right]^{T},

再将 B 的第一列分别与后面 4 列位置互换得到 4 个新的正交表, 分别为

B_2=\left[\begin{array}{cccccccccccccccc} 0123012301230123\\ 0000222233331111\\ 0123103223013210\\ 0123230132101032\\ 0123321010322301\\ \end{array}\right]^{T}, \ \ B_3=\left[\begin{array}{cccccccccccccccc} 0123103223013210\\ 0123012301230123\\ 0000222233331111\\ 0123230132101032\\ 0123321010322301\\ \end{array}\right]^{T}, B_4=\left[\begin{array}{cccccccccccccccc} 0123230132101032\\ 0123012301230123\\ 0123103223013210\\ 0000222233331111\\ 0123321010322301\\ \end{array}\right]^{T}, \ \ B_5=\left[\begin{array}{cccccccccccccccc} 0123321010322301\\ 0123012301230123\\ 0123103223013210\\ 0123230132101032\\ 0000222233331111\\ \end{array}\right]^{T}.

上述 B_1,B_2,B_3,B_4,B_5 任意 m 个组合在一起可以得到一个最优正交表 L_{16m}(4^{5}), 其中只重复 m 个全 0 行, 如

{\begin{align*} A_{1}&=L_{80}(4^{5})\\ &=\left[\begin{array}{cccccccccccccccc} 0000111122223333&0123012301230123&0123103223013210&0123230132101032&0123321010322301\\ 0123012301230123&0000222233331111&0123012301230123&0123012301230123&0123012301230123\\ 0123103223013210&0123103223013210&0000222233331111&0123103223013210&0123103223013210\\ 0123230132101032&0123230132101032&0123230132101032&0000222233331111&0123230132101032\\ 0123321010322301&0123321010322301&0123321010322301&0123321010322301&0000222233331111\\ \end{array}\right]^{T}, \end{align*}}

其中只重复 5 个全 0 行.

同理, 根据推论 4.2, 若将 B_1 的第一列 (1,2,3) 作水平变换变为 (3,1,2), 则最终可以得到一个新的最优正交表

{\begin{align*} A_{2}&=L_{80}(4^{5})\\ &=\left[\begin{array}{cccccccccccccccc} 0000111122223333&0123012301230123&0123103223013210&0123230132101032&0123321010322301\\ 0123012301230123&0000333311112222&0123012301230123&0123012301230123&0123012301230123\\ 0123103223013210&0123103223013210&0000333311112222&0123103223013210&0123103223013210\\ 0123230132101032&0123230132101032&0123230132101032&0000333311112222&0123230132101032\\ 0123321010322301&0123321010322301&0123321010322301&0123321010322301&0000333311112222\\ \end{array}\right]^{T}, \end{align*}}

其中只重复 5 个全 0 行, 根据推论 4.2 可知 A_{1},A_{2} 有 32 个重复行, 分别为第 1-20 行, 第 33-36 行, 第 49-52 行, 第 65-68 行.

定理4.2 存在只重复 2 个全 0 行的最优正交表 L_{2s^{n}}(s^{\frac{s^{n}-1}{s-1}}), 其中 n\geq3,s\geq3, s 为素数或素数幂.

B_1 是引理 2.3 中 n 个独立列构造的正交表. 我们将对 B_1 进行列变换, 得到一个新的正交表 B_2, 使

L_{2s^{n}}(s^{\frac{s^{n}-1}{s-1}})=\left(\begin{array}{c} B_1\\ B_2\\ \end{array}\right)

为期望的正交表.

因为 B_1n 个独立列构造, (a_1,a_2,\cdots,a_n)=F_{s}^{n} 为全排列, 所以 (a_1,a_2,\cdots,a_n) 中存在 (1,1,\cdots,1),(2,2,\cdots,2),\cdots,(s-1,s-1,\cdots,s-1), 而 (a_1,a_2,\cdots,a_n) 并上 (a _1+a_2), 得到的 (n+1) 列中一定不包含上述水平全部相同的行. 类似于定理 3.1 的证明, 对 B_1 中这 (n+1) 列进行列变换 f=[1234\cdots(n-1)n(n+1)] 得到 B_2, 并且 B_2B_1 的第 i 行在变换之后不同, 2\leq i\leq s^{n}.

又因为 B_1 的汉明距离为 s^{n-1}, 且 s\geq3,n\geq3, 所以去掉 (n+1) 列之后的汉明距离 s^{n-1}-n-1 仍大于 1, 所以 B_1 的另外 s^{\frac{s^{n}-1}{s-1}}-n-1 列保持不变且每行都不相同, 因此变换之后的 B_2 的第 j 行与 B_1 的第 i 行不同, i\neq j,2\leq i,j\leq s^{n}.

综上, 进行该列变换之后两个正交表 B_1B_2 除了全 0 行外没有重复行, 由定义 2.3 可知上述 L_{2s^{n}}(s^{\frac{s^{n}-1}{s-1}}) 满足最优正交表的条件, 且只重复两个全 0 行.

注4.2 B_1 中选择不包含 (1,1,\cdots,1),(2,2,\cdots,2),\cdots,(s-1,s-1,\cdots,s-1)(n+1) 列时, 只需选择 (a_1,a_2,\cdots,a_n) 并上任意两个独立列生成的列即可, 如: (a_1,a_2,\cdots,a_n) 并上 (a_2+a_3) 等.

注4.3 根据定理 3.1 可知, 为了得到 B_2, 对上述 B_1 中选定的列进行的列变换也可以有同注 3.2 的 n! 种变换, f 取其中之一就能保证 B_2B_1 的的第 i 行在变换之后不同.

例4.2 最优正交表 L_{54}(3^{13}) 的构造, 其中只重复两个全 0 行.

根据定理 4.2 可以得到正交表 B_1B_2

B_1=L_{27}(3^{13})=\left[\begin{array}{ccccccccccccccccccccccccccc} 000000000111111111222222222\\ 000111222000111222000111222\\ 012012012012012012012012012\\ 000111222111222000222000111\\ 000111222222000111111222000\\ 012012012120120120201201201\\ 012012012201201201120120120\\ 012120201012120201012120201\\ 012201120012201120012201120\\ 012120201120201012201012120\\ 012201120201120012120012201\\ 012120201201012120120201012\\ 012201120120012201201120012\\ \end{array}\right]^{T},
B_2=\left[\begin{array}{cccccccccccccccc} 000111222000111222000111222\\ 012012012012012012012012012\\ 000111222111222000222000111\\ 000000000111111111222222222\\ 000111222222000111111222000\\ 012012012120120120201201201\\ 012012012201201201120120120\\ 012120201012120201012120201\\ 012201120012201120012201120\\ 012120201120201012201012120\\ 012201120201120012120012201\\ 012120201201012120120201012\\ 012201120120012201201120012\\ \end{array}\right]^{T}\!\!\!,
L_{54}(3^{13})\!=\!\left(\begin{array}{c} B_1\\ B_2\\ \end{array}\right)\!=\!\left[\begin{array}{ccccccccccccccccccccccccccc} 000000000111111111222222222&000111222000111222000111222\\ 000111222000111222000111222&012012012012012012012012012\\ 012012012012012012012012012&000111222111222000222000111\\ 000111222111222000222000111&000000000111111111222222222\\ 000111222222000111111222000&000111222222000111111222000\\ 012012012120120120201201201&012012012120120120201201201\\ 012012012201201201120120120&012012012201201201120120120\\ 012120201012120201012120201&012120201012120201012120201\\ 012201120012201120012201120&012201120012201120012201120\\ 012120201120201012201012120&012120201120201012201012120\\ 012201120201120012120012201&012201120201120012120012201\\ 012120201201012120120201012&012120201201012120120201012\\ 012201120120012201201120012&012201120120012201201120012\\ \end{array}\right]^{T}\!\!\!,

就是只重复 2 个全 0 行的最优正交表.

5 具有重复行的 s 水平 m-最优正交表的构造,其中 {s\geq2}

本节研究 m-最优正交表的构造, 在前两节构造的最优正交表的基础之上通过删除一些列, 获得 m-最优正交表. 虽然删列这个操作很简单, 但是具体删除多少列, 删除哪些列, 才能得到 m-最优正交表, 是我们本节重点讨论的问题. 首先给出下面的定理.

定理5.1 对于 n\geq4, 且 (n+1) 为奇素数; m=1,2,\cdots,n+1, 存在只重复 m 个全 0 行的 m-最优正交表 L_{2^{n}m}(2^{2^{n}-1-x}), 其中 0\leq x<\frac{2^{n}}{m+1}, x 为整数.

根据定理 3.2 我们可以得到最优正交表

L_{2^{n}m}(2^{2^{n}-1})=\left(\begin{array}{c} B_1\\ \vdots\\ B_m\\ \end{array}\right).

其中 B_2,\cdots,B_m 是由 B_1 中不包含全 1 行的 (n+1) 列进行列变换得到的, 而剩余 2^{n}-n-2 列保持不变, 因此如果删除的 x 列不包含进行列变换的 (n+1) 列, 那么 B_iB_j 删除这 x 列后的第 a 行也不同, i\neq j1\leq i,j\leq m; 2\leq a\leq2^{n}.

由定义 2.3 可知, 一个正交表 L_{\lambda s^{2}}(s^{k}) 满足 m-最优正交表的条件为 m=\left\lfloor\frac{\lambda s^{2}}{k(s-1)+1}\right\rfloor, 因此从最优正交表 L_{2^{n}m}(2^{2^{n}-1}) 中删除 x 列得到的正交表 L_{2^{n}m}(2^{2^{n}-1-x}) 是一个 m-最优正交表需要满足 m=\left\lfloor\frac{2^{n}m}{2^{n}-x}\right\rfloor, 即

m\leq\frac{2^{n}m}{2^{n}-x}<m+1,\quad 0\leq x< \frac{2^{n}(m+1)-2^{n}m}{m+1},\quad 0\leq x< \frac{2^{n}}{m+1}.

由引理 2.2 可知 B_1,\cdots,B_m 的汉明距离为 2^{n-1}. 当 n\geq4(n+1) 为奇素数时, 在删除 x 列后的汉明距离总是大于等于 1, 即 2^{n-1}-x\geq1, 因此删除 x 列后的 B_i 的第 a 行与 B_j 的第 b 行不同, i\neq j1\leq i,j\leq m; a\neq b2\leq a,b\leq2^{n}.

综上, 最优正交表 L_{2^{n}m}(2^{2^{n}-1}) 删除 x 列(不包含进行列变换的 (n+1) 列)之后得到的 L_{2^{n}m}(2^{2^{n}-1-x}) 是期望的正交表.

注5.1 原正交表有 2^{n}-1 列, 可以删除的列总共有 2^{n}-n-2 列, 而当 m=1 时, \frac{2^{n}}{m+1}=2^{n-1} 最大, 但是 n\geq4, 2^{n}-n-2 总是大于 2^{n-1}, 因此不存在列不够删除的情况.

例5.1 在定理 5.1 中让 n=4 我们可以获得 m-最优正交表 L_{16m}(2^{15-x}), 其中 m=1,2,3,4,5; 0\leq x<\frac{16}{m+1}x 为整数, 如 m=1,x=7m=5,x=2 时可以得到如下 m-最优正交表 L_{16}(2^{8}),L_{80}(2^{13})

L_{16}(2^{8})=\left[\begin{array}{cccccccccccccccc} 0000000011111111\\ 0000111100001111\\ 0000111111110000\\ 0011001100110011\\ 0011001111001100\\ 0011110000111100\\ 0011110011000011\\ 0101010101010101\\ \end{array}\right]^{T},

L_{80}(2^{13})= { \left[\begin{array}{cccccccccccccccc} 0000000011111111&0000111100001111&0000111111110000&0011001100110011&0101010101010101\\ 0000111100001111&0000111111110000&0011001100110011&0101010101010101&0000000011111111\\ 0000111111110000&0011001100110011&0101010101010101&0000000011111111&0000111100001111\\ 0011001100110011&0101010101010101&0000000011111111&0000111100001111&0000111111110000\\ 0011001111001100&0011001111001100&0011001111001100&0011001111001100&0011001111001100\\ 0011110000111100&0011110000111100&0011110000111100&0011110000111100&0011110000111100\\ 0011110011000011&0011110011000011&0011110011000011&0011110011000011&0011110011000011\\ 0101010101010101&0000000011111111&0000111100001111&0000111111110000&0011001100110011\\ 0101010110101010&0101010110101010&0101010110101010&0101010110101010&0101010110101010\\ 0101101001011010&0101101001011010&0101101001011010&0101101001011010&0101101001011010\\ 0101101010100101&0101101010100101&0101101010100101&0101101010100101&0101101010100101\\ 0110011001100110&0110011001100110&0110011001100110&0110011001100110&0110011001100110\\ 0110011010011001&0110011010011001&0110011010011001&0110011010011001&0110011010011001\\ \end{array}\right]^{T}.}

引理5.1 s\geq3,n\geq3 时, \frac{s^{n}-1}{s-1}-n-1 总是大于 \frac{s^{n}}{3(s-1)}.

首先我们对不等式 \frac{s^{n}-1}{s-1}-n-1>\frac{s^{n}}{3(s-1)} 进行化简得到 2s^{n}-3ns-3s+3n>0. 令 f(s,n)=2s^{n}-3ns-3s+3n>0.

我们计算 fsn 的偏导数

\begin{array}{l} \frac{\partial f}{\partial s}=\frac{\partial 2 s^{n}-3 n s-3 s+3 n}{\partial s}=2 n s^{n-1}-3 n-3,\\ \frac{\partial f}{\partial n}=\frac{\partial 2 s^{n}-3 n s-3 s+3 n}{\partial n}=2 s^{n} \ln s-3 s+3. \end{array}

先分析 \frac{\partial f}{\partial s}=2ns^{n-1}-3n-3. 当 s\geq3,n\geq3 时, 2ns^{n-1} 是一个正数且迅速增长, 而 -3n-3 是一个负数, 且对于大的 sn, 其增长速度远远小于 2ns^{n-1}. 因此 \frac{\partial f}{\partial s}s\geq3,n\geq3 的范围内是正数.

再分析 \frac{\partial f}{\partial n}=2s^{n}\ln s-3s+3. 当 s\geq3,n\geq3 时, 2s^{n}\ln s 是一个正数且迅速增长, 而 -3s+3 是一个负数, 且对于大的 sn, 其增长速度远远小于 2s^{n}\ln s. 因此 \frac{\partial f}{\partial n}s\geq3,n\geq3 的范围内也是正数.

因为 \frac{\partial f}{\partial n}>0\frac{\partial f}{\partial s}>0, 这表明 f(s,n) 随着 sn 的增长而单调递增. 另外我们还需要检查边界条件来确保 f(s,n)s\geq3,n\geq3 时为正

\begin{align*} f(3,3) =2\cdot3^{3}-3\cdot3\cdot3-3\cdot3+3\cdot3 =27. \end{align*}

由于 f(s,n)s\geq3,n\geq3 的范围内单调递增且在边界条件下为正, 可以推断 f(s,n)>0s\geq3,n\geq3 时总是成立, 因此不等式 2s^{n}-3ns-3s+3n>0 总是成立, 即 \frac{s^{n}-1}{s-1}-n-1>\frac{s^{n}}{3(s-1)} 得证.

定理5.2 s\geq3 是素数或素数幂, n\geq3, 存在只重复 2 个全 0 行的 m-最优正交表 L_{2s^{n}}(s^{\frac{s^{n}-1}{s-1}-x}), 其中 0\leq x<\frac{s^{n}}{3(s-1)}, x 为整数.

根据定理 4.2 我们可以得到最优正交表

L_{2s^{n}}(s^{\frac{s^{n}-1}{s-1}})=\left(\begin{array}{c} B_1\\ B_2\\ \end{array}\right)

其中 B_2 是由 B_1 进行列变换得到的, 如果删除的 x 列不包含进行列变换的 (n+1) 列, 那么 B_1B_2 删除这 x 列后的第 i 行仍不同, 2\leq i\leq s^{n}.

由定义 2.3 可知, 从最优正交表 L_{2s^{n}}(s^{\frac{s^{n}-1}{s-1}}) 中删除 x 列得到的正交表 L_{2s^{n}}(s^{\frac{s^{n}-1}{s-1}-x}) 是一个 m- 最优正交表需要满足

2=\left\lfloor\frac{2s^{n}}{(\frac{s^{n}-1}{s-1}-x)(s-1)+1}\right\rfloor,

\begin{align*} & 2\leq\frac{2s^{n}}{(\frac{s^{n}-1}{s-1}-x)(s-1)+1}<3,\\ & 2\leq\frac{2s^{n}}{s^{n}-(s-1)x}<3,\\ & 0\leq x<\frac{s^{n}}{3(s-1)}. \end{align*}

由引理 2.2 可知 B_1,B_2 的汉明距离为 s^{n-1}, 所以在删除 x 列后的汉明距离总是大于 1, 即 s^{n-1}-x>1, 因此删除 x 列后的 B_1 的第 i 行与 B_2 的第 j 行不同, i\neq j2\leq i,j\leq s^{n}.

最后, 根据引理 5.1 可知 B_1,B_2 中可以删除的 \frac{s^{n}-1}{s-1}-n-1 列总是大于 \frac{s^{n}}{3(s-1)}, 因此不会存在列不够删除的情况. 最优正交表 L_{2s^{n}}(s^{\frac{s^{n}-1}{s-1}}) 删除 x 列(不包含进行列变换的 (n+1) 列)之后得到的 L_{2s^{n}}(s^{\frac{s^{n}-1}{s-1}-x}) 是期望的正交表.

注5.2 由于定理 4.1 中构造最优正交表的方法是先进行水平变换再进行所有列的置换, 我们无法确定在删除的 x 列后是否只有全 0 行重复, 因此不讨论此时对应的 m-最优正交表.

例5.2 s=3,n=3 时, 由定理 5.2 可以得到 m-最优正交表

L_{54}(3^{9})\!=\!\left(\begin{array}{c} B_1'\\ B_2'\\ \end{array}\right)\!=\!\left[\begin{array}{ccccccccccccccccccccccccccc} 000000000111111111222222222&000111222000111222000111222\\ 000111222000111222000111222&012012012012012012012012012\\ 012012012012012012012012012&000111222111222000222000111\\ 000111222111222000222000111&000000000111111111222222222\\ 000111222222000111111222000&000111222222000111111222000\\ 012012012120120120201201201&012012012120120120201201201\\ 012012012201201201120120120&012012012201201201120120120\\ 012120201012120201012120201&012120201012120201012120201\\ 012201120012201120012201120&012201120012201120012201120\\ \end{array}\right]^{T}\!\!\!,

其中只重复 2 个全 0 行.

6 总结

本文主要研究各种具有重复行的强度 2 的最优正交表的构造, 通过水平变换和列变换得到具有重复行的 2 水平的最优正交表、具有重复行的 s 水平的最优正交表和具有重复行的 m-最优正交表的构造方法, 为试验设计提供了大量优秀的正交表. 由于本文构造的均为强度 2 的最优正交表, 今后, 可以基于这些构造方法, 尝试构造具有重复行的更高强度的最优正交表.

参考文献

庞善起, 张应山.

正交表的乘法

数学物理学报, 2007, 27A(3): 568-576

[本文引用: 1]

Zhang Y S.

Multiplication of orthogonal arrays

Acta Math Sci, 2007, 27A(3): 568-576

[本文引用: 1]

Hedayat A S, Sloane N J A, Stufken J. Orthogonal Arrays:Theory and Applications. New York: Springer, 1999

[本文引用: 3]

Rao C R.

Factorial experiments derivable from combinatorial arrangements of arrays

J R Stat Soc, 1947, 9(1): 128-139

[本文引用: 1]

罗纯. 广义差集矩阵理论和正交表构造. 北京: 科学出版社, 2015

[本文引用: 1]

Luo C. Generalized Difference Set Matrix Theory and Orthogonal Array Construction. Beijing: Science Press, 2015

[本文引用: 1]

茆诗松, 周纪芗, 陈颖. 试验设计. 北京: 中国统计出版社, 2004

[本文引用: 1]

Zhou J X, Chen Y. Experiment Design. Beijing: China Statistics Press, 2004

[本文引用: 1]

Arnaud L, Cerf N J.

Exploring pure quantum states with maximally mixed reductions

Phys Rev A, 2013, 87(1): Article 012319

[本文引用: 1]

Kim H J, Lee Y.

t-CIS codes over GF(p) and orthogonal arrays

Discrete Appl Math, 2017, 217(3): 601-612

[本文引用: 1]

聂嘉乐, 余旌胡.

一类受随机扰动的动态优化问题的环境检测与响应

数学物理学报, 2022, 42A(5): 1560-1574

[本文引用: 1]

Nie J L, Yu J H.

Environmental detection and response to a kind of dynamic optimization problem subjected to random disturbance

Acta Math Sci, 2022, 42A(5): 1560-1574

[本文引用: 1]

Du J, Fu S J, Qu L J, et al.

New constructions of q-variable 1-resilient rotation symmetric functions over F_{p}

Sci China Inf Sci, 2016, 59: Article 079102

[本文引用: 1]

Aggarwal M L, Budhraja V.

On construction of some new symmetric and asymmetric orthogonal arrays

JDMSC, 2002, 5(3): 215-225

[本文引用: 1]

Bierbrauer J.

Introduction to Coding Theory

Boca Raton, FL: Chapman and Hall/CRC, 2016

[本文引用: 1]

Stinson D R. Combinatorial Designs. New York: Springer, 2004

[本文引用: 1]

Antonio J D S, Andrea L, Fabio Z.

Symplectic duality between complex domains

Monatsh Math, 2010, 160(4): 403-428

[本文引用: 1]

管乾清, 开晓山, 朱士信.

F_{4^{m}} 上厄米特自正交常循环码

电子学报, 2017, 45(6): 1469-1474

DOI:10.3969/j.issn.0372-2112.2017.06.027      [本文引用: 1]

有限域上常循环码具有丰富的代数结构,其编译码电路容易实现,因而在信息传输实践中具有重要的应用.该文研究了一类有限域上任意长度的厄米特自正交常循环码的结构,给出了此类有限域上厄米特自正交常循环码的生成多项式与存在条件,确立了此类有限域上厄米特自正交常循环码的计数公式,并且利用此类有限域上偶长度的厄米特自正交常循环码构造了最优的量子码.

Kai X S, Zhu S X.

Hermitian self-orthogonal constacyclic codes over F_{4^{m}}

Acta Ecol Sin, 2017, 45(6): 1469-1474

DOI:10.3969/j.issn.0372-2112.2017.06.027      [本文引用: 1]

有限域上常循环码具有丰富的代数结构,其编译码电路容易实现,因而在信息传输实践中具有重要的应用.该文研究了一类有限域上任意长度的厄米特自正交常循环码的结构,给出了此类有限域上厄米特自正交常循环码的生成多项式与存在条件,确立了此类有限域上厄米特自正交常循环码的计数公式,并且利用此类有限域上偶长度的厄米特自正交常循环码构造了最优的量子码.

冯克勤, 陈豪. 量子纠错码. 北京: 科学出版社, 2010

[本文引用: 1]

Chen H. Quantum Error Correction Codes. Beijing: Science Press, 2010

[本文引用: 1]

Qin H, Li D.

Connection between uniformity and orthogonality for symmetrical factorial designs

J Stat Plan Infer, 2006, 136(8): 2770-2782

[本文引用: 1]

Sun F S, Liu M Q, Lin D K J.

Construction of orthogonal Latin hypercube designs

Biometrika, 2009, 96(4): 971-974

[本文引用: 1]

Mukerjee R, Sun F S, Tang B.

Nearly orthogonal arrays mappable into fully orthogonal arrays

Biometrika, 2014, 101(4): 957-963

[本文引用: 1]

Pang S Q, Wang X N, Wang J, et al.

Construction and count of 1-resilient rotation symmetric boolean functions

Inf Sci, 2018, 450: 336-342

[本文引用: 1]

Fang K T, Li R Z, Sudjianto A.

Design and Modeling for Computer Experiments

Boca Raton, FL: Chapman and Hall/CRC, 2005

[本文引用: 1]

Ji L J, Yin J X.

Constructions of new orthogonal arrays and covering arrays of strength three

J Comb Theory Ser A, 2010, 117(3): 236-247

[本文引用: 1]

Goyeneche D, Bielawski J, Zyczkowski K.

Multipartite entanglement in heterogeneous systems

Phys Rev A, 2016, 94(1): Article 012346

[本文引用: 1]

Bouwmeester D, Pan J W, Mattle K, et al.

Experimental quantum teleportation

Nature, 1997, 390: 575-579

[本文引用: 1]

Goyeneche D, Raissi Z, Martino S D, et al. Entanglement and quantum combinatorial designs. Phys Rev A, 2018, 97(6): Article 062326

[本文引用: 1]

Goyeneche D, Zyczkowski K.

Genuinely multipartite entangled states and orthogonal arrays

Phys Rev A, 2014, 90(2): Article 022316

[本文引用: 1]

Jozsa R, Linden N.

On the role of entanglement in quantum-computational speed-up

Proc R Soc A, 2003, 459(2036): 2011-2032

[本文引用: 1]

Li M, Wang Z, Wang J, et al.

The norms of Bloch vectors and classification of four qudits quantum states

Europhysic Letters, 2019, 125(2): Article 20006

[本文引用: 1]

Wang Y, Fang K T.

A note on uniform distribution and experimental design

Kexue Tongbao, 1981, 26(6): 485-489

[本文引用: 1]

周永道, 刘幼妹.

非同构饱和正交设计的最小矩混杂优势准则

数学物理学报, 2009, 29A(5): 1145-1152

[本文引用: 1]

Zhou Y D, Liu Y M.

Minimum moment aberration majorization in non-isomorphic asymmetrical saturated designs

Acta Math Sci, 2009, 29A(5): 1145-1152

[本文引用: 1]

Schoen E D, Eendebak P T, Man M V M.

Complete enumeration of pure-level and mixed-level orthogonal arrays

J Comb Des, 2010, 18(2): 123-140

[本文引用: 1]

Pang S Q, Chen L Y.

Generalized Latin matrix and construction of orthogonal arrays

Acta Math Appl Sin, 2017, 33(4): 1083-1092

[本文引用: 1]

Colbourn C J, Stinson D R, Veitch S.

Constructions of optimal orthogonal arrays with repeated rows

Discrete Math, 2019, 342(9): 2455-2466

DOI:10.1016/j.disc.2019.05.021      [本文引用: 5]

We construct orthogonal arrays OA(lambda)(k, n) (of strength two) having a row that is repeated m times, where m is as large as possible. In particular, we consider OAs where the ratio mp. is as large as possible; these OAs are termed optimal. We provide constructions of optimal OAs for any k >= n +1, albeit with large lambda. We also study basic OAs; these are optimal OAs in which gcd(m, lambda) = 1. We construct a basic OA with n = 2 and k = 4t + 1, provided that a Hadamard matrix of order 8t + 4 exists. This completely solves the problem of constructing basic OAs with n = 2, modulo the Hadamard matrix conjecture. (C) 2019 Elsevier B.V.

Stinson D R.

Bounds for orthogonal arrays with repeated rows

Bull Inst Combin Appl, 2019, 85: 60-73

[本文引用: 3]

Plackett R L, Burman J P.

The design of optimum multifactorial experiments

Biometrika, 1946, 33(4): 305-325

[本文引用: 1]

Bird E M, Street D J.

D-optimal asymmetric orthogonal array plus p run designs

J Stat Plan Infer, 2016, 170: 64-76

[本文引用: 1]

Pang S Q, Zhang X, Lin X, et al.

Two and three-uniform states from irredundant orthogonal arrays

npj Quantum Information, 2019, 5(52): 1-10

[本文引用: 1]

Fisher R A.

An examination of the different possible solutions of a problem in incomplete blocks

Ann Eugen, 1940, 10: 52-75

[本文引用: 1]

Johnson S M.

A new upper bound for error-correcting codes

IRE Trans, 1962, 8: 203-207

[本文引用: 1]

Ray-Chaudhuri D K, Wilson R M.

On t-designs

Osaka J Math, 1975, 12(3): 737-744

[本文引用: 1]

Wilson R M.

Inequalities for t-designs

J Comb Theory Ser A, 1983, 34(3): 313-324

[本文引用: 1]

Wilson R M.

Incidence matrices of t-designs

Linear Algebra Appl, 1982, 46: 73-82

[本文引用: 1]

Cheng C S.

Orthogonal arrays with variable numbers of symbols

Ann Stat, 1980, 8(2): 447-453

[本文引用: 1]

Yu Z, Peng Z, Kristofer J.

Optimal compound orthogonal arrays and single arrays for robust parameter design experiments

Technometrics, 2007, 49(4): 440-453

[本文引用: 1]

El-Zanati S, Jordonh H, Seelinger G, et al.

The maximum size of a partial 3-spread in a finite vector space over GF(2

). Des Codes Cryptogr, 2010, 54(2): 101-107

[本文引用: 1]

Lin C Y.

Optimal blocked orthogonal arrays

J Stat Plan Infer, 2014, 145: 139-147

[本文引用: 1]

Du J, Chen Z Y, Chen G Z, et al.

Constructions of optimal 2-level orthogonal arrays with repeated rows

Adv Math Commun, 2025, 19(2): 379-396

[本文引用: 1]

Stanton R G, Sprott D A.

Block intersections in balanced incomplete block designs

Can Math Bull, 1964, 7, 539-548

[本文引用: 1]

Mann H B.

A note on balanced incomplete block designs

Ann Math Stat, 1969, 40: 679-680

[本文引用: 1]

Montgometry D C. Design and Analysis of Experiments. New York: Springer Verlag, 1976

[本文引用: 1]

Culus J F, Toulouse S.

How far from a worst solution a random solution of a k CSP instance can be?

Lecture Notes in Comput. Sci, 2018, 10979: 374-386

[本文引用: 1]

Mukerjee R, Qian P Z G, Wu C F J.

On the existence of nested orthogonal arrays

Discrete Math, 2008, 308(20): 4635-4642

[本文引用: 1]

Du J, Wen Q Y, Zhang J, et al.

New construction of symmetric orthogonal arrays of strength t

IEICE Trans Fundamentals, 2013, 96(9): 1901-1904

[本文引用: 1]

Mukerjee R, Wu C F J.

On the existence of saturated and nearly saturated asymmetrical orthogonal arrays

Ann Stat, 1995, 23(6): 2102-2115

[本文引用: 2]

Zhang Y L.

On schematic orthogonal arrays of strength two

Ars Comb, 2009, 91: 147-163

[本文引用: 1]

Pang S Q, Yan R, Li S.

Schematic saturated orthogonal arrays obtained by using the contractive replacement method

Commun Stat-Theor M, 2017, 46(18): 8913-8924

[本文引用: 1]

Boyvalenkov P, Marinova T, Stoyanova M.

Nonexistence of a few binary orthogonal arrays

Discrete Appl Math, 2017, 217(2): 144-150

[本文引用: 1]

Pang S Q, Zhang X, Zhang Q J.

The hamming distances of saturated asymmetrical orthogonal arrays with strength 2

Commun Stat-Theor M, 2020, 49(16): 3895-3910

[本文引用: 1]

杨子胥. 正交表的构造. 济南: 山东人民出版社, 1978

[本文引用: 1]

The Construction of Orthogonal Array. Jinan: Shandong People's Publishing House, 1978

[本文引用: 1]

庞善起. 正交表的构造方法及其应用. 成都: 电子科技大学出版社, 2004

[本文引用: 1]

The Construction Method and Application of Orthogonal Array. Chengdu: University of Electronic Science and Technology of China Press, 2004

[本文引用: 1]

/