矩阵的特征值 - 多项式

Question

我计算矩阵的特征值。我想编写一个可以帮助我的程序（一步一步）。我知道特征函数，但我想看一个多项式（像这里）a link 然后计算多项式元素。

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

我在问路 能不能像“charpoly(B,x)”一样在R中求解(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