如何使用矩陣的不同冪生成圖

Question

我有以下轉換矩陣。

n <-10
A<-matrix(0,n,n)
diag(A[-1,]) <-0.5 
diag(A[,-1]) <-0.5 
A[1,n]<-0.5
A[n,1]<-0.5

如何通過產生不同的A來得到如下圖並計算1-范數？

Answer 1

要計算矩陣的力量，你可以使用expm包或matrixcalc包：

A <- toeplitz(c(1,2,3)) # a square matrix
A
#       [,1] [,2] [,3]
# [1,]    1    2    3
# [2,]    2    1    2
# [3,]    3    2    1

library(expm)
A %^% 2
#       [,1] [,2] [,3]
# [1,]   14   10   10
# [2,]   10    9   10
# [3,]   10   10   14

library(matrixcalc)
matrix.power(A, 2)
#       [,1] [,2] [,3]
# [1,]   14   10   10
# [2,]   10    9   10
# [3,]   10   10   14

對於情節：

powers <- 0:8
Apowers <- lapply(powers, function(k) A %^% k)
norms <- sapply(Apowers, norm, type = "1")

plot(powers, norms)