불가능 이중 적분기 최적 제어와 더글라스‑라프슨 알고리즘

초록

본 논문은 제어 구간의 상·하한이 과도하게 제한되어 발생하는 이중 적분기 시스템의 불가능 최적 제어 문제를 다룬다. 제어와 상태 제약 집합 사이의 갭 함수를 L² 노름으로 최소화하는 ‘최적 근사’ 해를 구하고, 그 해가 최대 한 번의 전환을 갖는 bang‑bang 형태임을 분석한다. 두 개의 대수 방정식을 풀어 전환 시점과 갭 함수를 얻는 방법을 제시하고, 뉴턴·일반화 뉴턴 방법을 통한 수치 해법을 검증한다. 또한, 완화된 더글라스‑라프슨(DR) 알고리즘을 적용하여 불가능 문제에 대한 수치 실험을 수행하고, 기존 이산화‑후‑최적화 방식 대비 효율성을 입증한다.

상세 분석

논문은 먼저 일반적인 선형 시스템에 대한 최적 제어와 제약 집합 A(동역학·경계조건 만족)와 B(제어 상·하한) 사이의 거리(갭) 최소화 문제를 정리한다. A와 B가 교차하지 않을 때, 즉 문제 자체가 불가능할 때 ‘최적 근사 쌍(uA, uB)’을 정의하고, 갭 벡터 v = uA − uB 를 L²‑노름 최소화 문제(Pf)로 전환한다. 이때 라그랑주 승수 λ(t) 를 도입하면 v(t)=−Bᵀ(t)λ(t) 라는 닫힌 형태가 도출되고, λ는 adjoint 방정식 ˙λ=−Aᵀλ 를 만족한다. 중요한 결과는 제어가 구성요소별로 제어가능(component‑wise controllable)하면, u*B는 v의 부호에 따라 상·하한을 그대로 취하는 bang‑bang 형태가 최적임을 보인다.

특히 이중 적분기(¨y = u) 문제에 적용하면, 최적 제어는 상·하한 사이를 전환하는 단일 스위칭 bang‑bang 형태가 된다. 저자들은 이를 두 개의 다항식 방정식(스위칭 시점 τ와 갭 함수의 두 파라미터)으로 환원한다. 이 방정식은 비선형이지만, 뉴턴 방법과 일반화 뉴턴 방법(가중치 행렬을 이용한 변형)으로 빠르게 수렴한다는 실험적 증거를 제시한다. 또한 a→0(제어 상·하한이 매우 작아지는 경우)에서 스위칭 시점 τ가 a의 1/3 제곱에 비례한다는 점근적 분석을 제공해, 작은 제어 한계에서도 근사 해를 예측할 수 있다.

알고리즘적 측면에서는 기존 연구에서 제시된 DR 알고리즘을 ‘완화된’ 형태로 재구성한다. 여기서는 두 투영 연산(π_A, π_B)을 번갈아 적용하고, 과도한 제약으로 인해 발생하는 비교적 큰 갭을 점진적으로 감소시킨다. 실험에서는 동일한 초기·최종 상태와 제어 한계 a에 대해, DR 알고리즘이 Newton 기반 직접 해법보다 CPU 시간과 메모리 사용량에서 우수함을 보이며, 특히 고해상도 시간 이산화(수천 개 구간)에서도 수렴성을 유지한다.

전체적으로 논문은 (1) 불가능 선형 제어 문제를 최적 근사 프레임워크로 정형화, (2) 이중 적분기에 대한 명시적 bang‑bang 해와 스위칭 시점 계산법을 제공, (3) 수치 해법의 수렴 특성을 체계적으로 분석, (4) DR 알고리즘을 불가능 문제에 적용해 실용적인 계산 도구로 확장한다는 네 가지 주요 기여를 한다.