[논문 리뷰] Support Recovery for Sparse Deconvolution of Positive Measures

Quentin Denoyelle, Vincent Duval외 1명|arXiv (Cornell University)|2015. 06. 27.

Sparse and Compressive Sensing Techniques참고 문헌 20인용 수 9

한 줄 요약

이 논문은 BLASSO 볼록 최적화 프레임워크를 사용하여 양의 라돈 측도의 스파스 디컨볼루션에서 지지의 복원을 연구한다. 노이즈와 정규화 파rameter가 $ t^{2N-1} $ 보다 더 빠르게 감소할 경우, 해는 원래의 스파이크 위치와 강도를 정확하게 복원하며, 노이즈가 사라질수록 진짜 파rameter로 수렴함을 확립한다.

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 $ 를 결정하는 것.
• 노이즈 및 정규화 파rameter가 0으로 수렴함에 따라 복원된 강도와 위치가 진짜 값으로 수렴하는 수렴 속도를 확립하는 것.
• 신호 대 노이즈 비율과 최소 간격이 해의 유일성과 정확성 보장에 미치는 역할을 규명하는 것.

제안 방법

• 전체 변동(norm)을 사용한 라돈 측도 공간 위의 볼록 최적화로 스파스 디컨볼루션 문제를 설정하며, 이는 연속적 상황에서 $ \ell^1 $ 정규화와 동치이다.
• BLASSO(Basis Pursuit in the L1 sense for measures) 프로그램을 적용: $ |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} $ 보다 더 빠르게 감소함에 따라 진짜 측도로 수렴함을 보인다.
• 노이즈 $ w/\lambda $, $ w/t^{2N-1} $, $ \lambda/t^{2N-1} $ 가 충분히 작을 경우, 해는 정확히 $ N $ 개의 스파이크를 가지며 안정적인 복원을 보장한다.

실험 결과

연구 질문

• RQ1 노이즈 수준 $ w $, 정규화 파rameter $ \lambda $, 최소 스파이크 간격 $ t $ 가 어떤 조건을 만족할 때 BLASSO 프로그램이 스파스한 양의 측도의 정확한 지지를 복원하는가?
• RQ2 노이즈와 정규화 파arameter가 $ t^{2N-1} $ 에 비해 얼마나 빠르게 감소해야 스파이크 위치와 강도를 정확히 복원할 수 있는가?
• RQ3 노이즈 $ w \to 0 $ 이고 정규화 파arameter $ \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 $ 가 노이즈 및 정규화 파arameter에 비해 충분히 크면 해는 유일하며 원래 측도를 정확히 복원한다.
• 이중 증명자 $ \eta_{\lambda_n,t_n} $ 가 $ C^{2N} $ 노름에서 $ \eta_W $ 로 수렴함을 통해 최적성 조건이 극한에서 만족됨을 보인다.
• 분석은 $ X = \mathbb{R} $ 와 $ X = \mathbb{T} $ 의 경우 모두 성립하며, 실수선과 토러스 상에서 복원 프레임워크의 강건성을 보여준다.