[논문 리뷰] Aliasing error of the exp$(β\sqrt{1-z^2})$ kernel in the nonuniform fast Fourier transform

Alex H. Barnett|arXiv (Cornell University)|2020. 01. 26.

Numerical methods in inverse problems참고 문헌 34인용 수 13

한 줄 요약

이 논문은 비균일 빠른 푸리에 변환(NUFFT)에서 반원의 지수(ES) 커널의 앨리어징 오차를 철저히 분석하여, 커널 폭 $ w $ 에 대해 지수 수렴을 달성하고, 카이머–베셀 및 프로라트 구형파동함수(PSWF) 커널의 최적 수렴률과 游적하게 수렴함을 증명한다. 분석은 최강하강법, 윤곽선 적분 및 위상이 부여된 색인 합을 사용하여 타입 1 및 2 NUFFT의 오차를 근사하고, ES 커널이 신호 처리 및 영상 응용 분야에서 스펙트럼 정확도를 위해 거의 최적의 정확도를 갖는 고성능 대안으로서의 위치를 확립한다.

ABSTRACT

The most popular algorithm for the nonuniform fast Fourier transform (NUFFT) uses the dilation of a kernel $ϕ$ to spread (or interpolate) between given nonuniform points and a uniform upsampled grid, combined with an FFT and diagonal scaling (deconvolution) in frequency space. The high performance of the recent FINUFFT library is in part due to its use of a new "exponential of semicircle" kernel $ϕ(z)=e^{β\sqrt{1-z^2}}$, for $z\in[-1,1]$, zero otherwise, whose Fourier transform $\hatϕ$ is unknown analytically. We place this kernel on a rigorous footing by proving an aliasing error estimate which bounds the error of the one-dimensional NUFFT of types 1 and 2 in exact arithmetic. Asymptotically in the kernel width measured in upsampled grid points, the error is shown to decrease with an exponential rate arbitrarily close to that of the popular Kaiser--Bessel kernel. This requires controlling a conditionally-convergent sum over the tails of $\hatϕ$, using steepest descent, other classical estimates on contour integrals, and a phased sinc sum. We also draw new connections between the above kernel, Kaiser--Bessel, and prolate spheroidal wavefunctions of order zero, which all appear to share an optimal exponential convergence rate.

연구 동기 및 목표

• 비균일 빠른 푸리에 변환(NUFFT)에서 반원의 지수(ES) 커널을 사용할 경우 앨리어징 오차를 엄밀히 유계화하는 것.
• 기존의 커널들인 카이머–베셀 및 프로라트 구형파동함수(PSWF)와 비교하여 ES 커널의 游적 수렴률을 규명하는 것.
• FINUFFT 라이브러리에서 ES 커널의 고성능을 이론적으로 기반화하기 위해 커널 폭 $ w $ 에 따라 오차가 지수적으로 감소함을 증명하는 것.
• ES 커널, 카이머–베셀, PSWF 간의 관계를 탐색하며, 특히 그들의 푸리에 꼬리 행동과 수렴률 측면에서의 유사성을 분석하는 것.
• 지수 형태 $ e^{\beta\sqrt{1-z^2}} $ 가 NUFFT 맥락에서 최적의 $ \mathcal{O}(e^{-cw}) $ 수렴률을 갖는 모든 커널의 근본적 특성인지를 탐구하는 것.

제안 방법

• 타입 1 및 2 NUFFT에 대한 앨리어징 오차 유계를 도출하기 위해, $ |z| \leq 1 $ 에서 $ \phi_{\text{ES},\beta}(z) = e^{\beta(\sqrt{1-z^2}-1)} $ 로 정의된 ES 커널을 사용한다.
• 조건수 수렴하는 푸리에 변환 꼬리에 대한 합을 제어하기 위해 최강하강법과 고전적 윤곽선 적분 추정을 적용한다.
• 핵심 주파수 영역에서 커널의 푸리에 변환 행동을 분석하기 위해 위상이 부여된 색인 합을 사용한다.
• ES 커널의 푸리에 변환 꼬리 감쇠를 카이머–베셀 및 PSWF 커널과 비교하여, 지수 감쇠율 측면에서 游적 동등함을 보여준다.
• 엄밀한 점근적 분석을 통해 오차가 $ \mathcal{O}(e^{-c w}) $ 로 감소함을 증명하며, $ c = \pi\sqrt{1 - 1/\sigma} $ 로 주어지며, $ \sigma > 1 $ 일 때 최적 수렴률과 일치함을 보인다.
• 점근 전개를 통해 ES 커널과 PSWF 간의 연결 고리를 설정하며, 가우시안 극한 및 WKBJ 근사도 포함한다.

실험 결과

연구 질문

• RQ1 NUFFT에서 반원의 지수 커널이 이론적 최적 수렴률에 도달하는 앨리어징 오차 감쇠율을 갖는가?
• RQ2 ES 커널의 푸리에 변환 꼬리는 카이머–베셀 및 프로라트 구형파동함수 커널과 비교해 어떻게 다른가?
• RQ3 최강하강법과 윤곽선 적분 기법을 사용하여 ES 커널의 앨리어징 오차를 엄밀히 유계화할 수 있는가?
• RQ4 지수 형태 $ e^{\beta\sqrt{1-z^2}} $ 는 NUFFT 맥락에서 최적의 $ \mathcal{O}(e^{-cw}) $ 수렴률을 갖는 모든 커널의 특징인가?
• RQ5 다른 대수적 구조를 지닌다 해도, ES, KB, PSWF 커널 간의 경험적 성능 유사성은 무엇에 기인하는가?

주요 결과

• NUFFT에서 ES 커널의 앨리어징 오차는 커널 폭 $ w $ 에 따라 지수적으로 감소하며, 수렴률 $ c = \pi\sqrt{1 - 1/\sigma} $ 을 보이며, 카이머–베셀 및 PSWF 커널의 최적 수렴률과 일치한다.
• 오차 유계는 커널의 푸리에 변환 꼬리에 대한 조건수 수렴합에 최강하강법과 윤곽선 적분 기법을 적용하여 확립된다.
• ES 커널의 푸리에 변환은 $ e^{-\beta} $ 정도의 꼬리 감쇠를 보이며, 이는 카이머–베셀 및 PSWF 커널과 동일하여 고성능을 설명한다.
• 논문은 ES 커널이 이론적 최적에 근접한 점근적 수렴을 달성함을 증명하며, NUFFT 커널에 대해 추측된 최선의 수렴률과 일치함을 확인한다.
• 경험적 비교 결과, 하이브리드 형태 $ \psi_{\text{SlepH}}(z) \propto I_0(\beta\sqrt{1-z^2})(1 + \sqrt{1-z^2})^{-1/2} $ 는 KB 커널보다 PSWF를 더 정확히 근사하며, 상대 오차는 $ < 0.2/\beta $ 이다.
• PSWF의 푸리에 변환의 $ [-\beta, \beta] $ 대역 외부에서의 $ L^2 $-노름은 $ \sim 4\sqrt{\pi\beta} e^{-2\beta} $ 로 감쇠하며, 이는 $ e^{-\beta} $ 의 지수 감쇠율을 확인하며, ES 및 KB 커널과 일치한다.