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

Resumo Técnico: Determinação da Função de Rolle na Aproximação Interpolatória de Hermite

Enunciado do Problema
O artigo aborda o desafio de melhorar a precisão das aproximações interpolatórias de Hermite. Enquanto o polinômio de Hermite padrão $H_{2n+1}(x)$ interpola uma função $f(x)$ e sua derivada $f'(x)$ em $n+1$ nós, o erro pontual $\Delta(x | H_{2n+1})$ é governado por um termo envolvendo a derivada de ordem $(2n+2)$ de $f$ avaliada em um ponto desconhecido $\xi(x)$, conhecido como função de Rolle. O problema central é que $\xi(x)$ é geralmente desconhecido, impedindo o cálculo ou correção direta do termo de erro. O autor busca um método para determinar $\xi(x)$ numericamente, a fim de facilitar uma aproximação mais precisa de $f(x)$.

Metodologia
A metodologia proposta estende uma técnica previamente aplicada à interpolação de Lagrange para o caso de Hermite. A abordagem envolve as seguintes etapas:

1. Derivação de uma Equação Diferencial: A partir da fórmula de erro padrão para interpolação de Hermite, o autor diferencia o termo de erro em relação a $x$. Ao aplicar a regra da cadeia e isolar a derivada da função de Rolle, $\frac{d\xi}{dx}$, uma equação diferencial ordinária (EDO) de primeira ordem é derivada. Esta EDO relaciona $\frac{d\xi}{dx}$ às funções conhecidas $f(x)$, $H_{2n+1}(x)$ e suas derivadas, bem como ao desconhecido $\xi(x)$.
2. Formulação do Problema de Valor Inicial: Para resolver a EDO, é necessário um valor inicial $\xi_z = \xi(x_z)$ em um ponto $x_z$ próximo a (mas não em) um nó de interpolação. Como o erro é zero nos próprios nós, $\xi$ não pode ser determinado neles. O autor seleciona um ponto $x_z$ próximo a um nó e resolve numericamente a equação algébrica de erro para encontrar um $\xi_z$ inicial válido.
3. Solução Numérica: A EDO derivada é resolvida numericamente usando um método de Runge-Kutta (RK) de sétima ordem. O artigo observa que podem existir múltiplas soluções matemáticas para o valor inicial; a correta é selecionada com base na restrição de que $\xi(x)$ deve estar estritamente dentro do intervalo definido pelos nós de interpolação ($x_0 < \xi(x) < x_n$).
4. Correção do Erro: Uma vez que $\xi(x)$ é determinado, o termo $f^{(2n+2)}(\xi(x))$ na fórmula de erro torna-se conhecido. O autor aproxima este termo usando um polinômio $H_\xi(x)$ (via ajuste de mínimos quadrados ou splines cúbicos). Este polinômio é multiplicado pelo fator nodal $\prod (x-x_k)^2$ e dividido por $(2n+2)!$ para formar um polinômio de erro $E(x)$.
5. Aproximação Refinada: A aproximação final, aprimorada, é construída adicionando o polinômio de erro ao polinômio de Hermite original: $f(x) \approx H_{2n+1}(x) + E(x)$.

Principais Resultados
A metodologia foi testada na função $f(x) = e^x \sin x$ sobre os nós $\{0, 3\pi/2\}$ (onde $n=1$).

• Determinação da Função de Rolle: A solução numérica gerou com sucesso uma função de Rolle $\xi(x)$ que satisfaz as restrições de intervalo. A solução numérica da EDO foi verificada comparando os lados esquerdo e direito da equação de erro, mostrando-os essencialmente indistinguíveis.
• Melhoria da Precisão:
  • O polinômio de Hermite original $H_3(x)$ apresentou um erro máximo de aproximação de 7,04.
  • Ao adicionar um polinômio de erro $E(x)$ derivado de um ajuste de mínimos quadrados de grau 11, o erro máximo foi reduzido para $7,3 \times 10^{-8}$.
  • Usando um spline cúbico para gerar a aproximação do termo derivado, o erro máximo foi reduzido ainda mais para aproximadamente $10^{-12}$.
• Preservação das Propriedades Interpolatórias: O artigo verifica que a aproximação corrigida $H_{2n+1}(x) + E(x)$ mantém as propriedades interpolatórias originais de Hermite, correspondendo tanto aos valores da função quanto aos valores das derivadas nos nós.
• Integração Numérica: O artigo demonstra brevemente que integrar o polinômio corrigido sobre o intervalo produz um resultado significativamente mais preciso (erro de $\sim 3 \times 10^{-12}$) em comparação com a integração do polinômio de Hermite original (erro de $\sim 0,7$).

Significância e Alegações
O artigo alega que determinar a função de Rolle via um problema de valor inicial permite a aproximação explícita do erro pontual na interpolação de Hermite. A principal significância reside na capacidade de construir uma nova aproximação polinomial que é ordens de magnitude mais precisa do que o polinômio de Hermite padrão, mantendo estritamente as condições de interpolação originais.

O autor observa que a precisão do resultado final é atualmente limitada pela precisão do método RK usado para resolver $\xi(x)$, sugerindo que trabalhos futuros poderiam focar no controle de erro dentro do solver RK. O artigo conclui que esta técnica, construindo sobre trabalhos anteriores em interpolação de Lagrange, demonstra potencial substancial para melhorar a qualidade da aproximação numérica.