MB536 - EP - 2022 - 1 - Solucionario A Publicar

SOLUCIONARIO DEL EXAMEN PARCIAL DE METODOS NUMERICOS
(MB536)
PARTE I
Responda a las siguientes preguntas
𝑥
Si lim𝑥→0 𝑒 𝑥 −1=1 es la solución analítica.
𝑥
1. (0.5P) ¿Calcular usando la calculadora 𝑦 = 𝑒 𝑥 −1 para 𝑥 = 10−15. ¿Cuál es el
error porcentual cometido? ¿Por qué el error es grande?

Solución
x=1e-15
format long
y1=x/(exp(x)-1)
y1=0.95
(1-y1)/1 *100
4.9% debido al error de cancelación
2. (0.5P) Sustituya 𝑒 𝑥 usando la serie de Taylor en torno a 𝑥 = 0 hasta el término

𝑥
cuadrático. Luego, calcule la forma aproximada de 𝑒 𝑥 −1.¿ Cuál es el error
porcentual cometido?
Solución
y2=1/(1+x/2) →1
error= (1-y2)*100 %error percentual -→0%
3. Supongamos que estamos usando una representación de coma flotante de 13 bits

donde hay 1 bit de signo y 4 bits de exponente.
a) (0.5P) ¿La mantisa tiene cuántos bits?
Rpta: 8 bits
b) (0.5P) ¿Cuál es el mínimo exponente (L) y máximo exponente (U) con signo?
BIAS=2^(4-1)-1=7
L=0001-7 = -6
U=1110-7 =7
Rpta: L= -6; U=7
4. (1P) ¿Qué regla no es correcta en el método de Eliminación Gaussiana?. Seleccione
solo una de las opciones:
a) Multiplicar filas por una constante.
b) Intercambiar filas.
c) Sumar un número a una fila.
d) Sumar a una fila el múltiplo de otra fila
P á g i n a 1 | 10
(MB536)
(1P) Determinar el valor (V/F) de la siguiente proposición:
Si la matriz A de orden 4 admite una factorización LU por Crout entonces la traza

de la matriz U es igual a 4.
a) Verdadero
b) Falso.
−4 0 2
5. (1P) Sea la siguiente matriz: 𝐴 = [ 0 −4 7]
0 0 3
¿Es conveniente aplicar el método directo de la potencia? SIN REALIZAR
ITERACIONES. Justifique su respuesta.
Respuesta:
Por ser una matriz triangular: 𝜆1 = −4, 𝜆2 = −4, 𝜆3 = 3
No es conveniente ya que el valor propio dominante debe ser único.
Es decir, debe cumplir: |𝜆1 |>|𝜆2 |≥|𝜆3 |
1 2
6. (1P) Sea la siguiente matriz: 𝐴 = [ ]. ¿Es posible diagonalizar la matriz A?
0 1
Justifique su respuesta.
Respuesta:
𝜆1 = 1 es de multiplicidad 2.
𝑘1 = 2
Determinación del subespacio propio:
(𝐴 − 𝜆1 ) 𝑋 = 0
0 2 𝑥1 0
[ ] [𝑥 ] = [ ]
0 0 2 0
1
𝑋 = 𝑥1 [ ] 𝑥1 ≠ 0
0
Siendo la dimensión:
𝑑1 = 1
Como no satisface 𝑘1 = 𝑑1 , la matriz A no es diagonalizable.
P á g i n a 2 | 10
(MB536)
7. (1P) Complete el código que falta con respecto al test de la dominancia estricta de
la matriz A de orden n.
s=0;
for i=1:n
j=………………………………………………;
j(i)=……………………………………;
if abs(A(i,i))>sum(abs(A(i,j)))
s=s+1;
end
if n==s
disp(------------------------------------)
end
end
Respuesta:
n=4;
s=0;
for i=1:n
j=1:n;
j(i)=[];
if abs(A(i,i))>sum(abs(A(i,j)))
s=s+1;
end
if n==s
disp('A tiene diagonal estrictamente dominante')
end
end
P á g i n a 3 | 10
(MB536)
8. (1P) Sea el siguiente código de programación MATLAB, para analizar la convergencia del
método de iteración de punto fijo para la ecuación x3-3x-20=0:
g=____________________
dg=diff(g)
x0=3 % Donde |g’(x)| alcanza su máximo valor.
m=abs(double(subs(dg,x0)))
if m<1
disp('Punto Fijo convergente')
else
disp('Punto Fijo divergente')
end
% m=____________________
Complete el programa de tal manera que la condición lógica sea verdadera y muestre el
valor de m.
Respuesta(s)
((3*x+20)/x)^0.5 m=0.3574
(3*x+20)^(1/3) m=0.1059
P á g i n a 4 | 10
(MB536)
PARTE II
Problema 1
Para determinar las pérdidas en tuberías que transportan un flujo turbulento, en

ocasiones se emplea la siguiente ecuación para determinar el factor de fricción de
Darcy 𝑓 en cálculos computacionales:
𝑏
𝑓=𝑎+
𝑅𝑐
En donde:
∈ 0.225 ∈
𝑎 = 0.104 ( ) + 0.532 ( )
𝐷 𝐷
∈ 0.44
𝑏 = 88 ( )
𝐷
∈ 0.134
𝑐 = 1.62 ( )
𝐷
∈
Si 𝐷 = 0.001, es la rugosidad relativa calculado con una precisión de 1% y el número
de Reynolds R está entre 19944 y 19980. Estime:
a) (1.5P) El error absoluto esperado en las variables a, b y c.

b) (1.5P) La cota de error absoluto y relativo esperado en el cálculo del factor de
fricción.
c) (1.0P) El intervalo probable donde se encontrará el valor exacto de f.
Solución:
∈
a) Haciendo: 𝐷 = 𝑒𝐷
𝑒𝐷 = 0.001
𝜀𝑒𝐷 = 0.01 ∗ 0.001 = 10−5
𝑅 = 19968
𝜀𝑅 = 18
𝜕𝑎
𝜀𝑎 = | | 𝜀𝑒𝐷 = |5.4776| ∗ 0.00001 = 5.4776 ∗ 10−5
𝜕𝑒𝐷
𝜕𝑏
𝜀𝑏 = | | 𝜀𝑒𝐷 = |1.8533 ∗ 103 | ∗ 0.00001 = 0.0185
𝜕𝑒𝐷
𝜕𝑐
𝜀𝑐 = | | 𝜀𝑒𝐷 = |86.024| ∗ 0.00001 = 8.6024 ∗ 10−4
𝜕𝑒𝐷
P á g i n a 5 | 10
(MB536)
b)
𝜕𝑓 𝜕𝑓 𝜕𝑓 𝜕𝑓
𝜀𝑓 = |𝜕𝑎| 𝜀𝑎 + |𝜕𝑏| 𝜀𝑏 + | 𝜕𝑐 | 𝜀𝑐+|𝜕𝑅𝑒| 𝜀𝑅
𝜀𝑓 = |1|𝜀𝑎 + |0.0017|𝜀𝑏 + |−0.0017|𝜀𝑐+|−2.3507 ∗ 10−7 |𝜀𝑅
𝜀𝑓 = 1.5343 ∗ 10−4
c)
𝑎 = 0.0225
𝑏 = 4.2119
c=0.6420
𝑏
𝑓 = 𝑎 + 𝑅𝑐 =0.02982164
𝑓 − 𝜀𝑓 ≤ 𝐹 ≤ 𝑓 + 𝜀𝑓
0.02966821 ≤ 𝐹 ≤ 0.02997507
P á g i n a 6 | 10
(MB536)
Problema 2
En la Figura 1 se muestra tres tanques en serie en un proceso industrial. Los tanques

tienen sección transversal unitaria y flujos de
agua x,y, z respectivamente.
El flujo de entrada en el primer tanque es u, el

flujo de salida que conecta los tanques 1 y 2 es
6 (x-y), el flujo de salida en los tanques 2 y 3 es
5 (y-z), y el flujo de salida del tanque 3 es de 5z.
Suponiendo que el sistema se encuentra en
estado estacionario, se desea encontrar el nivel
de agua (x, y, z) en cada tanque.
Considere u=2 lt/min. Figura 1 Sistema de tanques industriales
Nota: Considere la ley de conservación de masa en estado estacionario.
Se pide:
a) (1P) Modele el sistema lineal a resolver.

b) (1P) Verifique que el método de Jacobi converge. Justifique. (NO REALICE
ITERACIONES).
c) (1P) Realice tres iteraciones usando el método de Jacobi.
d) (1P) ¿Cuál es el error porcentual en la tercera iteración?
Nota: Tome como valor exacto x=𝐴−1 𝑏. Use la norma infinita. x(0) = [0 0 0]T.
Solución:
a)
b) No tiene diagonal estrictamente dominante por lo que es

necesario el cálculo del radio espectral para determinar si el
método de Jacobi Converge.
P á g i n a 7 | 10
(MB536)
c) y d)
P á g i n a 8 | 10
(MB536)
Problema 3
En un experimento militar se ha lanzado un proyectil

descrito por una trayectoria parabólica representada en
la figura:
La siguiente función
𝑣0 cos(𝛼)(𝑣0 sin(𝛼)+√sin2 (𝛼)𝑣02 +2𝑔𝑑 )

𝑅(𝛼 ) = ,
𝑔
permite obtener el alcance R del proyectil. El alto mando de la fuerza militar desea
determinar el ángulo de lanzamiento 𝛼 (rad) que permita obtener un alcance de
R=162.87 m.
𝑚 m
Considere las condiciones iniciales: 𝑔 = 9.81 𝑠2 ; d = 100 m. ; v0 = 30 s
a) (1P) Localice la raíz o raíces de la ecuación con intervalos de longitud 0.1.

Seleccione el(los) intervalos que asegure la existencia de al menos una raíz según
el Teorema de Bolzano.
i. [0.2;0.3]
ii. [0.3;0.4]
iii. [0.4;0.5]
b) (2P) A partir del intervalo inicial obtenido en (a) aplique el método de la bisección
para aproximar la raíz realizando 02 iteraciones e indique el error en cada iteración.
c) (1P) Halle el número de iteraciones necesarios teóricas para aproximar la raíz en
el intervalo [0.2;0.5] con una tolerancia 10−5 .
Solución:
a)
𝑣0 cos(𝛼) (𝑣0 sin(𝛼) + √sin2(𝛼)𝑣02 + 2𝑔𝑑 )
𝑓(𝛼) = − 162.87
𝑔
𝛼 rad. 𝑓(𝛼)
0.2 -11.0189
0.3 -4.9608
0.4 -0.8985
0.5 0.7487
Existe una raíz entre 0.4 y 0.5
P á g i n a 9 | 10
(MB536)
b)
Iteración 0
Xi=0.4 f(0.4)=-0.8985
Xs=0.5 f(0.5)=0.7487
Xr=(0.4+0.5)/2=0.45 f(0.45)= 0.2515
eXr=(0.5-0.4)/2=0.05
f(0.4)*f(0.45)<0, la raíz está en el lado izquierdo del intervalo
Iteración 1
Xi=0.4 f(0.4)= -0.8985
Xs=0.45 f(0.45)= 0.2515
Xr=(0.4+0.45)/2=0.425 f(0.425)=-0.2448
eXr=(Xs-Xi)/2=0.025
f(0.425)*f(0.45)<0, la raíz está en el lado derecho del intervalo
Iteración 2
Xi=0.4 f(0.4)= -0.8985
Xs=0.45 f(0.45)= 0.2515
Xr=(0.4+0.45)/2=0.425 f(0.425)=-0.2448
eXr=(Xs-Xi)/2=0.025
f(0.425)*f(0.45)<0, la raíz está en el lado derecho del intervalo
i Xi Xr Xs f(Xi) f(Xr) f(Xs) Err

0 0.4 0.5 0.45 - + + 0.025
1 0.4 0.45 0.425 - + - 0.125
2 0.425 0.45 0.4375 - + - 0.0625
La raiz aproximada es Xr=0.4375.

c)
|0.5 − 0.2|
< 𝑇𝑂𝐿
2𝑁+1
0.3
< 0.00001
2𝑁+1
𝑁 > 13.8727
Se requieren como mínimo 14 iteraciones para alcanzar la precisión deseada.
P á g i n a 10 | 10

MB536 - EP - 2022 - 1 - Solucionario A Publicar

Cargado por

Copyright:

Formatos disponibles

MB536 - EP - 2022 - 1 - Solucionario A Publicar

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

MB536 - EP - 2022 - 1 - Solucionario A Publicar

Cargado por

Copyright:

Formatos disponibles

SOLUCIONARIO DEL EXAMEN PARCIAL DE METODOS NUMERICOS

error porcentual cometido? ¿Por qué el error es grande?

2. (0.5P) Sustituya 𝑒 𝑥 usando la serie de Taylor en torno a 𝑥 = 0 hasta el término

3. Supongamos que estamos usando una representación de coma flotante de 13 bits

(1P) Determinar el valor (V/F) de la siguiente proposición:

Si la matriz A de orden 4 admite una factorización LU por Crout entonces la traza

Por ser una matriz triangular: 𝜆1 = −4, 𝜆2 = −4, 𝜆3 = 3

No es conveniente ya que el valor propio dominante debe ser único.

Es decir, debe cumplir: |𝜆1 |>|𝜆2 |≥|𝜆3 |

Determinación del subespacio propio:

Como no satisface 𝑘1 = 𝑑1 , la matriz A no es diagonalizable.

x0=3 % Donde |g’(x)| alcanza su máximo valor.

disp('Punto Fijo convergente')

disp('Punto Fijo divergente')

Para determinar las pérdidas en tuberías que transportan un flujo turbulento, en

a) (1.5P) El error absoluto esperado en las variables a, b y c.

𝜀𝑒𝐷 = 0.01 ∗ 0.001 = 10−5

𝜀𝑓 = |1|𝜀𝑎 + |0.0017|𝜀𝑏 + |−0.0017|𝜀𝑐+|−2.3507 ∗ 10−7 |𝜀𝑅

En la Figura 1 se muestra tres tanques en serie en un proceso industrial. Los tanques

El flujo de entrada en el primer tanque es u, el

Nota: Considere la ley de conservación de masa en estado estacionario.

a) (1P) Modele el sistema lineal a resolver.

b) No tiene diagonal estrictamente dominante por lo que es

En un experimento militar se ha lanzado un proyectil

𝑣0 cos(𝛼)(𝑣0 sin(𝛼)+√sin2 (𝛼)𝑣02 +2𝑔𝑑 )

a) (1P) Localice la raíz o raíces de la ecuación con intervalos de longitud 0.1.

i Xi Xr Xs f(Xi) f(Xr) f(Xs) Err

La raiz aproximada es Xr=0.4375.

También podría gustarte