固有频率问题

lnly

clc
clear
%syms w m l a b j1 j2 k1 k2 k3 k4 k5 k6
w=1000;m=1.5;l=219;a=109;b=110;;k1=2.5;k2=0.01;k3=249.33;k4=203.24;k5=1.25e+3;k6=1.01e+3;
D=49.1
j1=pi*D^4/32;
j2=pi*D^4/64;
M=[m*b/l,m*a/l,0,0;0,0,m*b/l,m*a/l;-j2/l,j2/l,0,0;0,0,-j2/l,j2/l];
C=[k3*k2,k4*k2,0,0;0,0,k3*k2,k4*k2;-a*k3*k2,b*k4*k2,w*j1/l,-w*j1/l;-w*j1/l,w*j1/l,-a*k3*k2,b*k4*k2];
K=[(k5+k3*k1),(k6+k4*k1),0,0;0,0,(k5+k3*k1),(k6+k4*k1);-a*(k3*k1+k5),b*(k6+k4*k1),0,0;0,0,-a*(k3*k1+k5),b*(k6+k4*k1)];
I=eye(4);
mm=zeros(4);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%转换为标准特征值问题
G=[mm,I;-inv(M)*K,-inv(M)*C];
cri=eig(G)
计算的结果：
以为mm单位
cri =
  1.0e+003 *
  -0.0002 + 2.0001i
  -0.0002 - 2.0001i
  -0.0000 + 0.0001i
  -0.0000 - 0.0001i
  -0.0015 + 0.0475i
  -0.0015 - 0.0475i
  -0.0015 + 0.0475i
  -0.0015 - 0.0475i
以为m单位
cri =
  1.0e+005 *
  -0.9534 + 2.9411i
  -0.9534 - 2.9411i
  -0.9469 + 2.9211i
  -0.9469 - 2.9211i
  -0.0000 + 0.0122i
  -0.0000 - 0.0122i
  -0.0000 + 0.0122i
  -0.0000 - 0.0122i

这是一个求转子固有频率的程序，为什么输入的参数分别以m作单位和以mm作单位后的计算的特征值差别怎那么大啊？而且以为mm单位时，算出的一阶频率特别小，不知道是什么原因？

[ 本帖最后由 mjhzhjg 于 2007-8-1 18:16 编辑 ]

hector1982

m为米，mm为毫米。所以计算的结果不一样啊！因为对应的结构也不一样。