[论文解读] A bound on the minimum of a real positive polynomial over the standard simplex

Saugata Basu, Richard Leroy+1 more|ArXiv.org|Feb 19, 2009

Advanced Optimization Algorithms Research参考文献 5被引用 12

一句话总结

本文利用实代数几何技术，推导出在标准单纯形上对正多元多项式最小值的显式、非退化的下界。通过将一元结果推广至结式和西尔维斯特矩阵，并对结式多项式应用柯西不等式，作者建立了仅依赖于维度 $k$、次数 $d$ 和系数位大小 $ au$ 的可计算下界，无需非退化性假设。

ABSTRACT

We consider the problem of bounding away from 0 the minimum value m taken by a polynomial P of Z[X_1,...,X_k] over the standard simplex, assuming that m>0. Recent algorithmic developments in real algebraic geometry enable us to obtain a positive lower bound on m in terms of the dimension k, the degree d and the bitsize of the coefficients of P. The bound is explicit, and obtained without any extra assumption on P, in contrast with previous results reported in the literature.

研究动机与目标

• 为标准单纯形 $\Delta \subset \mathbb{R}^k$ 上的正多项式 $P \in \mathbb{Z}[X_1,\dots,X_k]$ 的最小值 $m > 0$ 提供一个显式下界。
• 消除先前方法所需关于 $P$ 的非退化性或特殊假设。
• 推导出仅依赖于 $k$、$d$ 和 $ au$（系数位大小）的下界 $m_{k,d,\tau}$。
• 通过算法实代数几何方法，提供一个显式、可计算且普适的下界，优于以往结果。

提出的方法

• 通过分析结式 $R(Z) = \operatorname{Res}_T(P(T) - Z, P'(T))$ 推广一元情形，以定位最小值出现的临界点。
• 利用西尔维斯特矩阵结构和哈达玛行列式界，控制结式多项式 $R(Z)$ 的系数。
• 对结式多项式应用柯西不等式，推导出 $m$（最小正根）的下界。
• 通过在单纯形的每个面上使用有理一元表示（RUR），将一元方法推广至多元多项式。
• 控制 RUR 多项式的次数和系数位大小 ($d_u, \tau_u, d_{P,u}, \tau_{P,u}$)，以控制结式的大小。
• 通过 $D, \rho, \rho'$ 的递归界构造最终下界 $m_{k,d,\tau} = \left[2^{\rho'+1}\sqrt{dD+1}\right]^{-D}\left[2^{\rho}\sqrt{D+1}\right]^{-dD}$。

实验结果

研究问题

• RQ1 能否在不依赖非退化性假设的前提下，为标准单纯形上正多项式的最小值导出一个普适的显式下界？
• RQ2 正多项式在单纯形上的最小值如何与其次数 $d$、维度 $k$ 和系数位大小 $ au$ 相关？
• RQ3 使用算法实代数几何时，可获得的 $m = \min_{\Delta} P$ 的最紧显式下界是什么？
• RQ4 与 Canny 及 de Loera-Santos 的先前结果相比，所提方法在指数形式和参数依赖性方面有何定量差异？

主要发现

• 本文建立了完全显式的下界 $m_{k,d,\tau} = \left[2^{\rho'+1}\sqrt{dD+1}\right]^{-D}\left[2^{\rho}\sqrt{D+1}\right]^{-dD}$，其中 $D = 2d(2d-1)^{k-1}$，$\rho$ 和 $\rho'$ 通过 $k,d,\tau$ 的递归表达式定义。
• 该下界对所有在单纯形上的正多项式均有效，无需非退化性或横截性条件，与 Canny 的间隙定理或洛瓦谢维奇型界不同。
• 该下界优于以往结果，提供了对 $k$、$d$ 和 $ au$ 的显式依赖关系，并给出定量估计 $\frac{1}{m_{k,d,\tau}} \leq \left(2^{\tau}\right)^{2^{k+3}d^{k+1}k}2^{2^{k+7}d^{k+2}k^{2}}$。
• 与 de Loera 和 Santos 的一般估计相比，该下界更紧，因为其指数部分是显式给出的，而非包含未知的普适常数。
• 该方法通过将系统 $\partial P_\sigma / \partial Y_i = 0$ 替换为单个方程 $Q = \sum_i (\partial P_\sigma / \partial Y_i)^2 = 0$ 实现，从而覆盖退化情形。
• 一元情形下得到的下界为 $m > \frac{3^{d/2}}{2^{(2d-1)\tau}(d+1)^{2d-1/2}}$，略优于 [BCR] 中近期结果。