Diferenciación Numérica

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 4

Diferenciación Numérica

Las fórmulas de derivación numérica aparecen en el desarrollo de algoritmos para la


solución de problemas de contorno en ecuaciones diferenciales ordinarias (y en ecuaciones
en derivadas parciales). En general, podemos obtener aproximaciones numéricas de la
derivada en un punto derivando alguna función interpolante, por ejemplo un polinomio de
Lagrange, algún trazador cúbico, etc. Sin embargo, en la práctica pequeños errores en los
datos pueden producir resultados no aproximados en las derivadas. Se puede llegar a
demostrar la fórmula de los puntos igualmente espaciados a partir de derivar un polinomio.

Figura1 Aproximación a la primera derivada (a) hacia adelante (b) hacia atrás (c)
central.

Aproximación a la primera derivada usando dos puntos usando polinomio interpolante de


grado 1 (diferencias finitas)

a) Hacia adelante

' f ( x i +h )−f ( xi )
f ( x i )=
h

b) Hacia atrás

' f ( x i )−f ( x i −h )
f ( x i )=
h

c) Central

' f ( x i +h )−f ( xi −h )
f ( x i )=
2h

Aproximación a la primera derivada usando tres puntos (polinomio interpolante de grado


2 - diferencias finitas)
El Polinomio de Newton igualmente espaciado:
2
Δf 0 Δ f0
p2 ( x )=f 0 + ( x− x0 )+ 2
( x− x0 )( x−x 1 )
1! h 2!h

( x−x 0 )
s=
Haciendo h

Δf 0 Δ 2 f 0
p2 ( x )=f 0 + s+ s (s−1)
1! 2!

'
f ( x i )=
dp 2 ( x )
ds
ds dp 2 ( x ) 1
dx = dx h = h
1
[ Δf 0 Δ2 f 0
1!
+
2!
(2 s−1 ) ]
si xi fi f 'i Error Fórmula de cinco
puntos central:

−f 2 +4 f 1 −3 f 0 2
0 x0 f0 h ' ''
f '0= f (ξ )
2h 3

1 x1 f1 −f 0 +f 2 2
f '1 =
h ' ''
f (ξ )
2h 6

2 x2 f2 3 f 2 −4 f 1 +f 0 2
f '2=
h ' ''
f (ξ )
2h 3

f i−2−8 f i−1 + 8 f i+1−f i+2 h4 ( 5)


f 'i = f (ξ )
12 h + 30

Ejercicio 1

Encontrar la pendiente de la recta que es tangente al gráfico de la función


f (x )=√ x 2 +2 x −√3 x en el punto x=3.85 completando la siguiente tabla. Registre los
cálculos usando 4 cifras decimales

h -0.02 -0.01 -0.001 0 0.001 0.01 0.02


x+h 3.83 3.84 3.849 3.85 3.851 3.86 3.87
f(x) 4.37 4.37 4.37 4.37 4.37 4.37 4.37
f(x+h) 4.3519 4.3625 4.3720 4.37 4.3742 4.3837 4.3943
f(x+h)-f(x) 1.0588 1.0587 1.0586 1.058 1.0586 1.0585 1.0584
h 6

Ejercicio 2
Usando tres puntos calcule la derivada en los siguientes puntos:
2
f (x )=x ln( x )+1

si xi fi f 'i Derivada Exacta Error

0 1 1.000 −f 2 +4 f 1 −3 f 0 1 0.0238
0 f '0=
2(0 . 2) =0.9762

1 1.2 1.262 '


−f 0 +f 2 1.6376 0.0111
5 f 1= =
2( 0. 2) 1.6487

2 1.4 1.659 3 f 2 −4 f 1 +f 0 2.3421 0.0208


5 f '2=
2(0 . 2) =2.3213

Ejercicio 3
Sea f (x )=−0.1 x 4 −0.15x 3 −0.5 x 2−0.25 x+1.2, aproxime la primera derivada en
x i=0 . 5, usando h=0. 25 . Para realizar esta tarea debe usar la fórmula de dos, tres y cinco
puntos.
Solución
x i=0 . 5 , x i−1 =x i−h=0 . 25 , x i−2 =x i−2 h=0, x i+1=x i +h=0 . 75 , x i+2=x i +2 h=1
Nodos: f (x i )=0 . 925 , f ( x i−1 )=1 .1035 , f ( xi−2 )=1. 2, f ( x i+1 )=0. 6363 , f ( x i+2 )=0 .2 .
Usando dos puntos:
' f ( x i +1 )−f ( x i )
f ( x i )≈ =−1 .155 , δ =26 . 5
Hacia adelante: h Valor exacto =
−0. 9125
' f ( x i )−f ( xi −1 )
f ( x i )≈ =−0 .714 , δ =21. 7
Hacia atrás: h
Nota: Con dos puntos no se puede usar la fórmula central, necesitamos tres puntos.
Usando tres puntos:
' −f ( x i+2 )+4 f ( x i+1 )−3 f ( x i )
f ( x i )≈ =−0 . 8594 , δ =5 .82
Hacia adelante: 2h

' 3 f ( x i )−4 f ( x i−1 )+ f ( x i−2 )


f ( x i )≈ =−0 .8781 , δ=3 .77
Hacia atrás: 2h
' −f ( x i−1 )+ f ( x i +1 )
f ( x i )≈ =−0 . 9344 , δ =2. 4
Central: 2h
Usando cinco puntos centrales
' f ( x i−2 )−8 f ( x i−1 )+ 8 f ( x i +1 )−f ( x i +2 )
f ( x i )≈ =−0 . 9125 ,δ =0
12 h

RMGJ

Lima, Mayo 2016

También podría gustarte