SOSTOOLS v4.00 — 고속 SOS 최적화와 새로운 dpvar 구조 소개

초록

SOSTOOLS v4.00은 다항식 최적화와 SOS(합의 제곱) 프로그램을 위한 MATLAB 툴박스로, 파싱 단계의 메모리·시간 복잡도를 크게 낮춘 새로운 dpvar 구조를 도입했다. 파싱 시간이 전체 SDP 해결 시간의 10 % 이하로 감소했으며, MOSEK·SeDuMi·SDPT3·CSDP·SDPNAL·SDPNAL+·SDPA 등 주요 SDP 솔버와 연동한다. 또한 다항식 다중곱(multiploly)과 맞춤형 결정 변수 선언을 지원해 대규모 비선형 안정성 분석 등 실용적인 응용에 적합하다.

상세 분석

SOSTOOLS v4.00은 기존 버전에서 가장 큰 병목으로 지적되던 파싱 단계의 연산량과 메모리 사용을 근본적으로 재설계하였다. 핵심은 dpvar 라는 새로운 다항식 구조체를 도입한 점이다. 기존에는 sossosvar, sospolyvar, sosmatrixvar 등이 MATLAB Symbolic Math Toolbox에 의존해 다항식과 결정 변수를 동시에 관리했으며, 이는 수천 개 이상의 결정 변수가 등장하는 대규모 SOS 문제에서 메모리 폭증과 속도 저하를 초래했다. dpvar 구조는 **결정 변수(스칼라 SDP 변수)**와 **독립 변수(다항식의 입력 변수)**를 명확히 분리한다. 내부적으로는 결정 변수들을 벡터 형태로 저장하고, 다항식 계수를 효율적인 인덱스 매핑을 통해 접근한다. 이 설계 덕분에 파싱 복잡도는 결정 변수 수에 대해 거의 선형으로 증가한다. 실제 벤치마크에서 파싱 시간이 전체 솔버 실행 시간의 10 % 이하로 감소했으며, 특히 차수 4 이상의 다변량 다항식(10 변수 이상)으로 구성된 비선형 안정성 검증 문제에서 5배~10배 이상의 속도 향상이 보고되었다.

또한, v4.00은 멀티폴리(multiploly) 툴박스와 완전 호환된다. 멀티폴리는 심볼릭 연산 없이 다항식 연산을 수행하도록 설계된 경량 객체이며, dpvar와 결합하면 파싱 단계에서 심볼릭 연산을 완전히 배제할 수 있다. 이는 메모리 사용량을 크게 줄이고, MATLAB R2021a 이상의 최신 버전과도 원활히 동작한다.

솔버 인터페이스 측면에서, 기존에 지원하던 SeDuMi, SDPT3, CSDP, SDPNAL, SDPNAL+, SDPA 외에 MOSEK를 새롭게 추가하였다. MOSEK는 대규모 SDP에 특화된 상용 솔버로, 특히 고밀도 행렬과 복합 제약조건을 포함하는 문제에서 뛰어난 성능을 보인다. 사용자는 sospsimplify 함수를 통해 SOS 제약을 사전 정리하고, sosquadvar 으로 사용자 정의 결정 변수를 선언할 수 있어, 기존 sossosvar · sospolyvar · sosmatrixvar 의 한계를 넘어서는 유연성을 제공한다.

응용 분야는 크게 네 가지로 정리된다. 첫째, Lyapunov 함수 탐색을 통한 비선형 시스템 안정성 증명; 둘째, 전역·제한 최적화에서 SOS를 이용한 다항식 경계 추정; 셋째, 행렬 양정성 검증 및 구조적 특이값 상한 계산; 넷째, MAX‑CUT, Chebyshev 다항식, 확률적 경계 등 전통적인 조합 최적화와 확률론적 문제에 SOS 이완을 적용한다. 각 장에서는 MATLAB 코드 예제와 함께 실제 데이터에 적용한 결과를 제시하며, 특히 대규모 문제(수천 개 변수, 수만 개 제약)에서도 파싱이 병목이 되지 않음을 강조한다.

전반적으로 SOSTOOLS v4.00은 파싱 효율성, 확장성, 솔버 호환성 세 축을 균형 있게 강화함으로써, 연구자와 엔지니어가 SOS 기반 설계·분석을 손쉽게 수행하도록 만든다. 향후 SDP 솔버와 하드웨어 가속기의 발전과 맞물려, SOS 방법론이 제어·신호·통신·기계·전력 시스템 등 다양한 분야에 보다 널리 적용될 기반을 제공한다.