BAB 2 Example 69

Example 69—-Di R, jika kita ingin menghitung invers dari matriks yang bisa kita gunakan fungsi inv() dari paket matlib. Misalnya, jika kita memiliki matriks

library(matlib)
A <- matrix(c(1,-2,-1,2,3,2,3,-2,1), nrow = 3, ncol = 3)
inv(A)

##             [,1]       [,2]       [,3]
## [1,]  0.58333333  0.3333333 -1.0833333
## [2,]  0.33333333  0.3333333 -0.3333333
## [3,] -0.08333333 -0.3333333  0.5833333

Kemudian, R kembali

inv(A)

##             [,1]       [,2]       [,3]
## [1,]  0.58333333  0.3333333 -1.0833333
## [2,]  0.33333333  0.3333333 -0.3333333
## [3,] -0.08333333 -0.3333333  0.5833333

Jika kita ingin melihat keluarannya dalam bentuk bilangan rasional, maka kita bisa menggunakan fungsi fractions() dari paket MASS

library(MASS)
fractions(inv(A))

##      [,1]   [,2]   [,3]  
## [1,]   7/12    1/3 -13/12
## [2,]    1/3    1/3   -1/3
## [3,]  -1/12   -1/3   7/12

Kemudian R outputs

 fractions(inv(A))

##      [,1]   [,2]   [,3]  
## [1,]   7/12    1/3 -13/12
## [2,]    1/3    1/3   -1/3
## [3,]  -1/12   -1/3   7/12

Dalam R kita dapat menggunakan fungsi solve() untuk mencari invers dari a matriks, bukan fungsi inv(). Misalnya dengan matriks

Kemudian kita bisa

library(matlib)
A <- matrix(c(1,-2,-1,2,3,2,3,-2,1), nrow = 3, ncol = 3)
solve(A)

##             [,1]       [,2]       [,3]
## [1,]  0.58333333  0.3333333 -1.0833333
## [2,]  0.33333333  0.3333333 -0.3333333
## [3,] -0.08333333 -0.3333333  0.5833333

Kemudian kami memperoleh hasil yang sama dengan fungsi inv() :

 solve(A)

##             [,1]       [,2]       [,3]
## [1,]  0.58333333  0.3333333 -1.0833333
## [2,]  0.33333333  0.3333333 -0.3333333
## [3,] -0.08333333 -0.3333333  0.5833333