[논문 리뷰] On the Chebyshev approximation of a function with two variables

Ernest Scheiber|arXiv (Cornell University)|2015. 04. 18.

Advanced Computational Techniques in Science and Engineering참고 문헌 4인용 수 7

한 줄 요약

이 논문은 이변수 함수에 대한 균일 수렴을 보장하는 이차원 이산 푸리에 변환(2D-DFT)을 이용한 이변수 체비셰프 근사 다항식을 구성하는 방법을 제시한다. 이 방법은 체비셰프 급수 전개와 라그랑주-체비셰프 보간을 활용하여 테스트 함수에 대한 수치적 검증을 통해 높은 정확도를 달성하며, 표준 수학 소프트웨어에서 사용 가능한 간단한 구현을 통해 체브펑의 저질서 가우스 소거법 대안을 제공한다.

ABSTRACT

There is presented an approach to find an approximation polynomial of a function with two variables based on the two dimensional discrete Fourier transform. The approximation polynomial is expressed through Chebyshev polynomials. There is given an uniform convergence result.

연구 동기 및 목표

• 이변수 함수 근사에 대해 체브펑의 저질서 가우스 소거법에 대한 대안을 개발하는 것.
• 일반적인 수학 소프트웨어에서 이용 가능한 표준 2D-DFT 도구를 활용하여 이변수 함수의 고정밀 근사 구현을 가능하게 하는 것.
• 부드러운 조건 하에서 라그랑주-체비셰프 보간 다항식이 진짜 함수로의 균일 수렴을 보장하는 것.
• 체비셰프 근사 기반의 적분 및 편미분 계산을 위한 실용적인 수치 방법을 제공하는 것.
• 오차 및 적분 비교를 통한 수치 실험을 통해 방법의 타당성을 검증하는 것.

제안 방법

• 함수 값이 텐서 격자상에서 샘플링된 점들에 대해 2차원 이산 푸리에 변환(2D-DFT)을 이용하여 이변수 체비셰프 근사 다항식을 구성하는 방법.
• 체비셰프 계수는 체비셰프-레지온드르 노드에서 샘플된 함수 값의 복소수 2D-DFT를 통해 계산되며, 주기성과 트라프제졸드 근사법을 활용한다.
• 근사 다항식은 체비셰프 다항식의 곱 $ T_n(x)T_m(y) $ 의 이重합으로 표현되며, 계수는 DFT 출력에서 유도된다.
• 체비셰프 노드에서 함수와 일치하는 다항식을 구성하기 위해 라그랑주-체비셰프 보간을 사용하여 유일성과 수렴성을 보장한다.
• 함수가 연속적인 이阶 편미분을 가진다는 가정 하에 균일 수렴이 증명된다.
• 수치적 적분 및 편미분 계산은 체비셰프 급수 계수에서 유도된 닫힌 형태의 표현식을 통해 수행된다.

실험 결과

연구 질문

• RQ1 2D 이산 푸리에 변환은 이변수 함수 근사에 대한 체비셰프 계수를 효과적으로 계산하는 데 사용될 수 있는가?
• RQ2 부드러운 조건 하에서 도출된 라그랑주-체비셰프 보간 다항식이 원래 함수로의 균일 수렴을 보장하는가?
• RQ3 오차 및 계산 비용 측면에서 2D-DFT 기반 방법은 기존의 체브펑의 저질서 근사 방법과 비교해 어떤가?
• RQ4 표준 수학 소프트웨어에서 표준 DFT 루틴을 사용하여 이 방법을 효율적으로 구현할 수 있는가?
• RQ5 절단 및 반올림 오차는 절단된 체비셰프 급수를 사용한 수치적 적분 및 도함수 계산에 어떤 영향을 미치는가?

주요 결과

• 함수 $ f(x,y) = \cos(xy) $ 에 대해 $ 50 \times 50 $ 균일 격자에서 최대 절대 오차는 $ 8.2141 \times 10^{-5} $였으며, 잔차 지표는 $ 3.97247 \times 10^{-10} $였다.
• 함수 $ g(x,y) = \cos(10xy^2) + e^{-x^2} $ 에 대해 최대 절대 오차는 $ 2.98594 \times 10^{-13} $였고, 잔차 지표는 0이었으며, 이는 높은 수치 정밀도를 의미한다.
• 이 방법을 통해 계산한 $ f(x,y) = \cos(xy) $ 의 수치적 적분 결과는 $ 3.784330902 $였으며, 매트랩의 결과 $ 4\mathrm{Si}(1) \approx 3.78433228147 $ 와 매우 유사했다.
• 두 번째 테스트 함수에 대해 이 방법은 머신 정밀도를 달성하였으며, 적분 결과가 매트랩의 결과와 보고된 자릿수 내에서 정확히 일치하였다.
• 함수가 연속적인 이阶 편미분을 가진다는 조건 하에 다항식 $ L_{n,m} $ 이 진짜 함수 $ f $ 로의 균일 수렴이 보장된다.
• 이 방법은 체브펑의 가우스 소거법 기반 접근법에 대한 실용적인 대안을 제공하며, 표준 소프트웨어에서 더 단순한 구현이 가능하지만, 고정밀도일 경우 큰 메모리 요구량을 가질 수 있다.