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"))