热线电话:13121318867

登录
2019-02-28 阅读量: 4529
如何使用不同的矩阵幂来生成图

我有以下转换矩阵。

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来获得如下的情节?

在此输入图像描述

0.0000
3
关注作者
收藏
评论(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 回复