我正在寻找一种快速算法，用于从 MATLAB 中的稀疏刚度和质量矩阵计算特征值和特征向量。eig(K, M)不适用于稀疏矩阵。