关于 matlab 求解结构频率特征值的问题

江山美人

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

江山美人

请相关专业的大神指导一下，迫在眉睫！不胜感激！