собственные значения матрицы - полином

Question

Я рассчитываю собственные значения матрицы. Я хочу написать программу, которая поможет мне (шаг за шагом). Я знаю собственную функцию, но я хотел бы видеть многочлен (как здесь) ссылка а затем рассчитать полиномиальные элементы.

Я могу ввести матрицу, но я не знаю, как ввести новую матрицу после вычитания «а» из главной диагонали, не говоря уже о том, чтобы определить функцию, которую даст мне многочлен (определитель матрицы).

Я прошу указания Можете ли вы сделать то же самое как «charpoly (B, x)» и решить (polyB) в R

Answer 1

Используйте charpoly и roots из пакета Pracma:

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