Universidad Nacional de Ingeniería 19-11-05

Facultad de Ingeniería Mecánica P.A. 2005-2

Área de Ciencias Básicas y Humanidades

SOLUCIONARIO DE LA PRIMERA PRÁCTICA CALIFICADA

CALCULO NUMERICO (PARTE B)

APELLIDOS Y NOMBRE SECCION NOTA

Marque la alternativa que considere correcta o escriba su respuesta según el caso:

Pregunta 1

Cuál es la salida de S en la siguiente instrucción:

>> A = [1:2:5;4 7 8; ones(1,3)] ; S = sum(diag(A([3,2],[1:2:length(A)])))

a) 8 b) 15 c) 12 d.) 9 e) N.A

Pregunta 2

if any(any(tril(A,-1)));error('La matriz ……………………………………………………….');end

Completar lo que falta en el mensaje

La matriz A no es triangular superior

Pregunta 3

Si A es una matriz cuadrada, la instrucción: sum(diag(A(:,end:-1:1))), evalúa:

a) La suma de elementos de la matriz A

b) La suma de elementos de la diagonal principal de A
c) La suma de elementos de la diagonal secundaria de A
d) La suma de elementos de la ultima fila de A
e) N.A.

Pregunta 4

Si A =[ 1 1 1 1; 1 2 3 4; 1 3 6 10; 1 4 10 20]
for i=1:3
A(i+1,:)=A(i,:)
end

El nuevo valor de A será:

1 1 1 1
1 1 1 1
A
1 1 1 1
 
1 1 1 1
Pregunta 5

¿Cuál es el valor de S?

x = [ 1 2 -3 -4 5 ];
sum = 0;
k = 1;
while x(k)>=0 & k < length(x)
sum = sum + x(k);
k = k + 1;
end
S= sum

a) 3 b) 5 c) 1 d) 8 e) N.A

Pregunta 6
x  sin x
Escriba una función propia que dado el vector x calcule el valor de f  x  
tan x
y grafique de color magenta usando un ícono de ‘+’ en cada punto.

function [y]=f(x)

y=(x-sin(x))./tan(x);
plot(x,y,’m+’);

Pregunta 7

En número decimal 0.1 en un sistema de punto flotante F(2,10,-7,8) se expresara como:

0.000110011001100…=0.1100110011x2-3= 0.09997558593750

Pregunta 8

Hallar la raíz menor en modulo de la ecuación

x 2  40 x  0.25  0

Utilizando aritmética de 4 dígitos y comparar con el resultado obtenido utilizando aritmética exacta.
Calcular el error relativo y asegurarse de comprender de donde viene la perdida de dígitos
significativos.
a) 20.01% b) 20.55 % c) 16.15% d) 20.95% e) N.A.

La pérdida de cifras significativas es por

Efecto de cancelación al restar dos números parecidos

Pregunta 9

En que intervalo debe estar x* para que aproxime a 1000 con 4 cifras decimales exactas?
1000±0.5*10-4
Pregunta 10

Encuentre todos los valores de a tales que el sistema

x1  x2  x3  2
x1  2 x2  x3  3
x1  x2  (a 2  5) x3  a
tenga:

a) Solución única, a 2

b) infinitas soluciones, a=2

c) ninguna solución. a = -2

Pregunta 11

En la primera parte de la Eliminación Gaussiana la matriz A quedó reducida a la matriz B

 25 c 1 25 5 1 
   
A= 64 a 1 , B  0  4.8  1.56 . Cuál es el determinante de A?
   
144 b 1  0 0 0.7 

| A | = -84

Pregunta 12

Del ejemplo anterior dar los valores a ,b ,c

c = 5...........................a = 8.................................b = 12.............................

Pregunta 13

3 1
Sea la matriz A    , con respecto a la factorizacion de Doolite, cuál es la alternativa
 1 4 
correcta?

 3 1   1 0 3 1   1 0
a) U   b) L   c) U   d) L   1 / 3 1  e) N.A.
0  13 / 3 
1 / 3 1  0  13 / 3  

Los Profesores
Universidad Nacional de Ingeniería 19-11-05
Facultad de Ingeniería Mecánica P.A. 2005-2
Área de Ciencias Básicas y Humanidades

SOLUCIONARIO DE LA PRIMERA PRÁCTICA CALIFICADA

CALCULO NUMERICO (PARTE B)

APELLIDOS Y NOMBRE SECCION NOTA

Marque la alternativa que considere correcta o escriba su respuesta según el caso:

Pregunta 1

Si:
>> A=[2 3 4 6;6 2 3 2; 2 6 1 1; 2 3 2 4]
>> B=A(diag(A),3:-1:1)
Hallar : B(2,3)+B(3,1)

a) 8 b) 5 c) 10 d) 4 e) N.A

Pregunta 2

Si
>>A=zeros(2);
>>A(:)=8:-2:2
>>diag(A,-1)
El valor es:
a) 8 b) 4 c) 6 d)2 e) N.A.

Pregunta 3

Cual es el valor de z :

» x=(1:3)'
» M=[x.^0 x.^1 x.^2]
» A=M(2:3,2:3)
» b=M(1:2:3,2)
» z=A\b

a) [-0.5000 0.5000]’ b) [-0.5000 -0.5000]’ c) [0.5000 -0.5000]’ d) [0.5000 0.5000]’ e)

N.A.

Pregunta 4

Las instrucciones:

j=1;
N=length(v)
for i=N:-1:1
 x(j)=v(i);
j=j+1;
end

equivale a:
a) x=v(N:1)
b) x=v(end:-1:1)
c) x=v
d) v=x(end:-1:1)
e) N.A.

Pregunta 5

Cuál es la salida del siguiente segmento de programa

i=1;
while (i<4)
j=1;
while (j>i)
disp(j);
end
i=i+1;
end
disp(i);

a) 1234 b) 123 c) 234 d)34 e)4

Pregunta 6

Escriba un programa script que use un generador de números aleatorios 'rand' para determinar la
cantidad de números aleatorios necesarios para que la suma de estos números por lo menos sea
igual a 20. Su resultado deberá presentar formato y un mensaje de salida.

s=0; i=0;
while s<20
r=rand;
s=s+r;
i=i+1;
end
disp(‘Cantidad de Números aleatorios : ‘)
disk(i)

Pregunta 7
  0 .1d 2 d 3 d 53 x 2 e
En MATLAB un número se almacenará como x  , donde d k  0 o 1 ,
 1021  e  1024 . Para obtener el OVERFLOW en sistema decimal podemos usar la instrucción:

a) sum(2.^(970:1022))
b) sum(2.^(971:1023))
c) sum(2.^(972:1024))
d) sum(2.^(973:1025))
e) N.A.
Pregunta 8

En MATLAB el realmax es:

a) Es la precision de la maquina
b) El mayor valor que se puede almacenar (overflow)
c) realmax + 1 da inf
d) Es el epsilon de la maquina
e) N.A.

Pregunta 9

Sean x = 0.43574628 y y = 0.43574781. Si usamos aritmética (de punto flotante) decimal con
redondeo a 6 dígitos,
Hallar :
x y  fl ( fl ( x)  fl ( y ))

0.435746 - 0.435748 = 0.000002

Pregunta 10

Dado el siguiente sistema

x1  (1 / 2) x 2  (1 / 3) x3  1
(1 / 2) x1  (1 / 3) x 2  (1 / 4) x3  0
(1 / 3) x1  (1 / 4) x 2  (1 / 5) x3  0

Determinar la suma de todos los elementos pivotes en la solución del sistema por el método de
eliminación Gaussiana.

1 + 1/12 = 13/12

Pregunta 11

 1 10
Si A    , Es posible la reducción de Cholesky? Justifique.
10 12
.
No porque A no es definida positiva.

Pregunta 12

3 1 
Sea la matriz A    , con respecto a la factorización de Crout, cuál es la alternativa correcta?
1  2

1 5 / 3 3 0  1 1 / 3 2 0 
a) U   b) L   c) U   d) L    e) N.A.
0 1  
1  5 / 3 0 1  
1 5 / 3
Pregunta 13

En la descomposición de Cholesky, el sistema de ecuaciones Ax=b, la matriz A se descompone en:

a. L*UT
b. L*LT
c. U*L
d. L T*L

Los Profesores
Universidad Nacional de Ingeniería 19-11-05
Facultad de Ingeniería Mecánica P.A. 2005-2
Área de Ciencias Básicas y Humanidades

SOLUCIONARIO DE LA PRIMERA PRÁCTICA CALIFICADA

CALCULO NUMERICO (PARTE B)

APELLIDOS Y NOMBRE SECCION NOTA

Marque la alternativa que considere correcta o escriba su respuesta según el caso:

Pregunta 1

x 1
Para graficar la función f ( x )  e ln( x ) x  [1, 2] , realizamos la siguiente instrucción
en MatLab. (Seleccione todas las que aplique)

a. f = inline('exp(x+1).ln(x)'), x=[1:0.01:2]; plot(x,f(x));

b. f = inline('exp(x+1).*log(x)'), x=[1:0.01:2]; plot(x,f(x));
c. f = inline('exp(x+1).*ln(x)'), x=[1:0.01:2]; plot(x,f(x));
d. f = inline('exp(x+1).log(x)'), x=[1:0.01:2]; plot(x,f(x));
e. f = inline('e.^(x+1).*log(x)'), x=[1:0.01:2]; plot(x,f(x));

Pregunta 2

Si A es una matriz cuadrada de hilbert de orden 3 (aij= 1/(i+j-1)), las siguientes instrucciones:
>>format rat
>>B=tril(triu(A,-1), 1)
>>norm(B,inf)
Dará como resultado:
a) 2/3 b) 13/12 c) 3/2 d) ½ e) N.A.

Pregunta 3

Cual es el resultado de :

»ceil(-0.5)+fix(-0.5)+floor(-0.5)+round(-0.5)

a) –1 b) –1.5 c) –2 d) -2.5 e) N.A.

Pregunta 4

Escriba una función propia que ordene un vector en forma descendente:

function [v]=ordena(x)
n=length(x);
for i=1:n-1
for j=i+1:n
if x(i)<x(j)
temp=x(i);
x(i)=x(j)
x(j)=temp;
end
end
end
v=x;

o también,

function [v]=ordena(x)
v=-sort(-x);

Pregunta 5

Indique la salida para el siguiente programa

i=1;
j=0;
while (i <4)
j= j+1;
i=i+1;
end
disp([i j])

a) 4 5 b) 5 3 c) 5 6 d) 4 6 e) N.A [4 3]

Pregunta 6

Dada la función

function y=Fibonacci(n)
if (n<=2)
y=1;
else
y=Fibonacci(n-1)+Fibonacci(n-2);
end

Cuál es la salida del programa cuando n=5? (Seleccione la respuesta correcta)

a) 2 b) 3 c) 4 d) 5 e) N.A
Pregunta 7

En MATLAB el realmin es:

a) El menor valor que se puede almacenar

b) El menor valor positivo que se puede almacenar
c) El menor valor positivo normalizado que se puede almacenar
d) Es el underflow
e) N.A.

Pregunta 8

En número decimal 0.75 en un sistema de punto flotante F(2,8,-3,4) se expresara como:

0.11000000x2 0

Pregunta 9

Encontrar las raíces de la ecuación cuadrática

x 2  400.2 x  80  0

Usando la fórmula usual y aritmética decimal con redondeo a 4 dígitos. Dichas raíces satisfacen la
ecuación. Explique.

Solución exacta : 0.2 y 400

Aritmética decimal con 4 dígitos :

400.2  400.2 2  4 x80 400.2  160160  320

x12    399.995, 0.204999
2 2

P(0.2)  0 P( 400)  0 P(399.995)  1.99897 P (0.204999)  1.99858

Se observa que el polinomio esta mal condicionado.

Pregunta 10

Dada la siguiente descomposición LU de Doolittle de la matriz A .

 1 0 0 4 1 0 
L   1 / 2 1 0
 U  0 5 / 2 1 

 1 / 2 1 / 5 1  0 0 4 / 5

resolver el sistema de ecuaciones Ax = b siendo

b  1  2 7 
T

x=[5/4 -4 17/2]T