데이터 기반 ν‑갭 메트릭 추정 및 위상 검증 방법

본 논문은 제어 시스템 설계 시 모델 불일치를 정량적으로 평가하는 ν‑갭(ν‑gap) 메트릭을, 시스템 식별 없이 입력‑출력 데이터만으로 추정하는 새로운 프레임워크를 제시한다. 기존 ν‑갭은 두 시스템의 전달함수를 완전히 알아야 계산 가능했으며, 특히 인덱스 조건 C (함수 f₁, f₂의 외부 영점 개수 일치)을 만족해야 의미 있는 거리값을 제공한다. 이러한 제한을 극복하고자 저자는 두 가지 주요 기술을 결합한다. 첫 번째는 ν‑갭 자체를 파워 이터레이션 형태로 근사하는 방법이다. 입력 신호 u(k) 의 DFT U(e^{jω})와 실험으로 얻은 출력 스펙트럼 Y(e^{jω}) (실제 시스템 G)와 Y₀(e^{jω}) (명목 모델 G₀)를 이용해 ν‑갭을 주파수‑도함수 상한으로 표현한다. 이 식은 H∞ 노름과 동일하게 “최대 증폭”을 찾는 문제이므로, 파워 이터레이션과 유사하게 입력 스펙트럼을 비율 |Y−Y₀| / (|U|²+|Y|²)(|U|²+|Y₀|²) 에 비례하도록 재설계한다. 구체적인 업데이트 식(6)은 이전 입력 Uₙ 에 위 비율을 곱하고 정규화해 새로운 입력 Ũₙ₊₁을 만든다. 반복이 진행될수록 입력 에너지가 ν‑갭이 가장 크게 나타나는 주파수에 집중되며, 최종 비율값이 ν‑갭의 근사값이 된다. 두 번째는 인덱스 조건 C 을 데이터 기반으로 검증하는 winding‑number 계산이다. 기존에는 Nyquist 플롯을 눈으로 확인했지만, 저자는 Welch‑형 평균을 차용해 f₁, f₂를 추정한다. (7)식은 각 반복에서 얻은 Y, Y₀, U 의 교차 스펙트럼을 이용해 f₂(e^{jω})≈1+∑Y·U* / (max|U|²+ε₀) 와 같은 형태로 계산한다. f₁은 동일 방식으로 구한다. 이후 (8)식에서 연속적인 위상 차이를 누적해 winding‑number wno(f_i) 를 구한다. 두 함수의 winding‑number가 일치하면 C 조건이 만족된 것으로 판단하고, 불일치 시 알고리즘을 조기에 종료한다. 알고리즘 1은 위 두 과정을 통합해 ν‑갭 추정과 조건 검증을 동시에 수행한다. 알고리즘 2는 winding‑number 계산을 상세히 기술한다. 실험에서는 두 개의 고압 가스 터빈(HDGT) 모델과 교과서 수준의 합성 예제를 사용했다. HDGT 사례에서는 실제 시스템을 시뮬레이션으로 생성하고, 명목 모델 G₀와 비교해 ν‑갭을 추정했다. 결과는 MATLAB gapmetric 함수와 평균 오차 <0.02로 거의 일치했으며, 조건 위반 시(예: 위상 반전이 있는 경우) winding‑number 검증이 정확히 1≠0을 반환해 실패를 알렸다. 교과서 예제에서는 f₁, f₂의 영점 구조를 수식적으로 분석해 조건이 성립하지 않음을 보였고, 데이터 기반 검증 역시 동일한 결과를 도출했다. 논문의 주요 기여는 다음과 같다. ① 식별‑프리 방식으로 ν‑갭을 추정하는 파워 이터레이션 기반 알고리즘을 제시, 입력 설계와 주파수 스펙트럼 재활용을 통해 실시간 추정이 가능하도록 함. ② Welch‑형 평균을 이용한 데이터 기반 winding‑number 계산으로 인덱스 조건을 자동 검증, 기존의 복잡한 Nyquist 분석을 대체. ③ 시뮬레이션을 통해 MATLAB gapmetric과의 정량적 일치를 입증하고, 산업 현장(HDGT) 적용 가능성을 시연. 향후 연구에서는 다변량(MIMO) 시스템 확장, 잡음·비선형성에 대한 강인성 분석, 실시간 제어 루프 내에서의 온라인 ν‑갭 모니터링 등을 탐색할 여지가 있다.