DIFERENCIACIÓN NUMÉRICA

Elaboró: Norma Elva Espino Rojas

Si la aproximación es polinomial y con el criterio de ajuste exacto, la diferen-

ciación numérica consiste simplemente en diferenciar la fórmula del polinomio
interpolante que se utilizo. Sea en general
f (x) = pn (x) + Rn (x)
Es decir
dn f (x) dn pn (x)
dxn dxn

1 LAGRANGE:
Un polinomio de grado n tiene la forma
n
X n
Y x xj
pn (x) = f (xi )
i=0 j=0
xi xj
i6=j

Si es un polinomio lineal, se requieren dos puntos (x0 ; f (x0 )) ; (x1 ; f (x1 ))

p (x) = a0 (x x1 ) + a1 (x x0 )
x x1 x x0
p (x) = f (x0 ) + f (x1 )
x0 x1 x1 x0
Calculamos la derivada
f (x0 ) f (x1 )
p0 (x) = +
x0 x1 x1 x0
Pero si el polinomio es de grado dos
(x x1 ) (x x2 ) (x x0 ) (x x2 ) (x x0 ) (x x1 )
p (x) = f (x0 ) + f (x1 ) + f (x2 )
(x0 x1 ) (x0 x2 ) (x1 x0 ) (x1 x2 ) (x2 x0 ) (x2 x1 )
x2 (x1 + x2 ) x + x1 x2 x2 (x0 + x2 ) x + x0 x2
= f (x0 ) + f (x1 )
(x0 x1 ) (x0 x2 ) (x1 x0 ) (x1 x2 )
2
x (x0 + x1 ) x + x0 x1
+f (x2 )
(x2 x0 ) (x2 x1 )
La primer derivada del polinomio cuadrático es:
2x (x1 + x2 ) 2x (x0 + x2 ) 2x (x0 + x1 )
p0 (x) = f (x0 ) +f (x1 ) +f (x2 )
(x0 x1 ) (x0 x2 ) (x1 x0 ) (x1 x2 ) (x2 x0 ) (x2 x1 )
La segunda derivada
2f (x0 ) 2f (x1 ) 2f (x2 )
p00 (x) = + +
(x0 x1 ) (x0 x2 ) (x1 x0 ) (x1 x2 ) (x2 x0 ) (x2 x1 )

1
Ejemplo:
xi 1 0 1
f (xi ) 2 3 5
Calcular la primera y segunda deriva con los datos.
Solución: Como son tres pares ordenados en los datos, el polinomio será de
segundo grado. Sustituimos en las fórmulas encontradas anteriormente.

2x (x1 + x2 ) 2x (x0 + x2 ) 2x (x0 + x1 )

p0 (x) = f (x0 ) + f (x1 ) + f (x2 )
(x0 x1 ) (x0 x2 ) (x1 x0 ) (x1 x2 ) (x2 x0 ) (x2 x1 )
2x (0 + 1) 2x ( 1 + 1) 2x ( 1 + 0)
= (2) + (3) + (5)
( 1 0) ( 1 1) (0 ( 1)) (0 1) (1 ( 1)) (1 0)
2 (2x 1) (3) (2x) (5) (2x + 1)
= + +
2 ( 1) 2
5
= 2x 1 6x + 5x +
2
0 3
p (x) = x +
2
La segunda derivada:

2f (x0 ) 2f (x1 ) 2f (x2 )

p00 (x) = + +
(x0 x1 ) (x0 x2 ) (x1 x0 ) (x1 x2 ) (x2 x0 ) (x2 x1 )
2 (2) 2 (3) 2 (5)
p00 (x) = + +
( 1 0) ( 1 1) (0 ( 1)) (0 1) (1 ( 1)) (1 0)
4 6 10
= + +
2 ( 1) 2
= 2 6+5
p00 (x) = 1

¿Cuál es la derivada en x = 1:5?

3
p0 (1:5) = 1:5 +
2
= 3

2 NEWTON DIFERENCIAS DIVIDIDAS

Su polinomio tiene la forma

p (x) = b0 + b1 (x x0 ) + b2 (x x0 ) (x x1 ) + b3 (x x0 ) (x x1 ) (x x2 )
+ + bn (x x0 ) (x x1 ) (x xn 1)

2
donde

b0 = f (x0 )
b1 = f [x0 ; x1 ]
b2 = f [x0 ; x1 ; x2 ]
..
.
bn = f [x0 ; x1 ; : : : ; xn ]

Supongamos que nuestro polinomio es lineal

p (x) = b0 + b1 (x x0 )
= f (x0 ) + f [x0 ; x1 ] (x x0 )

Nuestra primera derivada es:

p0 (x) = f [x0 ; x1 ]

Si el polinomio es cuadrático

p (x) = b0 + b1 (x x0 ) + b2 (x x0 ) (x x1 )
p (x) = f (x0 ) + f [x0 ; x1 ] (x x0 ) + f [x0 ; x1 ; x2 ] (x x0 ) (x x1 )
p (x) = f (x0 ) + f [x0 ; x1 ] (x x0 ) + f [x0 ; x1 ; x2 ] x2 (x0 + x1 ) x + x0 x1

Su primer derivada es:

p0 (x) = f [x0 ; x1 ] + f [x0 ; x1 ; x2 ] [2x (x0 + x1 )]

Y su segunda derivada es:

p00 (x) = 2f [x0 ; x1 ; x2 ]

Ejemplo: Utilizamos la tabla del ejercicio de Lagrange

xi 1 0 1
f (xi ) 2 3 5

Solución: Necesitamos construir la tabla de diferencias divididas

xi f (xi ) Primera Diferencia Segunda Diferencia
1 2
f (x1 ) f (x0 ) 3 2
0 3 x1 x0 = f [x0 ; x1 ] = 0 ( 1) =1
f (x2 ) f (x1 ) 5 3 f [x1 ;x2 ] f [x0 ;x1 ] 2 1 1
1 5 x2 x1 = f [x1 ; x2 ] = 1 0 =2 x2 x0 = f [x0 ; x1 ; x2 ] = 1 ( 1) = 2

Para la primera derivada con el método de Newton Diferencias Divididas

p0 (x) = f [x0 ; x1 ]

3
 Usando la tabla entonces es:
p0 (x) = 1
Usando un polinomio de segundo grado sus primera es:
p0 (x) =
f [x0 ; x1 ] + f [x0 ; x1 ; x2 ] [2x (x0 + x1 )]
1
= 1 + (2x ( 1 + 0))
2
1
= 1 + (2x + 1)
2
0 3
p (x) = x +
2
Su segunda derivada:
p00 (x) =
2f [x0 ; x1 ; x2 ]
1
= 2
2
00
p (x) = 1
¿Cuál es la derivada en x = 1:5?
3
p0 (1:5) = 1:5 +
2
= 3

3 NEWTON DIFERENCIAS FINITAS

El polinomio de grado n con diferencias progresiva es:
42 f (x0 )
Pn (x) = 40 f (x0 ) + 4f (x0 ) s + s (s 1) +
2!
n
4 f (x0 )
+ s (s 1) s (s n + 1)
n!
donde
x x0
s =
h
x0 xn
h =
n
Recordemos que dichos polinomios sólo pueden ser aplicados si tienen sus vari-
ables independientes el mismo tamaño de paso, es decir, tienen un tamaño de h
de…nido.
Sea polinomio de grado uno
P (x) = 40 f (x0 ) + 4f (x0 ) s
x x0
P (x) = 40 f (x0 ) + 4f (x0 )
h

4
Calculemos su primer derivada

4f (x0 )
p0 (x) =
h
Si el polinomio es cuadrático

42 f (x0 )
P2 (x) = 40 f (x0 ) + 4f (x0 ) s + s (s 1)
2!
x x0 42 f (x0 ) x x0 x x0
P2 (x) = 40 f (x0 ) + 4f (x0 ) + 1
h 2! h h
x x0 42 f (x0 ) x x0 x x0 h
P2 (x) = 40 f (x0 ) + 4f (x0 ) +
h 2! h h
2
x x0 4 f (x0 ) 2
= 40 f (x0 ) + 4f (x0 ) + x 2x0 x xh + x20 + x0 h
h 2!h2

Calculamos su primer derivada:

4f (x0 ) 42 f (x0 )
P20 (x) = + [2x (2x0 + h)]
h 2!h2
Su segunda derivada es:

2 42 f (x0 ) 42 f (x0 )
P200 (x) = =
2!h2 h2
Ejemplo: Usamos la tabla de Lagrange, tienen el mismo tamaño de paso los
datos independientes
xi 1 0 1
f (xi ) 2 3 5

Solución: Construimos la tabla de diferencias

xi f (xi ) Primera Diferencia Segunda Diferencia
1 2
0 3 f (x1 ) f (x0 ) = f [x0 ; x1 ] = 3 2=1

1 5 f (x2 ) f (x1 ) = f [x1 ; x2 ] = 5 3=2 f [x1 ; x2 ] f [x0 ; x1 ] = f [x0 ; x1 ; x2 ] = 2 1=1

Calculamos la primera y segunda derivada con el polinomio de grado dos

con h = 1
NOTA: El valor de h en este ejercicio se calcula h = 1 (2 1) ; y elPvalor de n
es la cantidad de datos menos uno, es decir, tres datos, n = cantidad
de datos 1 = 3 1 = 2

5
 4f (x0 ) 42 f (x0 )
P20 (x) = + [2x (2x0 + h)]
h 2!h2
1 1
= + [2x (2 ( 1) + 1)]
1 2! (1)2
1
= 1 + (2x + 1)
2
3
P20 (x) = x+
2
42 f (x0 ) 1
P200 (x) = 2
= 2 =1
h (1)