owenqi 发表于 2022-2-8 12:27 
照着乘法的定义做就可以了。比如
谢谢!我按照您的方法试了一下,显示Error in A[i, ] : incorrect number of dimensions,想请教一下这是哪一步出了问题呢?
A <- matrix(data=1,6,4,2,5,7, nrow = 2, ncol = 3, byrow = FALSE, dimnames = NULL)
A
B <- matrix(data=4,11,13,6,8,3, nrow = 3, ncol = 2, nrow = 3, ncol = 2, byrow = FALSE, dimnames = NULL)
B
m <- nrow(A)
n <- ncol(B)
M <- matrix(NA, nrow = m, ncol =n )
for(i in 1:m)
for(j in 1:n)
M[i,j] <- sum(A[i,]*B[,j])
M
A %*% B