计算R中矩阵的前k个幂

Question

如果我想使用k包中的k = 10计算某个矩阵A的前matrix.power次幂（例如matrixcalc），将我真的需要写A2 <- matrix.power(A, 2)，A3 <- matrix.power(A, 3)，...，A10 <- matrix.power(A, 10)吗？对我来说似乎有点乏味，但是我尝试避免循环，因为我理解R推荐这样做，因为它效率不高。我可以使用apply函数做些魔术吗？基本上，我需要从矩阵A到k的所有幂的某些特定项的总和。

Answer 1

一个选项是lapply

library(matrixcalc)
lst1 <- lapply(2:100, matrix.power, x = A)

，然后用Reduce

求和

Reduce(`+`, lst1)