[논문 리뷰] On the stability and the uniform propagation of chaos properties of Ensemble Kalman-Bucy filters

Pierre Del Moral, Julian Tugaut|arXiv (Cornell University)|2016. 05. 30.

Meteorological Phenomena and Simulations참고 문헌 65인용 수 63

한 줄 요약

이 논문은 새로운 기능 부등식을 도입하고 포스터-라이아푸노프 기법, 커플링, 스펙트럼 분석을 결합함으로써, 앙상블 칼만-부시 필터에 대한 최초의 균일 안정성 및 혼합의 전파 결과를 확립한다. 이는 임의의 시간 구간 동안 균일한 L²-평균 오차 추정치를 증명하고, 고차원적이고 악조건인 시스템에서 필터의 수렴성과 안정성에 필요한 필수 및 충분 조건을 규명한다.

ABSTRACT

The Ensemble Kalman filter is a sophisticated and powerful data assimilation method for filtering high dimensional problems arising in fluid mechanics and geophysical sciences. This Monte Carlo method can be interpreted as a mean-field McKean-Vlasov type particle interpretation of the Kalman-Bucy diffusions. In contrast to more conventional particle filters and nonlinear Markov processes these models are designed in terms of a diffusion process with a diffusion matrix that depends on particle covariance matrices. Besides some recent advances on the stability of nonlinear Langevin type diffusions with drift interactions, the long-time behaviour of models with interacting diffusion matrices and conditional distribution interaction functions has never been discussed in the literature. One of the main contributions of the article is to initiate the study of this new class of models The article presents a series of new functional inequalities to quantify the stability of these nonlinear diffusion processes. In the same vein, despite some recent contributions on the convergence of the Ensemble Kalman filter when the number of sample tends to infinity very little is known on stability and the long-time behaviour of these mean-field interacting type particle filters. The second contribution of this article is to provide uniform propagation of chaos properties as well as Lp-mean error estimates w.r.t. to the time horizon. Our regularity condition is also shown to be sufficient and necessary for the uniform convergence of the Ensemble Kalman filter. The stochastic analysis developed in this article is based on an original combination of functional inequalities and Foster-Lyapunov techniques with coupling, martingale techniques, random matrices and spectral analysis theory.

연구 동기 및 목표

• 고차원적이고 악조건인 시스템에서 앙상블 칼만-부시 필터의 장기적 안정성과 수렴성에 대한 이론적 이해 부족을 해결하기 위해.
• 임의의 시간 구간 동안 표본 평균과 공분산 행렬에 대한 균일한 L²-평균 오차 추정치를 수립하기 위해.
• 공분산에 의존하는 공분산 행렬을 가진 비선형 확산 과정의 안정성을 분석하기 위해, 이전에 깊이 있게 연구되지 않은 클래스에 대해.
• 입자 수가 무한으로 갈 때 앙상블 칼만-필터(EnKF)의 균일 수렴성에 필요한 필수 및 충분한 정칙 조건을 유도하기 위해.
• 제닝스 유형의 부등식과 스펙트럼 분석을 통해 칼만-부시 필터 및 관련 리카티 방정식의 안정성을 정량화하기 위해.

제안 방법

• 공분산에 의존하는 공분산 행렬을 가진 비선형 매킨-블라소프 유형 확산 과정의 안정성을 정량화하기 위한 새로운 기능 부등식의 클래스를 개발한다.
• 포스터-라이아푸노프 기법과 커플링 방법을 적용하여 워샤르스탄 거리와 상대 엔트로피에서의 지수 안정성을 분석한다.
• 마링게일 기법, 랜덤 매트릭스 이론, 스펙트럼 분석을 조합하여 상호작용 입자 시스템의 장기적 행동을 연구한다.
• 칼만-부시 확산의 평균장 입자 해석을 사용하여 EnKF를 비선형 상호작용 입자 시스템으로 모델링한다.
• 표본 통계량(평균과 공분산)이 시간에 따라 진짜 값으로 수렴하는 것을 분석함으로써 균일한 L²-평균 오차 추정치를 도출한다.
• 시스템이 안정성을 유지할 수 있는 허용 가능한 변동 행렬 Q의 집합을 규명하고, 불안정성이 발생하는 발산 집합을 특정한다.

실험 결과

연구 질문

• RQ1 앙상블 칼만-부시 필터가 임의의 시간 구간 동안 균일 안정성을 보일 조건은 무엇인가?
• RQ2 EnKF 입자 시스템에서 혼합의 전파 성질이 시간에 대해 균일하게 어떻게 정량화될 수 있는가?
• RQ3 공분산에 의존하는 공분산 행렬을 가진 비선형 확산 과정의 안정성을 지배하는 기능 부등식은 무엇인가?
• RQ4 시스템의 정칙성과 입자 수가 무한으로 갈 때 EnKF의 균일 수렴성 사이의 정확한 관계는 무엇인가?
• RQ5 변동 행렬 Q에 대해 어떤 조건에서 칼만-부시 관측기에서 지수 발산이 발생하는가?

주요 결과

• 논문은 임의의 시간 구간 동안 앙상블 칼만-부시 필터의 표본 평균과 표본 공분산 행렬에 대해 최초로 균일한 L²-평균 오차 추정치를 확립한다.
• 입자 수가 무한으로 갈 때 EnKF의 균일 수렴성에 필요한 필수 및 충분한 정칙 조건을 규명한다.
• 변동 행렬 Q의 발산 집합은 두 영역의 합집합으로 명시적으로 특성화된다: 하나는 Q₁,₂ > 1.5Q₁,₁ − 1.8로 정의되고, 다른 하나는 Q₁,₁ ∈ ]−8.7, −4.7[로 정의된다.
• A = [[1,2],[1,3]]인 예제 시스템의 경우, 안정적인 Q의 부분집합은 Q₁,₂ > 1.5Q₁,₁ − 1.8 및 Q₁,₁ ∈ ]−4.7, ∞[로 정의된다.
• 분석 결과, 관측기 행렬의 고유값이 서로 다른 부호를 가지거나, 복소 고유값의 실수부가 양수가 되면 불안정성이 발생함을 밝혀냈다.
• 허용 가능한 변동 행렬 Q의 집합은 P와 Q를 포함하는 행렬식 및 양의 제약 조건의 조합으로 정의된다.