通过质量、刚度矩阵计算反求模态频率需要对提取出的矩阵进行什么处理？

目前按照网上的教程已经提取出模型的质量、刚度矩阵并使用命令流将二进制文件.stiff转换成十进制文件，通过python实现了将质量、刚度的十进制文件转换成常规刚度，但是在验证两个矩阵的过程中计算出的模态频率与ANSYS计算得出的模态频率对应不上，有无大佬知道计算时需要对矩阵进行什么处理，以下是我应用的公式与一些处理事项：

1、使用Lanczos方法求解广义特征值问题

2、将质量矩阵进行归一化处理

3、使用如图公式计算特征值，K是刚度矩阵，M是归一化后的质量矩阵， φ是提取刚度矩阵过程中的.mapping文件。