如何计算矩阵A和B的联合特征值,定义为R中方程det(lambda*A-B(=0的根?
在Matlab中,函数[V,D] = eig(A,B)
接受一个或两个矩阵作为输入参数(例如:如何在Matlab中找到两个矩阵的联合特征值?(。R函数e=eigen(A)
似乎不具有相同的特征。有没有其他方法可以计算R中的联合特征值?
您可以使用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