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

Sintesi Tecnica: Determinazione della Funzione di Rolle nell'Approssimazione Interpolatoria di Hermite

Enunciazione del Problema
Il lavoro affronta la sfida di migliorare la precisione delle approssimazioni interpolatorie di Hermite. Mentre il polinomio di Hermite standard $H_{2n+1}(x)$ interpola una funzione $f(x)$ e la sua derivata $f'(x)$ in $n+1$ nodi, l'errore puntuale $\Delta(x | H_{2n+1})$ è governato da un termine che coinvolge la derivata $(2n+2)$-esima di $f$ valutata in un punto incognito $\xi(x)$, noto come funzione di Rolle. Il problema fondamentale è che $\xi(x)$ è generalmente sconosciuto, impedendo il calcolo o la correzione diretta del termine di errore. L'autore cerca un metodo per determinare $\xi(x)$ numericamente per facilitare un'approssimazione più accurata di $f(x)$.

Metodologia
La metodologia proposta estende una tecnica precedentemente applicata all'interpolazione di Lagrange al caso di Hermite. L'approccio prevede i seguenti passaggi:

1. Derivazione di un'Equazione Differenziale: Partendo dalla formula standard dell'errore per l'interpolazione di Hermite, l'autore deriva il termine di errore rispetto a $x$. Applicando la regola della catena e isolando la derivata della funzione di Rolle, $\frac{d\xi}{dx}$, viene derivata un'equazione differenziale ordinaria (ODE) del primo ordine. Questa ODE mette in relazione $\frac{d\xi}{dx}$ con le funzioni note $f(x)$, $H_{2n+1}(x)$ e le loro derivate, nonché con l'incognita $\xi(x)$.
2. Formulazione del Problema ai Valori Iniziali: Per risolvere l'ODE, è richiesto un valore iniziale $\xi_z = \xi(x_z)$ in un punto $x_z$ vicino (ma non coincidente) a un nodo di interpolazione. Poiché l'errore è nullo negli stessi nodi, $\xi$ non può essere determinato in tali punti. L'autore seleziona un punto $x_z$ vicino a un nodo e risolve numericamente l'equazione algebrica dell'errore per trovare un $\xi_z$ iniziale valido.
3. Soluzione Numerica: L'ODE derivata viene risolta numericamente utilizzando un metodo di Runge-Kutta (RK) del settimo ordine. Il lavoro nota che possono esistere multiple soluzioni matematiche per il valore iniziale; quella corretta viene selezionata in base al vincolo che $\xi(x)$ deve giacere strettamente all'interno dell'intervallo definito dai nodi di interpolazione ($x_0 < \xi(x) < x_n$).
4. Correzione dell'Errore: Una volta determinato $\xi(x)$, il termine $f^{(2n+2)}(\xi(x))$ nella formula dell'errore diventa noto. L'autore approssima questo termine utilizzando un polinomio $H_\xi(x)$ (mediante adattamento ai minimi quadrati o spline cubiche). Questo polinomio viene moltiplicato per il fattore nodale $\prod (x-x_k)^2$ e diviso per $(2n+2)!$ per formare un polinomio di errore $E(x)$.
5. Approssimazione Raffinata: L'approssimazione finale, migliorata, viene costruita aggiungendo il polinomio di errore al polinomio di Hermite originale: $f(x) \approx H_{2n+1}(x) + E(x)$.

Risultati Chiave
La metodologia è stata testata sulla funzione $f(x) = e^x \sin x$ sui nodi $\{0, 3\pi/2\}$ (dove $n=1$).

• Determinazione della Funzione di Rolle: La soluzione numerica ha generato con successo una funzione di Rolle $\xi(x)$ che soddisfaceva i vincoli di intervallo. La soluzione numerica dell'ODE è stata verificata confrontando il lato sinistro e il lato destro dell'equazione dell'errore, dimostrando che erano essenzialmente indistinguibili.
• Miglioramento della Precisione:
  • Il polinomio di Hermite originale $H_3(x)$ ha mostrato un errore massimo di approssimazione di 7.04.
  • Aggiungendo un polinomio di errore $E(x)$ derivato da un adattamento ai minimi quadrati di grado 11, l'errore massimo è stato ridotto a $7.3 \times 10^{-8}$.
  • Utilizzando una spline cubica per generare l'approssimazione del termine derivato, l'errore massimo è stato ulteriormente ridotto a circa $10^{-12}$.
• Preservazione delle Proprietà Interpolatorie: Il lavoro verifica che l'approssimazione corretta $H_{2n+1}(x) + E(x)$ mantiene le proprietà interpolatorie originali di Hermite, corrispondendo sia ai valori della funzione che ai valori delle derivate nei nodi.
• Integrazione Numerica: Il lavoro dimostra brevemente che integrare il polinomio corretto sull'intervallo produce un risultato significativamente più accurato (errore $\sim 3 \times 10^{-12}$) rispetto all'integrazione del polinomio di Hermite originale (errore $\sim 0.7$).

Significato e Affermazioni
Il lavoro afferma che la determinazione della funzione di Rolle tramite un problema ai valori iniziali consente l'approssimazione esplicita dell'errore puntuale nell'interpolazione di Hermite. Il significato principale risiede nella capacità di costruire una nuova approssimazione polinomiale che è ordini di grandezza più accurata del polinomio di Hermite standard, mantenendo rigorosamente le condizioni di interpolazione originali.

L'autore nota che la precisione del risultato finale è attualmente limitata dalla precisione del metodo RK utilizzato per risolvere $\xi(x)$, suggerendo che il lavoro futuro potrebbe concentrarsi sul controllo dell'errore all'interno del solver RK. Il lavoro conclude che questa tecnica, basandosi su lavori precedenti nell'interpolazione di Lagrange, dimostra un potenziale sostanziale per migliorare la qualità dell'approssimazione numerica.