R CAex 阿尔伯斯的示例矩阵与“困难”特征分解

说明

eigen() 似乎很困难的稀疏矩阵示例，其未缩放版本已发布到网络上，并附有 E-mail 到 R-help ( https://stat.ethz.ch/mailman/listinfo/r-help )，作者：开放大学 Casper J Albers ， 英国。

用法

data(CAex)

格式

这是一个 72 \times 72 对称矩阵，在五个波段中有 216 个非零条目，存储为类 dgCMatrix 的稀疏矩阵。

细节

历史记录 (2006-03-30)：在早期版本中R,eigen(CAex)陷入无限循环而eigen(CAex, EISPACK=TRUE)本来还好。

例子

data(CAex, package = "Matrix")
str(CAex) # of class "dgCMatrix"

image(CAex)# -> it's a simple band matrix with 5 bands
## and the eigen values are basically 1 (42 times) and 0 (30 x):
zapsmall(ev <- eigen(CAex, only.values=TRUE)$values)
## i.e., the matrix is symmetric, hence
sCA <- as(CAex, "symmetricMatrix")
## and
stopifnot(class(sCA) == "dsCMatrix",
          as(sCA, "matrix") == as(CAex, "matrix"))