From an Elementary Proof of Error Representation for Hermite Quadrature to a Rediscovery of Legendre Polynomials and Rodrigues Formula

Título: De uma Prova Elementar da Representação de Erro para Quadratura de Hermite até a Redescoberta de Polinômios de Legendre e a Fórmula de Rodrigues

Autores: Tan Bui-Thanh, Giancarlo Villatoro e C G Krishnanunni.

---

1. O Problema

A integração numérica é fundamental para a computação científica. Métodos clássicos, como as fórmulas de Newton-Cotes (ex: regra do trapézio, regra de Simpson), aproximam integrais usando apenas os valores da função em pontos selecionados. No entanto, em muitos cenários de engenharia e física, não apenas os valores da função $f(x)$, mas também as suas derivadas até a ordem $n-1$ estão disponíveis nos pontos extremos do intervalo $[a, b]$.

A extensão natural para aproveitar essa informação é a quadratura de Hermite, que utiliza interpolação de Hermite para aproximar a função antes de integrar. O problema central abordado no artigo é a representação e limitação do erro dessa quadratura.

• Limitação dos Métodos Clássicos: As expressões de erro tradicionais para quadratura de Hermite exigem a existência da derivada de ordem $2n$ da função ($f^{(2n)}$). Isso restringe a aplicabilidade a funções muito regulares (suaves).
• Objetivo: Desenvolver uma representação de erro exata que exija apenas a existência da derivada de ordem $n$ ($f^{(n)}$), permitindo a aplicação a uma classe mais ampla de funções (menos regulares), e explorar as conexões teóricas inesperadas que surgem dessa abordagem.

2. Metodologia

Os autores propõem uma abordagem baseada em integração por partes reversa (Reverse Integration by Parts - RIBP), uma técnica elementar que evita o uso de ferramentas mais complexas como diferenças divididas de Newton ou teoremas de valor médio integral para derivadas de alta ordem.

O processo metodológico segue os seguintes passos:

1. Formulação da Quadratura: Define-se a regra de quadratura de Hermite de dois pontos, integrando o polinômio interpolador $H_n(f; x)$ que coincide com $f$ e suas derivadas até ordem $n-1$ em $a$ e $b$.
2. Integração por Partes Reversa (RIBP): Aplica-se o processo de integração por partes $n$ vezes de forma recursiva. Isso gera uma fórmula que expressa a integral exata como uma soma de termos de fronteira (envolvendo derivadas de $f$ em $a$ e $b$) e um termo residual (integral) envolvendo a $n$-ésima derivada de $f$ multiplicada por um núcleo (kernel) polinomial $K_n(x)$.
3. Casamento de Coeficientes: Os parâmetros livres introduzidos na RIBP (constantes de integração arbitrárias) são ajustados para que os coeficientes das derivadas de fronteira na fórmula de RIBP coincidam exatamente com os pesos da quadratura de Hermite derivada anteriormente.
4. Identificação do Núcleo de Erro: Ao igualar os coeficientes, o núcleo residual da RIBP é identificado como o erro exato da quadratura.

3. Contribuições Chave e Resultados

A. Representação de Erro com Condições de Regularidade Mais Brandas

A principal contribuição é a derivação de uma fórmula de erro exata que requer apenas que a função pertença a $C^{(n)}[a, b]$ (possua derivada contínua de ordem $n$), em contraste com os métodos clássicos que exigem $C^{(2n)}[a, b]$.

• Fórmula do Erro: $E_n = \int_a^b (-1)^n f^{(n)}(x) K_n(x) \, dx$.
• Vantagem: Isso torna a fórmula válida para funções que possuem derivadas de ordem $n$, mas cujas derivadas de ordem superior (até $2n$) podem não existir ou não ser contínuas.

B. Redescoberta dos Polinômios de Legendre

Ao determinar os parâmetros que tornam o método consistente, os autores descobrem que o núcleo de erro $K_n(x)$ (o polinômio resultante da RIBP) é, a menos de uma constante e de um deslocamento, o Polinômio de Legendre de ordem $n$ (especificamente, o polinômio de Legendre deslocado e não normalizado).

• Relação Elegante: Estabelece-se que os polinômios de Legendre são precisamente os núcleos de erro para regras de quadratura interpolatória de Hermite.

C. Redescoberta da Fórmula de Rodrigues

Através da análise das propriedades do núcleo $K_n(x)$ e da aplicação da integração por partes $n$ vezes ao termo residual, os autores derivam a Fórmula de Rodrigues para os polinômios de Legendre:
$$K_n(x) = \frac{1}{(2n)!} \frac{d^n}{dx^n} \left[ (x-a)^n (x-b)^n \right]$$
Isso fornece uma nova perspectiva sobre a origem dessa fórmula clássica, ligando-a diretamente à teoria de erro de quadratura.

D. Limites de Erro (Bounds)

O artigo fornece limites rigorosos para o erro $E_n$ baseados em normas $L^\infty$ e $L^2$:

• Quando apenas $f^{(n)}$ é conhecida, o erro é limitado pela norma da derivada $n$-ésima vezes a norma do núcleo de Legendre.
• Quando derivadas de ordem superior ($f^{(n+k)}$) estão disponíveis, limites mais apertados são derivados utilizando propriedades de ortogonalidade e aproximação polinomial.

E. Prova Alternativa via Teorema do Núcleo de Peano

Para validar os resultados e cobrir o caso de funções com $2n$ derivadas, os autores apresentam uma prova alternativa utilizando o Teorema do Núcleo de Peano, confirmando que, sob condições de alta regularidade, o núcleo de Peano coincide com o núcleo derivado via RIBP.

F. Regra de Quadratura Composta

O artigo estende a regra de um único intervalo para uma regra composta sobre múltiplos subintervalos, útil para aplicações práticas onde a função pode variar significativamente.

4. Significado e Impacto

• Simplicidade e Elegância: A abordagem utiliza apenas cálculo elementar (integração por partes), tornando a derivação de erros de quadratura de alta ordem acessível e intuitiva, sem a necessidade de teorias de interpolação complexas.
• Generalização de Aplicabilidade: Ao reduzir a exigência de regularidade de $2n$ para $n$, o método permite a aplicação de quadratura de Hermite a uma classe muito mais ampla de funções, incluindo aquelas com singularidades ou descontinuidades em derivadas de ordem superior.
• Conexão Teórica Profunda: A redescoberta dos polinômios de Legendre e da fórmula de Rodrigues a partir de um problema de análise numérica (erro de quadratura) revela uma ligação fundamental e elegante entre interpolação, integração e polinômios ortogonais.
• Ferramenta Prática: A derivação de limites de erro explícitos e a regra composta fornecem ferramentas diretas para engenheiros e cientistas computacionais estimarem a precisão de suas integrações numéricas quando dados de derivadas estão disponíveis.

Em resumo, o artigo transforma um problema técnico de análise numérica em uma ferramenta pedagógica e teórica poderosa, unificando conceitos de quadratura, interpolação e polinômios clássicos sob uma única estrutura de prova elementar.