[논문 리뷰] Approximate construction of rational approximations and the effect of error autocorrection. Applications

Grigori Litvinov|ArXiv.org|2001. 01. 05.

Numerical Methods and Algorithms참고 문헌 18인용 수 16

한 줄 요약

이 논문은 유리 근사 방법에서 계수 오차가 분자와 분모에서 상쇄되어, 선형 시스템의 조건수(condition number)가 높음에도 불구하고 근사 정확도를 유지하는 오차 자동보정 현상(autocorrection)을 조사한다. 연구는 심각한 계수 변동(예: 상대 오차 10⁻⁴)이 있음에도 불구하고, cos(x), arctan(x), 기타 초월함수의 근사에서 상대 오차가 최소 0.46×10⁻¹³에 이르게 되는 높은 정확도를 유지함을 보여준다.

ABSTRACT

Several construction methods for rational approximations to functions of one real variable are described in the present paper; the computational results that characterize the comparative accuracy of these methods are presented; an effect of error autocorrection is considered. This effect occurs in efficient methods of rational approximation (e.g., Pade approximations, linear and nonlinear Pade-Chebyshev approximations) where very significant errors in the coefficients do not affect the accuracy of the approximation. The matter of import is that the errors in the numerator and the denominator of a fractional rational approximant compensate each other. This effect is related to the fact that the errors in the coefficients of a rational approximant are not distributed in an arbitrary way but form the coefficients of a new approximant to the approximated function. Understanding of the error autocorrection mechanism allows to decrease this error by varying the approximation procedure depending on the form of the approximant. Some applications are described in the paper. In particular, a method of implementation of basic calculations on decimal computers that uses the technique of rational approximations is described in the Appendix. To a considerable extent the paper is a survey and the exposition is as elementary as possible.

연구 동기 및 목표

• 유리근사가 계수 계산 과정에서 상당한 반올림 오차가 발생하더라도 정확도를 유지하는 이유를 이해하는 것.
• Padé 및 Padé-Chebyshev 근사에서 오차 자동상관의 메커니즘을 조사하는 것.
• 표준적인 불안정 문제 기법(예: 정규화)이 이 맥락에서 정확도를 떨어뜨리고 계산 비용을 증가시킨다는 것을 보여주는 것.
• 오차 자동보정이 추가 변환 없이도 근사 알고리즘의 단순화와 정당화를 가능하게 한다는 것을 보여주는 것.
• 한정된 정밀도를 가진 십진수 컴퓨터에서 정확한 유리근사를 구현하기 위한 실용적 방법을 제공하는 것.

제안 방법

• 실수 구간에서 함수의 유리근사식 R(x) = P_n(x)/Q_m(x)의 계수를 계산하기 위해 선형 대수 시스템을 사용한다.
• 편차에 민감한 악조건의 선형 시스템을 해결하는 방법을 적용한다.
• ICL–4–50 및 ES–1045 컴퓨터에서 컴퓨터 실험을 수행하여, 다양한 반올림 오차를 가진 플랫폼 간 근사식을 비교한다.
• 계수 편차의 구조를 분석하고, 이들이 동일한 함수에 대한 새로운 유리근사식을 형성함을 보여준다.
• 오차 자동보정을 활용하여 구조적 절차를 수정하고 안정성을 향상시키는 방법을 제안한다.
• 부록에서 십진수 컴퓨터에서 사용 가능한 Fortran 알고리즘 구현을 기술한다.

실험 결과

연구 질문

• RQ1 왜 계수 계산에 상당한 상대 오차가 존재함에도 불구하고 유리근사식이 높은 정확도를 유지하는가?
• RQ2 유리근사식의 분자와 분모에서 발생하는 오차가 전체 근사 품질 유지를 위해 어떻게 상호작용하는가?
• RQ3 왜 표준 정규화 기법이 이 맥락에서 정확도 향상에 실패하며, 계산 비용에 어떤 영향을 미치는가?
• RQ4 오차 자동보정 효과를 체계적으로 활용하여 유리근사 알고리즘의 단순화 또는 향상이 가능한가?
• RQ5 오차 자동보정은 계수 계산이 수치적으로 불안정한 알고리즘을 얼마나까지 정당화할 수 있는가?

주요 결과

• [-π/4, π/4]에서 cos(x)에 대해 m=4, n=6인 유리근사식은 조건수 10⁹과 계수 편차가 6–9번째 소수자리 수준임에도 불구하고 상대 오차 0.55×10⁻¹³을 달성한다.
• [-1,1]에서 arctan(x)에 대해 ICL–4–50 및 ES–1045 컴퓨터 간에 동일한 상대 오차 0.5×10⁻¹¹를 유지하며, 계수 편차가 10⁻⁴이고 조건수는 10¹¹이다.
• 편차된 계수들이 동일한 함수에 대한 새로운 근사식을 형성함에 따라 근사식의 오차는 계수 편차에 영향을 받지 않으며, 이는 오차 상쇄를 초래한다.
• 표준 구간 산술 및 정규화 방법은 이 효과를 고려하지 못하며, 정확도를 떨어뜨리고 최대 2/3의 유효 자릿수를 잃을 수 있다.
• 오차 자동보정 현상은 Padé 근사에 국한되지 않으며, 선형 및 비선형 Padé–Chebyshev 근사로까지 확장된다.
• 이 논문은 오차 자동보정 덕분에 계수 계산이 불안정하더라도 안정적이고 고정밀한 근사가 가능함을 보여주며, 이는 저정밀도 십진수 컴퓨터에서도 실용적으로 활용 가능하다는 것을 입증한다.