[논문 리뷰] Error Estimates with Explicit Constants for Sinc Quadrature and Sinc Indefinite Integration over Infinite Intervals

Tomoaki Okayama|arXiv (Cornell University)|2013. 02. 06.

Numerical Methods and Algorithms참고 문헌 3인용 수 4

한 줄 요약

이 논문은 반무한 및 무한 구간에서 Sinc 사다리꼴 적분과 Sinc 부정적분에 대해 이중지수(DE) 변환을 사용하여 명시적이고 계산 가능한 오차 한계를 제공한다. 새로운 DE 변환 $\psi_{\text{DE3$\ddag$}}(t) = \log(1 + \mathrm{e}^\pi \sinh t})$를 도입하여 수렴 속도를 향상시키고 명시적인 상수를 제공함으로써 이전에 다루지 못했던 케이스에 대해 검증 가능한 자동 적분을 가능하게 한다.

ABSTRACT

The Sinc quadrature and the Sinc indefinite integration are approximation formulas for definite integration and indefinite integration, respectively, which can be applied on any interval by using an appropriate variable transformation. Their convergence rates have been analyzed for typical cases including finite, semi-infinite, and infinite intervals. In addition, for verified automatic integration, more explicit error bounds that are computable have been recently given on a finite interval. In this paper, such explicit error bounds are given in the remaining cases on semi-infinite and infinite intervals.

연구 동기 및 목표

• 반무한 및 무한 구간에서 Sinc 사다리꼴 적분과 Sinc 부정적분에 대해 이전에는 존재하지 않았던 명시적이고 계산 가능한 오차 한계를 제공하는 것.
• 사례 3(반무한 구간에서 지수적으로 감소하는 피적분함수)의 수렴 속도와 오차 추정을 향상시키기 위해 새로운 DE 변환을 도입하는 것.
• 기존 변환보다 더 나은 오차 감소 속도를 달성하는 새로운 DE 변환 $\psi_{\text{DE3$\ddag$}}(t) = \log(1 + \mathrm{e}^{\pi \sinh t})$를 수립하는 것.
• 오차 추정에 명시적인 상수를 제공하여 검증 가능한 자동 적분을 가능하게 하는 것 — 이는 엄밀한 수치 계산에 필수적이다.

제안 방법

• 무한 및 반무한 구간에서 이중지수 변환을 사용한 Sinc 사다리꼴 적분에 대해 명시적 오차 한계를 유도한다.
• 사례 3를 위해 새로운 DE 변환 $\psi_{\text{DE3$\ddag$}}(t) = \log(1 + \mathrm{e}^{\pi \sinh t})$를 도입하여 역함수가 초월함수이지만 수렴성을 향상시킨다.
• 새로운 변환에 대해 $|\text{오차}| \leq C \mathrm{e}^{-2\pi d n / \log(4dn/\mu)}$ 형태의 오차 한계를 유도하고, 명시적인 $C$를 제공한다.
• 복소해석 기법을 적용하여 스트립 영역에서 변환된 피적분함수의 크기를 제한한다. 특히 해석적 계속 및 최대모듈러스 원리를 사용한다.
• 변환된 함수의 감소성과 해석성 특성을 기술하기 위해 함수 클래스 $\mathbf{L}_{L,R,2\alpha,2\beta}^{\text{DE}}(\mathscr{D}_d)$를 도입한다.
• 복소평면에서 로그함수와 쌍곡함수의 근사치를 정밀하게 분석하여 증명을 수행하며, 특히 주요 유리지수함수의 $\sup$ 노름을 활용한다.

실험 결과

연구 질문

• RQ1 반무한 및 무한 구간에서 Sinc 사다리꼴 적분에 대해 이전에 알려진 점근적 한계 외에 명시적이고 계산 가능한 오차 한계를 도출할 수 있는가?
• RQ2 새로운 DE 변환 $\psi_{\text{DE3$\ddag$}}(t) = \log(1 + \mathrm{e}^{\pi \sinh t})$가 사례 3에서 기존 DE 변환보다 수렴 속도가 빠른가?
• RQ3 새로운 변환에 대한 오차 추정 $|\text{오차}| \leq C \mathrm{e}^{-2\pi d n / \log(4dn/\mu)}$에서 명시적인 상수 $C$는 무엇인가?
• RQ4 새로운 변환이 $\psi_{\text{DE3}}$ 및 $\psi_{\text{DE3$\dagger$}}$보다 더 나은 수렴 속도를 갖는 것으로 증명될 수 있는가?
• RQ5 새로운 변환은 사례 3에 대해 Sugihara의 존재하지 않는 정리의 의미에서 최적인가?

주요 결과

• 새로운 DE 변환 $\psi_{\text{DE3$\ddag$}}(t) = \log(1 + \mathrm{e}^{\pi \sinh t})$는 $\psi_{\text{DE3}}$ 및 $\psi_{\text{DE3$\dagger$}}$보다 더 빠른 수렴 속도를 보이며, 오차 한계 $|\text{오차}| \leq C \mathrm{e}^{-2\pi d n / \log(4dn/\mu)}$를 갖는다.
• 오차 한계에 포함된 명시적인 상수 $C$가 도출되었고, 이를 통해 반무한 및 무한 구간에서 검증 가능한 자동 적분이 가능해진다.
• 새로운 변환의 오차 한계는 기존의 $\mathrm{e}^{-2\pi d n / \log(8dn/\mu)}$ 및 $\mathrm{e}^{-2\pi d n / \log(2\pi d n /\mu)}$보다 우월하다.
• 함수 $F(\zeta) = f(\psi_{\text{DE3$\ddag$}}(\zeta)) \psi_{\text{DE3$\ddag$}}'(\zeta)$는 $L = 2(\tilde{c}_d)^{1-\alpha}K$ 및 $R = 2(\mathrm{e}^{\pi/12})^{1-\alpha}K$를 갖는 클래스 $\mathbf{L}_{L,R,2\alpha,2\beta}^{\text{DE}}(\mathscr{D}_d)$에 속하며, 이는 복소 스트립 영역 내에서 해석성과 감소성을 보장한다.
• 증명은 복소 스트립 영역에서 표현식 $\left| \frac{1+\log(1+\mathrm{e}^{\pi \sinh \zeta})}{\log(1+\mathrm{e}^{\pi \sinh \zeta})} \cdot \frac{1}{1+\mathrm{e}^{-\pi \sinh \zeta}} \right| \leq \tilde{c}_d$ 의 정밀한 추정에 기반한다.
• 새로운 변환이 Sugihara의 존재하지 않는 정리의 의미에서 최적이며, 사례 3에 대해 이와 유사한 유형의 변환 중에서 가장 우수한 것으로 입증된다.