위상 오류가 있는 음장 추정에 최적 수송 중심 활용

초록

본 논문은 마이크 배열에서 발생하는 센서 위치·응답 불확실성으로 인한 위상 오류를 고려한 평면파 계수 추정 문제를 최적 수송(OT) 바리센터 프레임워크로 해결한다. 복소 계수를 비음수 측정으로 리프팅하고, 위상 교정 과정을 OT 거리 최소화 문제로 전환해 볼록 최적화 형태의 추정식을 도출한다. γ 파라미터를 도입해 전송 행렬의 희소성을 보장하고, 시뮬레이션을 통해 기존 LASSO·Tikhonov 등과 비교해 노이즈와 위상 교란이 존재할 때 NMSE가 현저히 낮은 것을 확인하였다.

상세 분석

이 연구는 기존의 평면파 기반 음장 재구성 방법이 센서 캘리브레이션 오류, 특히 위상 변동에 대해 편향된 추정값을 내거나 비볼록 최적화 문제에 빠지는 한계를 정확히 짚어낸다. 저자들은 복소 평면파 계수 αℓ 를 절대값과 위상으로 분리하고, 위상 교란 Δ(q)ℓ 를 각 마이크 q마다 독립적인 위상 이동으로 모델링한다. 여기서 핵심 아이디어는 복소 계수를 단순히 스칼라값으로 다루는 대신, 단위 원 위에 정의된 비음수 측정 μ(q)ℓ 로 ‘리프팅’하는 것이다. 각 μ(q)ℓ 은 Dirac 델타 형태로 질량을 αℓ 의 위상에서 Δ(q)ℓ 만큼 이동시킨 위치에 집중시킨다. 이렇게 하면 위상 교란은 측정의 지원(support) 이동으로 해석될 수 있다.

OT 프레임워크를 도입해 두 측정 μ와 ν 사이의 거리 ˜T(μ,ν)를 최소화하는 전송 계획 m을 구한다. 여기서 사용된 지상 비용 c(ψ1,ψ2)=‖e^{iψ1}−e^{iψ2}‖²+γ는 두 위상 간 유클리드 거리의 제곱에 상수 γ를 더한 형태이며, γ>0는 전체 질량에 대한 페널티 역할을 한다. γ가 존재하면 최적 전송 행렬이 희소해져 각 μℓ이 단일 Dirac 질량을 유지하도록 강제한다—즉, 원래의 평면파 모델이 유지된다.

바리센터 μ⁽⁰⁾는 Q개의 측정 μ(q)와의 평균 OT 거리를 최소화하는 측정으로 정의된다. 이는 위상 교란을 평균화한 ‘위상 평균’으로 해석될 수 있다. 최종 추정값 ˆα는 μ⁽⁰⁾의 첫 번째 푸리에 계수를 통해 복원된다. 전체 추정 문제는

min_{μ⁽⁰⁾,μ(q)} Σ_q T(μ⁽⁰⁾,μ(q)) + η‖g(q),∫e^{iψ}dμ(q)(ψ)−p̃(q)‖²

라는 형태의 볼록 최적화로 변환된다. 여기서 η는 데이터 적합도와 OT 비용 간의 트레이드오프를 조절한다.

실험에서는 2D 평면에 3개의 랜덤 평면파를 9개의 마이크로폰 배열에 배치하고, 위상 교란 σ_Δ와 잡음 SNR=15dB를 변동시켜 비교하였다. 제안 방법은 LASSO, LAD-LASSO, Tikhonov 등 전통적인 스파스 회귀 기법에 비해 NMSE가 2~5dB 정도 개선되었으며, 특히 위상 교란이 커질수록 그 차이가 두드러졌다. 또한 마이크 수와 주파수가 변해도 안정적인 복원을 보여, 실제 음향 제어·오디오 재현 시스템에 적용 가능함을 시사한다.

이 논문의 주요 공헌은 (1) 위상 오류를 측정의 지원 이동으로 모델링한 리프팅 기법, (2) OT 바리센터를 이용해 비볼록 문제를 볼록화한 추정 프레임워크, (3) γ 파라미터를 통한 희소성 보장 및 물리적 해석 가능성, (4) 시뮬레이션을 통한 실효성 검증이다. 향후 연구에서는 다중 주파수·3차원 배열, 실시간 구현을 위한 알고리즘 가속화, 그리고 비선형 위상 교란 모델 확장 등을 탐색할 여지가 있다.