啊啊啊啊啊吖

2019-02-28   阅读量: 2120

数据分析师 R语言

如何使用不同的矩阵幂来生成图

我有以下转换矩阵。

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 0 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 回复

推荐课程