[논문 리뷰] Phase Recovery, MaxCut and Complex Semidefinite Programming

Irène Waldspurger, Alexandre d’Aspremont외 1명|arXiv (Cornell University)|2012. 06. 01.

Advanced X-ray Imaging Techniques참고 문헌 35인용 수 54

한 줄 요약

이 논문은 복소수 위상 복원 문제를 위한 볼록 준화 프로그래밍 이론인 PhaseCut을 제안한다. 이는 MaxCut SDP를 모방하여 블록 좌표 강하법을 통해 효율적인 행렬-벡터 반복을 가능하게 한다. 노이즈가 없는 조건에서 PhaseLift와 수정된 PhaseCut 간의 등가성을 입증하고, 이전 방법들에 비해 노이즈 및 희박한 설정에서 더 뛰어난 안정성을 보인다.

ABSTRACT

Phase retrieval seeks to recover a signal x from the amplitude |Ax| of linear measurements. We cast the phase retrieval problem as a non-convex quadratic program over a complex phase vector and formulate a tractable relaxation (called PhaseCut) similar to the classical MaxCut semidefinite program. We solve this problem using a provably convergent block coordinate descent algorithm whose structure is similar to that of the original greedy algorithm in Gerchberg-Saxton, where each iteration is a matrix vector product. Numerical results show the performance of this approach over three different phase retrieval problems, in comparison with greedy phase retrieval algorithms and matrix completion formulations.

연구 동기 및 목표

• 복소수 벡터 공간에서의 비볼록 위상 복원 문제를 복소수 단위 토러스 위의 이차 프로그래밍으로 재정의함으로써 새로운 볼록 준화를 가능하게 한다.
• MaxCut 준화 프로그램에서 영감을 얻은 실용적인 볼록 준화인 PhaseCut을 개발하여 복소수 측정치에 특화한다.
• 기존 Gerchberg-Saxton 방법과 동일한 반복 복잡도를 갖는 수렴 보장이 가능한 블록 좌표 강하 알고리즘을 설계한다.
• 노이즈가 없는 경우 PhaseCut와 PhaseLift 간의 이론적 등가성 및 강도 조건을 확립하고, 노이즈 있는 환경에서의 안정성 비교를 수행한다.
• 실험적으로 PhaseCut가 그레디 알고리즘 및 행렬 완성 공식화보다 노이즈 및 희박성 조건에서 더 뛰어난 강건성(robustness)을 보임을 입증한다.

제안 방법

• 모든 $ i $ 에 대해 $ |u_i| = 1 $ 인 복소수 위상 벡터 $ u riangleq rac{Ax}{|Ax|} $ 위에서 비볼록 이차 프로그래밍으로 위상 복원 문제를 재정의한다.
• Rank-one 제약 조건 $ X = xx^* $ 를 사용하여 MaxCut SDP와 유사한 복소수 정수형 프로그래밍으로 문제를 레이저링하여 볼록 준화인 PhaseCut를 유도한다.
• 각 반복에서 행렬-벡터 곱셈을 수행하는 블록 좌표 강하법을 사용하여 저비용 반복과 수렴 보장을 확보한다.
• 노이즈가 없는 조건에서 $ A $ 가 단사적이고 $ b = |Ax| $ 가 0이 아닌 경우, PhaseLift와 이론적으로 등가가 되도록 하기 위해 수정된 PhaseCut 공식을 도입한다.
• 단일 요소( singleton )로 구성된 제약 조건 행렬의 구조를 활용하여, PhaseLift보다 문제 크기가 더 크더라도 효율적인 계산을 가능하게 한다.
• 진짜 Gram 행렬과 복원된 Gram 행렬 간의 차이에 대한 추적 기반 경계를 분석하여 안정성 분석을 수행하며, $ ext{Tr}(V_{PC}^ot) riangleq d_1(V_{PC}, ext{range}(A)) $ 를 사용한다.

실험 결과

연구 질문

• RQ1 위상 복원 문제는 복소수 단위 토러스 위의 비볼록 이차 프로그래밍으로 재정의될 수 있으며, 이를 통해 새로운 볼록 준화가 가능할까?
• RQ2 MaxCut 유사한 정수형 프로그래밍 준화인 PhaseCut는 위상 복원 문제에서 PhaseLift와 비교해 유사하거나 더 나은 성능을 내는가?
• RQ3 PhaseCut 준화가 강도를 가지는 조건은 무엇이며, 이는 PhaseLift의 강도 조건과 어떻게 관련이 있는가?
• RQ4 노이즈 및 희박성 조건에서 PhaseCut는 PhaseLift 및 그레디 알고리즘과 비교해 안정성이 뛰어난가?
• RQ5 PhaseCut에 대한 블록 좌표 강하 알고리즘은 Gerchberg-Saxton 방법과 유사한 계산 효율성을 유지하면서도 수렴 보장을 받을 수 있는가?

주요 결과

• PhaseCut는 MaxCut SDP와 유사한 구조를 지닌 위상 복원 문제의 볼록 준화이며, 블록 좌표 강하법에 의한 효율적 해법이 가능하다.
• 논문은 PhaseLift의 강도 조건이 수정된 PhaseCut의 강도 조건을 유도하고, 그 반대도 미약한 가정 하에 성립함을 증명하여 노이즈가 없는 경우 이론적 등가성을 확립한다.
• 노이즈가 없는 설정에서 $ A $ 가 단사적이고 $ |Ax| $ 가 0이 아닌 성분을 모두 포함할 경우, PhaseLift와 수정된 PhaseCut는 등가이다.
• 수치적 비교를 통해 PhaseCut는 노이즈가 있는 환경에서 PhaseLift보다 더 뛰어난 안정성을 보이며, 특히 $ b = |Ax| $ 가 희박할 경우 두드러진다.
• PhaseCut에 대한 블록 좌표 강하 알고리즘은 Gerchberg-Saxton 방법과 동일한 반복 복잡도(행렬-벡터 곱셈)를 가지며, 이는 효율성을 보장한다.
• 이론적 경계는 $ ext{Tr}(V_{PC}^ot) riangleq d_1(V_{PC}, ext{range}(A)) riangleq ext{Tr}((I - AA^ op)V_{PC}(I - AA^ op)) riangleq ext{Tr}(V_{PC}(I - AA^ op)) $ 를 포함하며, 이 값은 $ orm{b_{ ext{n,PC}}}_2^2 $ 로 유계이므로 안정성 보장을 이끌어낸다.