피드백 제어를 통한 근접 이웃 및 완전 그래프에서의 꼬임 상태 안정화

초록

본 논문은 동일한 진동자를 갖는 Kuramoto 모델에 선형·비선형 피드백을 적용하여, 근접 이웃 그래프와 완전 단순 그래프에서 존재하는 q‑twisted 상태의 존재와 안정성을 분석한다. 연속극한(CLT)에서 중심다양체 축소법을 이용해 선형 피드백 이득 b₁을 제어 매개변수로 삼아, b₁이 특정 임계값을 통과하면 꼬임 상태가 불안정해지고, 위상 지연 σ에 따라 변조형 혹은 진동형 새로운 꼬임 해가 생성되는 분기 현상을 보인다. 수치 실험을 통해 이론적 결과를 검증한다.

상세 분석

이 연구는 Kuramoto 모델(KM)의 “꼬임 상태”(q‑twisted state)를 피드백 제어를 통해 조절·안정화할 수 있음을 수학적으로 증명한다. 먼저, 동일 진동수 ω와 위상 지연 σ를 갖는 KM을 (1.1)식으로 정의하고, 선형(b₁)·비선형(b₃) 피드백 항을 추가한다. 그래프는 n개의 정점으로 이루어진 근접 이웃 그래프(κ≤½)와 κ=½일 때의 완전 단순 그래프를 포함한다. 그래프 가중치 행렬 Wₙ은 graphon W(x,y)으로 수렴하도록 가정하고, n→∞ 극한에서 연속극한(CLT) 방정식(1.6)을 도출한다.

핵심은 목표 궤도 ˆu(t,x)=2πqx+Ωt를 설정하고, 오차 v(t,x)=u(t,x)−ˆu(t,x)로 변환함으로써 비선형 항이 v³ 형태로 나타나는 (1.5) 형태의 동역학을 얻는 것이다. 여기서 q∈ℕ은 꼬임 차수를 의미한다.

연속극한에서의 선형화는 Fourier 모드 e^{2πiℓx}에 대해 고유값 λ_ℓ(b₁)=−b₁+Λ_ℓ(κ,σ,q) 로 표현된다. Λ_ℓ는 그래프온(graphon)과 위상 지연, 꼬임 차수에 의해 결정되는 실수 함수이며, ℓ=±q 모드가 가장 큰 실수부를 가진다. b₁을 감소시키면 λ_{±q}가 0을 통과하는 임계값 b₁^{c}가 존재하고, 이는 피드백 제어에 의한 피치포크(pitchfork) 혹은 호프(또는 서브크리티컬) 분기를 유발한다.

σ=0(위상 지연 없음)인 경우, 임계점에서 실수 고유값이 0이 되면서 대칭이 깨지는 피치포크 분기가 일어나며, 새로운 변조형 꼬임 해(진폭이 일정하고 위상이 공간적으로 변조된 형태)가 안정적으로 나타난다. 반면 σ≠0이면 복소 고유값 쌍이 실축을 가로질러 실부가 양수가 되면서 호프(Hopf) 분기가 발생하고, 진동 주파수를 갖는 시간 진동형 꼬임 해가 탄생한다.

비선형 피드백 b₃>0은 고차 항 v³을 제공해 포화 효과를 주어, 분기 후 생성되는 새로운 해가 안정적인 영역을 넓힌다. 특히 b₃가 충분히 크면, b₁이 임계값 이하일 때 변조·진동 해가 넓은 파라미터 구간에서 전역적으로 안정한다는 점이 강조된다.

수치 실험에서는 n=200 정도의 근접 이웃 그래프와 완전 그래프에 대해 (1.1)을 직접 시뮬레이션하고, b₁을 서서히 변화시켜 꼬임 상태의 전이 과정을 관찰한다. 결과는 이론적 예측과 일치하며, 특히 b₁이 임계값을 지나면 원래의 q‑twisted 상태가 급격히 붕괴하고, 변조형(σ=0) 혹은 진동형(σ≠0) 새로운 꼬임 패턴이 나타난다.

이 논문은 기존 연구에서 다루어지던 “동기화 제어”와는 달리, 비동기적인 꼬임 상태 자체를 목표로 삼아 피드백을 설계한다는 점에서 독창적이다. 또한, 연속극한을 이용해 무한대 네트워크의 거동을 정확히 분석하고, 중심다양체 축소법을 통해 저차원 정상형을 도출함으로써 복잡한 고차원 비선형 네트워크의 분기 구조를 명확히 파악한다. 이러한 방법론은 다른 유형의 비선형 네트워크(예: 비국소 결합, 이질적 진동수)에도 확장 가능성을 시사한다.