行列を参照行列n倍で掛ける。

Question

xは、すべての繰り返しの参照行列です。合計2つの行列があります。

matrix A 
1 4 1 4 
4 2 4 2 
2 3 2 3 
3 3 3 3 

matrix B 
1 4 1 4 
4 2 4 2 
2 3 2 3 
3 3 3 3 


matrix x 
4 1 4 3 
2 4 2 2 
3 2 3 5 
3 5 1 1 

ここで、行列AとBは同じです。

C <- B %*% x 
D <- C %*% x 
E <- D %*% x 
F <- E %*% x 

xは、すべての反復の参照行列です。 この乗算は15回実行されるので、ループや関数の書き込みに役立つ必要があります。

Answer 1

あなただけソトスにより示唆されるように、最後の行列は、我々は、再帰関数を使用することができますしたい場合：

A <- matrix(runif(9), 3) 
X <- matrix(runif(9), 3) 

repmult <- function(A,x,reps) 
{ 
    if(reps==0){ 
    return(A) 
    } 
    else 
    { 
    repmult(A%*%x,x,reps-1) 
    } 
} 

repmult(A,X,15) 

---

あなたにも、リスト内のすべての中間結果をしたい場合は、我々は機能を変更することができますthis SO questionの答え（あなたがその名前を変更したいかもしれないが）から：

Mpow <- function(A,x, n) { 
    L <- list(A%*%x) 
    if (n==1) return(L) 
    P <- A 
    for (i in 1:n) L[[i]] <- (P <- P %*% x) 
    return(L) 
} 

Mpow(A,X,3)