Fatal编程技术网
  • r/
  • R中矩阵和的计算

R中矩阵和的计算

r

R中矩阵和的计算,r,rstudio,R,Rstudio,我在考虑这个问题。如何在R中计数 A-平方矩阵，k-任何不带for的自然数？如果我正确解释了你的符号，可能在基R中类似于这样的东西 A <- matrix(c(1,2,3,4), nrow = 2) #example matrix k <- 10 B <- Reduce(`%*%`, (rep(list(A), k)), accumulate = TRUE) #list of A^(1:k) BB <- lapply

我在考虑这个问题。如何在R中计数

A-平方矩阵，k-任何不带for的自然数？

如果我正确解释了你的符号，可能在基R中类似于这样的东西

A <- matrix(c(1,2,3,4), nrow = 2)                         #example matrix

k <- 10

B <- Reduce(`%*%`, (rep(list(A), k)), accumulate = TRUE)  #list of A^(1:k)

BB <- lapply(1:k, function(k) B[[k]]/k)                   #list of A^(1:k)/k

Reduce(`+`, BB)                                           #sum of series BB

         [,1]    [,2]
[1,] 603684.8 1319741
[2,] 879827.1 1923425
你的意思是矩阵的轨迹，这可以使用下面的链接，它使用sumdiagA来计算方阵A的轨迹。