[논문 리뷰] On the backward stability of the second barycentric formula for interpolation

Walter F. Mascarenhas|arXiv (Cornell University)|2013. 10. 09.

Mathematical functions and polynomials인용 수 11

한 줄 요약

이 논문은 레바스 상수 Λx−,x+,x,w가 작을 조건 하에서 두 번째 바리센트 보간 공식의 역안정성을 입증한다. 노드와 가중치의 변동에 의한 역오차에 대해 엄밀한 상한과 하한을 제시하며, 수치적 평가가 약간 변형된 노드와 함수값에서의 정확한 평가와 동치로 유지됨을 보여주며, 오차 한계는 ϵn²log n (Salzer 가중치) 또는 ϵn log n (수치적 가중치) 비례로 증가한다.

ABSTRACT

We present a new stability analysis for the second barycentric formula for interpolation, showing that this formula is backward stable when the relevant Lebesgue constant is small.

연구 동기 및 목표

• 노드와 가중치에 대한 현실적인 변동 조건 하에서 두 번째 바리센트 보간 공식의 역안정성을 분석하기 위해.
• Higham(2004)이 악조건 하에서의 불안정성을 주장한 데 대해, 레바스 상수가 작을 경우 안정성이 유지됨을 보여주어 이의 모순을 해결하기 위해.
• 노드와 가중치의 변동을 동시에 고려하여 역오차에 대한 날카운 상한과 하한을 제공하기 위해.
• 실제 사례(예: 표준 가중치 계산 방식을 사용하는 체비셰프 노드)에서 공식이 여전히 역안정성을 유지함을 보여주기 위해.
• 편미분 이론 프레임워크를 통해 다항식 및 유리 함수 바리센트 보간의 분석을 통합하기 위해.

제안 방법

• 함수값 y에서 유리 보간을 생성하는 선형 연산자 Ix,w로 두 번째 바리센트 공식을 수식화하기 위해.
• Ix,w의 연산자 노름으로서 레바스 상수 Λx−,x+,x,w를 정의하여 오차의 최악의 경우 증폭 정도를 측정하기 위해.
• 노드의 변동에 대한 상대 오차 측정치 δjk, δ−j, δ+j 와 가중치의 변동에 대한 ζk 를 도입하여 수치적 부정확성을 정량화하기 위해.
• Higham(2002)의 오차 카운터 기법을 활용하여 노드 오차와 가중치 오차를 동시에 고려한 역오차 상한을 유도하기 위해.
• 편미분된 노드와 원래 노드 사이의 전단사 사상 χ를 정의하여 계산된 값을 근처 점에서의 정확한 보간값으로 매핑하기 위해.
• 역오차 분석 이론을 활용하여 계산된 값이 작은 상대 오차를 갖는 편미분된 노드와 함수값에서의 정확한 보간값과 동일함을 보여주기 위해.

실험 결과

연구 질문

• RQ1 노드와 가중치의 변동이 존재하는 상황에서 두 번째 바리센트 공식이 어떤 조건 하에 역안정성을 유지하는가?
• RQ2 노드와 가중치의 상대 오차가 동시에 보간 공식의 역오차에 어떤 영향을 미치는가?
• RQ3 Higham(2004)이 악조건 하에서의 불안정성을 주장했지만, 실제론 체비셰프 노드에서 공식이 안정성을 유지하는 이유는 무엇인가?
• RQ4 제2종 체비셰프 점을 사용한 다항식 보간에 대해 역오차의 날카운 상한은 무엇인가?
• RQ5 다른 가중치 계산 전략(Salzer의 닫힌 형태 vs. 수치적 평가)이 역오차에 어떤 영향을 미치는가?

주요 결과

• 레바스 상수 Λx−,x+,x,w 가 작을 경우, 노드와 가중치의 변동이 존재하더라도 두 번째 바리센트 공식은 역안정성을 유지한다.
• 제2종 체비셰프 노드의 경우, Salzer 가중치를 사용할 경우 역오차는 O(ϵn²log n) 이하로 제한되고, 수치적 가중치를 사용할 경우 O(ϵn log n) 이하로 제한된다.
• Section 3에서 제시된 하한을 통해, 라그랑주 보간의 경우 역오차 상한이 log n 요소를 제외하고 날카로움을 입증한다.
• 노드와 가중치의 변동을 동시에 고려한 분석을 통해, 한 가지를 무시할 경우 실제 역오차를 과소평가하게 된다.
• 계산된 값 fl(q(ˆx; ˆx, y, ˆw)) 는 근처 점 x 와 편미분된 함수값 ˜yk = yk(1 + αk) 에서의 정확한 보간값 q(x; x, ˜y, w) 와 동일하며, ∥α∥∞ 는 δ 와 레바스 상수의 함수로 유계이다.
• 노드의 편미션은 max{∥x − ˆx∥∞, |x− − ˆx−|, |x+ − ˆx+|} 이하로 제한되어, 계산된 점이 원래 점에 가까이 유지됨을 보장한다.