희소 신호 추적을 위한 지구 이동 거리 동적 정규화

본 논문은 시간에 따라 변하는 희소 신호를 추적하기 위한 새로운 동적 필터링 기법인 EMD‑DF(Earth Mover’s Distance Dynamic Filtering)를 제안한다. 기존의 Kalman 필터와 그 변형들은 ℓ2 혹은 ℓp 기반의 동적 정규화를 사용해 이전 추정값과 현재 추정값 사이의 차이를 최소화한다. 그러나 이러한 정규화는 신호 요소 간의 기하학적 관계, 즉 인접한 픽셀이나 주파수 인덱스 사이의 거리 정보를 반영하지 못한다. 결과적으로, 실제 신호가 작은 위치 이동만을 겪어도 큰 패널티가 부과되어 추적 성능이 저하된다. EMD는 두 히스토그램(또는 신호) 사이의 최소 운반 비용을 측정하는 최적 운송 거리로, 비용 행렬 R_{ij}에 의해 정의된 거리 정보를 자연스럽게 포함한다. 따라서 EMD를 동적 정규화 항으로 사용하면, 이전 프레임에서 현재 프레임으로의 질량(신호 에너지) 이동을 최소화하는 방향으로 추정이 이루어져, 인접한 위치 이동에 대해서는 낮은 비용, 먼 거리 이동에 대해서는 높은 비용을 부여한다. 전통적인 EMD는 흐름 변수 O(N²)를 필요로 하여 대규모 문제에 비현실적이다. 이를 해결하기 위해 저자들은 Beckmann 형태의 연속적인 EMD 모델을 도입한다. Beckmann 모델은 흐름 벡터 M을 최소화하는 형태로, 변수 수를 O(N)으로 줄인다. 또한, 신호의 총 질량이 서로 다를 경우를 다루기 위해 slack 변수 w와 v를 도입해 부분 운송(partial transport) 문제로 일반화한다. 변형된 최적화는 다음과 같다:  d_emd(x,y)=min_{M,w,v} ‖M‖_{2,1}  subject to div(M)+v−w=0, 0≤w≤x, 0≤v≤y, ‖w‖₁=‖v‖₁=min(‖x‖₁,‖y‖₁). 이 식을 기존의 BPDN‑DF 목적함수에 결합하면 최종 EMD‑DF 문제는  min_{x,M,u,v} ½‖yₙ−Ax‖₂² + λ‖x‖₁ + γ‖M‖_{2,1} − μu subject to 위의 제약을 만족한다. 여기서 u는 흐름에 의해 실제 이동된 질량을 나타내는 슬랙 변수이다. 이 최적화는 선형 제약을 가진 2차 프로그램이므로, CVX, Gurobi, Mosek 등 상용 솔버로 효율적으로 해결할 수 있다. 실험에서는 세 가지 시나리오를 통해 EMD‑DF의 장점을 검증한다. 첫 번째는 실제 랫트 해마 전기생리 데이터에 대한 시간‑주파수 스펙트럼 추정이다. 과잉완전 DFT 사전을 사용해 두 개의 주요 주파수를 추적했으며, EMD‑DF는 전통적인 STFT 스펙트로그램에 비해 주파수 해상도가 크게 향상되고, 5–10 Hz 대역의 누수가 사라졌다. 두 번째는 Kuramoto 모델 기반의 합성 여행 파동 데이터이다. 40 × 40 격자에서 생성된 파동 전선을 선형 측정 행렬 Φ와 잡음으로 관측했으며, EMD‑DF는 이전 프레임을 정확히 예측값으로 사용했을 때 거의 완벽한 복원을 달성했다. 반면 BPDN 및 BPDN‑DF는 예측 정보를 충분히 활용하지 못해 복원 품질이 크게 떨어졌다. 세 번째는 계산 효율성 평가이다. 일반적인 EMD와 Beckmann 기반 EMD‑DF를 비교했을 때, 상태 크기가 12 × 12(144)에서 48 × 48(2304)까지 증가함에 따라 Beckmann 방식은 실행 시간이 크게 감소했으며, 특히 대규모 문제에서 수십 배의 속도 향상을 보였다. 또한, 논문은 Sinkhorn 기반 최적 운송 정규화와의 차이점을 강조한다. Sinkhorn 방법은 엔트로피 정규화로 정확도와 속도 사이에 트레이드오프가 존재하고, 질량이 확산되어 희소성을 손상시킬 위험이 있다. 반면 Beckmann 기반 부분 운송은 정확한 최적 운송 거리를 유지하면서도 변수 차원을 O(N)으로 제한해, 희소 신호 추적에 적합한 정확도와 속도를 동시에 제공한다. 결론적으로, EMD‑DF는 희소 신호의 공간·주파수 구조를 보존하면서 실시간 추적이 가능한 효율적인 프레임워크를 제공한다. Beckmann 형태의 부분 운송 모델을 통해 대규모 문제에도 적용 가능하며, 신경과학을 비롯한 영상·주파수 분석 분야에서 기존 방법들을 능가하는 성능을 입증한다. 향후 연구에서는 비유클리드 거리, 비선형 동역학, 그리고 다중 타깃 상황에 대한 확장 가능성을 탐색할 수 있다.