[論文レビュー] An Iterative Method for Nonconvex Quadratically Constrained Quadratic Programs

Chuangchuang Sun, Ran Dai|arXiv (Cornell University)|Sep 8, 2016

Advanced Optimization Algorithms Research参考文献 23被引用数 3

ひとこと要約

本稿では、非凸な二次制約付き二次プログラム（QCQP）を解くための反復的ランク最小化（IRM）手法を提案する。各部分問題は凸な半定値計画（SDP）として定式化され、増大ラグランジュ法に基づく拡張Uzawaアルゴリズムを用いて解かれる。本手法は、保証された収束速度を伴う線形収束を達成し、商業的NLPソルバーよりも優れた目的関数値の低減を示す。

ABSTRACT

This paper examines the nonconvex quadratically constrained quadratic programming (QCQP) problems using an iterative method. One of the existing approaches for solving nonconvex QCQP problems relaxes the rank one constraint on the unknown matrix into semidefinite constraint to obtain the bound on the optimal value without finding the exact solution. By reconsidering the rank one matrix, an iterative rank minimization (IRM) method is proposed to gradually approach the rank one constraint. Each iteration of IRM is formulated as a convex problem with semidefinite constraints. An augmented Lagrangian method, named extended Uzawa algorithm, is developed to solve the subproblem at each iteration of IRM for improved scalability and computational efficiency. Simulation examples are presented using the proposed method and comparative results obtained from the other methods are provided and discussed.

研究の動機と目的

• 保証された収束性と向上した計算効率を備えた非凸QCQP問題の解法に取り組む。
• 行列緩和問題のランク1解に反復的に近づく手法を開発し、実行可能性と最適性を保証する。
• 効率的な部分問題ソルバを用いて、大規模QCQP問題におけるスケーラビリティと計算効率を向上させる。
• 反復的手法に理論的収束保証と線形収束速度を提供する。
• 宇宙機の姿勢制御のような実世界の最適制御問題（幾何的制約を伴う）に本手法を検証する。

提案手法

• IRM手法は、半正定値制約を満たす行列変数のランクを反復的に最小化することで、QCQPのランク1解に近づける。
• 各反復では、元のQCQP定式化から導かれた半正定値制約を有する凸SDP部分問題を解く。
• 増大ラグランジュ法に基づく拡張Uzawaアルゴリズムを用いて、各SDP部分問題を効率的かつスケーラブルに解く。
• 双対性理論とKarush-Kuhn-Tucker（KKT）条件を活用し、局所最適解への線形収束を確立する。
• 拡張ベクトルと正規化制約 $t^2 = 1$ を用いて、非同次QCQPを同次形式に再定式化する。
• アルゴリズムは、離散化と制約埋め込みを経て、最適宇宙機姿勢制御問題に適用される。

実験結果

リサーチクエスチョン

• RQ1 非凸QCQP問題に対して、収束が保証された反復的手法を開発できるか？
• RQ2 提案されたIRM手法は、局所最適解への線形収束速度を達成するか？
• RQ3 拡張Uzawaアルゴリズムは、IRM各反復で生じる大規模SDP部分問題を効果的に解けるか？
• RQ4 SNOPTなどの商業的NLPソルバーよりも、IRM手法の解の品質と計算コストはどのように異なるか？
• RQ5 複雑な幾何的制約を伴う実世界の最適制御問題に、IRM手法は成功裏に適用可能か？

主な発見

• IRM手法は、双対性とKKT解析によって裏付けられた理論的収束速度を伴い、局所最適解への線形収束を達成する。
• 拡張Uzawaアルゴリズムにより、各反復における大規模SDP部分問題の解法においてスケーラビリティと計算効率が確保される。
• 宇宙機の姿勢制御シミュレーションでは、IRM手法が目的関数値23.72を達成し、SNOPTの最良結果（31.79）と比較して25.39%低減した。
• IRM手法の収束には743.6秒を要したが、SNOPTは収束ケースの平均で8.4秒を要しており、解の品質と速度のトレードオフが示された。
• ランク1からの逸脱度を測る $ r_k $ の値は、数反復内でゼロに収束し、本手法がランク1解に近づける能力を裏付けた。
• 本手法は、3次元空間における角度制限やベクトルの向きに基づいて定義された禁止領域・必須領域を含む複雑な制約を効果的に処理した。