关于频响函数的问题

The_Old__Man

我以Euler-Bernoulli梁单元为基础 将一个量分成五个单元，然后根据有限元法集成整体质量[K]和刚度矩阵[M]
用matlab [gEigVector,gEigValue] = eigs(K, M, 15, 'sm')  求特征值及特征向量，特征值前三阶为接近零的小数

阶数        特征值                      频率(Hz)        圆频率^2(rad/s)  
1   -1.3871530e-005    0.0000000e+000   -1.3871530e-005
2   -2.5914562e-005    0.0000000e+000   -2.5914562e-005
3   -8.9737013e-005    0.0000000e+000   -8.9737013e-005
4    4.1952843e+008    3.2598741e+003    4.1952843e+008


然后将[K] [M] [C]=0.00001[K]  坐标模态化 为[k]， [m]，[c]

求其频响函数y=1/([k]-w^2*[m]+j*w*[c])

再将其坐标物理化为[H] 选择H（1，1）

得到其频率幅值图，为什么求解出来的是负值？ 是因为圆频率前三阶吗？

那位大师指点一下！谢谢啦