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}$ を分離することで、一階の常微分方程式（ODE）が導かれる。この ODE は、既知の関数 $f(x)$、$H_{2n+1}(x)$ およびそれらの導関数、ならびに未知の $\xi(x)$ と $\frac{d\xi}{dx}$ とを関連付ける。
2. 初期値問題の定式化：ODE を解くために、補間ノードに近く（ただしノードそのものではない）点 $x_z$ における初期値 $\xi_z = \xi(x_z)$ が必要である。ノード自体では誤差がゼロであるため、$\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)$（最小二乗法またはスプライン補間を用いて）で近似する。この多項式にノード因子 $\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}$ に減少した。
  • 導関数項の近似にスプライン補間を使用することで、最大誤差はさらに $10^{-12}$ 程度まで減少した。
• 補間性質の保持：論文は、補正された近似 $H_{2n+1}(x) + E(x)$ が元のエルミート補間性質を保持し、ノードにおける関数値と導関数値の両方に一致することを検証している。
• 数値積分：論文は、修正された多項式を区間上で積分すると、元のエルミート多項式を積分した場合（誤差 $\sim 0.7$）と比較して、はるかに正確な結果（誤差 $\sim 3 \times 10^{-12}$）が得られることを簡潔に示している。

意義と主張
本論文は、初期値問題を通じてロール関数を決定することで、エルミート補間における点ごとの誤差を明示的に近似できることを主張している。主な意義は、元の補間条件を厳密に維持しつつ、標準的なエルミート多項式よりも桁違いに正確な新しい多項式近似を構築できる点にある。

著者は、最終結果の精度が現時点では $\xi(x)$ を解くために使用された RK 法の精度によって制限されていることに言及しており、将来の研究は RK ソルバー内の誤差制御に焦点を当てることができると示唆している。論文は、この手法がラグランジュ補間に関する先行研究に立脚しており、数値近似の品質を向上させるための実質的な可能性を示していると結論付けている。