[论文解读] Global Solutions to Large-Scale Spherical Constrained Quadratic Minimization via Canonical Dual Approach

Yi Chen, David Yang Gao|arXiv (Cornell University)|Aug 20, 2013

Matrix Theory and Algorithms参考文献 20被引用 6

一句话总结

本文提出一种规范对偶方法，用于求解大规模球面约束二次最小化问题，特别针对标准方法失效的“困难情形”。通过将问题重新表述为无对偶间隙的一维规范对偶问题，该方法提供了识别困难情形的充分必要条件，并引入扰动方法以实现全局收敛。

ABSTRACT

This paper presents global optimal solutions to a nonconvex quadratic minimization problem over a sphere constraint. The problem is well-known as a trust region subproblem and has been studied extensively for decades. The main challenge is the so called 'hard case', i.e., the problem has multiple solutions on the boundary of the sphere. By canonical duality theory, this challenging problem is able to reformed as an one-dimensional canonical dual problem without duality gap. Sufficient and necessary conditions are obtained by the triality theory, which can be used to identify whether the problem is hard case or not. A perturbation method and the associated algorithms are proposed to solve this hard case problem. Theoretical results and methods are verified by large-size examples.

研究动机与目标

• 解决大规模球面约束二次最小化问题的全局解，特别是当在球面边界存在多个解的挑战性“困难情形”。
• 通过开发基于扰动的方法，克服困难情形中的数值困难，以恢复唯一解的存在性。
• 利用三枝性理论建立识别困难情形的充分必要条件。
• 通过规范对偶与扰动技术，提供一种适用于大规模问题的全局收敛算法。

提出的方法

• 通过规范对偶变换，将原始非凸二次最小化问题转化为无对偶间隙的规范对偶问题。
• 应用三枝性理论，基于对偶函数及其临界点的结构，推导出识别困难情形的充分必要条件。
• 通过在目标向量 f 中添加小参数 αi，引入扰动方法，确保 ∑αi² ≠ 0，以在扰动问题中恢复唯一解的存在性。
• 利用扰动问题近似原始困难情形问题的解，且在扰动参数满足特定界时保证收敛性。
• 使用对偶函数 Pd(σ) = ∑_{i=k+1}^n [f̂i² / (λi + σ)²] - r² 定位临界点 σ̄，其对应于最优对偶变量。
• 通过 x̄ = G_a(σ̄)^{-1}f 求解原始解，其中 G_a(σ) 为变换后的海森矩阵，确保全局最优性。

实验结果

研究问题

• RQ1 球面约束二次最小化问题在何种条件下具有唯一解，何时会变为“困难情形”？
• RQ2 规范对偶方法能否在无对偶间隙的情况下提供困难情形的全局解？
• RQ3 如何选择扰动参数，以确保扰动解收敛于困难情形的真实解？
• RQ4 何种理论条件可保证规范对偶函数中临界点的存在性与唯一性？

主要发现

• 规范对偶问题被重述为无对偶间隙的一维优化问题，从而实现全局解的恢复。
• 当 ∑_{i=1}^k f̂i² = 0 且 ∑_{i=k+1}^n f̂i² / (λi - λ1)² ≤ r² 时，识别出困难情形，此时对偶函数缺乏临界点。
• 针对困难情形，提出通过引入非零 αi 的扰动方法，以在对偶函数中恢复唯一临界点的存在性。
• 扰动条件 ∑αi² ≤ (λ2 - λ1)²(r² - ∑_{i=k+1}^n f̂i² / (λi - λ1)²)(1/√(2(1 - cos(ε/r)) - 1)^{-2} 确保扰动解在 ε 容差范围内收敛于真实解。
• 理论分析证实，对偶函数 Pd(σ) 的临界点 σ̄ 唯一，且对应于全局最优原始解 x̄ = G_a(σ̄)^{-1}f。
• 大规模问题的数值实验验证了所提方法的有效性与收敛性，即使在困难情形下亦成立。