[논문 리뷰] Recovering of dielectric constants of explosives via a globally strictly convex cost functional

Michael V. Klibanov, Nguyễn Trung Thành|arXiv (Cornell University)|2014. 08. 04.

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

한 줄 요약

이 논문은 라플라스 변환과 카르레만 가중치 함수를 사용한 비용 기능을 통해 경계 측정값에서 폭발물의 비유전율 계수를 복원하기 위한 전역 수렴성 있는 수치적 방법을 제안한다. 라거르 함수와 카르레만 가중치 함수를 이용해 엄격한 볼록 비용 기능를 구성함으로써, 기울기 기반 최적화의 전역 수렴성을 보장하여, 초기 추측값이 열악한 경우에도 신뢰할 수 있는 재구성을 가능하게 하며, 이는 비연속 및 연속 계수를 가진 1차원 수치 실험을 통해 검증되었다.

ABSTRACT

The inverse problem of estimating dielectric constants of explosives using boundary measurements of one component of the scattered electric field is addressed. It is formulated as a coefficient inverse problem for a hyperbolic differential equation. After applying the Laplace transform, a new cost functional is constructed and a variational problem is formulated. The key feature of this functional is the presence of the Carleman Weight Function for the Laplacian. The strict convexity of this functional on a bounded set in a Hilbert space of an arbitrary size is proven. This allows for establishing the global convergence of the gradient descent method. Some results of numerical experiments are presented.

연구 동기 및 목표

• 지속적 땅속 레이더를 통한 폭발물 탐지에서 비선형 비용 기능성으로 인해 국소 수렴성 문제를 겪는 기존 방법의 도전 과제를 해결하기 위해.
• 초기 추측값의 정확성에 의존하지 않고, 쌍곡 파동 방정식의 계수 추정을 위한 전역 수렴성 알고리즘 개발을 위해.
• 라플라스 도메인 설정에서 무한 적분의 절단을 $L_2(\underline{s}, \infty)$ 내의 정규직교 함수 기반의 급수 전개로 대체하기 위해.
• 힐버트 공간의 임의의 유계 집합에서 비용 기능의 엄격한 볼록성을 확립하여 기울기 하강법의 전역 수렴성을 보장하기 위해.
• 1차원 수치 실험을 통해 계산 가능성과 정확도 향상을 입증하며, 하이브리드 방법과 국소 수렴성 있는 해법 간의 비교를 수행하기 위해.

제안 방법

• 산란 전기장 성분의 경계 측정값을 사용하여, 비선형 편미분 방정식에 대한 계수 역문제로 문제를 설정한다.
• 시간에 의존하는 파동 방정식에 라플라스 변환을 적용하여, 라플라스 도메인 내의 비선형 타원형 방정식 시스템으로 변환한다.
• 라그레르 함수를 사용한 급수 전개를 통해 라플라스 도메인 내의 해를 표현하며, 이는 $L_2(\underline{s}, \infty)$ 내의 정규직교 기저를 형성한다.
• 카르레만 가중치 함수(CWF)를 포함한 최소 제곱 비용 기능을 구성하여 힐버트 공간 내의 임의의 유계 집합에서 엄격한 볼록성을 강제한다.
• 충분히 큰 매개수 $\lambda$에 대해 CWF가 비용 기능의 엄격한 볼록성을 보장함을 증명하여 기울기 방법의 전역 수렴성을 확보한다.
• 이중 단계 알고리즘을 적용한다: 첫 번째 단계에서는 전역 수렴성 방법을 통해 양호한 초기 근사값을 확보하고, 두 번째 단계에서는 국소 수렴성 방법(SQP)을 사용해 이를 정밀화한다.

실험 결과

연구 질문

• RQ1 폭발물의 전자기 영상에서 계수 역문제에 대해 전역 수렴성 방법을 구성할 수 있는가?
• RQ2 정규직교 기저를 사용한 급수 전개를 통해 라플라스 도메인 내의 무한 적분 절단을 피할 수 있는가?
• RQ3 비용 기능에 카르레만 가중치 함수를 포함함으로써, 힐버트 공간의 임의의 유계 집합에서 엄격한 볼록성이 보장되는가?
• RQ4 매우 열악한 초기 추측값(예: 배경 매질)에서 시작할 경우, 제안된 방법이 국소 수렴성 방법보다 우월한 성능을 보일 수 있는가?
• RQ5 1차원 수치 설정에서 비연속 및 연속 비유전율 프로파일의 재구성이 얼마나 정확한가?

주요 결과

• 카르레만 가중치 함수의 매개수 $\lambda$가 충분히 크면, 힐버트 공간의 임의의 유계 집합에서 비용 기능이 엄격한 볼록성을 띠며, 이는 기울기 하강법의 전역 수렴성을 보장한다.
• 1차원 수치 실험 결과, 전역 및 국소 최적화를 조합한 하이브리드 알고리즘이 조각별 상수 및 연속 비유전율 계수를 정확하게 재구성함을 보여주며, 큰 점프가 있는 경우(예: $c(x) = 1 + 14\chi[0.03,0.1]$)에도 성능이 우수하다.
• 균일한 배경 매질에서 초기화된 국소 수렴성 방법은 정확한 해에 수렴하지 못하는 반면, 제안된 하이브리드 방법은 모든 시험 케이스에서 성공적으로 수렴한다.
• 전역 수렴성 첫 번째 단계의 출력에서 시작할 경우, 국소 수렴성 방법의 재구성 정확도가 크게 향상됨을 예제 1–4를 통해 입증하였다.
• 예제 3에서 $c(x) = 1 - 0.5\chi[0.03,0.15]$로 배경 매질보다 낮은 비유전율 계수를 가진 타겟도 성공적으로 복원하였다.
• 알고리즘이 타겟 위치 및 계수 변화에 대해 강건하며, 측정점 근처에서 정확한 결과를 보이며 실질적인 GPR 데이터 사전 처리와 일치하는 경향을 보였다.