library(MASS)

Definisikan matriks

A <- matrix(1:25, nrow=5, ncol=5)
B <- matrix(26:50, nrow=5, ncol=5)

Cetak matriks A dan B

print(A)
##      [,1] [,2] [,3] [,4] [,5]
## [1,]    1    6   11   16   21
## [2,]    2    7   12   17   22
## [3,]    3    8   13   18   23
## [4,]    4    9   14   19   24
## [5,]    5   10   15   20   25
print(B)
##      [,1] [,2] [,3] [,4] [,5]
## [1,]   26   31   36   41   46
## [2,]   27   32   37   42   47
## [3,]   28   33   38   43   48
## [4,]   29   34   39   44   49
## [5,]   30   35   40   45   50

Operasi penjumlahan matriks

C <- A + B
print("Hasil penjumlahan matriks A dan B:")
## [1] "Hasil penjumlahan matriks A dan B:"
print(C)
##      [,1] [,2] [,3] [,4] [,5]
## [1,]   27   37   47   57   67
## [2,]   29   39   49   59   69
## [3,]   31   41   51   61   71
## [4,]   33   43   53   63   73
## [5,]   35   45   55   65   75

Operasi pengurangan matriks

D <- A - B
print("Hasil pengurangan matriks A dan B:")
## [1] "Hasil pengurangan matriks A dan B:"
print(D)
##      [,1] [,2] [,3] [,4] [,5]
## [1,]  -25  -25  -25  -25  -25
## [2,]  -25  -25  -25  -25  -25
## [3,]  -25  -25  -25  -25  -25
## [4,]  -25  -25  -25  -25  -25
## [5,]  -25  -25  -25  -25  -25

Operasi perkalian matriks

E <- A %*% t(B)
print("Hasil perkalian matriks A dan transpose matriks B:")
## [1] "Hasil perkalian matriks A dan transpose matriks B:"
print(E)
##      [,1] [,2] [,3] [,4] [,5]
## [1,] 2230 2285 2340 2395 2450
## [2,] 2410 2470 2530 2590 2650
## [3,] 2590 2655 2720 2785 2850
## [4,] 2770 2840 2910 2980 3050
## [5,] 2950 3025 3100 3175 3250

Operasi determinan matriks

detA <- det(A)
print("Determinan matriks A:")
## [1] "Determinan matriks A:"
print(detA)
## [1] 0

Operasi transpose matriks

tA <- t(A)
print("Transpose matriks A:")
## [1] "Transpose matriks A:"
print(tA)
##      [,1] [,2] [,3] [,4] [,5]
## [1,]    1    2    3    4    5
## [2,]    6    7    8    9   10
## [3,]   11   12   13   14   15
## [4,]   16   17   18   19   20
## [5,]   21   22   23   24   25

Operasi eigenvalues dan eigenvectors pada matriks A

eigA <- eigen(A)
eigA$values
## [1]  6.864208e+01+0.000000e+00i -3.642081e+00+0.000000e+00i
## [3]  4.257350e-15+0.000000e+00i -1.270981e-16+4.588876e-16i
## [5] -1.270981e-16-4.588876e-16i
eigA$vectors
##              [,1]           [,2]           [,3]                  [,4]
## [1,] 0.3800509+0i -0.76703416+0i  0.54621260+0i  0.1175132+0.0459634i
## [2,] 0.4124552+0i -0.48590617+0i -0.27228461+0i  0.4017692+0.0065072i
## [3,] 0.4448594+0i -0.20477817+0i -0.66418830+0i -0.7435988+0.0000000i
## [4,] 0.4772637+0i  0.07634982+0i -0.03961996+0i -0.1881630-0.2033750i
## [5,] 0.5096680+0i  0.35747782+0i  0.42988027+0i  0.4124793+0.1509044i
##                       [,5]
## [1,]  0.1175132-0.0459634i
## [2,]  0.4017692-0.0065072i
## [3,] -0.7435988+0.0000000i
## [4,] -0.1881630+0.2033750i
## [5,]  0.4124793-0.1509044i