[논문 리뷰] Solving an inverse elliptic coefficient problem by convex non-linear semidefinite programming

Bastian Harrach|arXiv (Cornell University)|2021. 05. 24.

Numerical methods in inverse problems참고 문헌 19인용 수 4

한 줄 요약

이 논문은 유한한 측정값을 가진 경우에 대해 역 타원형 계수 문제, 특히 로빈 전도 문제를 해결하기 위한 볼록 비선형 정준형 프로그래밍 설정을 제안한다. 뉴먼-디리클레 연산자의 단조성과 볼록성 성질을 활용하여, 역 문제를 국소 최솟값이 없는 유일하게 해가 존재하는 볼록 최적화 문제로 동치로 재구성함으로써, 로에브너 순서를 통한 명시적 오차 추정 및 측정 수 기준을 가능하게 한다.

ABSTRACT

Several applications in medical imaging and non-destructive material testing lead to inverse elliptic coefficient problems, where an unknown coefficient function in an elliptic PDE is to be determined from partial knowledge of its solutions. This is usually a highly non-linear ill-posed inverse problem, for which unique reconstructability results, stability estimates and global convergence of numerical methods are very hard to achieve. The aim of this note is to point out a new connection between inverse coefficient problems and semidefinite programming that may help addressing these challenges. We show that an inverse elliptic Robin transmission problem with finitely many measurements can be equivalently rewritten as a uniquely solvable convex non-linear semidefinite optimization problem. This allows to explicitly estimate the number of measurements that is required to achieve a desired resolution, to derive an error estimate for noisy data, and to overcome the problem of local minima that usually appears in optimization-based approaches for inverse coefficient problems.

연구 동기 및 목표

• 최적화 기반 역 계수 문제에서의 국소 최솟값 문제를 해결하기 위해.
• 유한한 측정값을 가진 역 타원형 계수 문제에 대해 전역 수렴성과 볼록 재구성된 최적화 문제를 제공하기 위해.
• 목표 해상도를 달성하기 위해 필요한 측정 수의 명시적 기준을 도출하기 위해.
• 노이즈가 있는 데이터에 대한 오차 추정을 역 계수 문제에서 유도하기 위해.
• 정준형 프로그래밍과 역 계수 문제를 연결하여 새로운 이론적 및 계산적 프레임워크를 제공하기 위해.

제안 방법

• 계수 벡터에서 대칭 행렬로 사상하는 행렬값 함수 F를 포함하는 유한 차원 역 문제로 역 로빈 전도 문제를 설정하기.
• 뉴먼-디리클레 연산자 Λ(γ)의 프레셰 미분 가능성, 볼록성 및 단조성 성질을 활용하여 F의 필요 조건을 확보하기.
• 로에브너 순서에 기반한 충분 조건 수립: 모든 j,k에 대해 F′(zj,k)dj ≺ 0이면, 역 문제는 볼록 최적화를 통해 유일하게 해가 존재한다.
• 역 문제를 ℓ¹-노름의 계수 벡터를 최소화하는 볼록 정준형 프로그래밍 문제로 재구성하기.
• 방향 도함수의 스펙트럼 노름 및 고유값 분석을 통해 안정성 및 오차 경계 유도하기.
• 측정 수를 증가시키며 기준을 반복적으로 적용하여, 방향 도함수 조건을 충족할 때까지 반복 적용하기.

실험 결과

연구 질문

• RQ1 유한한 측정값을 가진 비볼록 역 타원형 계수 문제를 볼록 최적화 문제로 재구성할 수 있는가?
• RQ2 주어진 해상도로 계수를 유일하게 재구성하기 위해 필요한 최소 측정 수는 얼마인가?
• RQ3 제안된 방법이 기존 최적화 기반 접근법에서 흔히 발생하는 국소 최솟값 문제를 제거할 수 있는가?
• RQ4 노이즈가 있는 측정값에 대해 역 계수 문제에서 오차 추정을 어떻게 도출할 수 있는가?
• RQ5 로에브너 순서와 단조성이 역 문제의 볼록 재구성에 어떻게 기여하는가?

주요 결과

• 유한한 측정값을 가진 역 로빈 전도 문제는 유일하게 해가 존재하는 볼록 비선형 정준형 프로그래밍 문제로 동치로 재구성될 수 있다.
• 방향 도함수의 성질에 기반한 유일한 해 존재를 위한 충분 조건이 도출되었으며, 이는 유한한 수의 정방향 해를 계산하여 수치적으로 검증할 수 있다.
• 측정 수는 조건을 충족할 때까지 단계적으로 증가시킴으로써 체계적으로 결정될 수 있다.
• 노이즈가 있는 데이터에 대해 ∥Ŷ − Yδ∥₂ ≤ δ일 경우, 재구성 오차는 ∥γ̂ − γδ∥∞ ≤ 2δ(n−1)/λ로 유계가 되며, 여기서 λ > 0은 방향 도함수의 최대 고유값들의 최소값이다.
• 재구성된 문제의 볼록성 덕분에 국소 최솟값을 피할 수 있으며, 특수하게 설계된 측정값이 필요하지 않다.
• 이 프레임워크는 정준형 프로그래밍과 역 계수 문제 간의 첫 번째 명시적 연결을 제공하며, 전역 수렴성과 안정성 분석의 새로운 길을 열어 놓았다.