Laboratorio 01

Matemáticas III
Matrices y
Sistema de Ecuaciones
Lineales
Semana 01
Hermes Pantoja Carhuavilca
([email protected])
Brigida Molina Carabaño
([email protected])
Rosulo Perez Cupe
([email protected])
Asistente: Angelo Pacheco
([email protected])
Profesores: H. Pantoja - B. Molina - R. Perez

Índice
1 Parte Teórica
2 Parte Práctica
Profesores: H. Pantoja -B. Molina - R. Perez Matemáticas III April 26, 2021 1 / 23
1
PARTE TEÓRICA
Matrices y Determinantes
Función Descripción
rank(A) Rango de la matriz A.
det(A) Determinante de la matriz A.
trace(A) Traza de la matriz A, suma de los elementos de la diagonal
inv(A) Calcula la inversa de la matriz A
A’ Transpuesta de la matriz A.
diag(A) Extrae la diagonal de la matriz A como vector columna
eye(n,n) Crea la matriz identidad de orden n
zeros(m,n) Crea la matriz nula de orden m × n
ones(m, n) Crea la matriz de orden m× n con todos sus elementos unos.
[n,m]=size(A) Devuelve el número de filas y columnas de la matriz A
Resolución de sistema de ecuaciones lineales
Función Descripción
linsolve (A, b) Calcula todas las soluciones de: A×X =b
A\b Calcula todas las soluciones de: A×X =b
inv(A)× b Calcula todas las soluciones de: A×X =b
Ejemplo
Example
Dado el Sistema de ecuaciones lineales
x + 2y − 3z = 5 (1)
−3x − y − z = −8 (2)
x −y +z =0 (3)
a) Construya matrices para que el sistema sea descrito por A × X = b

b) Resolver con una sola linea de código.
c) Determine el rango de la matriz A y de la matriz aumentada [A|P]
Solución
1 clc
2 % Ingresando la Matriz A
3 A=[1 2 -3;-3 -1 -1;1 -1 1]
4 % Ingresando vector columna b
5 b=[5 ;-8; 0]
6 %Resolviendo AX=b
7 X=linsolve(A,b)
8 X=A\b
9 X=inv(A)*b
10 % Matriz Aumentada Aa
11 Aa=[A b]
12 % Hallando los rangos
13 rA= rank(A)
14 rAb= rank(Aa)
Ejemplo
Example
Implemente una función en Matlab que determine si el sistema A × X = b tiene
solución única, infinitas soluciones o no tiene solución. (Sugerencia. Utilice el
teorema de Rouche-Frobenius)
Solución
1 function [X] = RoucheFrobenius(A,b)

2 [n,m]= size(A);
3 rA=rank(A);
4 rAb=rank([A b]);
5 if (rA ∼= rAb)
6 fprintf('El sistema no tiene solucion ');
7 elseif (rA==n)
8 fprintf('El sistema tiene solucion unica \n');
9 X=inv(A)*b;
10 fprintf (' x1= %.2f\n',X(1));
11 fprintf (' x2= %.2f\n',X(2));
12 fprintf (' x3= %.2f\n',X(3));
13 else
14 fprintf ('El sistema tiene infinitas soluciones ');
15 end
16 end
2
PARTE PRÁCTICA
Aplicación de matrices y SEL
El sector agropecuario holandés, junto con la totalidad del complejo agroalimentario (producción,
transformación, comercialización), arrojó una producción bruta de algo más de 36000 millones de
florines en 1992, cifra que tras descontar la adquisición de los medios de producción, las
amortizaciones y diversa tasas e impuestos se convirtió en 17400 millones de valor agregado bruto y
casi 13000 millones de valor agregado neto. Es por tal motivo que un granjero desea preparar una
fórmula alimenticia para engordar ganado para ello deberá de modelar un sistema de ecuaciones
lineales. Dispone de maíz, desperdicios, alfalfa y cebada, cada uno con ciertas unidades de
ingredientes nutritivos, de acuerdo con la tabla siguiente:
Table: Unidades de ingredientes nutritivos por kg de cada alimento disponible

ALIMENTO
Ingrediente Requerimiento diario
Maíz Desperdicio Alfalfa Cebada
nutritivo Unidades / kg
Carbohidrato 80 15 35 60 230
Proteína 28 72 57 25 180
Vitaminas 20 20 12 20 80
Celulosa 50 10 20 60 160
Costo $ 18 5 7 20 −
PREGUNTAS:
1. Determine un sistema de ecuaciones lineales con 4 incognitas de la forma
P × X = Q donde el vector columna X representa los kilogramos necesarios
de cada ingrediente (en el orden y numéricamente igual al presentado en la
tabla) para satisfacer en su totalidad el requerimiento diario. De como
respuesta la traza de la matriz P.
2. Determine el rango de la matriz P y el rango de la matriz ampliada [P|Q].
3. Determine si el sistema P × X = Q tiene solución única, infinitas soluciones o
no tiene solución. (Sugerencia. Utilice la función creada: RoucheFrobenius ).
Resolver el sistema de ecuaciones lineales si es compatible determinado.
4. Determine el costo de la mezcla.
5. Otro granjero decide optimizar la formula original agregando 2 alimentos de
tipos A y B, estos nuevos alimentos contiene: Proteina, Vitaminas, Calcio y
Sales Minerales.
El de tipo A sólo contiene:
Proteina, 50% mas que lo que contiene el maíz.
Calcio, una proporción equivalente al de carbohidrato que contiene la Alfalfa.
25 unidades por kg en sales minerales.
El de tipo B sólo contiene:
Vitaminas, un 5% mas que lo que contiene la cebada.
Calcio, un equivalente de Proteina que contiene la cebada.
30 unidades por kg en sales minerales.
Estos nuevos productos aplican nuevos requerimientos diarios de ingredientes
nutritivos.
Complete la siguiente tabla con los datos mencionados.
Continuación...
ALIMENTO
Ingrediente Requerimiento diario
Maíz Desperdicio Alfalfa Cebada A B
nutritivo Unidades / kg
Carbohidrato 80 15 35 60 - - 400
Proteína 28 72 57 25 ? - 660
Vitaminas 20 20 12 20 - ? 220
Celulosa 50 10 20 60 - - 260
Calcio - - - - ? ? 120
Sales Minerales - - - - ? ? 120
6. Utilizando los datos del item anterior, determine un sistema de 6 ecuaciones
lineales con 6 incógnitas de la forma M × Y = R donde el vector columna Y
representa los kilogramos necesarios de cada ingrediente (en el orden y
numéricamente igual al presentado en la tabla) para satisfacer en su totalidad
el requerimiento diario.
Halle la traza de la matriz M.
Resuelve el sistema utilizando comandos de Matlab.
7. Utilizando comandos de Matlab, construir una submatriz de M de orden 3 × 3
denominada P1, considerando la 1era, 3era y 4ta fila de M, además la 1era,
2da y 4ta columana de M. Luego, realice las siguientes operaciones:
Halle a traza de P1
Halle el determinante de P1
Halle la transpuesta de P1
Halle la diagonal de P1
Intercambie la primera fila con la 3era fila en P1.
Halle la matriz triángular inferior de P1
Halle la inversa de P1
Construye una matriz rectangular de orden 3 × 9 : S=[P1 O I].
Donde:
O: es la matriz cuadrada nula.
I: es la matriz identidad.
Continuación...
Construye las siguiente matriz con comandos de Matlab
1 0 0 0 1 2 3 4
 
 0 2 0 0 5 6 7 8 
 
 0 0 3 0 9 10 11 12 
 
 0 0 0 4 20 0 5 4 
T =  1 5 9 20 0 0 0 0


 
 2 6 10 0 0 0 0 0 
 
 3 7 11 5 0 0 0 0 
4 8 12 4 0 0 0 0
Luego extraer la siguiente submatriz
 
3 0 9 10
 0 4 20 0 
R=  9

20 0 0 
10 0 0 0
ALAB01
Entregable: 06 de Mayo hasta las 11:59 pm
PRUEBA DE ESTANQUEIDAD EN VÁLVULAS
Dentro de la etapa de la construcción de una embarcación marina se realiza el

proceso de instalación de tuberias interiores. Un sistema de tuberias importante es
la descarga sanitaria (flujo de agua de los baños, cocina, lavanderia etc.), para ello
se realiza las pruebas de estanqueidad a los accesorios como son las válvulas,
necesitando hallar el radio interno de las mangueras L1 , L2 y L3 mediante un
sistema de ecuaciones lineales.
Continuación...
Esta prueba consiste en bombear agua desde un recipiente cilindrico hacia las
válvulas cerradas y unidas mediante una brida (Ver figura 1). En un inicio del
bombeo:
En la mangera L4 , la velocidad del flujo es 6.93v .
En las mangeras L1 ,L2 ,y L3 la velocidad del flujo es v .
Este dato ayuda en aplicar la conservación de la materia:
Qe = Qs
Donde:
Qe : Flujo que entra (Qe = ve × AreaTransversal), ve : velocidad del flujo de entrada.
Qs : Flujo que sale (Qs = vs × AreaTransversal), vs : velocidad del flujo de salida.
Continuación...
El bombeo se sigue realizando hasta alcanzar la presión P = 0.5 Mpa
Fuerza
P=
AreaTransversal
que permance constante por el agua ocupada en todas las mangueras.
Luego se aplica un tiempo determinado para revisar posibles fugas de agua en las
partes externas de las válvulas debido a la presión.
Continuación...
En el inicio de la prueba de estanqueidad el agua se encuentra lleno hasta el punto
de descarga del bombeo manual (El llenado del tanque inicialmente está a una
altura H), luego se conecta las mangueras L1 , L2 , L3 y L4 , que al ser llenados de
agua cae el nivel desde una altura H=1m hasta h= 0.9844m del recipiente cilindrico
de radio interno 0.2m.
Por información del taller se tiene las especificaciones físicas de las mangueras Li
como su longitud (li ) y radio interno (ri ), i = 1, . . . , 4:
l1 = 0.9m y r1 .
l2 = 0.8m y r2 .
l3 = 0.6m y r3 .
l4 = 1.2m y r4 = 0.01m.
Además, se tiene como dato técnico que el doble del área transversal de la
manguera L1 más el triple del área transversal de la manguera L2 más el área
trasnversal de la manguera L3 es igual a 12.87x10−4 π.
Figura 1.
PREGUNTAS:
1 (4 pts) Determine un sistema de 3 ecuaciones lineales con 3 incógnitas de la
forma P.X = Q donde el vector columna X representa el r 2 de cada manguera
L1 , L2 y L3 .
2 (4 pts) Determine el rango de la matriz P y el rango de la matriz ampliada
[P|Q].
3 (4 pts) Implemente una función en Matlab que determine si el sistema tiene
solución única, infinitas soluciones ó no tiene solución. (Sug. Utilice el teorema
de Rouché-Frobenius )
4 (4 pts) Resolver el sistema P.X=Q utilizando la inverza de la matriz P como
método de solución.
5 (4 pts) Al realizar la prueba de estanqueidad las válvulas son selladas por una
brida que presenta 4 pernos los cuales soportan 140N cada uno. Determine
que válvulas soportan la prueba de estanqueidad.
Gracias

Laboratorio 01

Cargado por

Copyright:

Formatos disponibles

Laboratorio 01

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

Laboratorio 01

Cargado por

Copyright:

Formatos disponibles

Matemáticas III

Profesores: H. Pantoja - B. Molina - R. Perez

a) Construya matrices para que el sistema sea descrito por A × X = b

1 function [X] = RoucheFrobenius(A,b)

Table: Unidades de ingredientes nutritivos por kg de cada alimento disponible

Dentro de la etapa de la construcción de una embarcación marina se realiza el

El bombeo se sigue realizando hasta alcanzar la presión P = 0.5 Mpa

También podría gustarte