하이브리드 시스템을 위한 빠른 불확실성 정량화: 폴리노미얼 카오스와 웨이브렛·전송 이론 통합

초록

본 논문은 연속·이산 동작을 동시에 갖는 하이브리드 동역학 시스템에 대해, 전통적인 몬테카를로 방식보다 빠른 불확실성 정량화 기법을 제시한다. 폴리노미얼 카오스(gPC)를 기반으로 하되, 상태 전이와 리셋으로 인한 불연속성을 처리하기 위해 웨이브렛 기반 Wiener‑Haar 전개와 경계층 모델을 도입한다. 또한 확률 전송 방정식을 해석적 특성선 적분으로 풀어 하이퍼볼릭 PDE 형태로 불확실성을 전파하는 새로운 전송 이론 접근법을 제안한다. 스위칭 진동기와 점핑볼 예제로 성능을 검증한다.

상세 분석

이 연구는 하이브리드 동역학 시스템에서 파라미터 불확실성을 효율적으로 전파하기 위한 두 가지 핵심 기법을 결합한다. 첫 번째는 기존의 일반화 폴리노미얼 카오스(gPC) 프레임워크를 확장하여, 시스템이 여러 모드(연속·이산)로 전이될 때 각 모드에 대한 지시함수 1_R_i(x)를 도입함으로써 상태공간을 시간‑의존적 적분 영역으로 분할한다. 이때 계수 a_k(t)의 미분 방정식은 두 영역에 대한 기대값을 각각 적분해 얻으며, 적분 영역이 실시간으로 변하는 점이 계산 복잡도를 높인다.

두 번째는 이러한 불연속성을 보다 정밀히 포착하기 위해 Wiener‑Haar 웨이브렛 전개를 적용한다. 파라미터 λ의 누적분포함수(u(λ))에 대한 Haar 웨이브렛 ψ_{j,k}(u) 를 사용하면, 전개 계수는 구간별 상수값을 갖게 되어 적분이 구간별 사전계산된 λ 평균값으로 대체될 수 있다. 결과적으로 스위칭 경계 근처에서 급격히 변하는 확률밀도도 저차 웨이브렛으로 정확히 재현한다.

리셋 현상(예: 충돌 후 속도 반전)은 경계층(ε) 모델을 통해 연속적인 동역학으로 근사한다. 가드 조건 |g(y)|≥ε 영역에서는 원래 시스템을 그대로 사용하고, |g(y)|<ε 영역에서는 리셋 함수를 선형 보간식으로 대체한다. ε→0 한계에서 원 시스템과 동일한 동작을 보장하면서도, 미분 방정식 형태를 유지하므로 gPC 전개가 그대로 적용 가능하다.

마지막으로 전송 이론 기반 접근법을 제시한다. 확률 밀도 ρ(x,λ,t)를 상태‑파라미터 공간의 보존식으로 기술하고, 이를 Haar 전개에 삽입하면 계수들이 하이퍼볼릭 PDE 형태를 만든다. 특성선 적분을 이용해 각 계수를 시간에 따라 전파하면, 고차원 적분 없이도 정확한 순간통계(평균, 분산 등)를 얻을 수 있다. 이 방법은 특히 다중 모드가 복잡하게 얽힌 Zeno 현상이나 비선형 경계에서 유리하다.

실험에서는 스위칭 진동기와 점핑볼(바운싱 볼) 시스템을 대상으로, Monte‑Carlo와 Quasi‑Monte‑Carlo 결과와 비교했을 때, 제안 기법이 동일한 정확도에서 10‑100배 빠른 수렴 속도를 보였다. 특히 웨이브렛 전개는 높은 스위칭 비율에서도 계수 진동을 억제하고, 경계층 모델은 리셋에 따른 급격한 상태 변화를 부드럽게 처리한다. 전체적으로, 이 논문은 하이브리드 시스템의 불연속성을 정량화하는 새로운 수학적 도구와 구현 전략을 제공한다.