[论文解读] Compounding Doubly Affine Matrices

Adam Rogers, Ian Cameron+1 more|arXiv (Cornell University)|Nov 29, 2017

Matrix Theory and Algorithms参考文献 7被引用 3

一句话总结

本文提出了一种通用的解析方法，通过加权求和的左、右Kronecker积复合较小的ISDA矩阵，构建更大的整数序列双重仿射（ISDA）数组（如幻方、拉丁方和超立方体）。关键贡献在于推导出复合矩阵的特征值与奇异值的显式公式，实现精确的秩计算，并揭示了在不同维度间保持的谱特性，适用于任意维度下的高阶幻方与拉丁结构。

ABSTRACT

Weighted sums of left and right hand Kronecker products of Integer Sequence Doubly Affine (ISDA) as well as Generalized Arithmetic Progression Doubly Affine (GAPDA) arrays are used to generate larger ISDA arrays of multiplicative order (compound squares) from pairs of smaller ones. In two dimensions we find general expressions for the eigenvalues (EVs) and singular values (SVs) of the larger arrays in terms of the EVs and SVs of their constituent matrices, leading to a simple result for the rank of these highly singular compound matrices. Since the critical property of the smaller constituent matrices involves only identical row and column sums (often called semi-magic), the eigenvalue and singular value results can be applied to both magic squares and Latin squares. Additionally, the compounding process works in arbitrary dimensions due to the generality of the Kronecker product, providing a simple method to generate large order ISDA cubes and hypercubes. The first examples of compound magic squares are found in manuscripts that date back to the 10th century CE, and other representative applications are outlined through judicious examples.

研究动机与目标

• 开发一种通用且可解析处理的方法，利用Kronecker积从较小的种子矩阵构造更大的ISDA数组。
• 推导出二维复合ISDA矩阵的特征值与奇异值的闭式表达式。
• 将复合框架扩展至任意维度，实现高阶幻方立方体与超立方体的构造。
• 分析结构特性（如对角线和、全对角线性）在复合过程中的传播，尤其针对拉丁与幻方数组。
• 提供一个统一的谱框架，适用于半幻方、幻方与拉丁方，通过共享的特征值与奇异值特性。

提出的方法

• 该方法使用较小ISDA矩阵的左、右Kronecker积的加权和，从阶为$n$与$m$的种子矩阵生成阶为$mn$的更大复合矩阵。
• 复合矩阵的特征值被表示为构成矩阵特征值的乘积，其特定结构保持线和作为主导特征值。
• 复合矩阵的奇异值表示为种子矩阵奇异值的缩放乘积，从而实现完整奇异值谱的精确计算。
• 复合矩阵的秩被解析确定为$\text{rank}(\mathbf{C}) = \text{rank}(\mathbf{A}) + \text{rank}(\mathbf{B}) - 1$，并通过数值示例验证。
• 该方法利用Kronecker积在多维空间中的结合律与分配律，推广至$D$维数组。
• 该方法被应用于构造复合幻方（如12阶）、拉丁立方体（如8阶）与超立方体，提供了显式示例以验证其有效性。

实验结果

研究问题

• RQ1 复合ISDA矩阵的特征值与奇异值如何用其构成矩阵的谱特性表示？
• RQ2 通过Kronecker积复合形成的复合ISDA矩阵的精确秩是多少？其与种子矩阵秩的关系如何？
• RQ3 在复合过程中，全对角线性或对角线均匀性等特殊性质在何种条件下得以保持？
• RQ4 该复合方法能否推广至任意维度，以生成高阶幻方立方体与超立方体？
• RQ5 谱不变量如$R$-指数（奇异值四次方之和）在复合过程中如何表现？

主要发现

• 复合矩阵的特征值由种子矩阵特征值的乘积构成，主导特征值对应于线和，复共轭对保持不变。
• 复合矩阵的奇异值由$\sigma_{i}^{2} = (a_i b_j)^2$给出（针对缩放的Kronecker积），并推导出完整的奇异值谱表达式。
• 复合矩阵的秩精确为$\text{rank}(\mathbf{C}) = \text{rank}(\mathbf{A}) + \text{rank}(\mathbf{B}) - 1$，通过12阶幻方（秩为5）的数值验证。
• 由Lo Shu幻方（$n=3$）与$4$阶幻方复合而成的12阶幻方，其特征值为$\lambda_1 = 858$，$\lambda_2 = 4(2i\sqrt{6})$，$\lambda_3 = -4(2i\sqrt{6})$，奇异值满足$\sigma_i^2 = \{858^2, 9^2 \times 320, 9^2 \times 20, 4^2 \times 48, 4^2 \times 12\}$。
• 复合过程成功从两个$N=2$的拉丁立方体构造出有效的8阶拉丁立方体，保持了行、列与柱的和，同时构建了聚合与分散两种变体。
• 该方法普遍适用于ISDA数组，其中$k=1$（拉丁方）、$k=2$（幻方）、$k=3$（幻方立方体）及更高$k$，可实现任意维度下高阶超立方体的构造。