DIFERENCIACIÓN NUMERICA
CONCEPTO:
Técnica de análisis numérico para calcular una aproximación a la derivada de una
función en un punto utilizando los valores y propiedades de la misma.
APLICACIONES:
Este método se puede utilizar en dos casos:
Cuando estamos interesados en calcular el valor de alguna derivada en algún
punto prefijado, que a menudo se ha obtenido empíricamente.
Las fórmulas de derivación numérica se usan para obtener métodos numéricos
en la resolución de ecuaciones diferenciales ordinarias y en derivadas parciales.
FORMULAS:
PRIMERA DERIVADA:
PRIMERA DERIVADA
DIFERENCIAS 𝑓(𝑥0 + ℎ) − 𝑓(𝑥0 ) 0ℎ
𝑓′(𝑥0 ) =
FINITAS ℎ
PROGRESIVAS )
−3𝑓(𝑥0 + 4𝑓(𝑥0 + ℎ) − 𝑓(𝑥0 + 2ℎ) 0ℎ2
𝑓′(𝑥0 ) =
2ℎ
DIFERENCIAS 𝑓(𝑥0 + ℎ) − 𝑓(𝑥0 − ℎ) 0ℎ2
𝑓′(𝑥0 ) =
FINITAS 2ℎ
CENTRADAS 𝑓(𝑥0 − 2ℎ) − 8𝑓(𝑥0 − ℎ) + 8𝑓(𝑥0 + ℎ) − 𝑓(𝑥0 + 2ℎ) 0ℎ4
𝑓′(𝑥0 ) =
12ℎ
DIFERENCIAS 𝑓(𝑥0 ) − 𝑓(𝑥0 − ℎ) 0ℎ
𝑓′(𝑥0 ) =
FINITAS ℎ
REGRESIVAS 𝑓(𝑥0 − 2ℎ) − 4𝑓(𝑥0 − ℎ) + 3𝑓(𝑥0 ) 0ℎ2
𝑓′(𝑥0 ) =
2ℎ
SEGUNDA DERIVADA
SEGUNDA DERIVADA
DIFERENCIAS 𝑓(𝑥0 ) − 2𝑓(𝑥0 + ℎ) + 𝑓(𝑥0 + 2ℎ) 0ℎ
𝑓′′(𝑥0 ) =
FINITAS ℎ2
PROGRESIVAS 2𝑓(𝑥0 ) − 5𝑓(𝑥0 + ℎ) + 4𝑓(𝑥0 + 2ℎ) − 𝑓(𝑥0 + 3ℎ) 0ℎ2
𝑓′′(𝑥0 ) =
ℎ2
DIFERENCIAS 𝑓(𝑥0 − ℎ) − 2𝑓(𝑥0 ) + 𝑓(𝑥0 + ℎ) 0ℎ2
𝑓′′(𝑥0 ) =
FINITAS ℎ2
CENTRADAS −𝑓(𝑥0 − 2ℎ) + 16𝑓(𝑥0 − ℎ) − 30𝑓(𝑥0 ) + 16𝑓(𝑥0 + ℎ) − 𝑓(𝑥0 + 2ℎ) 0ℎ4
𝑓′′(𝑥0 ) =
12ℎ2
DIFERENCIAS 𝑓(𝑥0 − 2ℎ) − 2𝑓(𝑥0 − ℎ) + 𝑓(𝑥0 ) 0ℎ
𝑓′′(𝑥0 ) =
FINITAS ℎ2
REGRESIVAS −𝑓(𝑥0 − 3ℎ) + 4𝑓(𝑥0 − 2ℎ) − 5𝑓(𝑥0 − ℎ) + 2𝑓(𝑥0 ) 0ℎ2
𝑓′′(𝑥0 ) =
ℎ2
TERCERA DERIVADA
TERCERA DERIVADA
DIFERENCIAS −𝑓(𝑥0 ) + 3𝑓(𝑥0 + ℎ) − 3𝑓(𝑥0 + 2ℎ) + 𝑓(𝑥0 + 3ℎ) 0ℎ
𝑓′′′(𝑥0 ) =
FINITAS ℎ3
PROGRESIVAS −5𝑓(𝑥0 ) + 18𝑓(𝑥0 + ℎ) − 24𝑓(𝑥0 + 2ℎ) + 14𝑓(𝑥0 + 3ℎ) − 3𝑓(𝑥0 + 4ℎ) 0ℎ2
𝑓′′′(𝑥0 ) =
2ℎ3
DIFERENCIAS 0ℎ2
FINITAS −𝑓(𝑥0 − 2ℎ) + 2𝑓(𝑥0 − ℎ) − 2𝑓(𝑥0 + ℎ) − 𝑓(𝑥0 + 2ℎ)
𝑓′′′(𝑥0 ) =
CENTRADAS 8ℎ3
𝑓(𝑥0 − 3ℎ) − 8𝑓(𝑥0 − 2ℎ) + 13𝑓(𝑥0 − ℎ) − 13𝑓(𝑥0 + ℎ) + 8𝑓(𝑥0 + 2ℎ) − 𝑓(𝑥0 + 3ℎ) 0ℎ4
𝑓′′′(𝑥0 ) =
8ℎ3
DIFERENCIAS −𝑓(𝑥0 − 3ℎ) + 3𝑓(𝑥0 − 2ℎ) − 3𝑓(𝑥0 − ℎ) + 𝑓(𝑥0 ) 0ℎ
𝑓′′′(𝑥0 ) =
FINITAS ℎ3
REGRESIVAS 3𝑓(𝑥0 − 4ℎ) − 14𝑓(𝑥0 − 3ℎ) + 24𝑓(𝑥0 − 2ℎ) − 18𝑓(𝑥0 − ℎ) + 5𝑓(𝑥0 ) 0ℎ2
𝑓′′(𝑥0 ) =
2ℎ3
DEMOSTRACIÓN DE LA ALGUNAS FÓRMULAS:
A partir de la serie de Taylor.
ℎ 2 ′′ ℎ3
𝑓(𝑥 + ℎ) = 𝑓(𝑥) + ℎ𝑓 ′ (𝑥) + 2
𝑓 (𝑥) + 𝑓 ′′′ (𝑥) +
3!
⋯ Ecuación 1
ℎ2 ℎ3
ℎ𝑓 ′ (𝑥) = 𝑓(𝑥 + ℎ) − 𝑓(𝑥) − [ 𝑓 ′′ (𝑥) + 𝑓 ′′′ (𝑥) + ⋯ ]
2 3!
𝑓(𝑥 + ℎ) − 𝑓(𝑥) ℎ ℎ2
𝑓 ′ (𝑥) = − [ 𝑓 ′′ (𝑥) + 𝑓 ′′′ (𝑥) + ⋯ ]
ℎ 2 3!
Para 𝑥 = 𝑥0
𝑓(𝑥0 +ℎ)−𝑓(𝑥0 )
𝑓 ′ (𝑥0 ) ≈ PROGRESIVA
ℎ
Por Taylor:
ℎ 2 ′′ ℎ3
𝑓(𝑥 − ℎ) = 𝑓(𝑥) − ℎ𝑓 ′ (𝑥) + 2
𝑓 (𝑥) − 𝑓 ′′′ (𝑥) +
3!
⋯ Ecuación 2
ℎ2 ℎ3
ℎ𝑓 ′ (𝑥) = 𝑓(𝑥) − 𝑓(𝑥 − ℎ) + [ 𝑓 ′′ (𝑥) − 𝑓 ′′′ (𝑥) + ⋯ ]
2 3!
𝑓(𝑥) − 𝑓(𝑥 − ℎ) ℎ ℎ2
𝑓 ′ (𝑥) = + [ 𝑓 ′′ (𝑥) − 𝑓 ′′′ (𝑥) + ⋯ ]
ℎ 2 3!
Para 𝑥 = 𝑥0
𝑓(𝑥0 )−𝑓(𝑥0 +ℎ)
𝑓 ′ (𝑥0 ) ≈ REGRESIVA
ℎ
Restando la ecuación 1 con la 2:
ℎ 2 ′′ ℎ3
𝑓(𝑥 + ℎ) = 𝑓(𝑥) + ℎ𝑓 ′ (𝑥) + 2
𝑓 (𝑥) + 𝑓 ′′′ (𝑥) +
3!
⋯ Ecuación 1
ℎ 2 ′′ ℎ3
𝑓(𝑥 − ℎ) = 𝑓(𝑥) − ℎ𝑓 ′ (𝑥) + 2
𝑓 (𝑥) − 3! 𝑓 ′′′ (𝑥) + ⋯ Ecuación 2
2ℎ3
𝑓(𝑥 + ℎ) − 𝑓(𝑥 − ℎ) = 2ℎ𝑓 ′ (𝑥) + 𝑓 ′′′ (𝑥) + ⋯
3!
𝑓(𝑥+ℎ)−𝑓(𝑥−ℎ) ℎ2
= 𝑓′(𝑥) + 𝑓 ′′′ (𝑥) + ⋯
2ℎ 3!
𝑓(𝑥+ℎ)−𝑓(𝑥−ℎ)
𝑓′(𝑥) = + 0ℎ2
2ℎ
Para 𝑥 = 𝑥0
𝑓(𝑥0 +ℎ)−𝑓(𝑥0 −ℎ)
𝑓 ′ (𝑥0 ) ≈ CENTRAL
2ℎ
Considerando la formula progresiva de la primera derivada, se tiene la siguiente
fórmula para la segunda derivada:
𝑓(𝑥0 +ℎ)−𝑓(𝑥0 )
𝑓 ′ (𝑥0 ) ≈ PROGRESIVA
ℎ
𝑓′(𝑥0 +ℎ)−𝑓′(𝑥0 )
𝑓 ′′ (𝑥0 ) ≈
ℎ
𝑓(𝑥0 +2ℎ)−𝑓(𝑥0 +ℎ) 𝑓(𝑥0 +ℎ)−𝑓(𝑥0 )
−[ ]
′ ℎ ℎ
𝑓 ′(𝑥0 ) ≈
ℎ
𝑓(𝑥0 + 2ℎ) − 2𝑓(𝑥0 + ℎ) + 𝑓(𝑥0 )
𝑓 ′′ (𝑥0 ) =
ℎ2
Del mismo modo se pueden deducir las demás.
EJERCICIO:
Tenemos la siguiente ecuación:
𝑓(𝑥) = −0.1𝑥 4 − 0.15𝑥 3 − 0.5𝑥 2 − 0.25𝑥 + 1.2
Estimar la derivada en x=0.5 y un tamaño de paso: h= 0.25 usando diferenciación
numérica.
x f(x)
X0-2h 0 1.2
X0-h 0.25 1.103516
X0 0.5 0.925
X0+h 0.75 0.6363281
X0+2h 1 0.2
𝑓′(𝑥) = −0.4𝑥 3 − 0.45𝑥 2 − 1𝑥 1 − 0.25
𝑓′(0.5) = −0.4(0.5)3 − 0.45(0.5)2 − 1(0.5)1 − 0.25
𝑓 ′ (0.5) = −0.9125 Valor verdadero
En este caso aplicamos la fórmula de primera derivada centrada.
𝑓(𝑥0 − 2ℎ) − 8𝑓(𝑥0 − ℎ) + 8𝑓(𝑥0 + ℎ) − 𝑓(𝑥0 + 2ℎ)
𝑓′(𝑥0 ) =
12ℎ
1.2 − 8(1.103516) + 8(0.6363281) − 0.2
𝑓′(0.5) =
12(0.25)
𝑓 ′ (0.5) = −0.9125
𝑉𝑣 − 𝑉𝑎
𝐸𝑟𝑝 = | | 𝑥100%
𝑉𝑣
−.9125 − (−0.9125)
𝐸𝑟𝑝 = | | 𝑥100%
−0.9125
𝐸𝑟𝑝 = 0%
En el caso de segunda derivada seria lo siguiente:
𝑓(𝑥) = −0.1𝑥 4 − 0.15𝑥 3 − 0.5𝑥 2 − 0.25𝑥 + 1.2
𝑓′(𝑥) = −0.4𝑥 3 − 0.45𝑥 2 − 1𝑥 1 − 0.25
𝑓′′(𝑥) = −1.2𝑥 2 − 0.9𝑥 − 1
𝑓′′(0.5) = −1.2(0.5)2 − 0.9(0.5) − 1
𝑓 ′′ (0.5) = −1.75 Valor verdadero
En este caso aplicamos la fórmula de segunda derivada centrada.
x f(x)
X0-2h 0 1.2
X0-h 0.25 1.103516
X0 0.5 0.925
X0+h 0.75 0.6363281
X0+2h 1 0.2
−𝑓(𝑥0 − 2ℎ) + 16𝑓(𝑥0 − ℎ) − 30𝑓(𝑥0 ) + 16𝑓(𝑥0 + ℎ) − 𝑓(𝑥0 + 2ℎ)
𝑓′′(𝑥0 ) =
12ℎ2
−1.2 + 16(1.103516) − 30(0.925) + 16(0.6363281) − 0.2
𝑓′′(0.5) =
12(0.25)2
𝑓 ′′ (0.5) = −1.749992533
𝑉𝑣 − 𝑉𝑎
𝐸𝑟𝑝 = | | 𝑥100%
𝑉𝑣
−1.75 − (−1.749992533)
𝐸𝑟𝑝 = | | 𝑥100%
−1.75
𝐸𝑟𝑝 = 0,00044%
En el caso de tercera derivada seria lo siguiente:
𝑓(𝑥) = −0.1𝑥 4 − 0.15𝑥 3 − 0.5𝑥 2 − 0.25𝑥 + 1.2
𝑓′(𝑥) = −0.4𝑥 3 − 0.45𝑥 2 − 1𝑥 1 − 0.25
𝑓′′(𝑥) = −1.2𝑥 2 − 0.9𝑥 − 1
𝑓′′′(𝑥) = −2.4𝑥 − 0.9
𝑓′′′(0.5) = −2.4(0.5) − 0.9
𝑓 ′′′ (0.5) = −2.1 Valor verdadero
En este caso aplicamos la fórmula de tercera derivada centrada.
x f(x)
X0-3h -0.25 1.233203125
X0-2h 0 1.2
X0-h 0.25 1.103516
X0 0.5 0.925
X0+h 0.75 0.6363281
X0+2h 1 0.2
X0+3h 1.25 -0.430859375
𝑓(𝑥0 − 3ℎ) − 8𝑓(𝑥0 − 2ℎ) + 13𝑓(𝑥0 − ℎ) − 13𝑓(𝑥0 + ℎ) + 8𝑓(𝑥0 + 2ℎ) − 𝑓(𝑥0 + 3ℎ)
𝑓′′′(𝑥0 ) =
8ℎ3
1.233203125 − 8(1.2) + 13(1.103516) − 13(0.6363281) + 8(0.2) − (−0.430859375)
𝑓′′′(0.5) =
8(0.25)3
𝑓 ′′′ (0.5) = −2.0999584
𝑉𝑣 − 𝑉𝑎
𝐸𝑟𝑝 = | | 𝑥100%
𝑉𝑣
−2.1 − (−2.0999584)
𝐸𝑟𝑝 = | | 𝑥100%
−2.1
𝐸𝑟𝑝 = 0,0002%
