2019-02-28
阅读量:
4755
如何使用不同的矩阵幂来生成图
我有以下转换矩阵。
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-norms来获得如下的情节?






评论(1)


啊啊啊啊啊吖
2019-02-28
要计算矩阵的功效,您可以使用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
0.0000
0
0
推荐帖子
0条评论
0条评论
1条评论