可以使用MATLAB求解这个问题，
用Matlab 数值工具箱中库函数命令eig 来实现。
[ x ,d] = eig(ik ,im) ; % ik 和im 分别为结构的刚度矩阵和质量矩阵
d = sqrt (d)
for i = 1 :cn % cn 为结构的层数,即质点数
[ dl (i) ,j ] = min (d) ;
Xgd( : ,i) = x ( : ,j)
d(j) = max (d) + 1
end %以此循环对所求频率和振型进行排序
w = dl ; %所求自振频率
x = xgd ; %所求结构主振型
举个例子，可能会更清楚。
某三层钢筋混凝土结构,结构的各层特性参数为:第一层到第
三层质量m 分别为2 762 kg ,
2760 kg ,
2300 kg ,第一层到第三层
刚度k 分别为2. 485 ×104 N/ m , 1. 921 ×104 N/ m , 1. 522 ×
104 N/ m。地震波采用200gal El Centro 波,采样周期为0. 02 s。
