2014-03-25 21 views
0

我試圖做一個矩陣的形式,如從j=0我想提出三點矩陣的總和

求和來30[ (A^j%*%V%*%t(A)^j) ]。和AV33矩陣。因此,結果應該是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) 
} 

,但它無法正常工作。

你有什麼想法如何做矩陣求和如上?

+0

做什麼? –

+0

我試圖總結(A^0%*%V%*%A^0 + A^1%*%V%*%A^1 + A^2%*%V%*%A^2 +。 ... + A^30%*%V%*%A^30)。謝謝 – user3453816

+0

用'^ i'表示矩陣的功率,還是把矩陣的每個元素都加到第i個功率上?看到其他答案之一,我意識到我和我的前任一樣。 – TooTone

回答

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杜安的回答線的東西是你想要的。

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

這將返回求和矩陣對你來說,不是所有元素的總和爲整數要計算所有款項的總和

+0

這就是我想要的。非常感謝。 – user3453816

+0

正如託尼所說,你應該考慮矩陣的力量(以他的方式)或矩陣元素的力量(以我的方式)乾杯! –