如何创建稀疏对角矩阵？

Question

可能是一个非常愚蠢的问题，但我似乎找不到在 R 中创建稀疏对角矩阵的正确方法。我找到了函数：

• diag.spam()
• spdiags()

并将它们与下载的库Matrix和包一起使用spam，但 R 似乎无法识别这些函数。有人知道我需要下载的函数或库吗？

我需要它，因为我想创建大于 256 x 256 的对角矩阵。

Answer 1

包中的Diagonal()函数Matrix。（Matrix是一个“推荐”的包，这意味着它在你安装 R 时自动可用。）

library(Matrix)
m <- Diagonal(500)
image(m)

Diagonal(n)创建一个 nxn 单位矩阵。如果要创建具有指定对角线的对角矩阵x，请使用Diagonal(x=<your vector>)