R中矩阵和的计算

【问题标题】：Calculation of the sum of matrices in RR中矩阵和的计算
【发布时间】：2019-10-11 14:11:24
【问题描述】：

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

（A-方阵，k-任意自然数）没有“for”？

【问题讨论】：

你的意思是矩阵的迹吗，可以使用下面的链接rpubs.com/aaronsc32/matrix-trace，它使用sum(diag(A))来计算方阵A的迹。

【解决方案1】：

如果我正确地解释了你的符号，也许在基本 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

【讨论】：