求解一个结构的频率特征值

江山美人

     我在求解一个圆柱壳结构的频率特征值，先假定了位移函数，位移函数采用多项式求和的形式，理论利用能量法求得了刚度和质量矩阵，再用MATLAB编程求解频率特征值，但是遇到一个问题：      当项数取得不是很大的时候，随着项数的增大，频率特征值在慢慢收敛，但是仍未收敛到较精确的结果，而且项数再继续增大，频率特征值开始有点发散，不再收敛。           请问各位大侠，这是什么问题，根据理论，应该是项数越多越收敛才对，难道是因为malab程序计算过程中误差的原因吗,还是需要对得到的矩阵再处理一下吗？请大家指教，谢谢