帮忙看下程序，关于转子固有频率的

谁说我不会

clc
clear
%等截面轴参数
l=0.45;%轴的总长
d=0.01;
A=pi*d*d/4;
%轮盘参数
D=0.04;
h=0.02;
%盘轴材料参数（忽略轴的质量）
a=1;%截面剪切系数
u=0.3;%泊松比
rou=7800;%密度
E=2.1e11;%弹性模量
G=E/(2*(1+u));%切变模量
I=pi*(d^4)/64;%圆截面的惯性矩
K=2.1e7;%弹性支撑刚度
v=6*E*I/(a*G*A*l*l);%
mi=rou*pi*D^2/4*h;%轮盘的集质量
Jp=mi*D^2/8;    Jd=Jp/2;%转动惯量
Ji=Jp-Jd;
%参数的数组形式
L=[l l l l l l l l l];%轴的分段
M=[0 0 0 mi mi mi mi 0 0];
K=[K K K 0 0 0 0 K K];
v=[v v v v v v v v v];
J=[0 0 0 Ji Ji Ji Ji 0 0];
k=0;
for w=0:0.01:100;
    for i=1:9;
     T(:,:,i)=[1 1/K(i);
         -J(i)*w^2 1;];
    end
    H=T(:,:,1);
    for i2=2:9;
        H=T(:,:,i2)*H;
    end
    F=H(1,1)*H(2,2)-H(1,2)*H(2,1);
    if F*(-1)^k < 0 %求解临界转速
       k=k+1;
       wi(k)=w;
       w=wi(k);
       ni(k)=wi(k)*30/pi;
    end
end
ni;
wi;



为什么出不了结果？

ChaChing

出错提示!?

take a look
5.建议提问的网友分清 编程问题 和 专业问题 http://www.chinavib.com/thread-36746-1-1.html
6.提问的智慧！！！！(发帖前请认真阅读) http://www.chinavib.com/thread-21991-1-1.html
from http://home.vibunion.com/blog-63979-18250.html

hustxyong

你问的问题太泛了，具体点别人才好给你建议