若B可逆比较好计算,如下：
[V,D]=eig(A/B)
D是广义特征值构成的对角矩阵,V是相应特征向量.进一步用QR分解还可以得到正交基,不过你只要特征值,所以就不用了.