[论文解读] The Euclidean distance degree of an algebraic variety

Jan Draisma, Emil Horobeţ+3 more|arXiv (Cornell University)|Aug 31, 2013

Polynomial and algebraic computation参考文献 23被引用 20

一句话总结

本文引入了欧氏距离次数（ED degree）作为代数几何与优化中的基本不变量，用于度量从一般点到实代数簇的平方距离函数的临界点数量。通过拉格朗日乘数法与代数几何方法，开发了计算ED degree与平均ED degree的计算工具，应用于低秩矩阵、张量及几何建模问题，证明了对偶性不变性，并推导出关键情形（如Segre簇与秩一张量）的显式公式。

ABSTRACT

The nearest point map of a real algebraic variety with respect to Euclidean distance is an algebraic function. For instance, for varieties of low rank matrices, the Eckart-Young Theorem states that this map is given by the singular value decomposition. This article develops a theory of such nearest point maps from the perspective of computational algebraic geometry. The Euclidean distance degree of a variety is the number of critical points of the squared distance to a generic point outside the variety. Focusing on varieties seen in applications, we present numerous tools for exact computations.

研究动机与目标

• 将欧氏距离次数（ED degree）形式化为实代数簇上优化问题中代数复杂性的度量。
• 利用代数几何与拉格朗日乘数法，发展精确计算ED degree的计算技术。
• 在高斯抽样下研究平均ED degree（aED degree），将其与随机矩阵理论及统计模型相联系。
• 建立对偶性不变性：对射影簇有ED degree(X) = ED degree(X*)，并将其与极类相关联。
• 推导关键应用中的ED degree显式公式，包括低秩矩阵簇与张量的Segre嵌入。

提出的方法

• 将ED degree定义为从一般点到簇的平方欧氏距离函数的非奇异临界点数量。
• 使用拉格朗日乘数法，通过条件“向量(u−x)与切空间TxX正交”来刻画临界点。
• 引入ED对应关系EX，即满足x是du在X上的临界点的点对(x,u)的集合，从而实现平均ED degree的计算。
• 将ED对应关系提升至余切丛，并利用对偶性证明：对射影簇有ED degree(X) = ED degree(X*)。
• 对秩一张量，推导出平均ED degree的公式，表示为具有高斯分布独立条目的随机对称矩阵的行列式绝对值的期望。
• 应用这些工具，使用Macaulay2代码与符号计算，计算特定簇（如心脏线、低秩矩阵、Segre簇）的ED degree。

实验结果

研究问题

• RQ1 在欧氏空间中，寻找给定点到实代数簇的最近点的代数复杂性是什么？
• RQ2 如何对低秩矩阵填充与张量分解等应用中出现的簇精确计算ED degree？
• RQ3 簇的ED degree与其对偶簇的ED degree之间有何关系？ED degree在对偶变换下是否保持不变？
• RQ4 在高斯抽样下，平均ED degree（aED degree）的行为如何？能否通过随机矩阵期望表达？
• RQ5 能否为Segre簇及其他张量模型推导出闭式或高效公式？

主要发现

• 心脏线的ED degree为3，对一般外部点有三个实临界点，ED判别式（外摆线）将具有3个或1个实解的区域分隔开。
• 对于格式为$m_1 \times \cdots \times m_p$的秩一张量的Segre簇，ED degree由定理8.1给出：$\sum_{i=1}^p \binom{m_i}{2} + 1$，显式值如$2\times2\times2$为6，$2^5$为120。
• 在$2\times2\times2$张量格式下，平均ED degree约为4.287，显著低于其ED degree值6，表明大多数临界点为复数。
• 在标准高斯测度下，Segre簇的平均ED degree等于$\frac{\pi^{p/2}}{2^{m/2} \prod_{i=1}^p \Gamma(m_i/2)} \cdot \mathbb{E}[|\det(A)|]$，其中A为具有独立高斯条目的随机对称矩阵。
• 当$p \geq 3$时，平均ED degree在满足$m_p - 1 \geq \sum_{i=1}^{p-1} (m_i - 1)$的边界之外趋于稳定或略有下降，尽管目前尚无几何解释。
• 对偶性成立：对射影簇有$\mathrm{EDdegree}(X) = \mathrm{EDdegree}(X^*)$，且ED degree等于X的经典极类之和。