Lista04fisi MN 2013 II

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 5

Universidad Nacional Mayor de San Marcos

Facultad de Ingeniera de Sistemas e Informatica


E.A.P. de Ingeniera de Sistemas
Curso de Metodos Numericos 2013 II
Resoluci on numerica de un sistema de ecuaciones lineales
Lista de ejercicios 04
Resolucion numerica de un sistema de ecuaciones lineales
1
1. Emplee el metodo de eliminacion gaussiana para resolver los siguientes ejercicios:
(a)
_

_
7.6x
1
+ 0.5x
2
+ 2.4x
3
= 1.9
2.2x
1
+ 9.1x
2
+ 4.4x
3
= 9.7
1.3x
1
+ 0.2x
2
+ 5.8x
3
= 1.4
(b)
_

_
2x + 2y + 2z + 2w = 12
4x + 6y + 6z + 6w = 34
6x + 14y + 16z + 16w = 82
2x + 2y + 12z + 14w = 42
(c)
_

_
3x 0.1y 0.2z = 7.85
0.1x + 7y 0.3z = 19.3
0.3x 0.2y + 10z = 71.4
(d)
_

_
x + y + 3w = 4
2x + y z + 3w = 1
3x y z + 2w = 3
x + 2y + 3z w = 4
(e)
_

_
0.3x 0.2y + 10z = 71.4
0.1x + 7y 0.3z = 19.3
3x 0.1y 0.2z = 7.85
(f)
_

_
2x + 3y z = 5
4x + 4y 3z = 3
2x 3y + z = 1
(g)
_

_
2x + y + w = 1
3x +
1
2
y + z + w = 2
4x + 2y + 2z + w = 1
y + z + 2w = 0
2. En los problemas del item anterior, en caso de ser factible resuelva empleando el metodo
de descomposicion LU.
3. Obtenga el determinante de la matriz asociada a cada uno de los sistemas de ecuaciones
lineales en el caso de poderse resolverse mediante el item anterior.
1
Lista elaborada a partir de los enunciados de los textos: Metodos Numericos (Burden R. & Douglas
J.), Metodos Numericos con MatLab (Mathews J. & Fink K.), Metodos Numericos aplicados a la
Ingeniera (Nieve H. & Sanchez D.)
Prof. J. R. Avenda no Q. 1
Universidad Nacional Mayor de San Marcos
Facultad de Ingeniera de Sistemas e Informatica
E.A.P. de Ingeniera de Sistemas
Curso de Metodos Numericos 2013 II
Resoluci on numerica de un sistema de ecuaciones lineales
4. Considere la matriz A de orden n = 4, denida por:
a
i,j
=
_
i(n j + 1) , i j
a
j,i
, i > j
(a) Obtenga la descomposicion LU considerando los elementos L
i,i
= 1
(b) Resuelva el sistema de ecuaciones Ax = b donde b = (1, 2, 3, 4)
T
Sol: L =
_
_
_
_
1 0 0 0
3/4 1 0 0
1/2 2/3 1 0
1/4 1/3 1/2 1
_
_
_
_
, U =
_
_
_
_
4 3 1 1
0 15/4 10/4 5/4
0 0 10/3 5/3
0 0 0 5/2
_
_
_
_
, (0, 0, 0, 0, 1)
T
5. Use el metodo de Jacobi para resolver el primer item (en caso de ser posible), para
ello, aproxime la solucion con un error relativo por debajo de 10
3
.
6. Haga lo mismo que el item anterior, pero para ello emplee el metodo de Gauss-Seidel.
7. Considere el siguiente sistema de ecuaciones lineales:
_
_
_
_
8 9 1 0
1 8 0 2
1 0 4 1
1 2 3 5
_
_
_
_
x =
_
_
_
_
1
5
0
2
_
_
_
_
Resuelva empleando el esquema de Gauss Seidel, vericando que se cumpla las hipotesis
del metodo, considere un error relativo menor que 10
5
. Adjunte una tabla con las
aproximaciones y los respectivos errores. Sol: (1, 1, 0, 1)
T
8. La secuencia de Fibonacci es generada usando la ecuacion de diferencias:
f
j
= f
j1
+ f
j2
; j > 1 , f
0
= 0 , f
1
= 1
(a) Demuestre que f
n
f
n+2
f
2
n+1
= (1)
n+1
, n = 0, 1,
(b) Tiene solucion el siguiente sistema de ecuaciones? Justique.
_
f
n
f
n+1
f
n+1
f
n+2
__
x
1
x
2
_
=
_
f
n+2
f
n+3
_
(c) Si su respuesta en (b) es armativa, de la solucion.
9. El siguiente sistema tiene solucion ? Justique.
_
_
_
_
_
_
_
_
_
1
1
2
0 . . . .

1
2
1
1
2
. . . .
0
1
2
1 . . . .
. . .
.
.
.
. . .
. . . .
1
2
1
1
2
. . . . 0
1
2
1
_
_
_
_
_
_
_
_
_
_
_
_
_
_
_
_
_
_
P
1
P
2
P
3
.
.
.
p
n1
P
n
_
_
_
_
_
_
_
_
_
=
_
_
_
_
_
_
_
_
_
1
2
0
0
.
.
.
0
0
_
_
_
_
_
_
_
_
_
Prof. J. R. Avenda no Q. 2
Universidad Nacional Mayor de San Marcos
Facultad de Ingeniera de Sistemas e Informatica
E.A.P. de Ingeniera de Sistemas
Curso de Metodos Numericos 2013 II
Resoluci on numerica de un sistema de ecuaciones lineales
10. Sea Ax = b un sistema de ecuaciones lineales y x = Bx + d un sistema equivalente,
con:
B =
_
_
a 3 a a + 4
a 0 c + 3
0 4 c
_
_
(a) Determine a y c para que B sea la matriz D del metodo de Jacobi y obtener una
matriz A.
(b) Es convergente el metodo de Jacobi en este caso?.
11. Probar que la iteracion de Jacobi para la resolucion de un sistema Ax = b converge en
un n umero nito de pasos cuando la matriz A es triangular superior e invertible.
12. Se dispone de 3 marcas de fertilizantes que proporcionan Nitrogeno, Acido Fosf orico y
Potasio soluble al suelo. Una bolsa de cada marca proporciona las siguientes unidades
de cada nutriente:
Marca Nitrogeno Acido Fosforico Potasio
A 3 1 2
B 2 2 3
C 1 0 1
Para un crecimiento ideal, el suelo en cierta parte necesita 19 unidades de Nitrogeno,
16 de Acido Fosforico y 10 unidades de Potasio por Hectarea Cuantas bolsas de cada
fertilizante deben usarse por Hectarea para lograr un crecimiento ideal?. Emplee el
metodo de descomposicion LU para hallar la solucion.
13. Una fabrica de muebles de calidad tiene dos divisiones: un taller de maquinas donde se
fabrican las partes de los muebles y una division de ensamblado y terminado, donde las
partes se unen para formar el producto nal. Supongase que en el taller de maquinas
hay 12 empleados y 20 en la division de ensamblado y terminado, y que cada empleado
trabaja jornadas de 8 horas diarias. Supongase ademas que que solo se producen
dos artculos: sillas y mesas. Una silla requiere
384
17
horas de tiempo de taller de
maquina, ademas de
480
17
horas de tiempo de ensamblado y terminado. Mientras
que una mesa requiere
240
17
horas de tiempo de taller de maquinas y
640
17
horas de
tiempo de ensamblado y terminado. Suponiendo que hay una demanda ilimitada de
estos productos y que el fabricante desea mantener ocupados a todos los empleados,
cuantas sillas y cuantas mesas pueden producirse al da en esta planta?
Es posible plantear una sistema matricial?, en este caso obtenga la solucion por el
metodo de eliminacion gaussiana con sustitucion hacia adelante.
14. Una compa na de carga transporto tres tipos de ete en su transporte aereo ligero.
El espacio requerido por cada unidad de los tres tipos de carga eran de 5, 2 y 4 pies
c ubicos respectivamente. Cada unidad de los tres tipos de carga peso 2, 3 y 1 kilogramos
respectivamente, mientras que los valores unitarios de los tres tipos de carga fueron de
S/.10, S/.40 y S/.60 respectivamente. Determine el n umero de unidades de cada tipo
de carga transportada si el valor total de la carga fue de S/.13500, ocupo 1050 pies
c ubicos de espacio y peso 550 000 gramos.
Prof. J. R. Avenda no Q. 3
Universidad Nacional Mayor de San Marcos
Facultad de Ingeniera de Sistemas e Informatica
E.A.P. de Ingeniera de Sistemas
Curso de Metodos Numericos 2013 II
Resoluci on numerica de un sistema de ecuaciones lineales
(a) Formule el problema de forma clara, presentando la forma matricial. Sea explcito.
(b) Emplee el metodo de descomposicion LU para resolver el problema; para esto
considere los elementos en la diagonal de L igual a 2.
15. Tres especies de bacterias se alimentan con tres productos diferentes A, B y C. Una
bacteria de la primera especie consume 1.3 unidades de los productos A y B, ademas
2.3 unidades del producto C diariamente. Una bacteria de la segunda especie consume
1.1 unidades del producto A, 2.4 unidades del producto B y 3.7 unidades del producto
C cada da. Una bacteria de la tercera especie consume 8.1 unidades del producto A,
2.9 unidades del producto B y 5.1 del producto C cada da. Si se suministran cada
da 16000 unidades del producto A, 28000 unidades del B y 44000 unidades del C,
Cuantas bacterias de cada especie pueden mantenerse en este ambiente?
(a) Formule el problema de forma clara, presentando la forma matricial. Sea explcito.
(b) Resuelva empleando el metodo de Eliminacion gaussiana con sustitucion hacia
atras (se debe indicar las operaciones). Sol: A 2339, B 10128, A 224
16. El due no de un bar ha comprado refrescos, cerveza y vino por un importe de 500
Nuevos soles (sin impuestos). El valor del vino es 60 Nuevos soles menos que el de los
refrescos y de la cerveza conjuntamente. Teniendo en cuenta que los refrescos deben
pagar un IGV del 6%, por la cerveza el 12% y por el vino el 30%, lo que hace que la
factura total con impuestos sea de 592.4 Nuevos soles. Obtenga la cantidad invertida
en cada tipo de bebida empleando el metodo de descomposicion LU, considerando los
valores de la diagonal de L todos ellos iguales a 1.
17. Se desea determinar la parabola de ecuacion y = A + Bx + Cx
2
, la cual pasa por los
puntos (0,0), (2,4) y (0.5,1) Es posible emplear el metodo de Jacobi para hallar los
coecientes indicados? Si la respuesta es armativa, proceda a calcularlos con un error
relativo menor que 10
4
, caso contrario indicar que no es posible.
18. El ingeniero Juan Machuca supervisa la produccion de 3 tipos de componentes electricos.
Se requiere de 3 clases de materiales (metal, plastico y hule) para la produccion de es-
tas componentes.

El ha determinado que las cantidades necesarias para producir cada
componente estaran dadas por la siguiente tabla:
Componente Metal gr/comp. Plastico gr/comp Hule gr/comp
1 14.86 0.21 0.81
2 9.12 0.35 1.2
3 4.07 0.12 2.93
Total en un da 2.17 Kg 3.57 Kg 2.87 Kg
Cuantas componentes se puede producir por da? Plantee la resolucion y resuelva
aplicando el esquema iterativo adecuado.
19. Considere el siguiente sistema de ecuaciones lineales:
_

_
10x y + 2z = 6
x + 11y z + 3w = 25
12x 2y + 12z w = 5
3y z + 8w = 15
Prof. J. R. Avenda no Q. 4
Universidad Nacional Mayor de San Marcos
Facultad de Ingeniera de Sistemas e Informatica
E.A.P. de Ingeniera de Sistemas
Curso de Metodos Numericos 2013 II
Resoluci on numerica de un sistema de ecuaciones lineales
(a) Se puede resolver empleando un esquema iterativo?, justique su respuesta.
(E
3
E
3
E
1
)
(b) De ser armativa su respuesta en el item anterior, debera aproximar la solucion
usando el metodo de Jacobi con un error relativo menor que 10
5
; presente las
iteraciones obtenidas en su prueba escrita. Caso contrario omitir este item. (x
1, y 2, z 1, w 1)
20. En el graco las fuerzas sobre el puente verican las ecuaciones:
_

_
F
1
+

2
2
f
1
+ f
2
= 0

2
2
f
1
+

3
2
f
4
= 0
f
2
+ f
5
= 0

3
2
f
4
f
5
= 0
_

2
2
f
1
F
2
= 0

2
2
f
1
f
3
+
1
2
f
4
= 0
f
3
10000 = 0

1
2
f
4
F
3
= 0
(a) Es posible resolverlo mediante un esquema iterativo?, justique.
(b) De ser armativo el item anterior, emplee el metodo de Jacobi con un error menor
que 10
5
. (F
1
0, F
2
2366.0227, F
3
13659.95356, f
1
33460.22721,
f
2
23659.7336, f
3
10000, f
4
27320.24095, f
5
23659.7336)
Referencias:
1. Burden F. & Faires D. Metodos Numericos. International Thomson Editores. 2004.
Prof. J. R. Avenda no Q. 5

También podría gustarte