R中两个矩阵的联合特征值

Question

如何计算定义为方程根的矩阵 A 和 B 的联合特征值 det(lambda * A - B) = 0 在 R?

在 Matlab 中，函数[V,D] = eig(A,B) 接受一个或两个矩阵作为输入参数（例如：How can I find the joint eigenvalues of two matrices in MATLAB?）。 R 函数e=eigen(A) 似乎没有相同的功能。是否有其他方法可以计算 R 中的联合特征值？

Answer 1

您可以使用 geigen 包做到这一点。

library(geigen)

A <- toeplitz(c(2,1))
B <- toeplitz(c(4,3))

jointEigen <- geigen(B, A)

lambda <- jointEigen$values[1]
det(lambda*A - B)
# 0