matlab 振型矩阵,请教如何求出有阻尼系统的固有频率及相应振型矩阵

用矩阵形式表示方程：M*X''+C*X'+K*X=F(t)

方程实际上是六自由度系统的振动响应方程，就是说X是一个6×1的矩阵，M，C，K为质量，阻尼，刚度矩阵。而且C，K不是对角阵，各个自由度之间有耦合。下面是一些参量代码：

function dy=odefun(t,y)

%六自由度常微分非齐次方程

m=87.7;

I_x=2.57;I_y=4.94;I_z=3.04;

a=0.2;b=0.15;h=0.3;h_1=h/2;h_2=h/2;

k_x=1.5e5;k_y=8.75e5;k_z=1e5;

k_ox=5e4;k_oy=5e4;k_oz=5e4;

c_x=0.12*(2*sqrt(m*k_x));c_y=0.12*(2*sqrt(m*k_y));c_z=0.12*(2*sqrt(m*k_z));

c_ox=0.12*(2*sqrt(m*k_ox));c_oy=0.12*(2*sqrt(m*k_oy));c_oz=0.12*(2*sqrt(m*k_oz));

%系统参数

c_11=4*c_z+2*c_oz;

c_22=(4*c_x+2*c_ox)*b^2+(4*c_y+2*c_oy)*a^2;

c_33=4*c_x+2*c_ox;

c_44=4*c_x*h_1^2+2*c_oy*h_2^2+(4*c_z+2*c_oz)*b^2;

c_55=4*c_