Question

我正在尝试制作一个矩阵形式，例如

从j=0到30 [ (A^j%*%V%*%t(A)^j) ]的求和。 A和V 3为3矩阵。因此，结果应该是3乘3矩阵形式。

例如

A=matrix(1:9, ncol = 3)
V=matrix(1:9, ncol = 3)

我试过

for (i in 0:30) {
    AVA=sum(A^i%*%V%*%t(A)^i)
}

但它不起作用。

你知道怎么做如上所述的矩阵求和吗？

Answer 1

尝试

A=matrix(1:9, ncol = 3)
V=matrix(1:9, ncol = 3)
AVA=matrix(rep(0,9), ncol = 3)
Atopoweri = diag(3) # A to zeroth power
for (i in 0:n) {
    AVA = AVA + Atopoweri%*%V%*%t(Atopoweri)
    Atopoweri = Atopoweri %*% A
}

我假设^你的意思是矩阵的力量，但如果你的意思是元素的力量，那么@C Doan的回答就是你想要的东西。

Answer 2

AVA=matrix(0,ncol=3,nrow=3)
for(i in 0:30) AVA=AVA + A^i %*% V %*% t(A)^i

这将为您返回求和矩阵，而不是将所有元素求和为整数

我试图对三个矩阵求和

2 个答案: