精选优质文档-倾情为你奉上%求解转子系统前三个临界转速和主振型的传递矩阵法clcclear%等截面轴参数l1=0.12;d=0.04;A=pi*d*d/4;%轮盘参数D=0.5;h=0.025;%盘轴材料参数（忽略轴的质量）a=1;u=0.3;rou=7800;E=2.0e11;G=E/(2*(1+u);I=pi*(d4)/64;K1=2.0e7;v1=6*E*I/(a*G*A*l1*l1);mi=rou*pi*D2/4*h;%轮盘的集质量Jp=mi*D2/8; Jd=Jp/2;Ji=Jp-Jd;%参数的数组形式L=l1 l1 l1 l1 l1 l1 l1 l1 l1 l1 l1 l1 l1 0 0;M=0 mi mi mi mi mi mi 0 0 0 0 0 mi mi 0;K=K1 0 0 0 0 0 0 K1 0 0 0 K1 0 0 0;v=v1 v1 v1 v1 v1 v1 v1 v1 v1 v1 v1 v1 v1 0 0;J
