技术摘要：乘法复合算子的逆运算

问题陈述

本文探讨了 $k$ 阶乘法复合算子的逆问题。给定矩阵 $M \in \mathbb{R}^{\binom{n}{k} \times \binom{m}{k}}$，目标是确定是否存在矩阵 $A \in \mathbb{R}^{n \times m}$ 使得 $M = C_k(A)$，其中 $C_k(A)$ 表示 $A$ 的 $k$ 阶乘法复合。若此类矩阵存在（即 $M$ 是“可复合分解的”），本文旨在刻画所有此类原像 $A$ 的集合，并提供一种构造性算法以恢复它们。

作者指出，尽管复合矩阵在线性代数、动力系统和控制理论中已有深入研究（引用了 Schwarz、Li、Muldowney 等人的工作），但从 $C_k(A)$ 恢复 $A$ 这一特定的逆问题此前尚未被考虑。

方法与算法

提出的解决方案依赖于奇异值分解（SVD）和外代数的性质。该方法基于输入矩阵 $M$ 的非零奇异值互不相同的假设（这一条件可通过涉及可逆变换的通用预处理步骤来确保）。

核心算法将恢复过程分解为四个主要步骤：

1. 楔积分解（奇异向量恢复）：
   算法将 $M$ 的左奇异向量矩阵（$L$）和右奇异向量矩阵（$R$）的列视为 $k$ 阶外幂中的可分解向量。利用可分解向量 $z = u_1 \wedge \dots \wedge u_k$ 唯一确定子空间 $\text{span}\{u_1, \dots, u_k\}$ 的性质，作者构建了“楔积矩阵”$M_z$，其核对应于这些子空间。通过计算这些矩阵的核并寻找它们的一维两两交集，算法恢复了 $A$ 的单个奇异向量（符号除外）。此步骤在定理 2和算法 1中形式化。

2. 奇异值排序：
   由于 $M$ 的奇异值（$M=LSR^\top$ 中 $S$ 的对角元）不一定按与 $A$ 的复合奇异值相同的方式排序，算法利用恢复的奇异向量 $\hat{V}$ 计算 $C_k(\hat{\Sigma})$。这使得 $S$ 的对角元能够重新排序，以匹配复合奇异值的字典序。

3. 符号一致性恢复：
   楔积分解步骤恢复的奇异向量带有任意符号。由于左右奇异向量的符号必须一致才能重构 $A$，因此采用算法 2。它利用复合算子的乘法性质来对齐恢复的左右奇异向量的符号，确保 $C_k(\hat{V})$ 和 $C_k(\hat{W})$ 分别与 $L$ 和 $R$ 的符号匹配。

4. 奇异值恢复：
   一旦确定了奇异向量和有序的复合奇异值，原矩阵 $A$ 的奇异值即可通过求解线性系统来恢复。$A$ 的奇异值（$\sigma_i$）与 $C_k(\Sigma)$ 的对角元之间的关系是乘积关系（$\prod \sigma_i$）。取对数将此转化为线性系统 $y = Lx$，其中$x = \log(\sigma)$，如引理 8所示，该系统有唯一解。

主要贡献与结果

1. 唯一性刻画

本文根据 $M$ 的秩提供了原像集的完整刻画：

• 情形 $\text{rank}(M) \le 1$： 若 $k > \text{rank}(A)$，则 $C_k(A) = 0$。若 $k = \text{rank}(A)$ 且 $\text{rank}(M)=1$，则存在无穷多个矩阵 $A$ 使得 $C_k(A) = M$。具体而言，若 $A = UV^\top$ 是秩揭示分解，则任何满足 $\det(T)=1$ 的矩阵 $B = UTV^\top$ 都会产生相同的复合。
• 情形 $\text{rank}(M) > 1$： 若 $k < \text{rank}(A)$，则原像在整体符号意义下是唯一的。即，若 $C_k(A) = C_k(B)$，则 $A = B$（若 $k$ 为奇数）或 $A = \pm B$（若 $k$ 为偶数）。

2. 构造性算法与复杂度

作者提出了一种多项式时间算法来恢复矩阵 $A$。时间复杂度在定理 1中进行了分析：

• 对于 $k < \lceil \text{rank}(A)/2 \rceil$，复杂度为 $O(\max_{p \in \{m,n\}} [\binom{p}{k}^3 + \binom{p}{k}^2 \binom{r}{k}])$。
• 对于较大的 $k$，复杂度涉及一个额外的对数因子，与子空间的迭代交集有关。
• 该算法已在 MATLAB 中实现并公开可用。

3. 特殊情况

本文指出，对于 $k = n-1$ 且 $M$ 可逆的特定情况，存在闭式解，该解利用伴随矩阵的性质和 Sylvester-Franke 定理推导得出（推论 1）。然而，对于一般的 $n, m, k$，构造性算法是主要贡献。

意义与主张

本文声称解决了一个此前未被触及的线性代数基本逆问题。其意义在于：

• 理论完备性： 确立了逆复合存在的充要条件，并完整刻画了解集（唯一解与无穷多解）。
• 计算可行性： 提供了一种确定性的多项式时间算法来恢复矩阵，从理论存在性迈向了实际计算。
• 应用基础： 作者将这项工作定位为动力系统和控制理论应用的前提。具体而言，他们指出，确定一个系统是否能被变换为"$k$-正系统”（即其 $k$ 阶复合具有统一符号），需要能够将给定矩阵实现为复合的能力。这种逆能力使研究人员能够验证是否存在这样的坐标变换并构建它。

作者在展望未来工作时保持谦逊，明确指出：判断任意给定矩阵是否可复合分解仍是一个未解决的问题，而为不可分解矩阵寻找最佳近似是未来的研究方向。