假设A是一些方阵.如何在R中轻松取幂这个矩阵？
我已经尝试了两种方法:试用1使用for-loop hack和试验2更优雅但它仍然与A k简单相去甚远.
试验1
set.seed(10)
t(matrix(rnorm(16),ncol=4,nrow=4)) -> a
for(i in 1:2){a <- a %*% a}
试用2
a <- t(matrix(c(0,1,0,0,0,0,1,0,0,0,0,1,0,0,0,0),nrow=4))
i <- diag(4)
(function(n) {if (n<=1) a else (i+a) %*% Recall(n-1)})(10)