求助：为什么我用MATLAB计算的临界转速与论坛程序计算的差了很多

sogood

clc
clear
l1=0.075;l2=0.036;l3=0.068;l4=0.056;l5=0.02;l6=0.09;l7=0.09;l8=0.0925;l9=0.0925;l10=0.09;l11=0.09;
l12=0.035;l13=0.046;l14=0.071;l15=0.033;l16=0.046;
m1=3.305;m2=5.025;m3=5.22;m4=7.355;m5=6.3;m6=11.195;m7=17.5;m8=17.27;m9=17.065;m10=16.962;
m11=16.86;m12=11.28;m13=6;m14=6.82;m15=5.245;m16=4.88;m17=3.305;
I1=1.01e-5;I2=1.19e-5;I3=1.4e-5;I4=2.48e-5;I5=7.85e-5;I6=5.15e-5;I7=5.15e-5;I8=5.15e-5;I9=5.15e-5;
I10=5.15e-5;I11=5.15e-5;I12=4.6e-5;I13=2.48e-5;I14=1.4e-5;I15=1.19e-5;I16=1.01e-5;
v1=2.7648;v2=13.02;v3=3.9472;v4=7.7487;v5=108;v6=4.36;v7=4.36;
v8=4.12;v9=4.12;v10=4.36;v11=4.36;v12=27;v13=11.48;v14=3.62;v15=15.49;v16=7.34;
u=0.3;
E=2.0e11;
k1=1.0e9;
L=[l1 l2 l3 l4 l5 l6 l7 l8 l9 l10 l11 l12 l13 l14 l15 l16 0];
I=[I1 I2 I3 I4 I5 I6 I7 I8 I9 I10 I11 I12 I13 I14 I15 I16 I16];
M=[m1 m2 m3 m4 m5 m6 m7 m8 m9 m10 m11 m12 m13 m14 m15 m16 m17];
j1=22.1;j2=44.2;j3=44.8;j4=45.4;j5=44.8;j6=44.2;j7=22.1;
J=[0 0 0 0 0 j1 j2 j3 j4 j5 j6 j7 0 0 0 0 0];
K=[0 0 k1 0 0 0 0 0 0 0 0 0 0 0 k1 0 0];
V=[v1 v2 v3 v4 v5 v6 v7 v8 v9 v10 v11 v12 v13 v14 v15 v16 0];
%V=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
q1=3093;q2=6186;q3=6271.75;q4=6350;q5=6271.75;q6=6186;q7=3093;
Q=[0 0 0 0 0 q1 q2 q3 q4 q5 q6 q7 0 0 0 0 0];
k=0;
for w=0:0.01:10000;
for i=1:17;
T(:,:,i)=[1+(L(i)^3)*(1-V(i))*(M(i)*w^2-K(i))/(6*E*I(i)) L(i)+L(i)^2*J(i)*w^2/(2*E*I(i)) L(i)^2/(2*E*I(i)) L(i)^3*(1-V(i))/(6*E*I(i)) L(i)^3*(1-V(i))/(6*E*I(i))*Q(i);
   (L(i)^2)*(M(i)*w^2-K(i))/(2*E*I(i)) 1+L(i)*J(i)*w^2/(E*I(i)) L(i)/(E*I(i)) L(i)^2/(2*E*I(i)) L(i)^2/(2*E*I(i))*Q(i);
   L(i)*(M(i)*w^2-K(i)) J(i)*w^2 1 L(i) L(i)*Q(i);
   M(i)*w^2-K(i) 0 0 1 Q(i);
   0 0 0 0 1];
end
H=T(:,:,1);
for i2=2:17;
       H=T(:,:,i2)*H;
end
    F=H(3,1)*H(4,2)-H(3,2)*H(4,1);
    if F*(-1)^k < 0
       k=k+1;
       wi(k)=w;
       w=wi(k)
        ni(k)=wi(k)*30/pi
end
end

最后算出来的结果是 w =
  3.0927e+003

ni =
  2.9533e+004  （只算了w<10000的）
而论坛程序算出来前三阶为
CriticalSpeeds =
  1.0e+004 *


   0.3221
    0.6184
    0.8023

是我的程序错了吗，请高手指教