[論文レビュー] Support Recovery for Sparse Deconvolution of Positive Measures

Quentin Denoyelle, Vincent Duval+1 more|arXiv (Cornell University)|Jun 27, 2015

Sparse and Compressive Sensing Techniques参考文献 20被引用数 9

ひとこと要約

本稿は、正のラドン測度のスパースデコンボリューションにおけるサポート回復を、BLASSO凸最適化フレームワークを用いて研究する。ノイズと正則化パラメータが $ t^{2N-1} $ よりも速く減少する場合、解は元のスパイクの位置と振幅を一意に回復し、ノイズが消失するにつれて真のパラメータに収束することを確立する。

ABSTRACT

We study sparse spikes deconvolution over the space of Radon measures on $\mathbb{R}$ or $\mathbb{T}$ when the input measure is a finite sum of positive Dirac masses using the BLASSO convex program. We focus on the recovery properties of the support and the amplitudes of the initial measure in the presence of noise as a function of the minimum separation $t$ of the input measure (the minimum distance between two spikes). We show that when ${w}/λ$, ${w}/t^{2N-1}$ and $λ/t^{2N-1}$ are small enough (where $λ$ is the regularization parameter, $w$ the noise and $N$ the number of spikes), which corresponds roughly to a sufficient signal-to-noise ratio and a noise level small enough with respect to the minimum separation, there exists a unique solution to the BLASSO program with exactly the same number of spikes as the original measure. We show that the amplitudes and positions of the spikes of the solution both converge toward those of the input measure when the noise and the regularization parameter drops to zero faster than $t^{2N-1}$.

研究の動機と目的

• スパースな正の測度（ディラックデルタの重ね合わせ）の正確なサポート（スパイクの位置と振幅）がBLASSOプログラムによって回復される条件を分析すること。
• 安定かつ正確な回復を可能にするスパイク間の最小距離 $ t $ を特定すること。
• ノイズと正則化パラメータが 0 に近づく際の、回復された振幅および位置の収束速度を確立すること。
• 信号対ノイズ比と最小分離距離が解の一意性と正確性を保証する役割を特定すること。

提案手法

• 全 variation 範囲におけるラドン測度の空間上の凸最適化としてスパースデコンボリューション問題を定式化し、連続的設定下での $ \ell^1 $ 正則化に相当する。
• BLASSO（測度におけるL1意味での基礎 pursuit）プログラムを適用：$ |m|(X) $ を最小化し、制約条件 $ \Phi m = y_t + w $ を満たす。ここで $ y_t $ はスパース測度のノイズ付き観測値である。
• 核 $ \varphi $ と測定作用素 $ \Phi $ を用いて構築された双対証明 $ \eta_{\lambda,t} $ を用い、解が元のスパイク構造と一致することを保証する。
• 双対証明 $ \eta_W $ の再帰的構築を用い、それが $ (2N-1) $-非退化であり、$ \|\eta_W\|_{L^\infty} \leq 1 $ を満たすことを示し、これにより正確なサポート回復が保証されることを示す。
• 一次最適性条件を用いた漸近的挙動の分析により、$ \lambda, w \to 0 $ のとき解が $ t^{2N-1} $ より速く収束することを示す。
• 解が正確に $ N $ 個のスパイクを持つのは、$ w/\lambda $、$ w/t^{2N-1} $、および $ \lambda/t^{2N-1} $ が十分に小さい場合に限られ、これにより安定な回復が保証される。

実験結果

リサーチクエスチョン

• RQ1 ノイズレベル $ w $、正則化パラメータ $ \lambda $、最小スパイク間隔 $ t $ がどの程度の条件下で、BLASSOプログラムがスパース正の測度の正確なサポートを回復するか。
• RQ2 ノイズと正則化パラメータが $ t^{2N-1} $ に対してどの程度の速さで減少すれば、スパイクの位置と振幅の正確な回復が保証されるか。
• RQ3 $ w \to 0 $ および $ \lambda \to 0 $ のとき、回復された振幅と位置の漸近的挙動は何か。また、真の値への収束が成立する条件は何か。
• RQ4 ヘルミート多項式を用いて構築された双対証明 $ \eta_W $ が、解に正しい数のスパイクをもたらし、元の測度と一致することをどのように保証するか。
• RQ5 $ (2N-1) $-非退化性が双対証明に与える役割は何か。これにより一意性と正確な回復が保証される条件は何か。

主な発見

• $ w/\lambda $、$ w/t^{2N-1} $、$ \lambda/t^{2N-1} $ が十分に小さい場合、BLASSO解は正確に $ N $ 個のスパイクを持ち、元の測度のサポートと一致する。
• $ \lambda \to 0 $ および $ w \to 0 $ のとき、回復されたスパイクの振幅および位置は、$ t^{2N-1} $ よりも速い減少が保証される限り、真の値に収束する。
• ヘルミート級数の切り捨てとして構築された双対証明 $ \eta_W $ は、$ \|\eta_W\|_{L^\infty} \leq 1 $ および $ \eta_W^{(2N)}(0) < 0 $ を満たし、双対性条件を介して正確な回復を保証する。
• 信号対ノイズ比が高く、最小分離距離 $ t $ がノイズと正則化パラメータに対して十分に大きい場合、解は一意的であり、元の測度を正確に回復する。
• 双対証明 $ \eta_{\lambda_n,t_n} $ が $ C^{2N} $ 範囲で $ \eta_W $ に収束することにより、極限において最適性条件が満たされることを示す。
• 解析は $ X = \mathbb{R} $ および $ X = \mathbb{T} $ の両方に対して成り立つため、実直線およびトーラス上での回復フレームワークの頑健性が示される。