This paper presents a new interpretation of the Lyapunov and Riccati equations from the perspective of positive system theory. We show it is possible to construct positive systems related to these equations, and then certain conclusions -- such as the existence and uniqueness of solutions -- can be drawn from positive systems theory. Specifically, under standard observability assumptions, a strictly positive linear system can be constructed for Lyapunov equations, leading to exponential convergence in Hilbert metric to the Perron-Frobenius vector -- closely related to the solution of the Lyapunov equation. For algebraic Riccati equations, homogeneous strictly positive systems can be constructed, which exhibit more complex dynamical behaviors. While the existence and uniqueness of the solution can still be proven, only asymptotic convergence can be obtained.
- 논문 ID: 2511.17243
- 제목: Lyapunov and Riccati Equations from a Positive System Perspective
- 저자: Dongjun Wu (룬드 대학교), Yankai Lin (우한 공업 대학교)
- 분류: math.OC (최적화 및 제어)
- 제출 시간: 2025년 11월 24일 arXiv 제출
- 논문 링크: https://arxiv.org/abs/2511.17243
본 논문은 양의 시스템 이론(positive system theory)의 관점에서 Lyapunov 방정식과 Riccati 방정식에 대한 새로운 해석을 제시한다. 연구 결과에 따르면, 이들 방정식과 관련된 양의 시스템을 구성할 수 있으며, 양의 시스템 이론으로부터 해의 존재성과 유일성 등의 결론을 도출할 수 있다. 구체적으로, 표준 가관측성 가정 하에서 Lyapunov 방정식에 대해 순정 양의 선형 시스템을 구성할 수 있으며, Hilbert 메트릭 하에서 Perron-Frobenius 벡터로 지수 수렴하고, 이 벡터는 Lyapunov 방정식의 해와 밀접한 관련이 있다. 대수적 Riccati 방정식의 경우, 동차 순정 양의 시스템을 구성할 수 있으며, 시스템이 더 복잡한 동역학 거동을 나타내지만, 여전히 해의 존재성과 유일성을 증명할 수 있고, 수렴성은 지수 수렴에서 점근 수렴으로 감소한다.
Lyapunov 방정식과 대수적 Riccati 방정식은 제어 이론의 기초이며, 시스템 분석, 모델 축소, 안정화 및 최적 제어에 광범위하게 적용된다. 이들 방정식의 성질(해의 존재성, 유일성 및 수치 계산)은 광범위하게 연구되었지만, 전통적 방법에는 한계가 있다:
- Lyapunov 방정식은 명시적 해를 가지지만(연속시간에서는 적분 형태, 이산시간에서는 무한 급수), Riccati 방정식은 거의 명시적 해를 갖지 않음
- 두 종류의 방정식은 일반적으로 서로 다른 처리 방법이 필요하며, 통일된 프레임워크가 부족함
이 두 종류의 방정식은 선형 시스템 이론의 핵심이며, 해의 성질은 제어 시스템의 안정성 분석 및 최적 제어 설계에 직접 영향을 미친다. 통일된 이론 프레임워크를 구축하는 것은 이론적 이해를 심화시킬 뿐만 아니라 수치 알고리즘에 새로운 아이디어를 제공할 수 있다.
- Pates (2024)의 방법: Lyapunov 방정식을 컴팩트 볼록 집합 위의 고정점 문제로 보고, Brouwer 고정점 정리를 사용하여 해의 존재성을 증명한다. 그러나 이 기법은 Riccati 방정식으로 직접 확장될 수 없으며, Brouwer 정리는 유일성을 보장할 수 없다.
- 전통적 방법: Lyapunov 방정식과 Riccati 방정식을 각각 처리하며, 통일된 동역학 시스템 관점이 부족하다.
본 논문은 동역학 시스템을 구성하여 이 두 종류의 방정식을 통일적으로 처리하고, 양의 시스템 이론의 강력한 도구를 활용하여 해의 성질을 분석할 것을 제안한다. 핵심 혁신은 보조 변수(스칼라 선형 함수)를 도입하여 방정식을 동차화하여 양의 시스템 프레임워크에 적용 가능하게 하는 것이다.
- 새로운 이론적 연결 수립: 양의 시스템 이론과 고전 제어 방정식(Lyapunov 및 Riccati 방정식)을 처음으로 연결하여 완전히 새로운 이론적 관점을 제공한다.
- 동차화 기술: 스칼라값 선형 함수(예: tr(P))를 도입하여 방정식을 동차화하는 핵심 기술 단계를 제시하여, 비동차 방정식이 양의 시스템 프레임워크 하에서 통일적으로 처리될 수 있도록 한다.
- 양의 동역학 시스템 구성:
- Lyapunov 방정식에 대해 순정 양의 선형 시스템을 구성하고 지수 수렴성을 증명한다.
- 대수적 Riccati 방정식에 대해 동차 순정 양의 시스템을 구성하고 점근 수렴성을 증명한다.
- 통일된 해법 프레임워크: 방정식의 해를 적절히 구성된 양의 동역학 시스템의 정상 상태로 특성화하여, 동역학 시스템 관점에서의 해법 방법을 제공한다.
- 수렴성 분석: Hilbert 메트릭을 사용하여 동역학 시스템의 수렴 성질을 분석하고, 수치 알고리즘에 이론적 보장을 제공한다.
연속시간 Lyapunov 방정식:
A⊤P+PA+C⊤C=0
이산시간 Lyapunov 방정식:
A⊤QA−Q+C⊤C=0
이산시간 대수적 Riccati 방정식:
P=A⊤PA−A⊤PB(R+B⊤PB)−1B⊤PA+C⊤C
연속시간 대수적 Riccati 방정식:
A⊤P+PA−PBR−1B⊤P+C⊤C=0
여기서 P,Q∈S+n(반정부호 행렬 집합), A,B,C,R은 알려진 행렬이다.
원뿔(Cone) 정의: 본 논문은 주로 두 가지 원뿔을 사용한다:
- 양의 정사분면: R+n={x∈Rn:xi≥0,∀i}
- 반정부호 행렬 집합: S+n={X∈Rn×n:X⪰0}
Hilbert 메트릭: 원뿔 K의 x,y=0에 대해 다음과 같이 정의된다:
dH(x,y)=lnm(x/y)M(x/y)
여기서 M(x/y)=inf{λ:x≤λy}, m(x/y)=sup{μ:μy≤x}
핵심 정리 (정리 1 - 비선형 Perron-Frobenius 이론):
K를 정규 원뿔, F:V→V를 연속, 동차, 보존 순서 매핑이라 하자. 그러면 F(x∗)=ρ(F)x∗인 x∗∈K∖{0}이 존재한다. F가 강 보존 순서이면, x∗는 양의 스칼라 의미에서 유일하고, dH(Fk(x),x∗)→0이다. F가 또한 선형이면, 수렴은 지수적이다.
연속시간 경우:
단계 1: 선형화 시스템 구성
간단한 선형 방정식 Ap+c=0(A는 Metzler 행렬, c≥0)을 고려한다. 전통적 방법은 동역학 시스템 x˙=Ax+c를 구성하지만, 이는 동차가 아니다. "선형화" 기법을 통해:
x˙=(A+αc1⊤)x
여기서 α>0은 조절 매개변수, 1은 모두 1인 벡터이다.
단계 2: Lyapunov 방정식에 대한 양의 시스템 구성P˙=A⊤P+PA+αtr(P)C⊤C(5)
핵심 보조정리 (보조정리 1): (A,C)가 가관측이면, 시스템 (5)은 S+n에서 순정 양이다.
증명 개요: 해의 표현식은
P(t)=eA⊤tP0eAt+α∫0ttr(P(s))eA⊤(t−s)C⊤CeA(t−s)dsv=0이 존재하여 P(t)v=0이면, 모든 s∈[0,t]에 대해 CeAsv=0이 성립하는데, 이는 가관측성과 모순된다. 따라서 P(t)≻0이다.
단계 3: 매개변수 조절 및 해 구성
추론 1에 의해, 다음을 만족하는 유일한 P∗≻0과 λ∈R이 존재한다:
A⊤P∗+P∗A+αtr(P∗)C⊤C=λP∗
A의 고유벡터 v(고유값 η)에 대해,
λ=2Re(η)+αtr(P∗)v∗P∗v∥Cv∥2
- α가 충분히 작을 때, λ<0(A가 안정적이므로)
- α가 충분히 클 때, λ>0(∥Cv∥2=0이므로)
- 연속성에 의해, λ=0인 α>0이 존재한다.
이때 P#=αtr(P∗)P∗는 Lyapunov 방정식의 유일한 정부호 해이다.
수렴성: 시스템 (5)는 Hilbert 메트릭 하에서 지수 수렴한다:
dH(P(t),P#)≤ke−λtdH(P0,P#)
이산시간 경우:
반복 시스템을 구성한다:
Xk+1=F(Xk)=A⊤XkA+αtr(Xk)C⊤C
유사한 분석에 따르면, 가관측성과 Schur 안정성 가정 하에서, 유일한 해가 존재하고 시스템은 지수 수렴한다.
이산시간 경우:
단계 1: 동차 시스템 구성Xk+1=F(Xk)=A⊤S(Xk)A+αtr(Xk)C⊤C
여기서
S(X):=X−XB(αtr(X)R+B⊤XB)−1B⊤X
핵심 보조정리 (보조정리 3): 연산자 F는 연속, 동차, 보존 순서 및 오목이다. (A,C)가 가관측이면, Fn은 강 보존 순서이다.
증명 요점:
- 보존 순서성은 Schur 여인수의 단조성에서 나온다(보조정리 6).
- 오목성은 Schur 여인수의 알려진 성질이다.
- 강 보존 순서성은 재귀 논증을 통해: v⊤Xk+1v=0이면, Cv=0이고 v⊤A⊤XkAv=0이며, n번 재귀한 후
col(C,CA,…,CAn−1)v=0
가관측성에 의해 v=0이므로, Fn(X)∈IntS+n이다.
단계 2: 매개변수 조절
- α가 작을 때, A가 Schur 안정적이면(가제어성 가정으로 보장 가능), ρ(F)<1
- α→∞일 때, ρ(F)→∞(가관측성에 의해)
- 원뿔 스펙트럼 반경의 연속성에 의해, ρ(F)=1인 α>0이 존재한다.
이때 P=αtr(X∗)X∗는 Riccati 방정식의 유일한 정부호 해이다.
수렴성: Hilbert 메트릭 하에서 점근 수렴한다(지수적 아님):
dH(Xk,P)→0 as k→∞
연속시간 경우:
시스템을 구성한다:
P˙=αA⊤P+αPA−tr(P)1PBR−1B⊤P+α2tr(P)C⊤C
동등하게,
P˙=Aα(t)⊤P+PAα(t)+α2tr(P)C⊤C
여기서 Aα(t)=αA−2tr(P)1BR−1B⊤P(t)
핵심 보조정리 (보조정리 4): (A,C)가 가관측이면, 시스템은 순정 양이다.
유사한 매개변수 조절 논증에 따르면, α>0이 존재하여 시스템이 유일한 정상 상태 해 P를 가지며, Hilbert 메트릭 하에서 점근 수렴한다.
- 동차화 기법: tr(P)를 스칼라화 인수로 도입하여 원래의 비동차 방정식을 동차 동역학 시스템으로 변환하는 것이 양의 시스템 이론을 적용하는 핵심이다.
- 통일된 프레임워크: 양의 동역학 시스템을 구성함으로써 Lyapunov 방정식과 Riccati 방정식을 통일된 분석 프레임워크에 포함시킨다.
- 매개변수 조절 방법: 매개변수 α를 조절하여 원뿔 스펙트럼 반경을 특정 값(Lyapunov 방정식의 경우 0, Riccati 방정식의 경우 1)에 도달하게 하여, 원래 방정식의 해를 교묘하게 구성한다.
- Hilbert 메트릭 분석: 전통적인 유클리드 메트릭이 아닌 Hilbert 메트릭을 사용하여 수렴성을 분석하는 것이 양의 시스템의 기하학적 구조에 더 적합하다.
- 가관측성의 새로운 역할: 가관측성이 시스템의 순정 양성을 보장하는 데 있어 핵심적 역할을 하는 것을 드러내어, 고전적 가정에 새로운 기하학적 해석을 제공한다.
본 논문은 주로 이론 작업이며, 수치 예제를 통해 이론 결과를 검증한다:
예제 1 (그림 1): 간단한 선형 방정식 Ap+c=0에 대해, 행렬 A+αc1⊤의 최대 고유값이 α에 따라 변하는 곡선을 보여주며, 고유값이 0인 α가 존재함을 검증한다.
예제 2 (그림 2): 이산시간 Lyapunov 방정식에 대해, 연산자 L(X)=A⊤XA+αtr(X)C⊤C의 원뿔 스펙트럼 반경이 α에 따라 변하는 곡선을 보여주며, A가 Schur 안정적이고 (A,C)가 가관측일 때 스펙트럼 반경이 1인 α가 존재함을 보여준다.
- 매개변수 연속성: 그림 1과 2는 스펙트럼 반경/고유값이 매개변수 α에 대해 연속임을 명확히 보여주며, 이론의 중간값 정리 적용을 지지한다.
- 수렴 거동: 논문이 상세한 수렴 곡선을 제공하지는 않지만, 이론은 Lyapunov 방정식의 지수 수렴과 Riccati 방정식의 점근 수렴을 보장한다.
정리적 결과:
- 사실 1 (연속시간 Lyapunov 방정식): (A,C)가 가관측이고 A가 안정적이면, 유일한 정부호 해 P#이 존재하며, 구성된 동역학 시스템은 Hilbert 메트릭 하에서 P#로 지수 수렴한다.
- 사실 2 (이산시간 Lyapunov 방정식): A가 Schur 안정적이고 (A,C)가 가관측이면, 유일한 정부호 해 Q가 존재하며, 시스템은 지수 수렴한다.
- 사실 3 (이산시간 Riccati 방정식): (A,B)가 가제어이고 (A,C)가 가관측이면, 유일한 정부호 해 P가 존재하며, 시스템은 Hilbert 메트릭 하에서 점근 수렴한다.
- 사실 4 (연속시간 Riccati 방정식): (A,B)가 가제어이고 (A,C)가 가관측이면, 유일한 정부호 해 P가 존재하며, 시스템은 점근 수렴한다.
- 가관측성과 가제어성의 쌍대 작용 (비고 5):
- 가관측성은 원뿔 스펙트럼 반경을 증가시키는 경향이 있다.
- 가제어성은 원뿔 스펙트럼 반경을 감소시키는 경향이 있다.
- 둘 다 필수적이며, 함께 해의 존재성과 유일성을 보장한다.
- 수렴성 차이:
- Lyapunov 방정식: 선형 시스템 → 지수 수렴
- Riccati 방정식: 비선형 동차 시스템 → 점근 수렴(지수적 아님)
- 확장성: 방법은 더 일반적인 Lyapunov 방정식으로 확장될 수 있다:
- 확률 시스템: A⊤P+PA+∑i=1rAi⊤PAi+C⊤C=0
- 특이 시스템: A⊤PE+E⊤PA+C⊤C=0
양의 제어 시스템의 Bellman 방정식:
p=s+A⊤p−E⊤∣r+B⊤p∣
가정 H1(안정화 이득의 존재)과 H2(가관측성 유사 조건) 하에서, 시스템을 구성함으로써
pk+1=α(1⊤pk)s+A⊤pk−E⊤∣α(1⊤pk)r+B⊤pk∣
해의 존재성과 유일성을 증명하여, 양의 시스템 최적 제어에서 방법의 응용을 보여준다.
- 문헌 9,11,18,1,14: 이들 방정식의 성질, 수치 알고리즘 및 응용을 체계적으로 연구한다.
- 전통적 방법: 명시적 해(Lyapunov) vs 반복 알고리즘(Riccati)
- Pates 19: Lyapunov 방정식을 컴팩트 볼록 집합 위의 고정점 문제로 보고, Brouwer 고정점 정리를 사용한다.
- 한계: 유일성을 보장할 수 없으며, Riccati 방정식으로 확장하기 어렵다.
- Rantzer 및 Valcher 21: 양의 시스템 이론의 종합
- Lemmens 및 Nussbaum 16: 비선형 Perron-Frobenius 이론
- 응용: 합의 문제 24,22, 최적 제어 20,6
- Bertsekas 5, Lee 15, Lai 및 Xiong 12: Riccati 방정식의 값 반복 수렴성
- 본 논문의 기여: 동차화를 통해 값 반복을 양의 시스템 이론 프레임워크에 적용 가능하게 한다.
- Rantzer 20, Li 및 Rantzer 17: 양의 제어 시스템의 명시적 Bellman 방정식
- 본 논문의 확장: 양의 시스템 이론을 사용하여 통일적으로 분석한다.
- 이론적 통일: Lyapunov 방정식과 대수적 Riccati 방정식을 양의 시스템 이론 프레임워크에 성공적으로 포함시켜 통일된 분석 방법을 제공한다.
- 해의 성질: 표준 가정(가관측성, 가제어성, 안정성) 하에서 해의 존재성, 유일성 및 수렴성을 증명한다.
- 수렴 특성:
- Lyapunov 방정식: 지수 수렴(선형 시스템)
- Riccati 방정식: 점근 수렴(비선형 동차 시스템)
- 방법론적 기여: 동차화 기술은 대수 방정식 처리를 위한 새로운 도구를 제공한다.
- 수렴 속도: Riccati 방정식의 경우, 점근 수렴만 보장되며 지수 수렴이 아니므로, 수치 구현에서 느린 수렴 속도를 초래할 수 있다.
- 매개변수 선택: 이론은 적절한 α의 존재를 보장하지만, 논문은 α를 효율적으로 선택하는 방법을 제공하지 않는다.
- 수치 실험 부족: 논문은 주로 이론 작업이며, 방법의 실제 효율성과 기존 알고리즘과의 비교를 검증하는 대규모 수치 실험이 부족하다.
- 확장성: 시변 시스템, 비선형 시스템 등 더 일반적인 경우로의 확장은 아직 탐구되지 않았다.
- 오차 한계: 논문은 반복 알고리즘의 오차 한계와 수렴 속도의 구체적 상수를 정량화하지 않는다.
논문이 명시한 연구 방향:
- 시변 방정식: 시변 Lyapunov 및 Riccati 방정식으로 확장
- 오차 한계 정량화: 제시된 반복 방법에 대한 정확한 오차 한계 수립
- 알고리즘 구현: 매개변수 α의 적응적 선택 전략을 포함한 효율적인 수치 알고리즘 개발
- 응용 확장: 모델 예측 제어, 강건 제어 등 분야에서의 응용 탐구
- 이론적 혁신성이 강함:
- 양의 시스템 이론과 고전 제어 방정식의 깊은 연결을 처음으로 수립한다.
- 동차화 기술은 보편적이며, 다른 대수 방정식 연구에 영감을 줄 수 있다.
- 가관측성, 가제어성 등 고전적 가정에 새로운 기하학적 해석을 제공한다.
- 방법의 통일성:
- 동일한 프레임워크로 Lyapunov 및 Riccati 방정식을 처리한다.
- 연속시간 및 이산시간 시스템의 통일적 처리
- 더 일반적인 방정식 형태로 확장 가능하다.
- 수학적 엄밀성:
- 증명이 완전하고 논리가 명확하다.
- 양의 시스템 이론의 깊은 결과(Perron-Frobenius 이론, Hilbert 메트릭)를 충분히 활용한다.
- 모든 핵심 보조정리는 상세한 증명을 가진다.
- 이론적 통찰이 깊음:
- 가관측성과 가제어성이 스펙트럼 반경 조절에서 하는 쌍대 작용을 드러낸다.
- Lyapunov 및 Riccati 방정식의 수렴성 차이의 본질적 원인을 설명한다.
- 동역학 시스템 관점에서 대수 방정식을 이해하기 위한 새로운 패러다임을 제공한다.
- 글쓰기가 명확함:
- 간단한 예제(선형 방정식)에서 복잡한 경우로 점진적으로 진행한다.
- 구조가 명확하고 논리가 단계적이다.
- 기호 체계가 통일되어 이해하기 쉽다.
- 실용성 검증 필요:
- 기존 수치 알고리즘(Bartels-Stewart 알고리즘, Newton 반복)과의 효율성 비교가 부족하다.
- 매개변수 α의 실제 선택 전략을 제공하지 않는다.
- 수렴 상수 k,λ의 추정 방법을 제시하지 않는다.
- 수치 실험 부족:
- 두 개의 간단한 수치 예제만 제시된다.
- 고차원 시스템, 병태 시스템 등 복잡한 경우의 테스트가 부족하다.
- 이론적 수렴 속도와 실제 수렴 속도의 일치성을 검증하지 않는다.
- Riccati 방정식의 수렴성:
- 점근 수렴만 보장되며, 수렴 속도의 정량적 추정이 불가능하다.
- 실제 응용에서는 지수 수렴이 더 이상적이다.
- 가정 조건:
- 표준적인 가관측성, 가제어성 및 안정성 가정에 의존한다.
- 가정이 만족되지 않을 때의 경우를 논의하지 않는다.
- 계산 복잡도:
- 알고리즘의 계산 복잡도를 분석하지 않는다.
- 각 반복에서 행렬 대각합을 계산해야 하므로 계산 비용이 증가할 수 있다.
- 일반화 가능성:
- 비동차 방정식, 제약 최적화 문제 등으로의 일반화가 명확하지 않다.
- 시변, 확률, 비선형 경우의 확장은 추가 연구가 필요하다.
- 이론적 기여:
- 제어 이론과 양의 시스템 이론 사이에 새로운 다리를 구축한다.
- 다른 대수 방정식의 동역학 시스템 방법 연구에 영감을 줄 수 있다.
- 동차화 기술은 방법론적 가치를 가진다.
- 학술적 가치:
- Lyapunov 및 Riccati 방정식에 새로운 이론적 관점을 제공한다.
- 가관측성, 가제어성 등 고전적 개념에 대한 이해를 심화시킨다.
- 양의 시스템 이론 응용에 관한 후속 연구를 자극할 수 있다.
- 실용적 가치:
- 단기적으로는 실용적 가치가 제한적이다(추가 알고리즘 개발 필요).
- 장기적으로는 대규모 시스템, 분산 계산에 새로운 아이디어를 제공할 수 있다.
- 양의 시스템 최적 제어에 직접 응용된다(섹션 4.3).
- 재현 가능성:
- 이론 증명이 완전하여 검증하기 쉽다.
- 알고리즘 설명이 명확하지만 구현 세부사항이 부족하다.
- 저자가 코드를 공개하여 연구를 촉진할 것을 권장한다.
- 이론 연구:
- 대수 방정식의 동역학 시스템 성질 연구
- 양의 시스템 이론의 새로운 응용 탐구
- 제어 시스템의 기하학적 구조 분석
- 특정 응용:
- 양의 시스템의 최적 제어(예: 전염병 모델, 교통 네트워크)
- 대규모 희소 시스템(양의 시스템 구조가 계산 이점을 가져올 수 있음)
- 분산 최적화(Hilbert 메트릭이 일치성 분석에 적합)
- 교육 목적:
- 양의 시스템 이론의 강력함을 보여준다.
- 대수 방정식과 동역학 시스템의 연결을 설명한다.
- 고급 제어 이론 과정의 보충 자료로 사용된다.
- 16 Lemmens & Nussbaum (2012): Nonlinear Perron-Frobenius Theory - 본 논문의 주요 이론 도구 출처
- 19 Pates (2024): Lyapunov equations: a (fixed) point of view - 직접적인 비교 연구
- 21 Rantzer & Valcher (2018): A tutorial on positive systems and large scale control - 양의 시스템 이론 종합
- 20 Rantzer (2022): Explicit solution to bellman equation for positive systems - 양의 시스템 최적 제어
- 8 Bushell (1973): Hilbert's metric and positive contraction mappings - Hilbert 메트릭의 고전 참고 자료
종합 평가: 이것은 이론적 혁신성이 매우 강한 연구로, 양의 시스템 이론을 고전 제어 방정식 분석에 성공적으로 적용하여 완전히 새로운 이론적 관점과 통일된 프레임워크를 제공한다. 동차화 기술은 방법론적 가치를 가지며 다른 분야의 연구에 영감을 줄 수 있다. 주요 부족점은 실용성 검증이 필요하며, 알고리즘 구현 및 수치 실험 측면에서 보완이 필요하다. 이론 연구자와 양의 시스템 이론에 관심 있는 학자들에게 깊이 있게 읽을 가치가 있는 논문이다.