비선형 제어 시스템을 위한 순차 이차 SOS 프로그래밍

초록

본 논문은 비선형 제어 시스템에서 발생하는 비볼록 합동제곱(SOS) 최적화 문제를 해결하기 위해, 필터 라인서치와 결합된 순차 이차 프로그래밍 알고리즘을 제안한다. 제안 방법은 각 반복에서 이차형 서브문제를 풀어 빠른 수렴을 유도하고, 전역 수렴성을 보장하기 위해 필터 기반 라인서치를 적용한다. 수치 실험을 통해 기존 좌표 하강법이나 선형 SOS 방법에 비해 반복 횟수와 계산 시간이 크게 감소함을 확인했으며, 구현 코드를 오픈소스로 제공한다.

상세 분석

논문은 비선형 시스템 분석·제어 설계에서 다루어지는 지역적 안정성, 도달 가능 집합, 제약조건을 포함한 제어 설계 문제들을 다항식 형태로 모델링하고, 이를 SOS 제약식으로 변환한다는 점에서 시작한다. 기존의 비볼록 SOS 문제 해결 방법은 좌표 하강법, 이분법, 혹은 혼합형 접근법에 의존했으며, 이들 방법은 초기값 의존성, 수렴 보장 부재, 반복 횟수 폭증 등의 한계를 가지고 있었다. 저자는 이러한 문제점을 극복하기 위해 두 가지 핵심 아이디어를 도입한다. 첫째, 각 반복 단계에서 선형화가 아닌 이차화된 서브문제(Quadratic SOS)를 풀어, 라그랑지안의 헤시안 근사(H_k)를 활용해 뉴턴‑유사 업데이트를 수행한다. 이는 전통적인 순차 이차 프로그래밍(SQP)과 유사하지만, SOS 제약을 다루기 위해서는 SOS 원뿔에 대한 투영 연산이 필요하고, 이를 효율적으로 구현하기 위해 SDP 혹은 2차원 원뿔 제약 형태로 변환한다. 둘째, 전역 수렴성을 확보하기 위해 필터 라인서치(filter line‑search) 전략을 적용한다. 필터는 목적 함수값과 제약 위반 정도를 동시에 관리하며, 페널티 파라미터에 민감하지 않은 복구 단계(restoration phase)를 제공한다. 이로써 비볼록 SOS 문제에서도 지역 최소점에 머무르지 않고, 넓은 초기값 영역에서 수렴을 기대할 수 있다.

수학적 분석 부분에서는 KKT 조건을 일반화한 변분 방정식 φ(ξ,λ)+N(ξ,λ)∋0 형태로 표현하고, Φ(ξ_k,H_k)를 φ의 1차 근사로 두어 집합값 뉴턴 업데이트식 (6)을 도출한다. 강한 정규성(strong regularity)과 Lipschitz 연속성을 가정한 뒤, 정리 2를 통해 ℓ_k와 ζ의 관계에 따라 선형, 초선형, 이차 수렴을 보인다. 특히 H_k를 정확한 2차 미분(헤시안)으로 선택하면 ℓ_k≡0이 되어 이차 수렴을 달성한다는 점이 강조된다.

알고리즘 구현 측면에서는 헤시안 근사 방법으로 BFGS와 제한된 메모리 BFGS(L‑BFGS)를 제시하고, SOS 제약 위반 검사는 기존 SDP 솔버의 프리듀얼 변수 값을 이용해 효율적으로 수행한다. 필터 라인서치에서는 현재 후보와 복구 후보를 각각 (f,‖c‖) 쌍으로 저장하고, 새로운 후보가 두 기준을 동시에 개선하면 수용한다. 이 과정은 페널티 파라미터 튜닝 없이도 안정적인 스텝 길이 α_k를 결정한다.

실험에서는 지역 안정성 추정, 도달 가능 집합 근사, 제어 장벽 함수 설계 등 네 가지 대표적인 비선형 제어 문제에 대해 기존 좌표 하강법, 순차 선형 SOS, 그리고 최신 비선형 SDP 솔버와 비교한다. 결과는 평균 반복 횟수가 3050% 감소하고, 전체 실행 시간이 4060% 단축됨을 보여준다. 특히 고차 차원(예: 6차 다항식, 10개 상태 변수) 문제에서도 메모리 사용량이 크게 증가하지 않아 실시간 적용 가능성을 시사한다. 마지막으로 저자는 MATLAB 기반 오픈소스 패키지를 공개하여 연구 재현성을 확보하고, 향후 비볼록 SOS 문제에 대한 전역 수렴 증명과 대규모 분산 구현을 향한 연구 방향을 제시한다.