[논문 리뷰] Stable Phase Retrieval in Infinite Dimensions

Rima Alaifari, Ingrid Daubechies외 2명|arXiv (Cornell University)|2016. 08. 31.

Advanced X-ray Imaging Techniques참고 문헌 32인용 수 10

한 줄 요약

이 논문은 무한차원 힐버트 공간에서 안정적인 위상 복원을 위한 새로운 개념을 제안한다. 전역 위상이 아닌 국소적 위상 모호성을 允가함으로써, 서로소인 주파수 또는 시간 영역에서 서로 다른 위상 인자로 신호를 복원한다. 저자들은 이 유연한 개념 하에서 가보 프레임과 코시 웨이브릿 프레임에 대해 안정적인 복원이 가능하다는 것을 증명하며, 무한차원에서의 고전적 위상 복원의 잘 알려진 불안정성 문제를 해결한다.

ABSTRACT

The problem of phase retrieval is to determine a signal $f\in \mathcal{H}$, with $\mathcal{H}$ a Hilbert space, from intensity measurements $|F(ω)|$, where $F(ω):=\langle f , φ_ω angle$ are measurements of $f$ with respect to a measurement system $(φ_ω)_{ω\in Ω}\subset \mathcal{H}$. Although phase retrieval is always stable in the finite dimensional setting whenever it is possible (i.e. injectivity implies stability for the inverse problem), the situation is drastically different if $\mathcal{H}$ is infinite-dimensional: in that case phase retrieval is never uniformly stable [8, 4]; moreover the stability deteriorates severely in the dimension of the problem [8]. On the other hand, all empirically observed instabilities are of a certain type: they occur whenever the function $|F|$ of intensity measurements is concentrated on disjoint sets $D_j\subset Ω$, i.e., when $F= \sum_{j=1}^k F_j$ where each $F_j$ is concentrated on $D_j$ (and $k \geq 2$). Motivated by these considerations we propose a new paradigm for stable phase retrieval by considering the problem of reconstructing $F$ up to a phase factor that is not global, but that can be different for each of the subsets $D_j$, i.e., recovering $F$ up to the equivalence $$ F \sim \sum_{j=1}^k e^{i α_j} F_j.$$ We present concrete applications (for example in audio processing) where this new notion of stability is natural and meaningful and show that in this setting stable phase retrieval can actually be achieved, for instance if the measurement system is a Gabor frame or a frame of Cauchy wavelets.

연구 동기 및 목표

• 무한차원 힐버트 공간에서 단사성은 안정성을 보장하지 못하는 위상 복원의 근본적 불안정성 문제를 다루기.
• 강도 측정치가 도메인의 서로소 부분집합에 집중될 경우 발생하는 경험적 불안정성의 원인을 규명하기.
• 전역 위상 인자를 요구하는 대신, 서로소 성분 각각에 대해 독립적인 위상 인자를 允가하는 새로운 안정성 프레임워크를 도입하기.
• 가보 및 코시 웨이브릿 프레임과 같은 구조적 측정 시스템에서 이 새로운 개념 하에서 안정적인 위상 복원이 가능하다는 것을 보여주기.
• 소음 존재 조건과 무한차원 설정 하에서도 새로운 동치관계 하에서의 안정성에 대한 이론적 보장을 제공하기.

제안 방법

• 힐버트 공간 $ \mathcal{H} $ 내 신호 $ f \in \mathcal{H} $ 를 강도 측정치 $ |\langle f, \varphi_\omega \rangle| $ 로 복원하는 위상 복원을 정의하며, $ (\varphi_\omega)_{\omega \in \Omega} $ 는 $ \mathcal{H} $ 에서의 프레임이다.
• 새로운 동치관계 도입: $ F \sim \sum_{j=1}^k e^{i\alpha_j} F_j $, 여기서 각 $ F_j $ 는 서로소 집합 $ D_j $ 에서 지지되며, 국소적 위상 자유도를 허용한다.
• 복소해석학 도구, 특히 해석적 파oincaré 부등식과 경계 정리(Trace theorems)를 사용하여 유사한 강도 크기를 가진 신호 간의 차이를 추정한다.
• $ |F| $ 의 그래디언트의 유계 $ L^p $-노름을 활용하여 신호 변동을 제어하며, 핵심 항등식 $ \|F'\|_{L^p(D_+)} = \|\nabla|F|\|_{L^p(D_+)} $ 을 활용한다.
• 보간법과 커팅 함수 $ \eta $ 를 사용하여 안정성 한계를 확립하며, 도메인 $ D_+ $ 와 $ D_0 $ 에서의 해석함수 이론과 경계 추정을 적용한다.
• 재구성 오차가 $ ||F_1| - |F_2|| $ 의 $ W^{1,p} $-노름에 의해 유계임을 증명하며, 이는 서로소 지지 영역에서의 강도 차이를 정량화한다.

실험 결과

연구 질문

• RQ1 왜 고전적 위상 복원은 단사성은 보장하지만 무한차원 힐버트 공간에서는 항상 불안정한가?
• RQ2 강도 측정치가 서로소 주파수 또는 시간 영역에 집중될 경우 발생하는 관측된 경험적 불안정성의 원인은 무엇인가?
• RQ3 전역 위상 인자를 완화하여 서로소 성분 각각에 대해 독립적인 위상 인자를 允가함으로써 위상 복원을 안정화시킬 수 있는가?
• RQ4 가보 프레임, 코시 웨이브릿 등과 같은 측정 시스템에서 이 새로운 국소적 위상 동치 하에서 안정적인 위상 복원이 성립하는 조건은 무엇인가?
• RQ5 새로운 개념 하에서 무한차원 설정에서 정량적 안정성 한계를 어떻게 유도할 수 있는가?

주요 결과

• 고전적 위상 복원은 무한차원 힐버트 공간에서조차 단일한 안정성 보장이 없으며, 차원이 증가함에 따라 심각한 열화가 발생한다.
• 불안정성의 주요 원인은 신호 $ F $ 가 $ F = \sum_{j=1}^k F_j $ 로 분해되며, $ \text{supp}(F_j) \subset D_j $ 를 만족하는 서로소 집합 $ D_j $ 에 강도가 집중될 경우이다.
• 각 $ D_j $ 에서 다른 위상 인자 $ e^{i\alpha_j} $ 를 允가함으로써, 새로운 동치관계 $ F \sim \sum_{j=1}^k e^{i\alpha_j} F_j $ 하에서 문제는 안정해진다.
• 이 새로운 프레임워크 하에서 가보 프레임과 코시 웨이브릿 프레임에 대해 안정적인 위상 복원이 달성되며, 안정성 상수는 측정 시스템의 기하학적 및 해석적 성질에 따라 달라진다.
• 재구성 오차는 $ W^{1,p} $-세미노름으로 유계이며, 이는 서로소 지지 영역에서의 강도 차이가 작을수록 신호 차이도 작다는 것을 보여준다.
• 증명은 복소해석학 도구에 기반한다: 해석적 파oincaré 부등식, 경계 정리, 커팅 함수 추정. 이는 최종적으로 $ \frac{\Delta}{\delta^2} \||F_1| - |F_2|\|_{W^{1,p}} $ 비례하는 bound 로 이어진다.