Tarea 5

DIVISIÓN DE CIENCIAS BÁSICAS E INGENIERÍA (CBI)
DEPARTAMENTO DE PROCESOS E HIDRAÚLICA

ÁREA DE RECURSOS ENERGETICOS
LICENCIATURA EN INGENIERÍA EN ENERGÍA
UEA: MÉTODOS NUMÉRICOS APLICADOS A LA INGENIERÍA
ESTUDIANTE:
ALVARADO SALAZAR RAFAEL
GARCÍA HERNÁNDEZ GERARDO IVAN
MARTÍNEZ CAMPOS JUAN JOSÉ
NÁJERA CANELA ANA LILIA
PROFESOR: ALEJANDRO VÁZQUEZ RODRÍGUEZ

RESUMEN DE RESULTADOS
1.- Valor exacto de la a) -1.033663127
derivada b) tabla de datos elaborada en Excel
Aproximación por
diferencias finitas
2.- 8.41
∇𝑓̌(1,1) = ( )
11.38
3.- 13 17
− −
36 36
a)𝐻(𝑓(1,1)) = ( 17 25)
− −
36 36
b)
1) 2) 3)
-0.3595 -0.4695 -0.3627 -0.4750 -0.3611 -0.4722
-0.4702 -0.6903 -0.4742 -0.6987 -0.4722 -0.6945
4.- Matriz Jacobiana
diferencia finitas en el
punto 𝑥 𝑇 = (1,1)
1.- El objetivo de este problema es comparar la aproximación por diferencias finitas de primer
orden centrada de la primera derivada de una función con el valor real de la derivada. Esto se hará
para
𝑓(𝑥) = ⅇ −2𝑥 − 𝑥
a) Use el calculo para determinar el valor correcto de la derivada en x=2
Para nuestra primera derivada se resuelve como
Por la regla de la suma, la derivada de 𝑓(𝑥) = ⅇ −2𝑥 − 𝑥 respecto de x es

ⅆ −2𝑥 ⅆ
[ⅇ ]+ [−𝑥]
ⅆ𝑥 ⅆ𝑥
ⅆ
Evaluamos [ⅇ −2𝑥 ] con regla de la cadena
ⅆ𝑥
ⅆ ⅆ
ⅇ −2𝑥 [−2𝑥] + [−𝑥]
ⅆ𝑥 ⅆ𝑥
ⅆ ⅆ
ⅇ −2𝑥 (−2 [𝑥]) + [−𝑥]
ⅆ𝑥 ⅆ𝑥
ⅆ
ⅇ −2𝑥 (−2𝑥1) + [−𝑥]
ⅆ𝑥
ⅆ
ⅇ −2𝑥 (−2) + [−𝑥]
ⅆ𝑥
ⅆ
−2ⅇ −2𝑥 + [−𝑥]
ⅆ𝑥
ⅆ
En donde ⅆ𝑥 [−𝑥] = 𝑥
ⅆ
−2ⅇ −2𝑥 − [𝑥]
ⅆ𝑥
−2ⅇ −2𝑥 − 1𝑥1
−2ⅇ −2𝑥 − 1
−1 − 2ⅇ −2𝑥
Por lo tanto
𝑓 ′ (𝑥) = −1 − 2ⅇ −2𝑥
Si tenemos un valor real x=2
𝑓 ′ (2) = −1 − 2ⅇ −2(2)
= −1 − 2ⅇ 4
2
= −1 − = −1.03663127
ⅇ4
b) Para evaluar la aproximación por diferencias finitas use el valor de h=0.001, y determine el error
como el valor absoluto de la diferencia de la derivada exacta y la derivada numérica.
Esto lo haremos por medio de Excel.
Tenemos nuestra función
Evaluando x=2 obtenemos
Ahora bien, se nos menciona que tenemos h=0.001 y se lo sumamos a nuestra función
Posteriormente restamos nuestra h
Para poder evaluar la aproximación para finitas centrales
En primer lugar ocuparemos nuestra primera derivada de nuestra función la cual es:
En donde evaluaremos en nuestro punto x=2, obteniendo como resultado:

Ahora bien, para nuestra finita adelante(numerica1)
ocupamos la formula
Hacia atrás (numérica 2)
Y centrada (numérica 3)
Por lo que nos queda
Y por último sacamos el error restando cada valor a el valor original que -1.033663
2. Problema
Dada la función
𝑓(𝑥, 𝑦) = 3𝑥𝑦 2 + 2ⅇ 𝑥𝑦
El gradiente analítico es:
3𝑦 2 + 2ⅇ 𝑥𝑦
∇𝑓(𝑥, 𝑦) = ( )
6𝑥𝑦 + 2𝑥ⅇ 𝑥𝑦
Evaluamos en 𝑋 𝑇 = (1,1) y obtenemos:
8.4365
∇𝑓(1,1) = ( )
11.4365
El 1% de variación de las variables está dando por:
𝛿1 = 𝛿2 = 0.01
Por lo tanto, usamos la siguiente ecuación:
𝑓(𝑥) − 𝑓(𝑥 − 𝛿𝑖 ⅇ𝑖 )
𝑔̃𝑖 =
𝛿𝑖
Entonces:
1
𝑥=( )
1
1 1 0.99
𝑥 − 𝛿1 ⅇ1 = ( ) − 0.01 ( ) = ( )
1 0 1.00
1 0 1.00
𝑥 − 𝛿2 ⅇ2 = ( ) − 0.01 ( ) = ( )
1 1 0.99
𝑓(1) = 3(1)(1)2 + 2ⅇ (1)(1) = 8.4365
𝑓(𝑥 − 𝛿1 ⅇ1 ) = 3(0.99)(1)2 + 2ⅇ (0.99)(1) = 8.3524

𝑓(𝑥 − 𝛿2 ⅇ2 ) = 3(1)(0.99)2 + 2ⅇ (1)(0.99) = 8.3227
𝑓(𝑥) − 𝑓(𝑥 − 𝛿1 ⅇ1 ) 8.4365 − 8.3524
𝑔
̃1 = = = 8.41
𝛿1 0.01
𝑓(𝑥) − 𝑓(𝑥 − 𝛿2 ⅇ2 ) 8.4365 − 8.3227
𝑔
̃2 = = = 11.38
𝛿2 0.01
Entonces el gradiente numérico que obtenemos es:
8.41
∇𝑓̌(1,1) = ( )
11.38
El error es la diferencia de entre gradientes:
ⅇ𝑟𝑟𝑜𝑟 = |∇𝑓(𝑥) − ∇𝑓̌(𝑥)| = 0.0538
Se nota que hay una mejor precisión en la primera que en la segunda.

3. Encuentre el valor de la matriz Hessiana numérica (por diferencias finitas
hacia adelante) en el punto xT = (1,1) para la función siguiente
𝑓(𝑥, 𝑦) = 𝐿𝑛(𝑥 2 + 3𝑥𝑦 + 2𝑦 2 )
Usando el 1% de variación en las variables, compare con el valor exacto de la
matriz hessiana analítica.
Resolviendo el problema de manera analítica se tiene que
𝜕2𝑓 𝜕2𝑓
𝜕𝑥 2 𝜕𝑥𝜕𝑦
𝐻(𝑓(𝑥, 𝑦)) =
𝜕2𝑓 𝜕2𝑓
(𝜕𝑦𝜕𝑥 𝜕𝑦 2 )
Obteniendo las derivadas parciales
𝜕𝑓 2𝑥 + 3𝑦
= 2
𝜕𝑥 𝑥 + 3𝑥𝑦 + 2𝑦 2
𝜕 2 𝑓 −2𝑥 2 − 6𝑥𝑦 − 5𝑦 2
=
𝜕𝑥 2 (𝑥 2 + 3𝑥𝑦 + 2𝑦 2 )2
𝜕2𝑓 −3𝑥 2 − 8𝑥𝑦 − 6𝑦 2
=
𝜕𝑥𝜕𝑦 (𝑥 2 + 3𝑥𝑦 + 2𝑦 2 )2
𝜕𝑦 3𝑥 + 4𝑦
= 2
𝜕𝑥 𝑥 + 3𝑥𝑦 + 2𝑦 2
𝜕2𝑓 −3𝑥 2 − 8𝑥𝑦 − 6𝑦 2
=
𝜕𝑦𝜕𝑥 (𝑥 2 + 3𝑥𝑦 + 2𝑦 2 )2
𝜕 2 𝑓 −5𝑥 2 − 12𝑥𝑦 − 8𝑦 2
=
𝜕𝑦 2 (𝑥 2 + 3𝑥𝑦 + 2𝑦 2 )2
Por lo que la matriz hessiana queda de la siguiente forma
−2𝑥 2 − 6𝑥𝑦 − 5𝑦 2 −3𝑥 2 − 8𝑥𝑦 − 6𝑦 2

(𝑥 2 + 3𝑥𝑦 + 2𝑦 2 )2 (𝑥 2 + 3𝑥𝑦 + 2𝑦 2 )2
𝐻(𝑓(𝑥, 𝑦)) =
−3𝑥 2 − 8𝑥𝑦 − 6𝑦 2 −5𝑥 2 − 12𝑥𝑦 − 8𝑦 2
2 2 2 (𝑥 2 + 3𝑥𝑦 + 2𝑦 2 )2 )
((𝑥 + 3𝑥𝑦 + 2𝑦 )
Evaluando la matriz en el punto xT = (1,1)
13 17
− −
𝐻(𝑓(1,1)) = ( 36 36)
17 25
− −
36 36
Que aproximadamente es
−0.36 −0.47
𝐻(𝑓(1,1)) = ( )
−0.47 −0.69
Usando diferencias finitas
Obtenemos el vector gradiente de f(x,y)
2𝑥 + 3𝑦
𝑥2 + 3𝑥𝑦 + 2𝑦 2
∇𝑓(𝑥, 𝑦) = ( )
3𝑥 + 4𝑦
𝑥 2 + 3𝑥𝑦 + 2𝑦 2
Y la matriz hessiana
−2𝑥 2 − 6𝑥𝑦 − 5𝑦 2 −3𝑥 2 − 8𝑥𝑦 − 6𝑦 2

(𝑥 2 + 3𝑥𝑦 + 2𝑦 2 )2 (𝑥 2 + 3𝑥𝑦 + 2𝑦 2 )2
𝑯(𝑓) =
−3𝑥 2 − 8𝑥𝑦 − 6𝑦 2 −5𝑥 2 − 12𝑥𝑦 − 8𝑦 2
2 2 2 (𝑥 2 + 3𝑥𝑦 + 2𝑦 2 )2 )
((𝑥 + 3𝑥𝑦 + 2𝑦 )
Se calculan las columnas de la matriz hessiana de la función por las tres
aproximaciones. Al 1% de cambio en las variables se tiene δ1=0.01 y δ2=0.01.
En la siguiente tabla se resumen los resultados de la aproximación hacia
adelante, atrás y centrada respectivamente.
Adelante Atras Centrada
-0.3595 -0.4695 -0.3627 -0.4750 -0.3611 -0.4722
-0.4702 -0.6903 -0.4742 -0.6987 -0.4722 -0.6945
Podemos notar que los resultados son muy próximos a los obtenidos de
manera analítica. En el archivo Excel adjunto se muestran los detalles del
calculo.
4. Contamos con nuestras ecuaciones:
𝑦 = −𝑥 2 + 𝑥 + .75
𝑥 2 = 𝑦 + 5𝑦𝑥
Por lo cual procedemos a calcular nuestra matriz Jacobiana:
1 − 2𝑥 −1
𝑱(𝑥, 𝑦) = ( )
−2𝑥 + 5𝑦 1 + 5𝑥
Sacando la matriz 𝑥 𝑇 = (1,1):

1
𝐽𝑓 (𝑥) = ( )
1
El ejercicio solo nos pide con diferencia finitas hacia delante, pero al ser buenos ingenieros
procederemos a hacer los tres métodos:
Teniendo como resultado nuestras funciones:

Por el ultimo obtenemos nuestras comparaciones de Matrices Numéricas con la Matriz Jacobiana
Analítica: así como sus respectivos errores:

Tarea 5

Cargado por

Copyright:

Formatos disponibles

Tarea 5

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Tarea 5

Cargado por

Copyright:

Formatos disponibles

DIVISIÓN DE CIENCIAS BÁSICAS E INGENIERÍA (CBI)

DEPARTAMENTO DE PROCESOS E HIDRAÚLICA

UEA: MÉTODOS NUMÉRICOS APLICADOS A LA INGENIERÍA

PROFESOR: ALEJANDRO VÁZQUEZ RODRÍGUEZ

Para nuestra primera derivada se resuelve como

Por la regla de la suma, la derivada de 𝑓(𝑥) = ⅇ −2𝑥 − 𝑥 respecto de x es

Esto lo haremos por medio de Excel.

Tenemos nuestra función

Evaluando x=2 obtenemos

Posteriormente restamos nuestra h

Para poder evaluar la aproximación para finitas centrales

En donde evaluaremos en nuestro punto x=2, obteniendo como resultado:

Hacia atrás (numérica 2)

Por lo que nos queda

𝑓(𝑥 − 𝛿1 ⅇ1 ) = 3(0.99)(1)2 + 2ⅇ (0.99)(1) = 8.3524

ⅇ𝑟𝑟𝑜𝑟 = |∇𝑓(𝑥) − ∇𝑓̌(𝑥)| = 0.0538

Se nota que hay una mejor precisión en la primera que en la segunda.

−2𝑥 2 − 6𝑥𝑦 − 5𝑦 2 −3𝑥 2 − 8𝑥𝑦 − 6𝑦 2

−2𝑥 2 − 6𝑥𝑦 − 5𝑦 2 −3𝑥 2 − 8𝑥𝑦 − 6𝑦 2

4. Contamos con nuestras ecuaciones:

Sacando la matriz 𝑥 𝑇 = (1,1):

Teniendo como resultado nuestras funciones:

También podría gustarte

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.