Determining the Rolle function in Hermite interpolatory approximation by solving an appropriate differential equation

기술적 요약: 에르미트 보간 근사에서 롤 함수 결정

문제 제기
본 논문은 에르미트 보간 근사의 정확도를 향상시키는 과제를 다룹니다. 표준 에르미트 다항식 $H_{2n+1}(x)$는 $n+1$개의 노드에서 함수 $f(x)$와 그 도함수 $f'(x)$를 보간하지만, 점별 오차 $\Delta(x | H_{2n+1})$는 미지의 점 $\xi(x)$에서 평가된 $f$의 $(2n+2)$차 도함수를 포함하는 항에 의해 결정되며, 이를 롤 함수라고 합니다. 핵심 문제는 일반적으로 $\xi(x)$가 알려져 있지 않아 오차 항을 직접 계산하거나 보정할 수 없다는 점입니다. 저자는 $f(x)$의 더 정확한 근사를 가능하게 하기 위해 $\xi(x)$를 수치적으로 결정하는 방법을 모색합니다.

방법론
제안된 방법론은 라그랑주 보간에 적용되었던 기법을 에르미트 사례로 확장한 것입니다. 이 접근법은 다음 단계들을 포함합니다:

1. 미분 방정식 유도: 에르미트 보간의 표준 오차 공식에서 시작하여 저자는 오차 항을 $x$에 대해 미분합니다. 연쇄 법칙을 적용하고 롤 함수의 미분 $\frac{d\xi}{dx}$를 분리함으로써 1계 상미분 방정식 (ODE) 을 유도합니다. 이 ODE 는 $\frac{d\xi}{dx}$를 알려진 함수들 $f(x)$, $H_{2n+1}(x)$ 및 그들의 도함수와 미지의 $\xi(x)$와 연관시킵니다.
2. 초기값 문제 공식화: ODE 를 풀기 위해 보간 노드와 가깝지만 (노드 자체는 아님) 노드 근처의 점 $x_z$에서 초기값 $\xi_z = \xi(x_z)$가 필요합니다. 노드 자체에서는 오차가 0 이므로 $\xi$를 결정할 수 없습니다. 저자는 노드 근처의 점 $x_z$를 선택하고 대수적 오차 방정식을 수치적으로 풀어 유효한 초기 $\xi_z$를 찾습니다.
3. 수치적 해법: 유도된 ODE 는 7 차 룽게 - 쿠타 (RK) 방법을 사용하여 수치적으로 풉니다. 논문은 초기값에 대해 여러 수학적 해가 존재할 수 있음을 지적하며, 올바른 해는 $\xi(x)$가 보간 노드로 정의된 구간 ($x_0 < \xi(x) < x_n$) 내에 엄격히 위치해야 한다는 제약 조건에 기반하여 선택됩니다.
4. 오차 보정: $\xi(x)$가 결정되면 오차 공식의 항 $f^{(2n+2)}(\xi(x))$가 알려지게 됩니다. 저자는 이 항을 다항식 $H_\xi(x)$ (최소제곱 적합 또는 3 차 스플라인을 통해) 로 근사합니다. 이 다항식은 노드 인자 $\prod (x-x_k)^2$와 곱해지고 $(2n+2)!$로 나누어져 오차 다항식 $E(x)$를 형성합니다.
5. 정제된 근사: 최종적으로 개선된 근사는 오차 다항식을 원래 에르미트 다항식에 더함으로써 구성됩니다: $f(x) \approx H_{2n+1}(x) + E(x)$.

주요 결과
이 방법론은 노드 $\{0, 3\pi/2\}$ (여기서 $n=1$) 에 걸쳐 있는 함수 $f(x) = e^x \sin x$에 대해 테스트되었습니다.

• 롤 함수 결정: 수치 해법은 구간 제약을 만족하는 롤 함수 $\xi(x)$를 성공적으로 생성했습니다. ODE 의 수치 해는 오차 방정식의 좌변과 우변을 비교하여 검증되었으며, 두 변은 본질적으로 구별할 수 없음을 보여주었습니다.
• 정확도 향상:
  • 원래 에르미트 다항식 $H_3(x)$는 최대 근사 오차 7.04를 보였습니다.
  • 11 차 최소제곱 적합에서 유도된 오차 다항식 $E(x)$를 더함으로써 최대 오차는 $7.3 \times 10^{-8}$로 감소했습니다.
  • 도함수 항의 근사를 생성하기 위해 3 차 스플라인을 사용하면 최대 오차는 약 $10^{-12}$로 추가로 감소했습니다.
• 보간 성질 보존: 논문은 보정된 근사 $H_{2n+1}(x) + E(x)$가 원래 에르미트 보간 성질을 유지하여 노드에서 함수 값과 도함수 값 모두를 일치시킨다고 검증합니다.
• 수치 적분: 논문은 보정된 다항식을 구간에서 적분하면 원래 에르미트 다항식을 적분한 결과 ($\sim 0.7$ 오차) 에 비해 훨씬 더 정확한 결과 ($\sim 3 \times 10^{-12}$ 오차) 를 산출함을 간략히 시연합니다.

의의 및 주장
본 논문은 초기값 문제를 통해 롤 함수를 결정함으로써 에르미트 보간에서 점별 오차를 명시적으로 근사할 수 있다고 주장합니다. 주요 의의는 원래 보간 조건을 엄격히 유지하면서 표준 에르미트 다항식보다 수백만 배 더 정확한 새로운 다항식 근사를 구성할 수 있다는 점에 있습니다.

저자는 최종 결과의 정확도가 현재 $\xi(x)$를 풀기 위해 사용된 RK 방법의 정밀도에 의해 제한되며, 향후 연구는 RK 솔버 내의 오차 제어에 초점을 맞출 수 있음을 지적합니다. 논문은 라그랑주 보간에 대한 이전 작업에 기반한 이 기법이 수치 근사 품질을 향상시키는 데 상당한 잠재력을 보여준다고 결론지었습니다.