[논문 리뷰] A stability barrier for reconstructions from Fourier samples

Ben Adcock, Anders C. Hansen외 1명|arXiv (Cornell University)|2012. 10. 29.

Mathematical functions and polynomials참고 문헌 44인용 수 4

한 줄 요약

이 논문은 분석적이고 비주기적인 함수를 그 첫 $ m $개의 푸리에 계수로부터 복원하는 데 있어 기본적인 안정성 장벽을 확립한다: 어떤 안정적인 방법이라도 $ m $에 대해 최대 루트지수적 수렴 속도를 달성할 수 있으며, 지수적 수렴은 지수적 불안정성(ill-conditioning)을 암시한다. 저자들은 다항함수의 푸리에 급수를 제한하는 데 있어 최대 행동에 대한 새로운 추정을 통해 이를 증명하며, 기하급수적으로 수렴하여 기계 정밀도에 도달하는 푸리에 확장이 이 장벽을 피할 수 있는 안정적이고 초급수적인 대안임을 보여준다.

ABSTRACT

We prove that any stable method for resolving the Gibbs phenomenon - that is, recovering high-order accuracy from the first $m$ Fourier coefficients of an analytic and nonperiodic function - can converge at best root-exponentially fast in $m$. Any method with faster convergence must also be unstable, and in particular, exponential convergence implies exponential ill-conditioning. This result is analogous to a recent theorem of Platte, Trefethen & Kuijlaars concerning recovery from pointwise function values on an equispaced $m$-grid. The main step in our proof is an estimate for the maximal behaviour of a polynomial of degree $n$ with bounded $m$-term Fourier series, which is related to a conjecture of Hrycak & Groechenig. In the second part of the paper we discuss the implications of our main theorem to polynomial-based interpolation and least-squares approaches for overcoming the Gibbs phenomenon. Finally, we consider the use of so-called Fourier extensions as an attractive alternative for this problem. We present numerical results demonstrating rapid convergence in a stable manner.

연구 동기 및 목표

• 푸리에 계수로부터 갈버그 현상을 해소하는 데 있어 안정적인 방법의 수렴 속도에 대한 기본적인 한계를 설정하는 것.
• 기하급수적 수렴을 보이는 기존 방법들이 왜 불안정한지, 그리고 왜 안정적인 방법들이 일반적으로 느린 수렴 속도를 보이는지 설명하는 것.
• 다항함수 기반 방법들인 역다항함수 복원법과 최소제곱법의 이론적 한계를 분석하는 것.
• 분석적 함수를 푸리에 데이터로부터 복원하는 데 있어 안정적이고 고차수의 대안으로 푸리에 확장을 제안하고 평가하는 것.

제안 방법

• 다항함수의 차수 $ n $에 대해 유계 $ m $-항 푸리에 급수를 가진 다항함수의 최대 행동에 대한 새로운 추정을 유도하며, 코퍼스마이스와 리블린의 결과를 푸리에 설정으로 확장한다.
• 이 추정을 이용해 $ \rho > 1 $, $ \tau \in (\frac{1}{2}, 1] $ 인 경우 수렴 속도가 $ \rho^{-m^\tau} $인 모든 방법은 불안정성의 정도가 $ \rho^{m^{2\tau - 1}} $이 되며, 이는 루트지수적 장벽을 확립한다.
• 역다항함수 복원법(IPRM)을 분석하고, 지수적으로 불안정하다는 것을 보이며, 최소제곱법을 통해 안정화하면 수렴 속도가 루트지수적 수준으로 떨어진다는 것을 밝힌다.
• 푸리에 확장을 안정적인 대안으로 제안하며, 이는 함수를 더 큰 간격으로 확장하고 그 푸리에 급수를 계산하는 방식으로 근사치를 구성한다.
• 다양한 시험 함수들(다양한 특이성을 가진)에 대해 다항최소제곱법과 푸리에 확장을 비교하기 위한 수치 실험을 수행한다.
• 확장 간격의 길이를 제어하는 데 사용되는 매개변수 $ T $ 를 도입하여, 푸리에 확장 방법에서 조건수와 수렴 속도 사이의 균형을 맞춘다.

실험 결과

연구 질문

• RQ1 분석적이고 비주기적인 함수를 그 첫 $ m $개의 푸리에 계수로부터 복원하는 데 있어 안정적인 방법이 달성할 수 있는 최선의 수렴 속도는 무엇인가?
• RQ2 실제로 기하급수적 수렴을 보이는 방법들이 왜 종종 심각한 불안정성을 보이는가?
• RQ3 다항함수 기반 최소제곱법이 안정성을 유지하면서도 루트지수적 수렴을 초월하는 수렴 속도를 달성할 수 있는가?
• RQ4 푸리에 확장이 기계 정밀도에 도달할 수 있도록 함으로써 안정성 장벽을 어느 정도 회피할 수 있는가?
• RQ5 경계층, 실수 또는 복소수 특이성을 가진 함수와 같은 다양한 함수 유형은 복원 방법의 성능에 어떤 영향을 미치는가?

주요 결과

• 푸리에 계수 $ m $개로부터 갈버그 현상을 해소하는 데 있어 안정적인 방법은 $ m $에 대해 최대 루트지수적 수렴 속도를 달성할 수 있으며, 즉 어떤 $ \rho > 1 $에 대해 $ \rho^{-m^{1/2}} $의 속도로 수렴한다.
• 지수적 수렴은 반드시 지수적 불안정성(ill-conditioning)을 암시하므로, 이러한 방법들은 수치적으로 불안정하다.
• 역다항함수 복원법(IPRM)이 지수적으로 불안정하다는 것이 입증되었으며, 고정밀도 계산에는 실용적이지 않다.
• IPRM를 최소제곱법으로 안정화하면 수렴 속도가 루트지수적 수준으로 떨어지며, 이는 이론적 안정성 장벽과 일치한다.
• 푸리에 확장은 기하급수적(종종 초급수적) 수렴을 기계 정밀도까지 보이며, 복잡한 특이성이나 고주파 진동을 가진 함수에서 다항최소제곱법을 능가한다.
• 실제로 확장 매개변수 $ T $ 의 선택이 오차에 미치는 영향은 미미하며, 더 많은 자유도가 증가함에 따라 각 자유도당 수렴 속도가 느려지는 것이 상쇄되기 때문이다.