Coeficiente de diferencias finitas

Si ${\displaystyle f\left(x\right)}$ es una función continua y diferenciable, y ${\displaystyle h_{x}}$ es un intervalo finito alrededor de ${\displaystyle x_{0}}$, la aproximación de una derivada enésima de la función ${\displaystyle f}$ con un orden de precisión ${\displaystyle m}$, puede representarse de manera general por el método de diferencias finitas centrales como

${\displaystyle \displaystyle f^{(n)}\left(x_{0}\right)\approx \displaystyle {\frac {\sum _{i=-\left\lfloor {\frac {m+n-1}{2}}\right\rfloor }^{\left\lfloor {\frac {m+n-1}{2}}\right\rfloor }{c_{i}^{(n,m)}f\left(x_{0}+ih_{x}\right)}}{h_{x}^{n}}}+O\left(h_{x}^{m}\right)}$

donde los coeficientes ${\displaystyle c_{i}^{(n,m)}}$ son específicos para cada derivada dependiendo del orden de precisión. Además, ${\displaystyle \lfloor y\rfloor }$ implica que el valor está redondeado al valor entero menor que ${\displaystyle y}$.

La siguiente tabla contiene los respectivos coeficientes para el cálculo de derivadas de una función por diferencias finitas centrales, para varios órdenes de precisión:[1]

Por ejemplo, la tercera derivada con un orden de exactitude de segundo grado es

${\displaystyle \displaystyle f'''(x_{0})\approx \displaystyle {\frac {-{\frac {1}{2}}f(x_{0}-2h_{x})+f(x_{0}-h_{x})-f(x_{0}+h_{x})+{\frac {1}{2}}f(x_{0}+2h_{x})}{h_{x}^{3}}}+O\left(h_{x}^{2}\right)}$

donde ${\displaystyle h_{x}}$ representa el paso finito entre cada intervalo de diferencias finitas.

De manera análoga al método de diferencias finitas centrales, es posible escribir la aproximación de la enésima derivada por el método de diferencias finitas hacia adelante como

${\displaystyle \displaystyle f^{(n)}\left(x_{0}\right)\approx \displaystyle {\frac {\sum _{i=0}^{m+n-1}{a_{i}^{(n,m)}f\left(x_{0}+ih_{x}\right)}}{h_{x}^{n}}}+O\left(h_{x}^{m}\right)}$

donde ${\displaystyle a_{i}^{(n,m)}}$ son los coeficientes correspondientes para este método. Algunos de estos se muestran en la siguiente tabla :[1]

Por ejemplo, por el método de diferencias finitas hacia delante las derivadas primera y segunda con una precisión de tercer y segundo orden, respectivamente, son

${\displaystyle \displaystyle f'(x_{0})\approx \displaystyle {\frac {-{\frac {11}{6}}f(x_{0})+3f(x_{0}+h_{x})-{\frac {3}{2}}f(x_{0}+2h_{x})+{\frac {1}{3}}f(x_{0}+3h_{x})}{h_{x}}}+O\left(h_{x}^{3}\right)}$

${\displaystyle \displaystyle f''(x_{0})\approx \displaystyle {\frac {2f(x_{0})-5f(x_{0}+h_{x})+4f(x_{0}+2h_{x})-f(x_{0}+3h_{x})}{h_{x}^{2}}}+O\left(h_{x}^{2}\right)}$

Mientras que las mismas aproximaciones obtenidas por el método de diferencias finitas hacia atrás están dadas por

${\displaystyle \displaystyle f'(x_{0})\approx \displaystyle {\frac {{\frac {11}{6}}f(x_{0})-3f(x_{0}-h_{x})+{\frac {3}{2}}f(x_{0}-2h_{x})-{\frac {1}{3}}f(x_{0}-3h_{x})}{h_{x}}}+O\left(h_{x}^{3}\right)}$

${\displaystyle \displaystyle f''(x_{0})\approx \displaystyle {\frac {-2f(x_{0})+5f(x_{0}-h_{x})-4f(x_{0}-2h_{x})+f(x_{0}-3h_{x})}{h_{x}^{2}}}+O\left(h_{x}^{2}\right)}$

En general, conseguir los coeficientes de la aproximación hacia atrás es muy simple. Para las derivadas de orden par (${\displaystyle n}$ = 2, 4, 6...) son los mismos que para la aproximación hacia delante. Por otro lado, para las derivadas de orden impar (${\displaystyle n}$ = 1, 3, 5...) basta con cambiar el signo de los coeficientes listados en la tabla anterior.

La tabla siguiente ilustra esto de manera resumida:

• Método de diferencia finita
• Diferencia finita
• Sténcil de cinco puntos