[논문 리뷰] Series reversion in Calder\'on's problem

Henrik Garde, Nuutti Hyvönen|arXiv (Cornell University)|2021. 05. 07.

Numerical methods in inverse problems참고 문헌 60인용 수 11

한 줄 요약

이 논문은 전기임피던스.tom그래피(EIT)에서 Calderón의 역도전도 문제를 해결하기 위한 새로운 시리즈 역행렬 방법을 제안한다. 이 방법은 부분 경계 측정값을 이용해 알려진 도전도 행렬에 대한 덧셈형 변화를 고차수 수치적으로 재구성할 수 있도록 한다. 정방향 맵의 해석성과 그 테일러 급수의 역행렬을 임의의 차수까지 적용함으로써, 프레셰 도함수의 가역성 조건 하에서 수렴성을 확보하였으며, 계산 복잡도는 선형화된 문제를 푸는 것과 동일하다. 이는 연속 모델과 매끄러운 완전 전극 모델 둘 다에 적용 가능하다.

ABSTRACT

This work derives explicit series reversions for the solution of Calder\'on's problem. The governing elliptic partial differential equation is $ abla\cdot(A abla u)=0$ in a bounded Lipschitz domain and with a matrix-valued coefficient. The corresponding forward map sends $A$ to a projected version of a local Neumann-to-Dirichlet operator, allowing for the use of partial boundary data and finitely many measurements. It is first shown that the forward map is analytic, and subsequently reversions of its Taylor series up to specified orders lead to a family of numerical methods for solving the inverse problem with increasing accuracy. The convergence of these methods is shown under conditions that ensure the invertibility of the Fr\'echet derivative of the forward map. The introduced numerical methods are of the same computational complexity as solving the linearised inverse problem. The analogous results are also presented for the smoothened complete electrode model.

연구 동기 및 목표

• Calderón의 역문제에서 알려진 도전도 계수에 대한 덧셈형 변화를 임의의 고차수로 재구성할 수 있는 수치적 방법의 집합을 개발하기 위해.
• 정방향 맵의 프레셰 도함수의 가역성 조건을 만족시키는 조건 하에서 이러한 방법의 수렴성을 확립하기 위해.
• 이 접근법을 연속 모델(CM)과 매끄러운 완전 전극 모델(SCEM) 둘 다로 확장하여 실질적인 EIT 측정에의 적용 가능성을 보장하기 위해.
• 특히 두 차원에서 시리즈 역행렬을 위해 필요한 투영 연산자 Q를 체계적으로 구성하는 체계적 프레임워크를 제공하기 위해.
• 이 방법의 계산 비용이 메쉬 해상도에 관계없이 선형화된 역문제를 푸는 것과 동일한 비율로 증가함을 보여주기 위해.

제안 방법

• 도전도에서 투영된 뉴먼-디리클레 연산자로의 정방향 맵 A ↦ Λ(A)가 행렬값 L∞ 계수 공간에서 해석적임을 증명하였다.
• 정방향 맵의 테일러 전개의 명시적 시리즈 역행렬을 활용하여, 변화 B를 차수 K까지 재구성하며, 점근적 공식 B = Σ_{j=1}^K F_j + O(||B||^{K+1})를 도출하였다.
• 투영된 프레셰 도함수 PDΛ(A; ·)P가 변화 B를 포함하는 알려진 부분공간 W에서 단사적이며, 닫힌 보조 부분공간로 사상되어야 하며, 이는 상대적 정방향 맵의 가역성을 보장한다.
• 적절한 바나흐 공간의 연산자에서 Q를 투영하여 W와의 호환성을 확보함으로써, 시리즈 역행렬을 위한 해석적 역행렬을 구성할 수 있도록 하였다.
• 두 차원에서는 힐베르트-슈미트 연산자 공간에서 닫힌 부분공간 위로의 수직 투영으로 Q를 체계적으로 선택하였으며, 이는 등각 사상과 단위 원판 ND 연산자의 스펙트럼 성질을 활용한다.
• W가 유한차원일 경우 Q의 명시적 구성 없이도 효율적으로 구현되며, 첫 번째 도함수 역행렬 단계에 표준 정규화 기법을 적용한다.

실험 결과

연구 질문

• RQ1 정방향 맵이 해석적이며, 투영된 프레셰 도함수가 가역적이라면, 시리즈 역행렬을 통해 Calderón 문제의 정방향 맵을 임의의 차수까지 역행렬로 구할 수 있는가?
• RQ2 계수 공간과 측정 데이터에 어떤 조건이 성립해야 시리즈 역행렬 방법이 덧셈형 변화를 재구성하는 데 수렴하는가?
• RQ3 두 차원에서 시리즈 역행렬을 안정적으로 수행하기 위해 투영 Q를 어떻게 체계적으로 구성할 수 있는가?
• RQ4 이 방법의 계산 비용이 선형화된 역문제를 푸는 것과 비교해 어느 정도 비례하는가?
• RQ5 시리즈 역행렬 프레임워크는 매끄러운 완전 전극 모델로 확장될 수 있으며, 이 경우 수렴성과 복잡도 성질은 유지되는가?

주요 결과

• 정방향 맵 A ↦ Λ(A)는 L∞ 행렬값 계수 공간에서 해석적이며, 이는 테일러 급수 전개와 이후 시리즈 역행렬을 가능하게 한다.
• 투영된 프레셰 도함수의 가역성 조건 하에서, 변화 B에 대한 수치적 재구성 방법의 수렴 차수는 O(||B||^{K+1})이다.
• 모든 고정된 K ∈ ℕ에 대해, 이 방법의 계산 복잡도는 해상도에 관계없이 K에 의존하는 상수배의 선형화된 역문제 해결 비용으로 제한된다.
• 이 방법은 연속 모델과 매끄러운 완전 전극 모델 양자리에 적용 가능하며, 둘 다에서 동일한 수렴성 및 복잡도 결과를 보인다.
• 두 차원에서는 Q를 힐베르트-슈미트 연산자 공간에서의 수직 투영으로 체계적으로 구성할 수 있으며, 이는 호환성과 안정성을 보장한다.
• 이 방법의 모든 불안정한 단계는 첫 번째 도함수 PDΛ(A; ·)P의 역행렬을 포함하며, 이는 선형 역문제에 대한 기존 정규화 기법을 직접 적용할 수 있음을 의미한다.