Sol EP Teo3-4 MN 2023 1

Métodos Numéricos
Examen Parcial
Pregrado
2023-1
Instrucciones: Se permite el uso de calculadora, una hoja de formulario.

Profesores: Jimmy Mendoza - Rósulo Pérez.
Secciones: 3 - 4.
Duración 90 minutos
CONCEPTUALIZACIÓN
Determine si las siguientes proposiciones son Verdadero (V) o Falso (F) según corresponda.
Justifique claramente cada una de sus respuestas.
1. [2 ptos] Si f (x) = 2x3 + x2 − 1, entonces cualquier valor semilla garantizará la conver-
gencia del método de Newton en la determinación de la raı́z real de la ecuación.
Solución:
Sabemos que la forma iterativa del método de Newton, viene dada por:
f (x(r) )
x(r+1) = x(r) −
f ′ (x(r) )
Analizando lo que sucede con la derivada:
f ′ (x) = 6x2 + 2x = 0 → f ′ (x) = 0 si x = 0, −1/3
Por lo tanto, existen puntos para los cuales no está garantizada la convergencia del
método de Newton.
2. [2 ptos] En el sistema de ecuaciones no lineales:
3
x + 3x + y = 0
kx + ey − 2 = 0
Si el parámetro k toma el valor de 1, entonces el sistema tiene una solución única usando
el método de Newton con una aproximación inicial (x0 , y0 ) = (0, 0).
Solución:
Definiendo:
F (x, y) = x3 + 3x + y
G(x, y) = x + ey − 2
Calculando el Jacobiano para un punto en general:
2
3x + 3 1 3 1
J(x, y) = → J(0, 0) =
1 ey 1 1
Calculando el determinante de dicho Jacobiano, vemos es 2, por lo tanto el sistema tiene
solución única.
1
PROCEDIMENTAL
3. [3 ptos] Dado el número N = 35.3125 representado en la base 10:
a) [1 pto] Expresar al número N en el sistema de base binaria.

Solución:
Primero necesitamos expresar:
35 en binario: 100011
0.3125 en binario: 0101
Entonces N = 100011.01012
b) [1 pto] Exprese el resultado anterior en el formato x = ±(1.d1 d2 ...dp )β × β E e
indique el valor de la mantisa y el exponente entero.
Solución:
Del item anterior:
35.3125 = 100011.01012 = 1.0001101012 × 25
c) [1 pto] Halle la secuencia de bits que representa al número N en la notación de

punto flotante IEEE-754 de precisión simple.
Solución:
Del item anterior, tenemos que el exponente es 5, para adaptar al estándar IEEE -
754, tenemos exp = 5+127 = 132 = 100001002 y la mantisa 00011010100000000000000,
como el número es positivo, el primer bit es 0, finalmente la respuesta es:
01000010000011010100000000000000
4. [3 ptos] El coeficiente de sustentación es una medida (adimensional) de la capacidad

de una superficie aerodinámica para generar sustentación. La expresión que permite
estimar dicho coeficiente es:
L
CL = 1 2
2
ρV A
Si se conocen de manera exacta los siguientes datos: ρ = 1.2 kg/m3 , V = 50 m/s.
Sin embargo, los valores de fuerza de sustentación L = 500 N y A = 10 m2 presentan
una incertidumbre de ±1 N y ±0.1 m2 , respectivamente. Halle el valor del error relativo
porcentual en la determinación del coeficiente de sustentación.
2
Solución:
El coeficiente de sustentación depende de dos variables L y A, por lo tanto:
L
CL = CL (L, A) = 1
2
ρV 2 A
Encontrando el valor estimado del coeficiente:

500
CL = = 0.03333...
0.5 × 1.2 × 502 × 10
Para encontrar el error absoluto usamos la serie de Taylor:

∂CL 1
∂L = 1500A .

∂CL L
∂A = 1500A2 .

Entonces:
1 L
△≈ ×1+ × 0.1
1500A 1500A2
Reemplazando:
△ ≈ 4 × 10−4
Para encontrar el error relativo porcentual:
4 × 10−4
er % = × 100 % = 1.20012001.... %
0.033333
APLICACIONES
5. [10 ptos] En la siguiente figura se muestra una barra, la cual se encuentra en contacto
con una pared a 120°C.
3
La ecuación diferencial que describe la variación de temperatura en la barra es:
d2 T

hc P
− (T − T∞ ) =0
dx2 Ti kA
Ti
El valor de cada parámetro del problema es: L = 0.50 m, P = 0.01 m, A = 2 × 10−5 m2 ,

T∞ = 20 ◦ C, k = 237 W/m.K, hc = 10 W/m2 .K. Considerando que es posible aproximar
la solución reemplazando la segunda derivada por una diferencia finita central:
d2 T

Ti−1 − 2Ti + Ti+1
=
2
dx Ti ∆x2
a) [2 ptos] Discretrice la ecuación diferencial y encuentre una ecuación en diferencias

que permita evaluar la temperatura en cada nodo. Expréselo en términos de hc , P ,
k, A y ∆x.
Solución:
De las igualdades lı́neas arriba, llegamos a:
Ti−1 − 2Ti + Ti+1 hc P
2 = (Ti − T∞ )
∆x kA
Lo cual es equivalente a:
Ti−1 − 2Ti + Ti+1 hc P hc P ∆x2

= (Ti − T∞ ) → Ti−1 − 2Ti + Ti+1 = (Ti − T∞ )
∆x2 kA kA
Ti−1 − 2Ti + Ti+1 = 0.211 (Ti − 20)
b) [2 ptos] Encuentre el sistema de ecuaciones lineales en la forma matricial Ax = b

para estimar las temperaturas en los nodos interiores Ti . Use los valores de los
parámetros del problema.
Solución:
Evaluando:
i = 1 → −T1 + T2 = 0.211(T1 − 20)
i = 2 → T1 − 2T2 + T3 = 0.211(T2 − 20)
i = 3 → T2 − 2T3 + T4 = 0.211(T3 − 20)
i = 4 → T3 − 2T4 + T5 = 0.211(T4 − 20)
Utilizando los datos brindados:
1) −1.211T1 + T2 = −4.22
2) T1 − 2.211T2 + T3 = −4.22
4
3) T2 − 2.211T3 + T4 = −4.22
4) T3 − 2.211T4 = −124.211
Convirtiendo a notación matricial:
    
−1.211 1 0 0 T1 −4.22
 1
 −2.211 1 0   T2  =  −4.22 
   
 0 1 −2.211 1   T3   −4.22 
0 0 1 −2.211 T4 −124.211
c) [2 ptos] Analice la convergencia del método iterativo de Gauss-Seidel (según el

criterio de la diagonal dominancia). Indique claramente su conclusión (converge o
diverge).
Solución:
La matriz de coeficientes claramente es diagonal estrictamente dominante, por lo
tanto el método de Gauss - Seidel converge.
d ) [2 ptos] Luego de descomponer la matriz A en D, L y U, halle la matriz de iteración
de Gauss-Seidel (Tgs ) para aplicar el proceso iterativo.
Solución:
Descomponiendo para encontrar Tgs , tenemos:
     
−1.211 0 0 0 0 −1 0 0 0 −1 0 0
 0 −2.211 0 0    0 0 −1 0    0 0 −1 0 
D= L= U = 
 0 0 −2.211 0  0 0 0 −1 0 0 0 −1
0 0 0 −2.211 0 0 0 0 0 0 0 0
Usando como dato, lo dado en la pizarra el dı́a del examen:

 
0 0.8257 0 0
0 0.3734 0.4522 0 
Tgs ≈ 
0

0.1689 0.2045 0.4522
0 0.0763 0.0925 0.2045
e) [2 ptos] Realice 02 iteraciones

  utilizando el método de Gauss-Seidel. Considere
30
50
como vector semilla x(0) = 
70.

90
Solución:
Efectuando:
x(k+1) = Tgs x(k) + cgs
5
Tenemos:  
44.7729
53.8185
x(1) =
66.9554

86.4615
 
47.9261
53.8677
x(2) =
65.3773

85.7477

Sol EP Teo3-4 MN 2023 1

Cargado por

Copyright:

Formatos disponibles

Sol EP Teo3-4 MN 2023 1

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

Sol EP Teo3-4 MN 2023 1

Cargado por

Copyright:

Formatos disponibles

Métodos Numéricos

Instrucciones: Se permite el uso de calculadora, una hoja de formulario.

3. [3 ptos] Dado el número N = 35.3125 representado en la base 10:

a) [1 pto] Expresar al número N en el sistema de base binaria.

35.3125 = 100011.01012 = 1.0001101012 × 25

c) [1 pto] Halle la secuencia de bits que representa al número N en la notación de

4. [3 ptos] El coeficiente de sustentación es una medida (adimensional) de la capacidad

Encontrando el valor estimado del coeficiente:

El valor de cada parámetro del problema es: L = 0.50 m, P = 0.01 m, A = 2 × 10−5 m2 ,

a) [2 ptos] Discretrice la ecuación diferencial y encuentre una ecuación en diferencias

Ti−1 − 2Ti + Ti+1 hc P hc P ∆x2

b) [2 ptos] Encuentre el sistema de ecuaciones lineales en la forma matricial Ax = b

c) [2 ptos] Analice la convergencia del método iterativo de Gauss-Seidel (según el

Usando como dato, lo dado en la pizarra el dı́a del examen:

e) [2 ptos] Realice 02 iteraciones

También podría gustarte