On Koopman Resolvents and Frequency Response of Nonlinear Systems

이 논문은 쿠퍼만 연산자 프레임워크를 활용하여 비선형 시스템의 주파수 응답을 유도하고, 이를 통해 고전적인 LTI 시스템의 접근법을 일반화한 보드 선도 작성 및 존재 조건을 제시합니다.

핵심

1. 문제: 왜 비선형 시스템은 분석하기 어려울까요?

**비선형 시스템 (Nonlinear System)**은 우리가 일상에서 마주치는 복잡한 현상들입니다. 예를 들어, 자동차의 핸들을 살짝 돌렸을 때 차가 부드럽게 움직이다가, 급하게 돌리면 갑자기 미끄러지거나 회전하는 것처럼, 원인과 결과가 비례하지 않는 경우를 말합니다.

• 기존의 방법 (선형 시스템): 과거의 공학자들은 시스템을 단순화해서 분석했습니다. "소리를 키우면 (입력) 소리도 비례해서 커진다 (출력)"는 식이죠. 이때는 **주파수 응답 (Frequency Response)**이라는 도구를 쓰면, "어떤 소리를 넣었을 때 시스템이 어떻게 반응하는지"를 그래프 (보드 선도) 로 쉽게 그릴 수 있었습니다.
• 문제점: 하지만 복잡한 비선형 시스템에서는 소리를 키운다고 해서 소리가 일정하게 커지지 않습니다. 새로운 주파수 (고조파) 가 튀어나오거나, 아예 다른 소리가 들리기도 합니다. 기존 도구로는 이 복잡한 춤을 따라가기 힘들었습니다.

2. 해결책: 거울을 통해 선형으로 바꾸다 (코프만 연산자)

이 논문은 **"복잡한 춤을 추는 비선형 시스템을, 거울에 비추면 단순한 선형 춤으로 바꿀 수 있다"**는 아이디어를 제시합니다.

• 비유:
  • 비선형 시스템: 혼란스러운 오케스트라가 제각기 다른 리듬으로 연주하는 상황입니다.
  • 코프만 연산자 (Koopman Operator): 이 혼란스러운 연주를 거대한 거울에 비추는 장치입니다.
  • 코프만 연산자의 마법: 이 거울을 통해 보면, 혼란스러운 연주자들이 사실은 **하나의 규칙적인 선율 (선형 시스템)**로 움직이고 있는 것처럼 보입니다. 비록 실제 세계에서는 복잡하지만, 이 '거울 세계'에서는 수학적으로 매우 단순하고 깔끔하게 다룰 수 있게 됩니다.

3. 이 논문의 핵심 기여: 새로운 주파수 응답 정의

이 논문은 그 '거울 세계'의 규칙을 이용해, 비선형 시스템의 주파수 응답을 새롭게 정의했습니다.

1. 입력을 주기적으로 바꾸기: 시스템에 일정한 리듬 (예: $u_0 e^{i\omega t}$) 으로 입력을 줍니다.
2. 거울 세계로 이동: 이 입력을 받은 시스템을 코프만 연산자의 거울 세계로 옮깁니다.
3. 주파수 응답 찾기: 거울 세계에서 시스템이 어떤 주파수로 반응하는지 분석합니다.
   • 기본 주파수: 입력된 리듬 그대로 반응하는 경우.
   • 고조파 (Harmonic): 입력 리듬의 2 배, 3 배 속도로 반응하는 경우 (비선형 시스템의 특징).
   • 서브하모닉 (Subharmonic): 입력 리듬의 절반, 1/3 속도로 반응하는 경우.

이론적으로 이 응답값은 복소수 함수로 표현되며, 이를 통해 기존 선형 시스템처럼 **보드 선도 (Bode Plot)**라는 그래프를 그릴 수 있게 됩니다. 이 그래프는 시스템이 특정 주파수에서 소리를 얼마나 키울지 (이득), 얼마나 지연시킬지 (위상) 를 보여줍니다.

4. 구체적인 예시: 2 차원 비선형 시스템

논문의 예시를 들어보면:

• 시스템 A ($x_2$): 입력 신호 $u$에 대해 선형처럼 반응합니다. 입력 주파수 $\omega$만 출력됩니다.
• 시스템 B ($x_1$): 시스템 A 의 출력 ($x_2$) 이 제곱 ($x_2^2$) 되는 비선형 항이 있습니다.
  • 여기서 신기한 일이 일어납니다. 입력이 $\omega$인데, 출력에는 $2\omega$ (2 배 주파수) 성분이 강하게 나타납니다.
  • 이 논문의 방법론을 사용하면, 이 **$2\omega$성분이 어디서 왔는지, 그 크기와 위상이 어떻게 되는지**를 정확히 계산해 낼 수 있습니다. 마치 오케스트라에서 베이스 기타 소리 ($\omega$) 가 들리는데, 갑자기 트럼펫 소리 ($2\omega$) 가 섞여 나오는 이유를 분석하는 것과 같습니다.

5. 언제 이 방법이 쓸모있을까요? (조건)

이 마법의 거울이 제대로 작동하려면 몇 가지 조건이 필요합니다.

1. 선형 시스템: 당연히 잘 작동합니다. 기존 방법과 결과가 일치합니다.
2. 전역적으로 안정적인 시스템: 시스템이 시간이 지나면 안정된 리듬 (주기적 해) 을 찾는 경우입니다.
3. 에르고드 (Ergodic) 시스템: 시스템이 특정 영역을 골고루 방문하는 복잡한 운동 (예: 카오스) 을 할 때도 적용 가능합니다.

6. 요약 및 결론

이 논문은 **"비선형 시스템도 거울 (코프만 연산자) 을 통해 선형 시스템처럼 분석할 수 있다"**는 것을 증명했습니다.

• 기존: 비선형 시스템은 분석하기 너무 복잡해서 주파수 응답을 그리기 힘들었다.
• 이제: 코프만 연산자를 이용해 시스템을 '선형화'된 거울 세계로 옮기면, 보드 선도를 그릴 수 있고, 고조파나 서브하모닉 같은 복잡한 반응도 수학적으로 정확히 계산할 수 있게 되었습니다.

이는 공학자들이 복잡한 비선형 시스템 (로봇, 전력망, 생체 시스템 등) 을 설계하고 제어할 때, 마치 선형 시스템을 다룰 때처럼 직관적이고 체계적인 도구를 쓸 수 있게 해주는 혁신적인 방법론입니다.

기술 요약

1. 연구 배경 및 문제 정의 (Problem)

• 배경: 주파수 응답 (Frequency Response) 은 선형 시불변 (LTI) 시스템의 분석 및 설계에 있어 가장 기본적이고 핵심적인 개념입니다. 그러나 비선형 시스템의 경우, 비선형성으로 인해 기존의 주파수 영역 해석이 어렵고, 기존 연구들은 주로 시간 영역 (Time-domain) 이나 상태 공간 (State-space) 기반의 현상론적 접근 (조화 평형법, 기술 함수, 볼테라 급수 등) 에 의존해 왔습니다.
• 문제점: 비선형 시스템에 대해 LTI 시스템과 유사한 엄밀한 이론에 기반한 주파수 응답을 정의하고, 이를 통해 시스템의 이득 (Gain) 과 위상 (Phase) 특성을 체계적으로 분석할 수 있는 방법이 부재했습니다.
• 목표: Koopman 연산자 (Koopman Operator) 프레임워크를 활용하여 비선형 시스템의 주파수 응답을 새로운 방식으로 공식화하고, 이를 통해 주파수 영역에서의 분석이 가능하도록 하는 것입니다.

2. 방법론 (Methodology)

이 논문은 Koopman 연산자 이론과 Resolvent (해결자) 이론을 결합하여 비선형 시스템의 주파수 응답을 유도합니다. 주요 방법론적 단계는 다음과 같습니다.

1. Skew-product Form (비선형 강제 시스템의 자동화):

   • 주기적인 입력 $u(t) = u_0 e^{i\omega t}$을 받는 비선형 시스템을 분석하기 위해, 입력을 상태 변수로 확장한 Skew-product 형태의 자동화 시스템 (Autonomous system) 을 구성합니다.
   • 이를 통해 비자율 시스템 (Non-autonomous) 을 동등한 자율 시스템으로 변환하여 Koopman 연산자 이론을 적용할 수 있게 합니다.

2. Koopman Resolvent 의 도입:

   • 변환된 시스템의 Koopman 생성자 (Generator) $L_{forced}$에 대해 Koopman Resolvent $R(s; L_{forced}) = (sI - L_{forced})^{-1}$를 정의합니다.
   • 출력 $y(t)$의 라플라스 변환 $\hat{y}(s)$가 Koopman Resolvent 의 작용으로 표현됨을 이용합니다 (선행 연구 [7] 기반).

3. 주파수 응답의 정의:

   • 기본 및 고조파 (Fundamental & Harmonic): 입력 주파수 $\omega$의 정수배 $n\omega$ 성분을 가진 정상 상태 주기 출력에 대한 응답 $H_n(\omega; g)$을 정의합니다.
   • 서브하모닉 (Subharmonic): 입력 주파수의 분수배 $\omega/n$ 성분을 가진 출력에 대한 응답 $H_{1/n}(\omega; g)$을 정의합니다.
   • 이 응답들은 **Koopman 모드 (Koopman Mode)**와 직접적으로 연결됩니다.

4. 스펙트럼 분석을 통한 유도:

   • 주파수 응답은 Koopman Resolvent 의 **극점 (Pole)**에서의 잔차 (Residue) 로 계산됩니다.
   • Specifically, $s = in\omega$ 또는 $s = i\omega/n$에서 Resolvent 가 1 차 극점을 가질 때, 해당 극점에서의 잔차를 통해 주파수 응답 함수를 구합니다.

3. 주요 기여 (Key Contributions)

1. 비선형 시스템 주파수 응답의 새로운 공식화:
   • Koopman Resolvent 이론에 기반한 비선형 시스템의 주파수 응답을 체계적으로 정의했습니다. 이는 LTI 시스템의 고전적 접근법을 비선형 영역으로 자연스럽게 확장한 것입니다.
2. 복소수 함수로서의 응답 및 보드 플롯 (Bode Plots) 작성 가능:
   • 유도된 주파수 응답은 구동 각주파수 $\omega$의 복소수 함수로, 이를 통해 **보드 플롯 (이득 및 위상 특성)**을 그릴 수 있음을 보였습니다. 이는 비선형 시스템의 주파수 특성을 직관적으로 시각화할 수 있는 길을 열었습니다.
3. 존재 조건 제시:
   • 주파수 응답이 존재하기 위한 충분 조건을 세 가지 동역학 클래스 (LTI 동역학, 전역적으로 안정한 동역학, 에르고드 동역학) 에 대해 제시했습니다.

4. 주요 결과 및 사례 (Results & Examples)

• 정리 1 및 2 (Theorems 1 & 2):
  • Koopman Resolvent 가 특정 주파수 ($in\omega$ 또는 $i\omega/n$) 에서 1 차 극점을 가진다면, 해당 주파수에서의 주파수 응답은 Resolvent 의 잔차 (Residue) 로 정확히 계산됨을 증명했습니다.
  • 이는 주파수 응답이 본질적으로 Koopman 모드임을 의미합니다.
• 수치 예시:
  • 1 차원 선형 시스템: 유도된 결과가 전통적인 선형 제어 이론의 주파수 응답 ($b/(i\omega - a)$) 과 일치함을 확인했습니다.
  • 2 차원 비선형 시스템: $x_2^2$ 항을 포함한 비선형 시스템을 분석했습니다.
    • 입력 $u$에 대해 출력 $x_2$는 1 차 지연 요소처럼 동작 ($H_1$만 존재).
    • 출력 $x_1$은 비선형성으로 인해 2 배 주파수 ($2\omega$) 성분이 발생하며, 이는 3 차 지연 요소처럼 동작 ($H_2$만 존재).
    • 이를 통해 보드 플롯을 작성하여 주파수별 이득과 위상 변화를 시각화할 수 있음을 시연했습니다.
• 선형 리프트 (Linear Lift) 와의 일치:
  • 비선형 시스템을 고차 선형 시스템으로 '리프트' (Lift) 하여 분석했을 때, Koopman Resolvent 를 통해 구한 주파수 응답과 동일한 결과를 얻음을 확인했습니다.

5. 의의 및 향후 전망 (Significance & Future Work)

• 이론적 의의:
  • 비선형 시스템의 주파수 응답을 **선형 연산자 (Koopman)**의 스펙트럼 성질로 해석할 수 있는 강력한 이론적 틀을 제공했습니다.
  • 기존에 시간 영역이나 상태 공간에서만 가능했던 비선형 시스템 분석을 주파수 영역으로 확장하여, 제어 공학의 직관적인 도구 (보드 플롯, 주파수 응답 함수) 를 비선형 시스템에 적용할 수 있게 했습니다.
• 실용적 가치:
  • **동역적 모드 분해 (Dynamic Mode Decomposition, DMD)**와 같은 수치 기법을 활용하여 실험 데이터나 시뮬레이션 데이터로부터 Koopman 모드를 추정하고, 이를 통해 비선형 시스템의 주파수 응답을 수치적으로 추정할 수 있는 길이 열렸습니다.
• 향후 연구 방향:
  • 주파수 응답 함수의 잘-정의됨 (Well-posedness), 연속성, 미분 가능성 분석.
  • 안정성 (Stability) 및 패시비티 (Passivity) 분석에 주파수 응답을 활용.
  • 비선형 시스템에 대한 입력 - 출력 관계 (Transfer Function) 의 정의 확장.
  • 연속 스펙트럼 (Continuous spectrum) 을 포함하는 경우 (예: 약한 혼합 시스템) 로의 이론 확장.

요약

이 논문은 Koopman 연산자 프레임워크를 활용하여 비선형 시스템의 주파수 응답을 엄밀하게 정의하고, 이를 **Koopman Resolvent 의 극점 (Pole)**과 Koopman 모드를 통해 계산하는 방법을 제시했습니다. 이를 통해 비선형 시스템에서도 LTI 시스템과 유사한 보드 플롯 작성이 가능해졌으며, 이는 비선형 시스템의 분석 및 합성을 위한 새로운 주파수 영역 기반의 강력한 도구가 될 것으로 기대됩니다.