[논문 리뷰] Simply Explicitly Invertible Approximations to 4 Decimals of Error Function and Normal Cumulative Distribution Function

Alessandro Soranzo, Emanuela Epure|arXiv (Cornell University)|2012. 01. 05.

Numerical Methods and Algorithms인용 수 12

한 줄 요약

이 논문은 오차 함수 erf(x)에 대한 새로운 명시적 역가능한 근사법을 제안하며, 모든 x ≥ 0 에서 절대 오차 <2.27×10⁻⁵ 이고 상대 오차 <1.21×10⁻⁴ 이며, Winitzki의 이전 방법보다 절대 오차를 약 5.5배 감소시켰다. 이 방법은 분수 함수 지수를 제곱근 표현식에 사용하여 이차방정식 해법을 통한 닫힌 형식의 역가능성을 가능하게 하며, 정규 누적분포함수 Φ(x), erfc(x), Q-함수에 대한 근사에도 확장 가능하며, 네 자리 소수 정밀도를 확보한다.

ABSTRACT

We improve the Modified Winitzki's Approximation of the error function $erf(x)\cong \sqrt{1-e^{-x^2\frac{\frac{4}π+0.147x^2}{1+0.147x^2}}}$ which has error $|\varepsilon (x)| &lt; 1.25 \cdot 10^{-4}$ $\forall x \ge 0$ till reaching 4 decimals of precision with $|\varepsilon (x)| &lt; 2.27 \cdot 10^{-5}$; also reducing slightly the relative error. Old formula and ours are both explicitly invertible, essentially solving a biquadratic equation, after obvious substitutions. Then we derive approximations to 4 decimals of normal cumulative distribution function $Φ(x)$, of erfc$(x)$ and of the $Q$ function (or cPhi).

연구 동기 및 목표

• Winitzki의 오차 함수 erf(x) 근사법의 정확도를 향상시켜 네 자리 소수 정밀도를 달성하기.
• 절대 오차를 1.25×10⁻⁴ 에서 2.27×10⁻⁵ 이하로 감소시키면서도 명시적 역가능성을 유지하기.
• 보조 오차 함수 erfc(x), 정규 누적분포함수 Φ(x), Q-함수(cPhi)에 대한 유사한 근사법을 네 자리 소수 정밀도로 유도하기.
• 모든 근사법이 대입 후 이차방정식으로 단순화되어 명시적 역가능성을 유지하도록 보장하기.
• 모든 도메인 x ≥ 0 에 걸쳐 오차 한계를 검증하며, x > 4 인 점 渐진적 행동도 포함하기.

제안 방법

• erf(x) 근사 향상을 위해 exp(-(a x² + b x⁴)/(1 + c x² + d x⁴)) 형태의 새로운 분수 함수 지수를 제안.
• erf(x) ≈ √(1 - exp(E(x))) 와 같은 변환을 사용하며, E(x) = (-1.2735457x² - 0.1487936x⁴)/(1 + 0.1480931x² + 0.0005160x⁴).
• 0 ≤ x ≤ 4 에서 분석적 비교를 통한 오차 한계 설정: Mathematica로 생성된 그래프 및 확대된 부분 도메인 사용.
• x > 4 인 경우, erf(x)와 근사값 η(x) 모두 1에 수렴하며, |1 - erf(x)| < 10⁻⁷ 이고 |1 - η(x)| < e⁻¹² 이므로, 병합 오차 < 10⁻⁵ 가 보장됨.
• Φ(x) 및 Q(x) 근사법을 erf(x) 공식에 스케일링을 적용하고 계수를 조정하여 유도: Φ(x) ≈ ½ + ½√(1 - exp(E'(x))), 여기서 E'(x)는 x⁴ 계수를 반으로 줄이고 분모 상수를 두 배로 조정.
• 모든 근사법이 대입 후 z = x² 를 통해 이차방정식 해법으로 단순화되어 명시적 역가능성을 유지하도록 보장.

실험 결과

연구 질문

• RQ1 Winitzki의 erf(x) 근사법은 네 자리 소수 정밀도를 달성하면서도 명시적 역가능성을 유지할 수 있는가?
• RQ2 분수 함수 지수의 계수 조정은 절대 오차를 5.5배 감소시키는 데 어떤 영향을 미치는가?
• RQ3 향상된 erf(x) 근사법은 체계적으로 Φ(x), erfc(x), Q(x)로 확장되어 유사한 정밀도를 달성할 수 있는가?
• RQ4 전체 도메인 x ≥ 0 에 걸쳐 오차를 엄격하게 한계를 설정하기 위한 분석적 및 수치적 기법은 무엇인가?
• RQ5 계산의 단순성과 역가능성을 유지하면서 상대 오차는 얼마나 줄일 수 있는가?

주요 결과

• 새로운 erf(x) 근사법은 모든 x ≥ 0 에서 절대 오차 <2.27×10⁻⁵ 를 달성하며, Winitzki의 1.25×10⁻⁴ 에 비해 5.5배 향상됨.
• 상대 오차는 <1.21×10⁻⁴ 으로 감소하여 도메인 전반에 걸쳐 일관된 향상이 이루어짐.
• Φ(x) 근사법은 절대 오차 <1.14×10⁻⁵ 를 달성하여 정규 누적분포함수의 네 자리 소수 정밀도를 가능하게 함.
• Q-함수 및 erfc(x) 근사법도 절대 오차 각각 <2.27×10⁻⁵ 와 <1.14×10⁻⁵ 를 달성하며, 유사한 정밀도 확보.
• 모든 근사법은 z = x² 대입 후 이차방정식 해법을 통해 명시적 역가능성이 보장됨.
• x > 4.125 인 영역에서는 erf(x) ≈ 1 이 새로운 공식보다 더 낮은 절대 오차를 보이지만, 새로운 공식은 여전히 전이 영역에서 유효하고 정확함.