[논문 리뷰] Finite Difference Weights Using The Modified Lagrange Interpolant

Burhan Sadiq, Divakar Viswanath|arXiv (Cornell University)|2011. 02. 15.

Numerical methods for differential equations인용 수 1

한 줄 요약

이 논문은 수정된 라그랑주 보간법을 사용하여 유한차분 계수를 최적화한 알고리즘을 제시한다. 기존 포른버그의 방법에 비해 산술 연산을 $4/(5m+5)$ 배 줄여 계산 비용을 감소시키면서도 동일한 정확도를 유지한다. 이는 고차수 도함수(최대 $m=16$)의 안정적 계산을 가능하게 하며, 유한차분 공식이 초수렴 정확도를 달성할 수 있는 조건을 밝혀내며, 격점이 실수일 경우 차수를 1 이상으로 높이는 것은 불가능하다는 것을 증명한다.

ABSTRACT

Let $z_{1},z_{2},...,z_{N}$ be a sequence of distinct grid points. A finite difference formula approximates the $m$-th derivative $f^{(m)}(0)$ as $\sum w_{k}f(z_{k})$, with $w_{k}$ being the weights. We derive an algorithm for finding the weights $w_{k}$ which is an improvement of an algorithm of Fornberg (\emph{Mathematics of Computation}, vol. 51 (1988), p. 699-706). This algorithm uses fewer arithmetic operations than that of Fornberg by a factor of $4/(5m+5)$ while being equally accurate. The algorithm that we derive computes finite difference weights accurately even when $m$, the order of the derivative, is as high as 16. In addition, the algorithm generalizes easily to the efficient computation of spectral differentiation matrices. The order of accuracy of the finite difference formula for $f^{(m)}(0)$ with grid points $hz_{k}$, $1\leq k\leq N$, is typically $\mathcal{O}(h^{N-m})$. However, the most commonly used finite difference formulas have an order of accuracy that is higher than the typical. For instance, the centered difference approximation $(f(h)-2f(0)+f(-h))/h^{2}$ to $f(0)$ has an order of accuracy equal to 2 not 1. Even unsymmetric finite difference formulas can exhibit such superconvergence or boosted order of accuracy, as shown by the explicit algebraic condition that we derive. If the grid points are real, we prove a basic result stating that the order of accuracy can never be boosted by more than 1.

연구 동기 및 목표

• 고차수 도함수의 유한차분 계수를 정확하게 계산하기 위한 계산 비용을 줄인 더 효율적인 알고리즘을 개발하는 것.
• 유한차분 공식이 초수렴 정확도(일반적인 것보다 높은 정확도)를 달성할 수 있는 조건을 분석하는 것.
• 스펙트럼 미분 행렬의 효율적 계산을 위한 방법을 일반화하는 것.
• 실수 격점이 있는 유한차분 스킴에서 정확도 향상의 이론적 한계를 설정하는 것.

제안 방법

• 알고리즘은 수정된 라그랑주 보간법의 형태를 사용하여 포른버그의 접근 방식보다 계산 효율성을 향상시킨다.
• 서로 다른 격점 $z_k$에서의 다항식 보간을 활용하여 $f^{(m)}(0)$를 $\sum w_k f(z_k)$를 통해 근사하는 계수 $w_k$를 유도한다.
• 정확도를 유지하면서도 산술 연산을 $4/(5m+5)$ 배 줄여 고차수 도함수 계산에 특히 효과적이다.
• 이 방법은 편미분 스펙트럼 방법에 사용되는 스펙트럼 미분 행렬의 구축으로 자연스럽게 일반화된다.
• 유한차분 공식이 초수렴 정확도를 달성할 수 있는지 여부를 판단하기 위한 명시적인 대수적 조건을 도출한다.
• 실수 격점이 있는 경우, 일반적인 $\mathcal{O}(h^{N-m})$ 비율을 초과해 정확도를 높일 수는 없다는 이론적 증명을 제공한다.

실험 결과

연구 질문

• RQ1 고차수 도함수의 정확한 유한차분 계수 계산을 위한 최소한의 계산 비용은 무엇인가?
• RQ2 유한차분 공식이 초수렴 정확도를 달성할 수 있는 조건은 무엇인가?
• RQ3 이 알고리즘이 스펙트럼 미분 행렬의 효율적 계산으로 확장될 수 있는가?
• RQ4 실수 격점이 있는 유한차분 스킴에서 정확도 향상의 이론적 최대치는 무엇인가?
• RQ5 수정된 라그랑주 보간법은 산술 연산 수와 안정성 측면에서 포른버그의 알고리즘에 비해 어떻게 향상되는가?

주요 결과

• 제안된 알고리즘은 포른버그의 방법에 비해 산술 연산 수를 $4/(5m+5)$ 배 줄여 효율성이 크게 향상된다.
• 이 방법은 $m=16$에 이르는 고차수 도함수의 계수를 정확하게 계산할 수 있다.
• 유한차분 공식은 초수렴 정확도를 달성할 수 있으며, 본 논문은 이러한 경우를 식별할 수 있는 명시적인 대수적 조건을 도출한다.
• 실수 격점이 있는 경우, 일반적인 $\mathcal{O}(h^{N-m})$ 비율을 초과해 정확도를 높일 수는 없다.
• 이 알고리즘은 스펙트럼 미분 행렬의 효율적 계산으로 자연스럽게 일반화되어 편미분 스펙트럼 방법에 넓은 응용을 가능하게 한다.