[논문 리뷰] Simplification Methods for Sum-of-Squares Programs

Peter Seiler, Qian Zheng외 1명|arXiv (Cornell University)|2013. 03. 04.

Advanced Optimization Algorithms Research참고 문헌 11인용 수 6

한 줄 요약

이 논문은 정방형 행렬의 영 대각성질을 기반으로 한 빠르고 계산 효율적인 단항식 축소 알고리즘을 제안한다. 이는 양의 준정적 행렬의 영 대각 성질을 이용하여, 해당 행렬의 대각 성분이 반드시 0이 되는 단항식을 식별하고 제거함으로써, 뉴턴 다면체 방법보다 항상 크지 않으며 종종 훨씬 더 작은 단항식 집합을 생성하고 계산 시간을 크게 줄인다. 이 방법은 또한 암묵적으로 0이 되어야 하는 변수를 제거하는 더 넓은 단순화 프레임워크로 확장되어 수치적 조건을 개선하고 솔버 성능을 향상시킨다.

ABSTRACT

A sum-of-squares is a polynomial that can be expressed as a sum of squares of other polynomials. Determining if a sum-of-squares decomposition exists for a given polynomial is equivalent to a linear matrix inequality feasibility problem. The computation required to solve the feasibility problem depends on the number of monomials used in the decomposition. The Newton polytope is a method to prune unnecessary monomials from the decomposition. This method requires the construction of a convex hull and this can be time consuming for polynomials with many terms. This paper presents a new algorithm for removing monomials based on a simple property of positive semidefinite matrices. It returns a set of monomials that is never larger than the set returned by the Newton polytope method and, for some polynomials, is a strictly smaller set. Moreover, the algorithm takes significantly less computation than the convex hull construction. This algorithm is then extended to a more general simplification method for sum-of-squares programming.

연구 동기 및 목표

• 분해 벡터 내 과도한 단항식 항으로 인해 발생하는 SOS 프로그래밍의 계산 병목 현상을 해결하기 위해.
• 비용이 많이 드는 볼록 껍질 계산이 필요한 뉴턴 다면체 방법에 비해 더 빠른 대안을 개발하기 위해.
• 행렬의 구조로 인해 반드시 0이 되어야 하는 단항식과 변수를 식별하고 제거하여 SOS 최적화 문제의 크기를 줄이기 위해.
• 기호 기반 단순화를 통해 수치적 조건을 향상시키고 솔버 시간을 줄임으로써 SOS 프로그래밍의 효율성을 높이기 위해.
• 영 대각 성질을 일반적인 단순화 프레임워크로 확장하여 약속이 있는 SOS 프로그램에 적용 가능하게 하기 위해.

제안 방법

• 영 대각 알고리즘은 SOS 공식화에서 등식 제약 조건에 의해 영이 되어야 하는 양의 준정적 행렬 Q의 대각 성분을 식별한다.
• 수학적 성질에 따르면, 양의 준정적 행렬의 대각 성분이 0이면 그 행과 열 전체가 0이 되어야 하므로, 이 성질을 이용해 단항식을 제거할 수 있다.
• 이 알고리즘은 MATLAB의 'find' 명령어를 활용하여 이러한 제약 조건을 효율적으로 탐지하며, 볼록 껍질 구축에 비해 최소한의 계산만을 요구한다.
• 이 방법은 단일 및 다중 변수를 포함하는 등식 제약 조건을 처리하여 부호 정보와 추가적인 영 제약 조건을 유도하도록 확장된다.
• 변수가 음이 아닌지, 음이 아닌지 또는 0인지 추적하기 위해 부호 벡터를 유지하며, 반복적 단순화를 가능하게 한다.
• 알고리즘은 반복적으로 단항식을 제거하고 제약 행렬 A의 해당 열을 0으로 설정함으로써 문제 크기를 줄이며, 타당성을 유지한다.

실험 결과

연구 질문

• RQ1 뉴턴 다면체 방법보다 더 빠르면서도 단항식 집합의 크기가 더 이상 크지 않은 단항식 축소 방법을 개발할 수 있는가?
• RQ2 양의 준정적 행렬의 영 대각 성질이 더 단순하고 효율적인 단항식 제거 전략을 가능하게 하는가?
• RQ3 영 대각 원리를 단항식 축소를 넘어서 전체 SOS 프로그램의 단순화로 일반화할 수 있는가?
• RQ4 부호 제약 조건에 기반한 변수 제거가 SOS 프로그램의 수치적 조건과 계산 시간을 얼마나 향상시킬 수 있는가?
• RQ5 제안된 단순화 방법은 뉴턴 다면체 및 SOSTOOLs와 같은 기존 방법들과 비교해 성능과 해의 크기 측면에서 어떻게 다른가?

주요 결과

• 영 대각 알고리즘은 뉴턴 다면체 방법보다 항상 크지 않은 단항식 집합을 반환하며, 일부 다항식에서는 엄청나게 더 작은 집합을 생성한다.
• 이 알고리즘은 볼록 껍질 구축을 피하기 때문에 뉴턴 다면체 방법보다 훨씬 더 적은 계산을 요구한다. 단지 선형 제약 조건 분석에 의존하기 때문이다.
• 이 방법은 SOSOPT 소프트웨어 패키지에 구현되어 실제 SOS 프로그래밍 문제에서 실용적인 효율성 향상을 보였다.
• 일반적인 단순화 절차는 암묵적으로 0이 되어야 하는 변수를 성공적으로 식별하고 제거하여 수치적 안정성을 향상시키고 솔버 시간을 줄였다.
• 반복적인 부호 추적 및 제약 조건 처리 논리는 하나 이상의 제약 조건이 있는 SOS 프로그램의 체계적인 단순화를 가능하게 한다.
• 이 방법은 문제 차원을 줄이면서도 타당성과 최적성을 유지하여 정확도를 잃지 않고 더 빠른 해답 시간을 제공한다.