Extrapolación de Richardson

El método de extrapolación de Richardson, desarrollado por Lewis Fry Richardson (1881-1953), permite construir a partir de una secuencia convergente otra secuencia más rápidamente convergente. Esta técnica se usa frecuentemente para mejorar los resultados de métodos numéricos a partir de una estimación previa, de igual forma mejora la precisión en el cálculo numérico de la derivada de una función, partiendo de la base de la serie de Taylor. Este proceso es especialmente utilizado para definir un método de integración: el método de Romberg.

Para una función variable en x, la primera derivada está definida por:

${displaystyle f'(x)=lim _{h o 0}{frac {f(x+h)-f(x)}{h}}}$

Una simple aproximación se tiene por la diferencia hacia adelante, de forma que:

${displaystyle f_{1}(x)={frac {f(x+h)-f(x)}{h}}}$

Esta aproximación está lejos del valor real, por tanto para hacer un análisis del error, expandimos en forma de serie de Taylor:

${displaystyle f(x+h)=f(x)+hf'(x)+{frac {h^{2}}{2}}f''(x)+{frac {h^{3}}{3!}}f'''(x)+cdots }$

Substrayendo f(x) de ambos lados y dividiendo por h, se tiene que:

${displaystyle f_{1}(x)={frac {f(x+h)-f(x)}{h}}=f'(x)+{frac {h}{2}}f''(x)+{frac {h^{2}}{3!}}f'''(x)+cdots =f'(x)+O(h)}$

Análogamente se derivan las demás fórmulas de aproximación, deduciendo por ejemplo, con diferencia hacia atrás o cambiando los valores de h; de esta forma se obtiene una expresión generalizada llamada extrapolación de Richardson:

Sea A, la respuesta exacta a la integral, y A(h) la estimación de A con orden ${displaystyle h^{k_{0}}}$. De tal forma que:

Donde:

Ahora bien: Usando tamaños de espaciamiento h y h/t, podemos aproximar a A como:

Multiplicando la última ecuación por ${displaystyle t^{k_{1}}}$

Sustrayendo (2) y (1), como se vio al inicio:

Despejando A:

De este modo, se ha obtenido una mejor aproximación de A al sustraer el término más grande en el error, ${displaystyle O(h^{k_{1}})}$. De igual manera se pueden remover más términos de error de modo que se obtengan mejores aproximaciones de A. Una relación de recurrencia general puede ser implementada en las aproximaciones al hacer:

${displaystyle A_{i+1}(h)approx {frac {t^{k_{i}}A_{i}{Big (}{frac {h}{t}}{Big )}-A_{i}(h)}{(t^{k_{i}}-1)}},}$ siendo ${displaystyle k_{i},}$ el orden del error

con:

Las aplicaciones más inmediatas de la Extrapolación de Richardson en los métodos numéricos son dos: derivación numérica mediante diferencias centradas y las fórmulas de Newton-Cotes para la integración numérica.

Dada una tabla equiespaciada de datos, el procedimiento de mejora de un resultado obtenido mediante derivación por diferencias centradas es el siguiente:

Igualando el valor real de la derivada en las dos estimaciones se obtiene:

siendo ${displaystyle R}$ el error de truncamiento (resto), que depende del espaciado elevado a una determinada potencia:

y siendo ${displaystyle n}$ el orden de error del método utilizado

Sustituyendo ${displaystyle k}$ en la diferencia del valor más preciso se tiene una nueva estimación de la derivada:

Esta nueva estimación del valor tiene una diferencia de orden eⁿ⁺² y el error del método se puede calcular aproximadamente mediante el valor del término corrector en valor absoluto:

Siguiendo este procedimiento se puede llegar a un resultado mejorado para la estimación de la derivada y también una estimación de su error. Además, el procedimiento puede aplicarse sucesivamente para así obtener resultados cada vez más precisos.