[논문 리뷰] Balanced truncation and singular perturbation approximation model order reduction for stochastically controlled linear systems

Martin Redmann, Melina A. Freitag|arXiv (Cornell University)|2016. 11. 16.

Model Reduction and Neural Networks참고 문헌 24인용 수 4

한 줄 요약

이 논문은 추가적인 Lévy 노이즈에 의해 구동되는 선형 스토크래틱 시스템으로 균형 절단(BT)과 단순화된 편미분 근사(SPA) 모델 차원 축소를 확장하며, 일반화된 라플라스 방정식을 통한 스토크래틱 접근 가능성과 관측 가능성 그람형을 도입한다. 두 방법에 대한 엄밀한 오차 한계를 유도하고, Lévy 노이즈가 있는 감쇠된 파동 방정식에서 수치적으로 검증하여 오차 추정치가 타이트하고 차원 축소가 효과적으로 이루어짐을 보여준다.

ABSTRACT

When solving linear stochastic differential equations numerically, usually a high order spatial discretisation is used. Balanced truncation (BT) and singular perturbation approximation (SPA) are well-known projection techniques in the deterministic framework which reduce the order of a control system and hence reduce computational complexity. This work considers both methods when the control is replaced by a noise term. We provide theoretical tools such as stochastic concepts for reachability and observability, which are necessary for balancing related model order reduction of linear stochastic differential equations with additive Lévy noise. Moreover, we derive error bounds for both BT and SPA and provide numerical results for a specific example which support the theory.

연구 동기 및 목표

• 제어 입력이 Lévy 노이즈로 대체되는 선형 스토크래틱 시스템에서 모델 차원 축소(MOR)를 위한 이론적 프레임워크를 개발하는 것.
• 추가적인 Lévy 노이즈를 갖는 시스템에서 균형 기반 MOR에 적합한 스토크래틱 접근 가능성과 관측 가능성의 개념을 도입하는 것.
• Lévy 노이즈를 갖는 스토크래틱 설정에서 균형 절단(BT)과 단순화된 편미분 근사(SPA)에 대한 사전 오차 한계를 유도하는 것.
• Lévy 노이즈가 있는 2차 감쇠된 파동 방정식에 이론을 적용하고 오차 한계를 수치적으로 검증하는 것.
• BT와 SPA가 실제로 Lévy 노이즈를 갖는 스토크래틱 시스템에서 작은 오차와 타이트한 오차 한계를 제공함을 보여주는 것.

제안 방법

• Lévy 노이즈의 특성, 특히 Lévy 과정의 2차 모멘트를 포함한 일반화된 라플라스 방정식을 통해 스토크래틱 접근 가능성과 관측 가능성 그람형을 정의한다.
• 스토크래틱 그람형의 곱을 사용하여 균형 변환을 구성함으로써 투영 기반 모델 차원 축소를 가능하게 한다.
• 무시된 한켈 특이값에 기반한 감소된 차원 시스템의 일반 오차 한계를 유도하며, 결정론적 오차 한계를 스토크래틱 케이스로 확장한다.
• 일반 오차 한계를 BT와 SPA에 특화하여, 오차가 무시된 특이값의 합에 의존함을 보여준다.
• 스펙트럼 갈레르킨 방법을 통해 공간 이산화하고, 이를 통해 유도된 고차원 시스템에 대해 MOR을 적용함으로써 2차 감쇠된 파동 방정식에 방법을 적용한다.
• 이토 미적분과 셈그룹 이론을 사용하여 감소된 시스템의 수렴성을 분석하고, 전체 시스템과 감소된 시스템 해 사이의 평균 제곱 오차에 대한 경계를 유도한다.

실험 결과

연구 질문

• RQ1 균형 절단과 단순화된 편미분 근사는 추가적인 Lévy 노이즈에 의해 구동되는 선형 스토크래틱 시스템으로 엄밀하게 확장될 수 있는가?
• RQ2 Lévy 노이즈를 갖는 시스템에서 균형을 정의하기 위해 필요한 접근 가능성과 관측 가능성 그람형의 스토크래틱 유사체는 무엇인가?
• RQ3 Lévy 노이즈를 갖는 스토크래틱 시스템의 맥락에서 BT와 SPA에 대한 사전 오차 한계는 무엇인가?
• RQ4 실제 응용에서 오차 한계는 어떻게 행동하며, 예를 들어 스토크래틱 PDE의 경우 오차 한계는 타이트한가?
• RQ5 제안된 MOR 기법은 Lévy 노이즈가 있는 2차 감쇠된 파동 방정식에 성공적으로 적용될 수 있으며, 정확성을 유지하는가?

주요 결과

• 논문은 일반화된 라플라스 방정식을 통한 스토크래틱 그람형 도입을 통해 균형 절단과 단순화된 편미분 근사를 추가적인 Lévy 노이즈를 갖는 선형 스토크래틱 시스템으로 성공적으로 확장한다.
• BT와 SPA에 대한 오차 한계가 유도되었으며, 이는 무시된 한켈 특이값의 합에 의존함을 보여주며, 결정론적 오차 한계를 스토크래틱 설정으로 일반화한다.
• 오차 한계는 타이트하며, Lévy 노이즈가 있는 2차 감쇠된 파동 방정식에서 수치적으로 검증되었으며, 실질적으로 작은 오차를 보여준다.
• 수치 결과는 BT와 SPA가 Lévy 노이즈 자극 조건 하에서도 시스템 행동의 손실이 최소화된 정확한 감소된 차원 모델을 제공함을 확인한다.
• 세미군 이론과 이토 등식을 통해 감소된 시스템이 전체 시스템으로 수렴하는 이론적 근거가 확립되었으며, 감소된 차원이 증가함에 따라 평균 제곱 오차는 0으로 수렴한다.
• 이 프레임워크는 Lévy 노이즈를 갖는 고차원 스토크래틱 시스템에 대한 실용적이고 이론적으로 타당한 MOR 접근법을 제공하며, 반복 샘플링이 필요한 몬테카를로 시뮬레이션에 특히 유용하다.