Universidad Mayor de “San Simón”

Facultad de Ciencias y Tecnología

Practica de Matlab

Universitario(a): Dorian Ojeda Quispe

Carrera: Ingeniería Industrial
Docente: Julio Medina Gamboa
Fecha: 20 de junio de 2023

Cochabamba - Bolivia

%Realice las operaciones indicadas con:

a=[1 3;2 b=[-2 0;1 4;-7 5];
IB
% Realizar las operaciones indicadas con:

a=[1 3;2 5;-1 2];

b=[-2 0;1 4;-7 5];
o=zeros(3,2);
p=ones(3,2);
c=[-1 1;4 6;-7 3];
% a)3a
res1=3*a

res1 =

3 9
6 15
3 6

% b) a + b
res2=a+b

res2 =

-1 3
3 9
-6 7

% c) c - a
res3=c-a

res3 =

-2 -2
2 1
-8 1

%d) a - c
res4=c-a

res4 =

-2 -2
2 1
-8 1

% e) 2c - 5a
res5=2*c-5*a

res5 =

-7 -13
-2 -13
-19 -4
% f) 0b
res6=0*b

res6 =

0 0
0 0
0 0

% g) b - a - 2c
res7=b-a-2*c

res7 =

-1 -5
-9 -13
6 -3

% h) 7c - b + 2a
res8=7*c-b+2*a

res8 =

-3 13
31 48
-40 20

% i) Calcular d tal que 2a + b - d = 0

d=2*a+b

d =

0 6
5 14
-5 9

% j) calcular tal que a + 2b - 3c +e = 0

e=-a-2*b+3*c

e =

0 0
8 5
-8 -3

% k) Calcular f tal que a + 2b + 3f = p

f=(1/3)*(p-a-2*b)

f =
 1.3333 -0.6667
-1.0000 -4.0000
4.6667 -3.6667

% Considere las siguientes matrices por bloques

c=[1 -1;2 0];
d=[2 4;4 5];
e=[1 1;-2 3];
f=[2 -3;5 0];
A=[c d;e f]

A =

1 -1 2 4
2 0 4 5
1 1 2 -3
-2 3 5 0

g=[1 4;2 -1];

h=[3;0];
j=[-3 2;0 1];
k=[1;2];
B=[g h;j k]

B =

1 4 3
2 -1 0
-3 2 1
0 1 2

% Comprobar que AB=([c d;e f]*[g h;j k])=[c*g+d*j c*h+d*k;e*g+f*j e*h+f*k]

% Solucion
X=A*B

X =

-7 13 13
-10 21 20
-3 4 -1
-11 -1 -1

X1=[c d;e f]*[g h;j k]

X1 =

-7 13 13
-10 21 20
-3 4 -1
-11 -1 -1
X2=[c*g+d*j c*h+d*k;e*g+f*j e*h+f*k]

X2 =

-7 13 13
-10 21 20
-3 4 -1
-11 -1 -1

% Determinar si el siguiente sistema tiene solucion:

% 4x -5y= 3
% 6x +7y= 4
% Solucion
A=[4 -5;6 7]

A =

4 -5
6 7

b=[3;4]

b =

3
4

%X=[x;y]
X=inv(A)*b

X =

0.7069
-0.0345

%x=0.7069 y=-0.345

% En los sistemas de ecuaciones dados, primero escribala de la forma Ax=b,

despues calcule A^(-1)x, por ultimo use la multiplicacion de matrices para
obtener el vector solucion.
% a) x-3y=4 2x+5y=7
A=[1 -3;2 5];
b=[4;7];
x=inv(A)*b

x =

3.7273
-0.0909

clear
% b) x-3y=4 2x+2y=3
A=[1 -3;2 2];
b=[4;3];
x=inv(A)*b

x =

2.1250
-0.6250

clear
% c) x+2y=3 2x+y-z=-1 3x+y+z=7
A=[1 2 0;2 1 -1;3 1 1];
b=[3;-1;7];
x=inv(A)*b

x =

0.7500
1.1250
3.6250

clear
% c) 2x+4z=7 -x+3y+z=-4 y+2z=5
A=[2 0 4;-1 3 1;0 1 2];
b=[7;-4;5];
x=inv(A)*b

x =

-6.8333
-5.3333
5.1667

%Una matriz A de n×n es normal si AA^(T)=A^(T)A. Comprobar que la siguiente

matriz es normal: [3 -1;1 3]
A=[3 -1;1 3]

A =

3 -1
1 3

A*[3 -1;1 3]

ans =

8 -6
6 8

[3 -1;1 3]*A
ans =

8 -6
6 8

% La matriz A es normal

% Una matriz cuadrada A de n×n es ortogonal si A^(-1)= A^(T). Decir si las

siguientes matrices son ortogonales o no:
% a) A=(1/3)*[1 -2 2;-2 1 2;-2 -2 1]
A=(1/3)*[1 -2 2;-2 1 2;-2 -2 1];

AT =

0.3333 -0.6667 -0.6667

-0.6667 0.3333 -0.6667
0.6667 0.6667 0.3333

inv(A)

ans =

0.7143 -0.2857 -0.8571

-0.2857 0.7143 -0.8571
0.8571 0.8571 -0.4286

% A no es ortogonal
% b) B=[3/5 -4/5 0;4/5 3/5 0;0 0 1]
B=[3/5 -4/5 0;4/5 3/5 0;0 0 1]

B=[3/5 -4/5 0;4/5 3/5 0;0 0 1];

BT=[3/5 4/5 0;-4/5 3/5 0;0 0 1]

BT =

0.6000 0.8000 0
-0.8000 0.6000 0
0 0 1.0000

inv(B)

ans =

0.6000 0.8000 0
-0.8000 0.6000 0
0 0 1.0000

% B es ortogonal
%Calcular (A^T)^-1 , (A^-1)^T y luego decir si son iguales o no
% a) A=[1 2;34]
% a) A=[1 2;3 4]
A=[1 2;3 4];
AT=[1 3;2 4];
inv(AT)

ans =

-2.0000 1.5000
1.0000 -0.5000

inv(A)

ans =

-2.0000 1.0000
1.5000 -0.5000

% (A^-1)T=
[-2 1.5;1 -0.5]

ans =

-2.0000 1.5000
1.0000 -0.5000

% Si son iguales
clear
% b) A=[2 0;6 3]
A=[2 0;6 3];
AT=[2 6;0 3];
inv(AT)

ans =

0.5000 -1.0000
0 0.3333

inv(A)

ans =

0.5000 0
-1.0000 0.3333

% (A^-1)^T=
[0.5 0;0 0.3333333333]

ans =
 0.5000 0
0 0.3333

[0.5 -1;0 0.3333333333]

ans =

0.5000 -1.0000
0 0.3333

% Si son iguales
% c) A=[2 1;3 2]
A=[2 1;3 2];
AT=[2 3;1 2];
inv(AT)

ans =

2 -3
-1 2

inv(A)

ans =

2.0000 -1.0000
-3.0000 2.0000

% (A^-1)^T=
[2 -3;-1 2]

ans =

2 -3
-1 2

% Si son lo mismo
clear

% d) A=[3 2 1;0 2 2;0 0 -1]

A=[3 2 1;0 2 2;0 0 -1];

AT=[3 0 0;2 2 0;1 2 -1];

inv(AT)

ans =

0.3333 0 0
-0.3333 0.5000 0
-0.3333 1.0000 -1.0000
inv(A)

ans =

0.3333 -0.3333 -0.3333

0 0.5000 1.0000
0 0 -1.0000

%(A^-1)^T=
[0.33333 0 0;-0.33333 0.5 0;-0.333333 1 -1]

ans =

0.3333 0 0
-0.3333 0.5000 0
-0.3333 1.0000 -1.0000

% Si son lo mismo
clear
% e) A=[1 1 1;0 2 3;5 5 1]
A=[1 1 1;0 2 3;5 5 1];
AT=[1 0 5;1 2 5;1 3 1];
inv(AT)

ans =

1.6250 -1.8750 1.2500

-0.5000 0.5000 0
-0.1250 0.3750 -0.2500

inv(A)

ans =

1.6250 -0.5000 -0.1250

-1.8750 0.5000 0.3750
1.2500 0 -0.2500

%(A^-1)^T=
[1.625 -1.875 1.25;-0.5 0.5 0;-0.125 0.375 -0.25]

ans =

1.6250 -1.8750 1.2500

-0.5000 0.5000 0
-0.1250 0.3750 -0.2500

Si son lo mismo

%Si A=[1 1 1;0 1 1;0 0 1], hallar A^2, A^3 y comprobar que: A^3-3A^2+3A=I
A=[1 1 1;0 1 1;0 0 1];
% A^2=
A*A

ans =

1 2 3
0 1 2
0 0 1

% A^3=
A*A*A

ans =

1 3 6
0 1 3
0 0 1

% A^3-3A^2+3A=

A*A*A-3*A*A+3*A

ans =

1 0 0
0 1 0
0 0 1

% Calcular
% a) A^-1X=A; A=[2 3;-1 2]
% X=
A=[2 3;-1 2];
% X=
A*A

ans =

1 12
-4 1

clear
% b) BXA=B; B=[4 -1;0 6], A=[2 3;-1 2]
B=[4 -1;0 6];
A=[2 3;-1 2];
% X=
inv(B)*B*inv(A)

ans =

0.2857 -0.4286
0.1429 0.2857
clear
% c) AX-3B=2X+A-B; A=[2 3;-1 2], B=[4 -1;0 6]
A=[2 3;-1 2];
B=[4 -1;0 6];
% X=

inv(A-2*eye(2))*(A+2*B)

ans =

1.0000 -14.0000
3.3333 0.3333

clear ans
% d) AXB+2BXB=CXB+CAB; C=[0 -1;1 0]
C=[0 -1;1 0];
% X=
inv(A+2*B-C)*(C*A*B)*inv(B)

ans =

0.0694 -0.2361
0.1528 0.1806

clear ans
% e) AX-3AB=X-3I
% X=

inv(A-eye(2))*(3*A*B-3*eye(2))

ans =

14.2500 -15.0000
2.2500 21.0000

clear
% f) XA+B=XC-D; A=[-8 1 3;-2 5 -2;4 -4 -2], B=[3 -7 6;1 4 -3;2 6 7], C=[3 -1
1;2 5 -3;-2 -3 -1], D=[-2 8 -5;-1 -4 3;-3 -7 -8]
A=[-8 1 3;-2 5 -2;4 -4 -2];
B=[3 -7 6;1 4 -3;2 6 7];
C=[3 -1 1;2 5 -3;-2 -3 -1];
D=[-2 8 -5;-1 -4 3;-3 -7 -8];
% X=
(-B-D)*inv(A-C)

ans =

-7.0000 0 -13.0000
0 0 0
7.0000 0 13.0000

% Si A=[0 -1;1 1], determine el primer entero positivo k tal que A^k=I
A=[0 -1;1 1];
A^2
ans =

-1 -1
1 0

A^3

ans =

-1 0
0 -1

A^4

ans =

0 1
-1 -1

A^5

ans =

1 1
-1 0

A^6

ans =

1 0
0 1

% k=6

% Encuentre una matriz X tal que X*[2 3;1 2]=[1 0;0 1]

% X=
eye(2)*inv([2 3;1 2])

ans =

2 -3
-1 2

% Sea A=[2 6;8 -6], encuentre un vector no nulo b=[x;y] tal que Ab=6b
A=[2 6;8 -6];
% b=
inv(6*eye(2)-A)*eye(2,1)
Warning: Matrix is singular to working precision.

ans =

NaN
NaN

6*eye(2)-A

ans =

4 -6
-8 12

inv(6*eye(2)-A)
Warning: Matrix is singular to working precision.

ans =

Inf Inf
Inf Inf

% La matriza resultante de 6*eye(2)-A no es invertible

% Realizar los calculos indicados para A=[1 -1 2;3 4 5;0 1 -1], B=[0 2 1;3 0
5;7 -6 0], C=[0 0 2;3 1 0;0 -2 4]
% a) 4C-2B+3A
A=[1 -1 2;3 4 5;0 1 -1];
B=[0 2 1;3 0 5;7 -6 0];
C=[0 0 2;3 1 0;0 -2 4];
% a) 4C-2B+3A
4*C-2*B+3*A

ans =

3 -7 12
15 16 5
-14 7 13

% b) 4c-2B+3A
4*C-2*B+3*A

ans =

3 -7 12
15 16 5
-14 7 13
% c) Calcular D tal que A+B+C+D es la matriz cero de 3×3
% D=
-A-B-C

ans =

-1 -1 -5
-9 -5 -10
-7 7 -3

% d) Calcular E tal que 3C-2B+8A-4E

% E=
(3*C-2*B+8*A)*(1/4)

ans =

2.0000 -3.0000 5.0000

6.7500 8.7500 7.5000
-3.5000 3.5000 1.0000

% e) Calcular F tal que A+2B+F-3C sea la matriz de 3×3 cuyos elementos son
todos uno.
% F=
ones(3)+3*C-2*B-A

ans =

0 -2 3
1 0 -14
-13 6 14