34.1 矩阵乘法

A <- matrix(c(1, 2, 2, 3), nrow = 2)
A
##      [,1] [,2]
## [1,]    1    2
## [2,]    2    3
B <- matrix(c(1, 2, 3, 4, 5, 6), nrow =2)
B
##      [,1] [,2] [,3]
## [1,]    1    3    5
## [2,]    2    4    6

通常的矩阵乘法也叫矩阵内积

A %*% B
##      [,1] [,2] [,3]
## [1,]    5   11   17
## [2,]    8   18   28
A ** 2
##      [,1] [,2]
## [1,]    1    4
## [2,]    4    9
A ^ 2
##      [,1] [,2]
## [1,]    1    4
## [2,]    4    9
A ** A
##      [,1] [,2]
## [1,]    1    4
## [2,]    4   27