[論文レビュー] Stable phase retrieval with low-redundancy frames

Bernhard G. Bodmann, Nathaniel Hammen|arXiv (Cornell University)|Feb 22, 2013

Advanced X-ray Imaging Techniques参考文献 19被引用数 14

ひとこと要約

本稿では、d次元の複素ヒルベルト空間における位相再構成を、ユニモジュラー定数を除き一意に回復可能とするために4d−4個のフレームベクトルで十分であることを確立している。また、8d−4個のベクトルでは、ノイズが存在する場合でも安定した回復が可能であり、フレーム係数が0から離れている限り、回復誤差はノイズレベルに比例する。この手法は、単位円上のディリクレ核と極化恒等式を用いた多項式再構成に基づいている。

ABSTRACT

We investigate the recovery of vectors from magnitudes of frame coefficients when the frames have a low redundancy, meaning a small number of frame vectors compared to the dimension of the Hilbert space. We first show that for vectors in d dimensions, 4d-4 suitably chosen frame vectors are sufficient to uniquely determine each signal, up to an overall unimodular constant, from the magnitudes of its frame coefficients. Then we discuss the effect of noise and show that 8d-4 frame vectors provide a stable recovery if part of the frame coefficients is bounded away from zero. In this regime, perturbing the magnitudes of the frame coefficients by noise that is sufficiently small results in a recovery error that is at most proportional to the noise level.

研究の動機と目的

• d次元ヒルベルト空間における位相再構成に必要なフレームベクトルの最小数を特定すること。
• 測定ノイズが存在する状況において、位相再構成がどの条件下で安定するかを確立すること。
• フレーム係数が0から離れている場合に、誤差がノイズレベルに比例することを保証する再構成手順を開発すること。
• 単位円上での多項式再構成技術を用いて、再構成プロセスの安定性を分析すること。
• 理論的境界を数値シミュレーションにより検証し、理論的ノイズ閾値を超えて線形誤差スケーリングが観察されることを示すこと。

提案手法

• 信号を単位円上での次数d−1以下の複素多項式として表現する。
• ノイズが付加された大きさ測定値を2(d−1)個の点で得、ディリクレ核補間を用いて|f(z)|²およびf′(z) overline{f(z)}を再構成する。
• 3つの回転位相シフト(ν^j, j=1,2,3)を用いた極化恒等式を適用し、摂動された大きさ|f(ω^l) + ν^j f′(ω^l)|²からf′(z) overline{f(z)}を抽出する。
• ニュートンの恒等式を用いて、再構成された|f(z)|²およびf′(z) overline{f(z)}から得られるモーメントから多項式係数を回復する。
• 一様ノルムの推定を用いて再構成誤差を評価し、測定ノイズがεで有界である場合に係数の誤差がO(ε)であることを示す。
• パーサバルの恒等式を用いてfのL²ノルムを推定し、信号をユニモジュラー定数を除き回復可能にする正規化を実施する。

実験結果

リサーチクエスチョン

• RQ1 d次元信号をユニモジュラー定数を除き一意に特定するためには、最小で何個のフレームベクトルが必要か？
• RQ2 大きさ測定値に小さなノイズが存在する場合、位相再構成を安定化させることは可能か？ どのような条件下で可能か？
• RQ3 大きさ測定値のノイズレベルに対して、再構成誤差はどのようにスケーリングされるか？
• RQ4 理論的安定性閾値と数値的シミュレーションでの実効性能の一致度はどの程度か？
• RQ5 フレーム係数が0から離れている場合、再構成手順はどの程度頑健に動作するか？

主な発見

• 4d−4個のフレームベクトルがあれば、任意のd次元信号をユニモジュラー定数を除き一意に特定可能である。
• 8d−4個のフレームベクトルでは、フレーム係数が0から離れている限り、安定した位相再構成が保証される。測定ノイズがεで有界であるとき、回復誤差はO(ε)である。
• 理論的安定性半径はε₀ = βm² / [(2d−1)d]であるが、数値的シミュレーションでは、εが50×ε₀に達するまで線形誤差スケーリングが継続することが観察された。
• 単位円上での|f(z)|²およびf′(z) overline{f(z)}の再構成誤差は(2d−1)εで有界であり、これにより安定な多項式再構成が可能である。
• L²ノルムの推定による正規化により、回復された多項式はユニモジュラー定数を除き正確であり、係数誤差はO(ε)である。
• 再構成手順は頑健である：数値実験により、理論的安定性閾値を超えて線形誤差スケーリングが確認され、理論的境界がタイトでない可能性が示唆された。