1.
Para una Matriz A2x2, usando MATLAB hallar:
Sea la matriz
>> format rational
A=[ 2, 4; 5, 6]
A=
2 4
5 6
a) La matriz inversa
>> inv(A)
ans =
-3/4 1/2
5/8 -1/4
b) La matriz traspuesta
>> A'
ans =
2 5
4 6
c) La matriz simétrica asociada
>> S=(A+A')/2
S=
2 9/2
9/2 6
d) La matriz anti simétrica asociada
>> N=(A-A')/2
N=
0 -1/2
1/2 0
e) El cuadrado de la matriz
>> A*A
ans =
24 32
40 56
f) La Factorización A = Pt L U
>> [L,U,P]=lu(A)
L=
1 0
2/5 1
U=
5 6
0 8/5
P=
0 1
1 0
g) El determinante de la matriz
>> det(A)
ans =
-8
h) El determinante de la matriz Traspuesta
>> det(A')
ans =
-8
i) El valor S=|A|−|A '|
>> s=det(A)-det(A')
s=
1/562949953421312
j) El determinante de la matriz inversa
>> det(inv(A))
ans =
-1/8
k) Polinomio característico
>> poly(A)
ans =
1 -8 -8
l) Autovalores
>> eig(A)
ans =
-881/980
8721/980
m,n) Autovectores
>> [ V,D]=eig(A)
V=
-417/515 -1265/2522
615/1048 -481/556
D=
-881/980 0
0 8721/980
o) Verificar A = P D P-1
>> V*D*inv(A)
ans =
-3309/992 1653/1117
-2049/464 191/115
50
p) La matriz A
>> V*(D^50)*inv(V)
ans =
8.672010722065207e+46 1.196560481402787e+47
1.4957006017534832e+47 2.063761553609307e+47
q) La matriz inversa usando el teorema de Cayley Hamilton
>> (A^2-3*A)/(-2)
ans =
-9 -10
-25/2 -19
[Link] Una Matriz A3x3, usando Matlab hallar
Sea la matriz
>> format rational
A=[ 3, 0, 1; 2, 0, -3; 1, -1, 0]
A=
3 0 1
2 0 -3
1 -1 0
a) La matriz inversa
>> inv(A)
ans =
3/11 1/11 0
3/11 1/11 -1
2/11 -3/11 0
b) La matriz traspuesta
>> A'
ans =
3 2 1
0 0 -1
1 -3 0
c) La matriz simétrica asociada
>> S=(A+A')/2
S=
3 1 1
1 0 -2
1 -2 0
d) La matriz anti simétrica asociada
>> N=(A-A')/2
N=
0 -1 0
1 0 -1
0 1 0
e) El cuadrado de la matriz
>> A*A
ans =
10 -1 3
3 3 2
1 0 4
f) La Factorización A = Pt L U
>> [L,U,P]=lu(A)
L=
1 0 0
1/3 1 0
2/3 0 1
U=
3 0 1
0 -1 -1/3
0 0 -11/3
P=
1 0 0
0 0 1
0 1 0
g) El determinante de la matriz
>> det(A)
ans =
-11
h) El determinante de la matriz Traspuesta
>> det(A')
ans =
-11
i) El valor S=|A|−|A '|
>> s=det(A)-det(A')
s=
0
j) El determinante de la matriz inversa
>> det(inv(A))
ans =
-1/11
k) Polinomio característico
>> poly(A)
ans =
1 -3 -4 11
l) Autovalores
>> eig(A)
ans =
-914/469
2422/765
747/419
m,n) Autovectores
>> [ V,D]=eig(A)
V=
369/3694 -4572/5117 -665/2354
-601/696 -593/1399 -60/67
-481/973 -871/5872 316/919
D=
-914/469 0 0
0 2422/765 0
0 0 747/419
o) Verificar A = P D P-1
>> V*D*inv(A)
ans =
-1497/1634 -345/2509 1355/479
-635/3218 361/774 310/231
413/1678 -334/2731 541/1152
50
p) La matriz A
>> V*(D^50)*inv(V)
ans =
1.0531211270482087e+25 -1.4993958743047412e+24 4.7471063494941303e+24
4.996025076074039e+24 -7.113160293908207e+23 4.7471063494941303e+24
4.996025076074039e+24 -2.4891872657990684e+23 7.880798449139212e+23
q) La matriz inversa usando el teorema de Cayley Hamilton
>> format rational
(A^2-3*A)/(-2)
ans =
-1/2 1/2 0
3/2 -3/2 -11/2
1 -3/2 -2
3. Para una matriz 4x4,, usando MATAB hallar
Sea la matriz
>> format rational
A=[ 3, 0, 2, 2; 2, 0, -2, 0; 1, -1, 0, 1; 0, 1, 0, 1]
A=
3 0 2 2
2 0 -2 0
1 -1 0 1
0 1 0 1
a) La matriz inversa
>> inv(A)
ans =
1/4 1/4 -1/4 -1/4
1/8 1/8 -5/8 3/8
1/4 -1/4 -1/4 -1/4
-1/8 -1/8 5/8 5/8
b) La matriz traspuesta
>> A'
ans =
3 2 1 0
0 0 -1 1
2 -2 0 0
2 0 1 1
c) La matriz simétrica asociada
>> S=(A+A')/2
S=
3 1 3/2 1
1 0 -3/2 1/2
3/2 -3/2 0 1/2
1 1/2 1/2 1
d) La matriz anti simétrica asociada
>> N=(A-A')/2
N=
0 -1 1/2 1
1 0 -1/2 -1/2
-1/2 1/2 0 1/2
-1 1/2 -1/2 0
e) El cuadrado de la matriz
>> A*A
ans =
11 0 6 10
4 2 4 2
1 1 4 3
2 1 -2 1
f) La Factorización A = Pt L U
>> [L,U,P]=lu(A)
L=
1 0 0 0
1/3 1 0 0
2/3 0 1 0
0 -1 1/5 1
U=
3 0 2 2
0 -1 -2/3 1/3
0 0 -10/3 -4/3
0 0 0 8/5
P=
1 0 0 0
0 0 1 0
0 1 0 0
0 0 0 1
g) El determinante de la matriz
>> det(A)
ans =
-16
h) El determinante de la matriz Traspuesta
>> det(A')
ans =
-16
i) El valor S=|A|−|A '|
>> s=det(A)-det(A')
s=
0
j) El determinante de la matriz inversa
>> det(inv(A))
ans =
-1/16
k) Polinomio característico
>> poly(A)
ans =
1 -4 -1 12 -16
l) Autovalores
>> eig(A)
ans =
-1828/945 + 0i
861/232 + 0i
1803/1622 + 1155/1159i
1803/1622 - 1155/1159i
m,n) Autovectores
>> [ V,D]=eig(A)
V=
-1074/7903 + 0i 1473/1648 + 0i 226/12731 + 697/1546i 226/12731 - 697/1546i
941/1251 + 0i 418/1079 + 0i 1802/3093 + 0i 1802/3093 + 0i
1540/2603 + 0i 1115/6373 + 0i -965/3153 + 201/1252i -965/3153 - 201/1252i
-1102/4299 + 0i 689/4822 + 0i 212/3279 - 1205/2087i 212/3279 + 1205/2087i
D=
-1828/945 + 0i 0 + 0i 0 + 0i 0 + 0i
0 + 0i 861/232 + 0i 0 + 0i 0 + 0i
0 + 0i 0 + 0i 1803/1622 + 1155/1159i 0 + 0i
0 + 0i 0 + 0i 0 + 0i 1803/1622 - 1155/1159i
o) Verificar A = P D P-1
>> V*D*inv(A)
ans =
6909/16193 + 938/4821i 356/555 - 253/3901i -23/10 - 1623/3575i 713/701 - 1623/3575i
-393/3812 + 371/1704i -549/1286 - 656/9039i -388/1329 - 4402/8665i 959/837 - 4402/8665i
-643/2404 - 970/20441i -37/2130 + 337/21305i -401/1305 + 255/2303i 538/1573 + 255/2303i
958/3533 - 464/2143i -52/991 + 277/3838i -1870/8793 + 1018/2015i 559/1760 + 1018/2015i
p) La matriz A50
>> V*(D^50)*inv(V)
ans =
2.3181848693186405e+28 2.7337656794522306e+27 1.1019631330313945e+28 2.11652098670043e+28
1.0047490581147114e+28 1.184870338811027e+27 4.776135133329127e+27 9.17343778752352e+27
4.537674915990142e+27 5.351143523550332e+26 2.1570110879778575e+27 4.142932825430858e+27
3.7059064286190605e+27 4.370263968117975e+26 1.7616249302853988e+27 3.3835216659082234e+2
q) La matriz inversa usando el teorema de Cayley Hamilton
>> format rational
(A^2-3*A)/(-2)
ans =
-1 0 0 -2
1 -1 -5 -1
1 -2 -2 0
-1 1 1 1