En matemáticas, para aproximar la derivada de una función con un orden arbitrario de exactitud, es posible de utilizar el método de diferencias finitas. Este puede ser central, hacia delante o hacia atrás.

Diferencia finita central

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

f ( n ) ( x 0 ) i = m n 1 2 m n 1 2 c i ( n , m ) f ( x 0 i h x ) h x n O ( h x m ) {\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 c i ( n , m ) {\displaystyle c_{i}^{(n,m)}} son específicos para cada derivada dependiendo del orden de precisión. Además, y {\displaystyle \lfloor y\rfloor } implica que el valor está redondeado al valor entero menor que y {\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 exactitud de segundo grado es

f ( x 0 ) 1 2 f ( x 0 2 h x ) f ( x 0 h x ) f ( x 0 h x ) 1 2 f ( x 0 2 h x ) h x 3 O ( h x 2 ) {\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 h x {\displaystyle h_{x}} representa el paso finito entre cada intervalo de diferencias finitas.

Diferencia finitas hacia delante y atrás

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

f ( n ) ( x 0 ) i = 0 m n 1 a i ( n , m ) f ( x 0 i h x ) h x n O ( h x m ) {\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 a i ( n , m ) {\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

f ( x 0 ) 11 6 f ( x 0 ) 3 f ( x 0 h x ) 3 2 f ( x 0 2 h x ) 1 3 f ( x 0 3 h x ) h x O ( h x 3 ) {\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)}
f ( x 0 ) 2 f ( x 0 ) 5 f ( x 0 h x ) 4 f ( x 0 2 h x ) f ( x 0 3 h x ) h x 2 O ( h x 2 ) {\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

f ( x 0 ) 11 6 f ( x 0 ) 3 f ( x 0 h x ) 3 2 f ( x 0 2 h x ) 1 3 f ( x 0 3 h x ) h x O ( h x 3 ) {\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)}
f ( x 0 ) 2 f ( x 0 ) 5 f ( x 0 h x ) 4 f ( x 0 2 h x ) f ( x 0 3 h x ) h x 2 O ( h x 2 ) {\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 ( n {\displaystyle n} = 2, 4, 6...) son los mismos que para la aproximación hacia delante. Por otro lado, para las derivadas de orden impar ( n {\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:

Véase también

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

Referencias


Diferencias finitas Diferencia finita Iteración

El Metodo de Diferencias Finitas PDF Diferencia finita Linealidad

Diferencias finitas

Diferencias Finitas

Diferencias Finitas PDF Análisis numérico Calor