Prueba Parcial 1, Métodos Numéricos 2015

UNIVERSIDAD AUSTRAL DE CHILE
FACULTAD DE CIENCIAS DE LA INGENIERIA
CENTRO DE DOCENCIA DE CIENCIAS BASICAS

PARA INGENIERIA.
BAIN053 METODOS
NUMERICOS
PARA INGENIERIA
Pauta Prueba Parcial 1
24 abril de 2015
I. En cada una de las siguientes preguntas de selecci
on m
ultiple, marque solo una alternativa.
1. Considere las matrices
1
A= 0
2
2
0 ;
1
0
1
0
A1
1/3
= 0
2/3
0 2/3
1
0 .
0 1/3
Al resolver el sistema Ax = b, en la solucion se comete un error relativo en norma 2 estrictamente mayor a

6 104 , debido a un error en el termino del lado derecho. Suponga que no hay errores en los coeficientes
de la matriz ni errores de redondeo. Entonces, siempre es cierto que:

b b > 2 104 kbk
2
2
Desarrollo:
1
Dado que los valores propios de A son 1, 1, 3, entonces los valores propios de la inversa son 1, 1, .
3
Luego,
||A||2 = 3 , ||A1 ||2 = 1 .
Por hip
otesis
6 104 <
||x x
||2
||x||2
y como K2 (A) = ||A||2 ||A1 ||2 = 3 se obtiene

6 104 <
Y por lo tanto, 6 104 < 3
||x x
||2
||b b||2
||b b||2
K2 (A)
=3
.
||x||2
||b||2
||b||2
||b b||2
, de donde se obtiene 2 104 ||b||2 < ||b b||2 .
||b||2
2. Se sabe que la soluci

on del sistema lineal, Ax = b, de orden 4, esta dada por:
x = (1, 1, 1, 1)t .
De igual modo se sabe que la soluci
on del sistema lineal perturbado (A + A) (x + x) = b, esta dado por:
x + x = (1.2, 0.9, 1.8, 1.5)t .
Considerando que se perturb
o un termino de la matriz A en 0.005, ||A|| = 45 y ||A1 || < 200, siempre
es cierto que:
4 103 K (A) < 9 103
Desarrollo : Como x
= x + x, de las hipotesis se obtiene

||x x
||
||x||
0.8
||A A||
0.005
103
=
=
= 0.8 ,
=
=
,
||x||
||x||
1
||A||
45
9
()
donde A = A + A es invertible dado que ||A1 A|| < 1. Esta u

ltima parte implica la desigualdad
||x x
||
||x||

K (A)
||A A||
||A||
||A A||
1 K (A)
||A||
de donde al reemplazar los valores obtenidos en () se llega a

0.8
K (A)
3
1 K (A)
10
9
103
.
9
As, al despejar K (A) se tiene directamente 4 103 K (A).

1
3. Considere tres metodos diferentes para resolver numericamente el P.V.I.

y 0 (x) = f (x, y(x)),
y(0) = 0,
x [0, 2]
y que el error |y(xk )yk | est

a acotado por 5h1 , 25h22 , 625h43 . El n
umero de iteraciones n1 , n2 , n3 necesarias
para que el error sea menor o igual a 1012 en cada uno de los metodos es, respectivamente:
n1 1013 , n2 107 y n3 104
Desarrollo : Dado que el intervalo es [0, 2], el tama
no del paso en cada caso es
para cada caso se tiene
2
, i = 1, 2, 3. Luego,
ni
625 16
25 4
52
1012 ,
1012 ,
1012 ,
n1
n22
n43
o equivalentemente
1013 n1 , 107 n2 , 104 n3 .

4. Dado el P.V.I.
y 0 (x) = 4xy,
y(0) = 1.
Considere el metodo numerico basado en el metodo de Euler:

y k+1 = yk + hf (xk , yk ),
yk+1 = yk + hf (xk+1 , y k+1 ),
k = 0, 1, . . .
k = 0, 1, . . .
Al utilizar este esquema con paso h = 0.5, el valor que falta en la siguiente tabla
xk
yk
0.0
1.0
0.5
2.0
1.0
es :
10.0
Desarrollo : El valor que falta es y2 , el cual de la definicion de yk se obtiene como

y2 = y1 + hf x2 , y1 + hf (x1 , y1 ) .
De la tabla se obtiene x0 = 0, x1 = 0.5, x2 = 1, y0 = 0, y1 = 1 y dado que h = 0.5, reemplazando se
concluye directamente que y2 = 10.

5. El perodo T [s] de un pendulo simple ideal satisface la igualdad
s
l
T = 2
,
g
donde l [m] denota el largo y g [m/s2 ] la aceleracion de gravedad. Si l se mide con un error relativo menor a 0.04, g se mide con un error relativo menor a 0.02 y se considera solo error de propagaci
on
en T , g y l, entonces para obtener un error relativo menor a 0.08 en el calculo de T siempre es cierto que:
utilizarse con un error relativo menor a 0.05
s
Desarrollo: Usando la f
ormula de propagacion de error para la funcion definida por T = 2
|ER (T)|
1
|T |
l
g

T
T
T

|

| + |l l| +
|g g|

g
l
1
1
|ER (
)| + |ER (l)| + |ER (
g )|
2
2
0.04 1
< |ER (
)| +
+
2
2
Imponiendo |ER (
)| +
0.04 0.02
+
< 0.08 se concluye
2
2
|ER (
)| < 0.05 .
6. Considere un computador imaginario que representa los n

umeros en aritmetica de punto flotante por
redondeo en la forma
0.d1 d2 d3 d4 d5 d6 d7 10m ;
d1 6= 0;
0 di 9,
i = 1, 2, 3, 4, 5, 6, 7;
9 m 9,
m Z.
Si a = 0.000001 , b = 0.00003 , c = 2000000 , x = 21236111115 1010 , y = 0.11512 101 , entonces el

valor de
x+y
(ac)b
obtenido por este computador es:
0.5458018 105
Desarrollo: Reemplazando, se obtiene
f l(x) = 0.2123611 101 , f l(y) = 0.1151200 101 .
Luego, f l(f l(x) + f l(y)) = 0.3274811 101 . Ademas
f l(f l(f l(a) f l(c)) f l(b)) = 0.6000000 104 .
As

fl
f l(f l(x) + f l(y))

f l(f l(f l(a) f l(c)) f l(b))

= fl
0.3274811 101
0.6000000 104
= 0.5458018 105 .
II. Resuelva el siguiente problema utilizando 4 decimales con redondeo.

El sistema de encendido de un autom
ovil esta representado por el circuito que se muestra en la figura
y est
a modelado por el P.V.I.
di(t) i(t)
d2 i(t)
+R
+
= 0,
L
2
dt
dt
C
i(0) = V0 , i0 (0) = 0,
R
siendo i(t) la corriente electrica del circuito en un instante de tiempo t y donde la tension electrica V0 representa
la batera y el alternador, el resistor R la resistencia, el inductor L la bobina de encendido, el capacitor C el
condensador, y el interruptor, paralelo al condensador, representa el encendido electronico del automovil.
Considere el metodo de Runge-Kutta de orden 3 descrito a continuacion:
h
(2K1 + 3K2 + 4K3 ) ;
9
yk+1
yk +
xk+1
xk + h;
k = 0, 1, 2, . . .
k = 0, 1, 2, . . .
donde:
K1
= f (xk , yk )
K2
= f

h
h
xk + , yk + K1
2
2
K3
= f

3
3
xk + h, yk + K2
4
4
para aproximar la soluci

on del P.V.I.

y 0 (x) =
y(x0 ) =
f (x, y(x))
y0
Para V0 = 12 [V ], R = 4 [], C = 1 [F ] y L = 8 [H]:

a) (20 puntos) Plantee el esquema numerico, basado en RK-3, para aproximar la solucion del P.V.I. asociado
al problema de la corriente electrica.
b) (10 puntos) Obtenga una aproximaci
on para la corriente electrica al cabo de 4 [s], con h = 4 [s].
c) (10 puntos) Sabiendo que la soluci
on analtica del P.V.I. es

t
t
i(t) = et/4 3 cos
+ 3 sin
4
4
determine el error relativo que se comete al estimar la corriente electrica a 4 [s] de iniciado el proceso.
Considere h = 4 [s].
Desarrollo:
1. Reemplazando los datos, el P.V.I. que modela el fenomeno planteado es
8 d i(t) + 4 di(t) + i(t) = 0,

dt2
dt
0
i(0) = 3, i (0) = 0,
Introduciendo la variable auxiliar z(t) = i0 (t)
0
i (t)
0
z (t)
i(0)
z(0)
se tiene el P.V.I. de primer orden

= z(t)
1
1
= i(t) z(t)
8
2
= 3
= 0
de donde
1
1
f2 (t, i, z) = i z.
8
2
Utilizando el metodo de Runge-Kutta de orden 3 se tiene que:
f1 (t, i, z) = z,
ik+1
= ik +
h
(2K1 + 3K2 + 4K3 ) ;
9
zk+1
= zk +
h
(2M1 + 3M2 + 4M3 ) ;
9
tk+1
= tk + h;
k = 0, 1, 2, . . .
k = 0, 1, 2, . . .
k = 0, 1, 2, . . .
donde:
K1 = f1 (tk , ik , zk ) ,
M1 = f2 (tk , ik , zk ) ;
k = 0, 1, 2, . . .

h
h
h
h
h
h
K2 = f1 tk + , ik + K1 , zk + M1 , M2 = f2 tk + , ik + K1 , zk + M1 ; k = 0, 1, 2, . . .
2
2
2
2
2
2

3
3
3
3
3
3
K3 = f1 tk + h, ik + K2 , zk + M2 , M3 = f2 tk + h, ik + K2 , zk + M2 ; k = 0, 1, 2, . . .
4
4
4
4
4
4

2. Para obtener una aproximaci
on de la corriente electrica al cabo de t = 4 [s], con h = 4 es necesario hacer
una iteraci
on la que se indica en la tabla:
k
0
1
tk
0
4
ik
3
2
zk
0
-0.8750
K1
0
K2
-0.75
K3
0
M1
-0.375
M2
0
M3
-0.3047
Se concluye que en t = 4 [s] una aproximacion para la corriente electrica es de 2 [A].

3. Sabiendo que la soluci
on analtica de la ecuacion anterior es

t
t
t/4
i(t) = e
3 cos
+ 3 sin
4
4
El error relativo en la corriente electrica es:
ER (i1 ) =
|i(4) i1 |
|1.5250 2|
=
= 0.3115
|i(4)|
|1.5250|

Prueba Parcial 1, Métodos Numéricos 2015

Cargado por

Copyright:

Formatos disponibles

Prueba Parcial 1, Métodos Numéricos 2015

Cargado por

Información del documento

Descripción 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

Prueba Parcial 1, Métodos Numéricos 2015

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD AUSTRAL DE CHILE

FACULTAD DE CIENCIAS DE LA INGENIERIA

CENTRO DE DOCENCIA DE CIENCIAS BASICAS

Al resolver el sistema Ax = b, en la solucion se comete un error relativo en norma 2 estrictamente mayor a

y como K2 (A) = ||A||2 ||A1 ||2 = 3 se obtiene

Y por lo tanto, 6 104 < 3

2. Se sabe que la soluci

donde A = A + A es invertible dado que ||A1 A|| < 1. Esta u

de donde al reemplazar los valores obtenidos en () se llega a

As, al despejar K (A) se tiene directamente 4 103 K (A).

3. Considere tres metodos diferentes para resolver numericamente el P.V.I.

y que el error |y(xk )yk | est

Considere el metodo numerico basado en el metodo de Euler:

6. Considere un computador imaginario que representa los n

Si a = 0.000001 , b = 0.00003 , c = 2000000 , x = 21236111115 1010 , y = 0.11512 101 , entonces el

f l(f l(x) + f l(y))

II. Resuelva el siguiente problema utilizando 4 decimales con redondeo.

para aproximar la soluci

Para V0 = 12 [V ], R = 4 [], C = 1 [F ] y L = 8 [H]:

8 d i(t) + 4 di(t) + i(t) = 0,

se tiene el P.V.I. de primer orden

Se concluye que en t = 4 [s] una aproximacion para la corriente electrica es de 2 [A].

También podría gustarte