矩阵的特征值-多项式

Question

我计算矩阵的特征值。我想写一个对我有帮助的程序（逐步）。我知道本征函数，但是我想看一个多项式（例如这里）a link 然后计算多项式元素。

我可以输入一个矩阵，但是在从主对角线减去“ a”后我不知道如何输入一个新矩阵，更不用说确定多项式（矩阵行列式）赋予我的函数了。

我在问路 您可以和R中的“ charpoly（B，x）”做同样的事情并解决（polyB）

Answer 1

使用pracma软件包中的charpoly和roots：

library(pracma)

m <- matrix(c(1,2,0,2,4,6,0,6,5), 3, byrow = TRUE)
cp <- charpoly(m)

cp
## [1]   1 -10 -11  36

roots(cp)
## [1] 10.71311 -2.22404  1.51093

# same (except for order)
eigen(m)$values
## [1] 10.71311  1.51093 -2.22404