Geometry, Conditioning, and Parity in One-Sided Directional Curvature Recovery, with a Controlled Comparison of Estimators

技术摘要：单侧方向曲率恢复中的几何、条件数与奇偶性

问题陈述
本文探讨了从沿射线发散的函数值采样中恢复海森矩阵（二阶导数信息）的基本极限。这种情况出现在导数自由优化、边界附近的有限差分近似以及逆问题中，其中域的补集是不可达的。在这些场景下，允许的方向受限于几何结构（例如，激活不等式约束、薄结构或边界），从而无法形成对称的模板。核心问题不在于如何设计特定的估计器，而在于在这些几何约束下理论上什么是可能的：哪些量是可识别的，条件数如何随受限访问而退化，以及可以达到什么样的精度阶数。

方法论与框架
本研究采用几何与分析方法，而非提出新的数值估计器。其核心方法论包括：

1. 测量模型： 使用定义在一条射线上的三个点构成的单侧方向探针：中心点 $x$，内部节点 $x + R h u$，以及端点 $x + h u$。该探针通过二阶差分来估计方向曲率 $u^T H u$。
2. 几何分析： 基于采样方向的分布分析海森矩阵的可识别性。作者利用对称 $3 \times 3$ 矩阵空间 ($\text{Sym}_3(\mathbb{R})$) 和球面谐波进行测量算子的分解。
3. 奇偶分解： 利用球面上齐次多项式的奇偶对称性（偶部/奇部）。探针展开被拆分为偶部（携带海森矩阵信息）和奇部（携带领先截断误差信息）。
4. 受控比较： 通过“公平协议”对八种估计器（包括二次回归、移动最小二乘法及各种方向方案）进行严格的基准测试。每种方法都获得了其自身的算子最优采样半径，并且共享相同的噪声实现和可行性约束。

主要贡献与结果

1. 可识别性界限（定理 4，命题 5）：

   • 评估底限： 在 $\mathbb{R}^3$ 中，在没有平滑先验的情况下，至少需要 10 个处于一般位置的点评估才能唯一确定海森矩阵。
   • 方向唯一性： 对于方向测量，唯一性的条件不仅是方向需张成 $\mathbb{R}^3$，而且是这些方向 不落在同一个二次锥面上。
   • 成本含义： 使用 $m$ 个方向的方向方案需要 $1 + 2m$ 次评估。由于唯一性要求 $m \ge 6$，因此最小成本为 13 次评估，超过了通用点集的 10 次评估底限。

2. 受限访问下的条件数（定理 10）：

   • 当允许的方向被限制在半角为 $\theta$ 的球面帽内时，恢复算子的奇异值按 $\{1, \theta, \theta, \theta^2, \theta^2, \theta^2\}$ 比例缩放。
   • 因此，条件数增长为 $\kappa(\theta) \asymp \theta^{-2}$。论文推导出了这种增长的精确常数（例如，对于表面均匀采样，$\kappa \theta^2 \to \sqrt{24} \approx 4.899$）。
   • 解释： 窄锥能很好地测量法向曲率，但对混合曲率和切向曲率的测量较差（分别为 $\theta$ 和 $\theta^2$ 阶）。平滑边界（半球，$\theta \approx 90^\circ$）是良性的；只有在角点、尖点或高度受限的激活集处才会出现条件数问题。

3. 奇偶分解与截断（命题 13, 14, 推论 15）：

   • 单侧探针展开为偶部（与海森矩阵相关）和奇部（与截断项相关）的球面谐波分量。
   • 无海森测量： 奇部投影可以独立于海森矩阵隔离出三阶方向导数（$D^3 f$），精度为 $O(h^2)$。这使得在无需预知海森矩阵的情况下测量截断项成为可能。
   • 可达阶数（命题 17）： 聚合的海森矩阵估计能够达到二阶精度，当且仅当 可行方向集具有 对跖对称性（即 $u \in U \implies -u \in U$）。如果集合是真正的单侧（没有对跖对），则无论使用多少个方向或权重方案，误差始终保持为 $O(h)$。

4. 噪声与试点设计（命题 19, 21, 22）：

   • 论文量化了探针的方差以及由此产生的可用步长（跨度）的下界。最优跨度按 $\sigma^{1/3}$ 比例缩放，这与标准有限差分理论一致。
   • 通过奇部投影提取截断常数需要一个“试点（pilot）”阶段。论文指出该试点过程非常昂贵，需要至少 18 个方向才能保持稳定，这可能超过了海森矩阵估计本身的成本。

5. 估计器的受控比较：

   • 在公平协议下，方向方案并未优于 标准的二次回归或 $\Lambda$-适定插值（$\Lambda$-poised interpolation）。
   • 在对称模板可行的情况下，中心差分最为精确。在不对称情况下，二次回归和移动最小二乘法的表现与方向方案持平或更优。
   • 方向方案未能获得优势的原因在于它们无法克服单侧访问固有的 $O(h)$ 截断误差（命题 17），并且受到与其他方法相同的几何条件数限制。

意义与主张
本文明确将贡献定位为对测量几何的 结构性界限，而非针对特定算法的改进。

• 对最优性的约束： 这些结果定义了“不可能”区域：没有估计器能在缺乏对跖对称性的情况下达到二阶精度，也没有估计器能绕过窄锥带来的 $\theta^{-2}$ 条件数惩罚。
• 未利用的结构： 识别出奇部球面谐波分量是截断项的无海森测量，这一发现被视为一种结构性洞察，尽管论文警告说提取该项在计算上是非常昂贵的。
• 方法论的严谨性： 论文的大部分篇幅致力于记录九个特定的实验伪影（例如参数传播误差、几何生成漏洞以及错误的公平性假设），这些伪影必须被剔除才能产生可靠的结果。作者认为，这种“负面结果”——即在公平条件下方向方案并不比回归法更具优势——比一个潜在的虚假正面结果更有价值。

谦逊态度与局限性
作者谨慎地区分了证明的定理与数值观察。他们并不声称单侧方向采样在边界附近具有优越性；事实上，其结果和比较表明事实并非如此。他们承认其奇偶分解结果可能是现有 Slepian 理论或受限 Veronese 嵌入文献的一个实例，并将其界定为“开放的位置问题”而非全新的发现。论文总结道，虽然单侧恢复的几何特性已得到良好刻画，但方向方案的实际效用受限于几何、截断与噪声之间的基本权衡。