我正在嘗試做一個嵌套點
M <- matrix(1:9,ncol=3)
x <- c(1,2,3)
m <- M
for (op in 1:1){
m <- m %*% M
}
z = x %*% m
結果實際上是 x.dot(M.dot(M)):
228, 516, 804
在 python 中,這個回圈可以通過以下方式減少:
x.dot(np.linalg.multi_dot([M]*2))
R有類似的東西嗎?
uj5u.com熱心網友回復:
library(expm)
x %*% (M %^% 2)
[,1] [,2] [,3]
[1,] 228 516 804
正如@akrun 評論的那樣,您還可以使用Reduce:
Reduce('%*%', rep(list(M), 2), init = x)
[,1] [,2] [,3]
[1,] 228 516 804
轉載請註明出處,本文鏈接:https://www.uj5u.com/yidong/462262.html
上一篇:識別值何時在r中重復
下一篇:R-資料框列中序列中數字的平均值
