[論文レビュー] Compounding Doubly Affine Matrices

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

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

ひとこと要約

本稿では、重み付き左および右のクロネッカー積を用いた小規模な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.

研究の動機と目的

• 小規模な種マトリクスを用いてクロネッカー積を用いてより大きなISDA配列を構築する一般的で解析的に取り扱いやすい手法の開発。
• 2次元における合成ISDA行列の固有値および特異値の閉形式表現の導出。
• 任意次元への合成フレームワークの拡張により、高次元の魔方陣立方体および超立方体の構築を可能にする。
• 特にラテン方陣および魔方陣において、構造的性質（例えば対角和、パンディアゴナル性）が合成過程でどのように伝搬されるかの分析。
• 半魔方陣、魔方陣、ラテン方陣に共通する固有値および特異値特性を統合的に適用可能なスペクトル的枠組みの提供。

提案手法

• 本手法は、サイズ$n$および$m$の種マトリクスから、サイズ$mn$のより大きな合成行列を、小規模なISDA行列の左および右のクロネッカー積の重み付き和を用いて生成する。
• 合成行列の固有値は、構成要素のマトリクスの固有値の積として得られ、線和が主固有値として保存される特定の構造を持つ。
• 合成行列の特異値は、種マトリクスの特異値のスケーリングされた積として表現され、全特異値スペクトルの正確な計算が可能になる。
• 合成行列のランクは解析的に$\text{rank}(\mathbf{C}) = \text{rank}(\mathbf{A}) + \text{rank}(\mathbf{B}) - 1$として決定され、数値例で検証されている。
• クロネッカー積の結合的および分配的性質を用いて、複数の次元にわたる$D$次元配列への一般化が可能である。
• 本手法は、例として12次魔方陣、8次ラテン立方体、超立方体の構築に適用され、検証のための明示的な例が提示されている。

実験結果

リサーチクエスチョン

• RQ1 合成ISDA行列の固有値および特異値は、その構成要素のマトリクスのスペクトル的性質を用いてどのように表現できるか？
• RQ2 クロネッカー積による合成によって得られるISDA行列の正確なランクは何か？また、種マトリクスのランクとどのように関係しているか？
• RQ3 パンディアゴナル性や対角一様性といった特別な性質が、合成過程で保存される条件は何か？
• RQ4 本合成手法は任意次元に一般化可能であり、高次元の魔方陣立方体および超立方体の生成を可能にするか？
• RQ5 特異値の4乗和である$R$-インデックスといったスペクトル不変量は、合成過程でどのように変化するか？

主な発見

• 合成行列の固有値は、種マトリクスの固有値の積として形成され、主固有値は線和に対応し、複素共役対も保存される。
• 合成行列の特異値は、$\sigma_{i}^{2} = (a_i b_j)^2$としてスケーリングされたクロネッカー積に対して与えられ、全スペクトルの明示的表現が導出された。
• 合成行列のランクは正確に$\text{rank}(\mathbf{C}) = \text{rank}(\mathbf{A}) + \text{rank}(\mathbf{B}) - 1$であり、12次魔方陣（ランク5）に対して数値的に検証された。
• 3次の洛書方陣（$n=3$）と4次魔方陣（$m=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\}$となる。
• 2次ラテン方陣（$N=2$）2つを用いた合成により、有効な8次ラテン立方体が構築され、行、列、ピラリット和が保存された。集約型および分散型の両方のバリエーションが構築された。
• 本手法は、$k=1$（ラテン）、$k=2$（魔方陣）、$k=3$（魔方陣立方体）、およびそれ以上の$k$を持つISDA配列に普遍的に適用可能であり、任意次元における高次元超立方体の構築を可能にする。