Scribd
Cargar
26 vistas11 páginas

Operaciones de Matrices y Análisis

Cargado por

Alexander Mel
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
26 vistas11 páginas

Operaciones de Matrices y Análisis

Cargado por

Alexander Mel
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

UNIVERSIDAD NACIONAL JORGE

BASADRE GROHOMAM

FACULTAD DE FAIN

ESCUELA PROFESIONAL
DE METALURGIA

Practica N° 04

Operación de matrices
Grupo: 02 (9-11)
Estudiantes: 2022-103023 Melchor Mamani Alexander Abel

Asignatura: Programación aplicada a procesos


Docente: Ing.Churata, Alexander

Tacna-Perú
2024
Practica Nº04

(Tema: Operaciones de matrices)

Alumno: Alexander Abel Melchor Mamani 2022-103023

1.Crea la siguiente matriz C:

c=[2 4 6 8 10;3 6 9 12 15;7 14 21 28 35]

c = 3×5
2 4 6 8 10
3 6 9 12 15
7 14 21 28 35

ua=c(:,3)

ua = 3×1
6
9
21

ub=c(2,:)

ub = 1×5
3 6 9 12 15

uc=[c(:,1),c(:,3),c(:,5)]

uc = 3×3
2 6 10
3 9 15
7 21 35

ud=[c(1,:),c(2,:)]

ud = 1×10
2 4 6 8 10 3 6 9 12 15

2.Cree la siguiente matriz A.

A=[1 2 3 4 5 6 7;2 4 6 8 10 12 14;21 18 15 12 9 6 3;5 10 15 20 25 30 35]

A = 4×7
1 2 3 4 5 6 7
2 4 6 8 10 12 14
21 18 15 12 9 6 3
5 10 15 20 25 30 35

B=A([1 3 4],[1 3 5 7])

B = 3×4
1 3 5 7
21 15 9 3
5 15 25 35

u=[A(3,:),A(:,5)',A(:,7)']

u = 1×15
21 18 15 12 9 6 3 5 10 9 25 7
14 ⋯

3.Utilizando las funciones Zeros, one y ey, cree los siguientes arrays.

a=zeros(5,5)

a = 5×5
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0

b=eye(4)

b = 4×4
1 0 0 0
0 1 0 0
0 0 1 0
0 0 0 1

c=ones(2)

c = 2×2
1 1
1 1

4.Utilizando el comando eye, cree un array. como el que se muestra màs abajo en la parte
izquierda. a continuacion, utilice los puntos para acceder a los elementos en el arrays,
cambiando el arrays para que sea como el de la parye derecha.

esme=eye(7)

esme = 7×7
1 0 0 0 0 0 0
0 1 0 0 0 0 0
0 0 1 0 0 0 0
0 0 0 1 0 0 0
0 0 0 0 1 0 0
0 0 0 0 0 1 0
0 0 0 0 0 0 1

esme(1:2,1:3)=2

esme = 7×7
2 2 2 0 0 0 0
2 2 2 0 0 0 0
0 0 1 0 0 0 0
0 0 0 1 0 0 0
0 0 0 0 1 0 0
0 0 0 0 0 1 0
0 0 0 0 0 0 1

esme(1:3,5:7)=5

esme = 7×7
2 2 2 0 5 5 5
2 2 2 0 5 5 5
0 0 1 0 5 5 5
0 0 0 1 0 0 0
0 0 0 0 1 0 0
0 0 0 0 0 1 0
0 0 0 0 0 0 1

esme(3,1:3)=3

esme = 7×7
2 2 2 0 5 5 5
2 2 2 0 5 5 5
3 3 3 0 5 5 5
0 0 0 1 0 0 0
0 0 0 0 1 0 0
0 0 0 0 0 1 0
0 0 0 0 0 0 1

esme(4,4)=6

esme = 7×7
2 2 2 0 5 5 5
2 2 2 0 5 5 5
3 3 3 0 5 5 5
0 0 0 6 0 0 0
0 0 0 0 1 0 0
0 0 0 0 0 1 0
0 0 0 0 0 0 1

esme(5,1:3)=4

esme = 7×7
2 2 2 0 5 5 5
2 2 2 0 5 5 5
3 3 3 0 5 5 5
0 0 0 6 0 0 0
4 4 4 0 1 0 0
0 0 0 0 0 1 0
0 0 0 0 0 0 1

esme(5,4:6)=7

esme = 7×7
2 2 2 0 5 5 5
2 2 2 0 5 5 5
3 3 3 0 5 5 5
0 0 0 6 0 0 0
4 4 4 7 7 7 0
0 0 0 0 0 1 0
0 0 0 0 0 0 1

esme(6,4:7)=9

esme = 7×7
2 2 2 0 5 5 5
2 2 2 0 5 5 5
3 3 3 0 5 5 5
0 0 0 6 0 0 0
4 4 4 7 7 7 0
0 0 0 9 9 9 9
0 0 0 0 0 0 1

esme(7,4:7)=9

esme = 7×7
2 2 2 0 5 5 5
2 2 2 0 5 5 5
3 3 3 0 5 5 5
0 0 0 6 0 0 0
4 4 4 7 7 7 0
0 0 0 9 9 9 9
0 0 0 9 9 9 9

5.Realizar operaciones matematicas (+-*/...) con arrays.

a1=[6 7 3 8 1 0;3 8 3 8 5 6;4 6 8 5 4 4;5 7 4 2 2 9]

a1 = 4×6
6 7 3 8 1 0
3 8 3 8 5 6
4 6 8 5 4 4
5 7 4 2 2 9

b1=[5 8 3 7 9 4;6 3 6 1 8 2;5 7 9 4 6 3;3 7 6 8 5 7]

b1 = 4×6
5 8 3 7 9 4
6 3 6 1 8 2
5 7 9 4 6 3
3 7 6 8 5 7

• Suma de las matrices

sumamos=a1+b1

sumamos = 4×6
11 15 6 15 10 4
9 11 9 9 13 8
9 13 17 9 10 7
8 14 10 10 7 16
• Multiplicaciòn elemento a elemento

multiplicacion_elemento = a1 .* b1

multiplicacion_elemento = 4×6
30 56 9 56 9 0
18 24 18 8 40 12
20 42 72 20 24 12
15 49 24 16 10 63

• Transpuestas de las matrìces

transpuesta_a1 = a1'

transpuesta_a1 = 6×4
6 3 4 5
7 8 6 7
3 3 8 4
8 8 5 2
1 5 4 2
0 6 4 9

transpuesta_b1 = b1'

transpuesta_b1 = 6×4
5 6 5 3
8 3 7 7
3 6 9 6
7 1 4 8
9 8 6 5
4 2 3 7

• Media

media_a1 = mean(a1(:))

media_a1 = 4.9167

media_b1 = mean(b1(:))

media_b1 = 5.5000

• Desviaciòn estàndar

desviacion_estandar_a1 = std(a1(:))

desviacion_estandar_a1 = 2.4480

desviacion_estandar_b1 = std(b1(:))

desviacion_estandar_b1 = 2.2458
• Maximo y minimo

maximo_a1 = max(a1(:))

maximo_a1 = 9

minimo_a1 = min(a1(:))

minimo_a1 = 0

• Suma de columna (individual)

suma_columnas_a1 = a1(:,1)+b1(:,1)

suma_columnas_a1 = 4×1
11
9
9
8

• Suma de filas (Individual

suma_filas_a1 =a1(2,:)+b1(2,:)

suma_filas_a1 = 1×6
9 11 9 9 13 8

6.Realizar operaciones con operadores relacionales y lògicos.

mayor_a1_b1 = a1 > b1

mayor_a1_b1 = 4×6 logical array


1 0 0 1 0 0
0 1 0 1 0 1
0 0 0 1 0 1
1 0 0 0 0 1

menor_a1_b1 = a1 < b1

menor_a1_b1 = 4×6 logical array


0 1 0 0 1 1
1 0 1 0 1 0
1 1 1 0 1 0
0 0 1 1 1 0

igual_a1_b1 = a1 == b1

igual_a1_b1 = 4×6 logical array


0 0 1 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 1 0 0 0 0

resultado_logico = (a1 > 5) & (b1 < 5)

resultado_logico = 4×6 logical array


0 0 0 0 0 0
0 1 0 1 0 1
0 0 0 0 0 0
0 0 0 0 0 0

resultado_logico2 = (a1 < 5) | (b1 > 5)

resultado_logico2 = 4×6 logical array


0 1 1 1 1 1
1 0 1 0 1 0
1 1 1 0 1 1
0 1 1 1 1 1

negacion_a1 = ~(a1 == 0)

negacion_a1 = 4×6 logical array


1 1 1 1 1 0
1 1 1 1 1 1
1 1 1 1 1 1
1 1 1 1 1 1

7.

TAREA PARA CASA:

1. ZARANDA SIZETEC

etapas = ["Under size"; "Over size";"Cabeza calculada"];


peso_TM_DIA = [132.93; 45.82];
suma_columna1 = sum(peso_TM_DIA);
peso_Porc=(100*(peso_TM_DIA))/suma_columna1 ;
suma_columna2 = sum(peso_Porc);
leyes_PorcSn = [21.41; 46.88];
finos_TM =(peso_TM_DIA .* leyes_PorcSn)/100;
suma_columna4 = sum(finos_TM);
suma_columna3 =(suma_columna4*100)/ suma_columna1;
recuperacion_Sn =(finos_TM *100)/ suma_columna4;
suma_columna5=sum(recuperacion_Sn);

peso_TM_DIA=[peso_TM_DIA;suma_columna1];
peso_Porc=[peso_Porc;suma_columna2];
leyes_PorcSn=[leyes_PorcSn;suma_columna3];
finos_TM=[finos_TM;suma_columna4];
recuperacion_Sn=[recuperacion_Sn;suma_columna5];
% Crear una tabla
tabla = table(etapas, peso_TM_DIA, peso_Porc, leyes_PorcSn, finos_TM,
recuperacion_Sn, ...
'VariableNames', {'Etapas', 'Peso_TM/DIA', 'Peso_%', 'Leyes_%Sn',
'Finos_TM', 'Recuperacion_Sn'})

tabla = 3×6 table


Etapas Peso_TM/DIA Peso_% Leyes_%Sn Finos_TM ⋯
1 "Under size" 132.9300 74.3664 21.4100 28.4603
2 "Over size" 45.8200 25.6336 46.8800 21.4804
3 "Cabeza calculada" 178.7500 100 27.9389 49.9407

2. Análisis Granulométrico Molienda Primaria

Malla=["5/8";"½ ";"3/8 ";"¼


";"4";"10";"14";"18";"30";"35";"40";"50";"70";"100";"140";"200";"270";"325";"400
";"<400"];
Peso_gr=[35.70;64.70;100.20;275.70;194.10;209.10;56.00;41.00;54.10;14.30;13.50;
26.00;24.10;19.90;8.70;19.50;13.80;4.80;3.80;68.10];
su_columna1=sum(Peso_gr);
peso_porctanje=(Peso_gr/su_columna1).*100;
su_columna2=sum(peso_porctanje);
acumulado=zeros(size(peso_porctanje));
acumulado(1)=peso_porctanje(1);
for i = 2:length(peso_porctanje)
acumulado(i)=acumulado(i-1) + peso_porctanje(i);
end
disp(acumulado);

2.8626
8.0507
16.0853
38.1926
53.7567
70.5236
75.0140
78.3017
82.6397
83.7864
84.8689
86.9537
88.8862
90.4819
91.1795
92.7432
93.8497
94.2346
94.5393
100.0000

pasante=100-acumulado;
tabla = table(Malla, Peso_gr, peso_porctanje, acumulado,pasante, ...
'VariableNames', {'Malla', 'Peso (gr)', 'Peso_%', '% Acumulado',
'Pasante'})

tabla = 20×5 table


Malla Peso (gr) Peso_% % Acumulado Pasante

1 "5/8" 35.7000 2.8626 2.8626 97.1374


2 "½ " 64.7000 5.1880 8.0507 91.9493
3 "3/8 " 100.2000 8.0346 16.0853 83.9147
4 "¼ " 275.7000 22.1073 38.1926 61.8074
5 "4" 194.1000 15.5641 53.7567 46.2433
6 "10" 209.1000 16.7669 70.5236 29.4764
7 "14" 56 4.4904 75.0140 24.9860
8 "18" 41 3.2876 78.3017 21.6983
9 "30" 54.1000 4.3381 82.6397 17.3603
10 "35" 14.3000 1.1467 83.7864 16.2136
11 "40" 13.5000 1.0825 84.8689 15.1311
12 "50" 26 2.0848 86.9537 13.0463
13 "70" 24.1000 1.9325 88.8862 11.1138
14 "100" 19.9000 1.5957 90.4819 9.5181
15 "140" 8.7000 0.6976 91.1795 8.8205
16 "200" 19.5000 1.5636 92.7432 7.2568
17 "270" 13.8000 1.1066 93.8497 6.1503
18 "325" 4.8000 0.3849 94.2346 5.7654
19 "400" 3.8000 0.3047 94.5393 5.4607
20 "<400" 68.1000 5.4607 100 -0

3.CALCULO DEL D50

Micras=[600;500;450;300;212;150;106;75;53;44;33;22];
Feed_P=[12.35;3.53;5.29;12.19;16.72;14.93;7.29;7.56;4.49;1.37;0.97;13.29];
UF_P=[7.24;3.44;5.80;15.42;22.99;20.84;9.29;8.20;3.11;0.66;0.38;2.63];
OF_P=[0.20;0.16;0.20;2.59;5.34;7.69;6.06;11.80;11.92;4.62;3.79;45.64];

Mineral_t_h=[39.6,33.4];
resta=(Mineral_t_h(1)-Mineral_t_h(2));
Mineral_t_h=[Mineral_t_h,resta];
Feed_t_h=[(Mineral_t_h(1).*Feed_P)/100];
UF_t_h=[(Mineral_t_h(2).*UF_P)/100];
OF_t_h=[(Mineral_t_h(3).*OF_P)/100];
UF_Rend=[UF_t_h ./ Feed_t_h .*100];

Solidos=[35.2,77.6,25.4]

Solidos = 1×3
35.2000 77.6000 25.4000

Agua_t_h=[(100-Solidos) .* Mineral_t_h ./ Solidos];

tabla = table(Micras, Feed_P, UF_P,OF_P,Feed_t_h,UF_t_h,OF_t_h,UF_Rend,...


'VariableNames', {'Micras', 'Feed %P', 'UF %P','F %P','Feed
t/h','UF t/h','OF t/h','UF Rend'})

tabla = 12×8 table


Micras Feed %P UF %P F %P Feed t/h UF t/h ⋯
1 600 12.3500 7.2400 0.2000 4.8906 2.4182
2 500 3.5300 3.4400 0.1600 1.3979 1.1490
3 450 5.2900 5.8000 0.2000 2.0948 1.9372
4 300 12.1900 15.4200 2.5900 4.8272 5.1503
5 212 16.7200 22.9900 5.3400 6.6211 7.6787
6 150 14.9300 20.8400 7.6900 5.9123 6.9606
7 106 7.2900 9.2900 6.0600 2.8868 3.1029
8 75 7.5600 8.2000 11.8000 2.9938 2.7388
9 53 4.4900 3.1100 11.9200 1.7780 1.0387
10 44 1.3700 0.6600 4.6200 0.5425 0.2204
11 33 0.9700 0.3800 3.7900 0.3841 0.1269
12 22 13.2900 2.6300 45.6400 5.2628 0.8784

tabla = table(Mineral_t_h, Solidos, Agua_t_h, ...


'VariableNames', {'Mineral t/h', '% Sólidos', 'Agua t/h'})

tabla = 1×3 table


Mineral t/h ⋯
1 39.6000 33.4000 6.2000

También podría gustarte