R语言中矩阵运算

 

1、加减（矩阵加减只能在同维度矩阵之间进行，即行列数相同）

a1 <- c(2, 3, 5, 1, 7, 4, 8, 4, 2)
a2 <- c(3, 1, 3, 1, 4, 8, 4, 3, 4)

m1 <- matrix(a1, nrow = 3, byrow = T)
m2 <- matrix(a2, nrow = 3, byrow = T)
m1
m2

m1 +  m2
m1 -  m2

 

 

2、矩阵相乘， 只有左侧矩阵的列数等于右侧矩阵的行数时才能进行， 相乘矩阵的行数等于左侧矩阵的行数，列数等于右侧矩阵的列数。

a1 <- c(2, 1, 3, 2, 4, 3)
a2 <- c(2, 2, 3, 1, 2, 1)

m1 <- matrix(a1, nrow = 2, byrow = T)    ##2行3列矩阵
m2 <- matrix(a2, nrow = 3, byrow = T)    ##3行2列矩阵
m1
m2
m1 * m2
m1 %*% m2      ## 矩阵之间的相乘应该是 %*%， 相乘的结果为左侧矩阵的行数， 右侧矩阵的列数