[논문 리뷰] Distributed source identification for wave equations: an observer-based approach

Marianne Chapouly, Mazyar Mirrahimi|arXiv (Cornell University)|2010. 05. 03.

Stability and Controllability of Differential Equations참고 문헌 4인용 수 5

한 줄 요약

이 논문은 초기 조건이 알려진 0인 파동 방정식에서 경계 노이먼 측정값만을 사용하여 분포형 소스 항 $ q(x) $ 를 특정하기 위한 관측기 기반의 오고-되돌아오기 반복 알고리즘을 제안한다. 소스 식별 문제를 결합된 파동-진동계에 대한 초기 상태 추정 문제로 변환함으로써, 이 방법은 관측 시간 $ T=2 $ 에서 진짜 소스로 수렴함을 보장하며, 리아푸노프 분석과 10%의 노이즈에 대한 수치 시뮬레이션을 통해 검증된다.

ABSTRACT

In this paper, we consider a wave equation on a bounded interval where the initial conditions are known (are zero) and we are rather interested in identifying an unknown source term $q(x)$ thanks to the measurement output $y$ which is the Neumann derivative on one of the boundaries. We use a back-and-forth iterative procedure and construct well-chosen observers which allow to retrieve $q$ from $y$ in the minimal observation time.

연구 동기 및 목표

• 초기 조건이 알려진 0인 파동 방정식에서 공간적으로 분포된 소스 항 $ q(x) $ 를 특정하는 역 문제를 해결하기 위해.
• 유한한 시간 간격 동안 경계 노이먼 측정값 $ y(t) = u_x(t,0) $ 로부터 $ q(x) $ 를 재구성하는 효율적인 추정 알고리즘을 개발하기 위해.
• 1차원 파동 방정식에서 관측 가능성을 위한 이론적 하한선인 최소 관측 시간 $ T=2 $ 에서 소스 재구성을 달성하기 위해.
• 무한차원 시스템의 성격을 고려하여 리아푸노프 기법과 라살의 불변성 원리를 사용하여 추정기의 수렴을 보장하기 위해.
• 노이즈가 있는 조건에서 수치적으로 방법을 검증하여 강건성과 빠른 수렴을 입증하기 위해.

제안 방법

• 동차 파동 방정식과 조화 진동자로 이루어진 계열 시스템을 도입하여 원래의 소스 식별 문제를 등가의 초기 상태 추정 문제로 변환하기 위해.
• 시간에 따라 전진하는 관측기와 후진하는 관측기를 각각 구성하고, 시간 간격 $ T=2 $ 에 걸쳐 반복적으로 적용하기 위해.
• 오차 동역학이 0으로 수렴함을 보이기 위해 리아푸노프 함수를 사용하여 수렴성을 증명하기 위해.
• 전진 관측기가 $ t=0 $ 에서 $ t=T $ 까지 상태를 추정하고, 후진 관측기가 $ t=T $ 에서 $ t=0 $ 으로 상태를 재구성하는 오고-되돌아오기 절차를 적용하기 위해, 이때 이득은 $ \gamma_1 = 1 $, $ \gamma_2 = 0.5 $ 이다.
• 원래 시스템과 변환된 시스템(2)-(3) 사이의 등가성에 기반하여, $ q(x) $ 가 파동 성분의 초기 조건가 된다.
• 유한 차분을 이용한 공간 및 시간 이산화 방법을 사용하여 수치 시뮬레이션을 수행하며, $ \Delta x = 0.05 $, $ \Delta t = 2.5 \times 10^{-4} $, 출력에 10%의 백색 노이즈를 적용하기 위해.

실험 결과

연구 질문

• RQ1 초기 조건이 0일 때, 경계 노이먼 측정값으로부터 파동 방정식의 분포형 소스 항 $ q(x) $ 를 특정할 수 있는가?
• RQ2 관측기 기반의 오고-되돌아오기 절차를 사용하여 최소 관측 시간 $ T=2 $ 에서 소스 재구성을 달성할 수 있는가?
• RQ3 리아푸노프 안정성 이론은 분포형 소스 식별을 위한 무한차원 관측기에 어떻게 적용될 수 있는가?
• RQ4 알려진 초기 조건은 원래 관측 불가능한 시스템을 관측 가능하게 만드는 데 어떤 역할을 하는가?
• RQ5 수치 시뮬레이션에서 관측기 기반 알고리즘이 측정 노이즈에 얼마나 강건한가?

주요 결과

• 제안된 오고-되돌아오기 관측기 알고리즘이 최소 관측 시간 $ T=2 $ 에서 진짜 소스 $ q(x) $ 로 수렴하며, 이는 $ (0,1) $ 에서 파동 방정식의 관측 가능성에 대한 이론적 하한선이다.
• 리아푸노프 함수는 단조 감소하며 0으로 수렴함을 보여, 오차 동역학의 점근적 안정성과 추정기의 수렴성을 증명한다.
• 10%의 백색 노이즈가 포함된 수치 시뮬레이션에서, 추정기는 50회 반복 이내에 $ q(x) = x - x^2 $ 에 수렴하며, $ \hat{q} $ 와 $ q $ 사이의 $ L^2 $-거리가 0으로 수렴함을 보였다.
• 식 (42)에 정의된 리아푸노프 함수는 반복 과정 동안 단조 감소하며, 추정 과정의 안정성을 확인한다.
• 이론적 분석에 따르면 오차 동역학의 유일한 불변 집합은 원점이며, 이는 오차가 0으로 수렴함을 의미하므로 $ q(x) $ 의 정확한 재구성을 보장한다.
• 초기 조건이 0임을 알고 있다면 전체 초기 상태가 알려지지 않았더라도 소스를 성공적으로 식별할 수 있으며, 이는 원래 시스템의 관측 불가능성을 깨뜨린다.