[논문 리뷰] Minimizing Polynomial Functions

Pablo A. Parrilo, Bernd Sturmfels|ArXiv.org|2001. 03. 26.

Advanced Optimization Algorithms Research참고 문헌 2인용 수 72

한 줄 요약

이 논문은 다항함수의 전역 최소화를 위한 정수형 프로그래밍 기반의 이완 방법을 제안하며, 제곱합과 포지티브스텔렌츠를 활용한다. 이 방법은 기존의 대수적 방법들인 그로버 기반, 결과식, 호모토피 방법들과 비교해 높은 차수의 다항함수에서 뛰어난 성능을 보이며, 차수가 고정된 경우 전역 최소값을 다항시간에 계산할 수 있음을 보여준다.

ABSTRACT

We compare algorithms for global optimization of polynomial functions in many variables. It is demonstrated that existing algebraic methods (Gröbner bases, resultants, homotopy methods) are dramatically outperformed by a relaxation technique, due to N.Z. Shor and the first author, which involves sums of squares and semidefinite programming. This opens up the possibility of using semidefinite programming relaxations arising from the Positivstellensatz for a wide range of computational problems in real algebraic geometry. This paper was presented at the Workshop on Algorithmic and Quantitative Aspects of Real Algebraic Geometry in Mathematics and Computer Science, held at DIMACS, Rutgers University, March 12-16, 2001.

연구 동기 및 목표

• 전역 다항 최적화를 위한 전통적인 대수적 방법과 새로운 이완 기법을 비교하고 대조한다.
• 정수형 프로그래밍 기반의 이완 기법이 그로버 기반, 결과식, 호모토피 방법보다 계산 효율성이 뛰어나다는 것을 입증한다.
• 포지티브스텔렌츠가 실대수기하 문제에서 비가능성 증명서로 사용될 수 있음을 확립한다.
• 다항함수의 전역 최소화가 차수가 고정된 경우 정수형 프로그래밍을 통해 다항시간에 해결될 수 있음을 보여준다.
• 실다항식 문제를 해결하기 위한 새로운 계산 패러다임을 제안한다: 포지티브스텔렌츠 + 정수형 프로그래밍.

제안 방법

• 전역 최소화 문제를 f(x) − λ 가 ℝ[x₁,…,xₙ] 내에서 제곱합(SOS)이 되는 가장 큰 λ 를 찾는 문제로 재구성한다.
• 정수형 프로그래밍(SDP)을 사용해 그러한 λ 를 계산하며, 이는 전역 최소값 f* 에 대한 하한을 제공한다.
• SDP의 이중성을 적용해 λ = f* 일 때 최적성을 확인하고, 이 때 최소화자 p* 를 복원한다.
• 포지티브스텔렌츠를 활용해 다항부등식 및 다항방정식 시스템의 비가능성 증명서를 SOS 항등식을 통해 구성한다.
• 차수 D 가 점차 증가하는 SDP 이완을 사용해 다항식 시스템의 비가능성 여부를 테스트하거나 제약 조건 하에 f 를 최소화한다.
• SOS 다항식을 양의 준정부호 행렬로 표현하고, 내부점 방법을 사용해 결과적인 SDP 를 해결한다.

실험 결과

연구 질문

• RQ1 제곱합 기반 정수형 프로그래밍 이완 기법이 전통적인 대수적 방법들인 그로버 기반, 결과식, 호모토피 방법보다 전역 다항 최적화에서 성능이 뛰어나다고 할 수 있는가?
• RQ2 제곱합 이완 기법이 다항함수의 정확한 전역 최소값 f* 를 얼마나 잘 도출하는가?
• RQ3 포지티브스텔렌츠는 다항방정식 및 부등식 시스템의 비가능성 증명서를 효과적으로 구성하는 데 사용될 수 있는가?
• RQ4 제곱합과 정수형 프로그래밍을 사용해 다항함수의 전역 최소값을 계산할 때의 계산 복잡도는 무엇이며, 특히 차수가 고정된 경우 어떻게 되는가?
• RQ5 제한된 차수 D 를 가진 SDP 이완 기법이 실대수기하 문제를 해결하는 데 실용적이고 확장 가능한 방법이 될 수 있는가?

주요 결과

• 제곱합 기반 정수형 프로그래밍 이완 기법은 계산 실험에서 기존의 대수적 방법들인 그로버 기반, 호모토피 방법보다 일관되게 뛰어난 성능을 보였다.
• 차수 고정 다항식의 경우, 정수형 프로그래밍을 통해 전역 최소값을 다항시간에 계산할 수 있다.
• 계산 실험에서 이완 기법은 거의 모든 테스트 케이스에서 정확한 전역 최소값(λ = f*)을 도출하였으며, 최대 n = 15개의 변수를 가진 문제에서도 성공하였다.
• SDP 이중성을 통한 최적성 증명서를 제공하며, λ = f* 일 때 최소화자 p* 를 복원할 수 있다.
• 포지티브스텔렌츠를 통해 SOS 항등식을 구성함으로써 다항식 시스템의 비가능성을 체계적으로 증명할 수 있으며, 이는 SDP 를 통해 계산할 수 있다.
• 차수 2의 SDP 이완 기법이 {x − y² + 3 ≥ 0, y + x² + 2 = 0} 시스템의 비가능성을 명시적인 SOS 항등식을 통해 성공적으로 증명하였다.