matlab程序怎么改,新手求教怎么修改程序

clear all

l1=0.06730;l2=0.03770;l3=0.04300;l4=0.03150;l5=0.05375;l6=0.05375;l7=0.05375;l8=0.05375;l9=0.02450;l10=0.01500;l11=0.04200;l12=0.03650;l13=0.03000;

D1=0.075;D2=0.075;D3=0.087;D4=0.092;D5=0.080;D6=0.080;D7=0.080;D8=0.080;D9=0.080;D10=0.073;D11=0.065;D12=0.065;D13=0.056;

d=0.028;

A1=pi*(D1.^2-d.^2)/4;A2=pi*(D2.^2-d.^2)/4;A3=pi*(D3.^2-d.^2)/4;A4=pi*(D4.^2-d.^2)/4;A5=pi*(D5.^2-d.^2)/4;A6=pi*(D6.^2-d.^2)/4;A7=pi*(D7.^2-d.^2)/4;A8=pi*(D8.^2-d.^2)/4;A9=pi*(D9.^2-d.^2)/4;A10=pi*(D10.^2-d.^2)/4;A11=pi*(D11.^2-d.^2)/4;A12=pi*(D12.^2-d.^2)/4;A13=pi*(D13.^2-d.^2)/4;

h=0.025;

a=2/3;

u=0.25;

rou=7800;

rouz=6010;

E=2.0e11;

G=E/(2*(1+u));

Qa=rouz*pi*((0.080/2).^2)*0.215;

Ma=0;

I1=pi*(D1.^4-d.^4)/64;I2=pi*(D2.^4-d.^4)/64;I3=pi*(D3.^4-d.^4)/64;I4=pi*(D4.^4-d.^4)/64;I5=pi*(D5.^4-d.^4)/64;I6=pi*(D6.^4-d.^4)/64;I7=pi*(D7.^4-d.^4)/64;I8=pi*(D8.^4-d.^4)/64;I9=pi*(D9.^4-d.^4)/64;I10=pi*(D10.^4-d.^4)/64;I11=pi*(D11.^4-d.^4)/64;I12=pi*(D12.^4-d.^4)/64;I13=pi*(D13.^4-d.^4)/64;I14=pi*(D13.^4-d.^4)/64;

Kt1=12*E*I2/0.08;Kt2=12*E*I12/0.08;

Kr1=48*E*I2/(0.08.^3);Kr2=48*E*I12/(0.08.^3);

v1=6*E.*I1./(a*G.*A1.*l1.^2);v2=6*E.*I2./(a*G.*A2.*l2.^2);v3=6*E.*I3./(a*G.*A3.*l3.^2);v4=6*E.*I4./(a*G.*A4.*l4.^2);v5=6*E.*I5./(a*G.*A5.*l5.^2);v6=6*E.*I6./(a*G.*A6.*l6.^2);v7=6*E.*I7./(a*G.*A7.*l7.^2);v8=6*E.*I8./(a*G.*A8.*l8.^2);v9=6*E.*I9./(a*G.*A9.*l9.^2);v10=6*E.*I10./(a*G.*A10.*l10.^2);v11=6*E.*I11./(a*G.*A11.*l11.^2);v12=6*E.*I12./(a*G.*A12.*l12.^2);v13=6*E.*I13./(a*G.*A13.*l13.^2);

m1=rou*pi*l1*(D1.^2-d.^2)/2;m2=rou*pi*l1*(D1.^2-d.^2)/2+rou*pi*l2*(D2.^2-d.^2)/2;m3=rou*pi*l2*(D2.^2-d.^2)/2+rou*pi*l3*(D3.^2-d.^2)/2;

m4=rou*pi*l3*(D3.^2-d.^2)/2+rou*pi*l4*(D4.^2-d.^2)/2;m5=rou*pi*l4*(D4.^2-d.^2)/2+rou*pi*l5*(D5.^2-d.^2)/2+Qa/2;m6=rou*pi*l5*(D5.^2-d.^2)/2+rou*pi*l6*(D6.^2-d.^2)/2+Qa;m7=rou*pi*l6*(D6.^2-d.^2)/2+rou*pi*l7*(D7.^2-d.^2)/2+Qa;m8=rou*pi*l7*(D7.^2-d.^2)/2+rou*pi*l8*(D8.^2-d.^2)/2+Qa;m9=rou*pi*l8*(D8.^2-d.^2)/2+rou*pi*l9*(D9.^2-d.^2)/2+Qa/2;m10=rou*pi*l9*(D9.^2-d.^2)/2+rou*pi*l10*(D10.^2-d.^2)/2;m11=rou*pi*l10*(D10.^2-d.^2)/2+rou*pi*l11*(D11.^2-d.^2)/2;m12=rou*pi*l11*(D11.^2-d.^2)/2+rou*pi*l12*(D12.^2-d.^2)/2;m13=rou*pi*l12*(D12.^2-d.^2)/2+rou*pi*l13*(D13.^2-d.^2)/2;m14=rou*pi*l13*(D13.^2-d.^2)/2;

Jp1=m1*h.^2/8;Jp2=m2*h.^2/8;Jp3=m3*h.^2/8;Jp4=m4*h.^2/8;Jp5=m5*h.^2/8;Jp6=m6*h.^2/8;Jp7=m7*h.^2/8;Jp8=m8*h.^2/8;Jp9=m9*h.^2/8;Jp10=m10*h.^2/8;Jp11=m11*h.^2/8;Jp12=m12*h.^2/8;Jp13=m13*h.^2/8;Jp14=m14*h.^2/8;

Jd1=Jp1/2;Jd2=Jp2/2;Jd3=Jp3/2;Jd4=Jp4/2;Jd5=Jp5/2;Jd6=Jp6/2;Jd7=Jp7/2;Jd8=Jp8/2;Jd9=Jp9/2;Jd10=Jp10/2;Jd11=Jp11/2;Jd12=Jp12/2;Jd13=Jp13/2;Jd14=Jp14/2;

Ji1=Jp1-Jd1;Ji2=Jp2-Jd2;Ji3=Jp3-Jd3;Ji4=Jp4-Jd4;Ji5=Jp5-Jd5;Ji6=Jp6-Jd6;Ji7=Jp7-Jd7;Ji8=Jp8-Jd8;Ji9=Jp9-Jd9;Ji10=Jp10-Jd10;Ji11=Jp11-Jd11;Ji12=Jp12-Jd12;Ji13=Jp13-Jd13;Ji14=Jp14-Jd14;

L=[l1 l2 l3 l4 l5 l6 l7 l8 l9 l10 l11 l12 l13 0];

M=[0 m2 m3 m4 m5 m6 m7 m8 m9 m10 m11 m12 m13 0];

Kt=[0 Kt1 0 0 0 0 0 0 0 0 0 Kt2 0 0];

Kr=[0 Kr1 0 0 0 0 0 0 0 0 0 Kr2 0 0];

v=[v1 v2 v3 v4 v5 v6 v7 v8 v9 v10 v11 v12 v13 0];

J=[0 Ji2 Ji3 Ji4 Ji5 Ji6 Ji7 Ji8 Ji9 Ji10 Ji11 Ji12 Ji13 0];

I=[I1 I2 I3 I4 I5 I6 I7 I8 I9 I10 I11 I12 I13 I14];

Q=[0 0 0 0 Qa/2 Qa Qa Qa Qa/2 0 0 0 0 0];

k=0;

Tit=['第一阶频率的振型图';'第二阶频率的振型图';'第三阶频率的振型图';'第四阶频率的振型图'];

for w=0:0.1:25000;

for i=1:14;

T(:,:,i)=[1+(L(i)^3)*(1-v(i))*(M(i)*w^2-Kr(i))/(6*E*I(i)) L(i)+L(i)^2*(J(i)*w^2+Kt(i))/(2*E*I(i)) L(i)^2/(2*E*I(i)) L(i)^3*(1-v(i))/(6*E*I(i)) L(i)^2*Ma/(2*E*I(i))+L(i)^3*(1-v(i))*Q(i)/(6*E*I(i));

(L(i)^2)*(M(i)*w^2-Kr(i))/(2*E*I(i)) 1+L(i)*(J(i)*w^2+Kt(i))/(E*I(i)) L(i)/(E*I(i)) L(i)^2/(2*E*I(i)) L(i)*Ma/(E*I(i))+L(i)^2*Q(i)/(2*E*I(i));

L(i)*(M(i)*w^2-Kr(i)) J(i)*w^2+Kt(i) 1 L(i) Ma+L(i)*Q(i);

M(i)*w^2-Kr(i) 0 0 1 Q(i);

0 0 0 0 1];

end

H=T(:,:,1);

for i2=2:14;

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);

f=wi(k)/(2*pi)

ni(k)=wi(k)*30/pi;

end

end

for i1=1:1;

w=wi(i1);

for j=1:14;

T(:,:,j)=[1+(L(j)^3)*(1-v(j))*(M(j)*w^2-Kr(j))/(6*E*I(j)) L(j)+L(j)^2*(J(j)*w^2+Kt(j))/(2*E*I(j)) L(j)^2/(2*E*I(j)) L(j)^3*(1-v(j))/(6*E*I(j)) L(j)^2*Ma/(2*E*I(j))+L(j)^3*(1-v(j))*Q(j)/(6*E*I(j));

(L(j)^2)*(M(j)*w^2-Kr(j))/(2*E*I(j)) 1+L(j)*(J(j)*w^2+Kt(j))/(E*I(j)) L(j)/(E*I(j)) L(j)^2/(2*E*I(j)) L(j)*Ma/(E*I(j))+L(j)^2*Q(j)/(2*E*I(j));

L(j)*(M(j)*w^2-Kr(j)) J(j)*w^2+Kt(j) 1 L(j) Ma+L(j)*Q(j);

M(j)*w^2-Kr(j) 0 0 1 Q(j);

0 0 0 0 1];

end

H=T(:,:,1);

for j=2:14;

H=T(:,:,j)*H;

end

b=-H(4,1)/H(4,2);

X(:,1)=([1 b 0 0 1]');

for n=2:15;

X(:,n)=T(:,:,n-1)*X(:,n-1);

end

% YY=X(2,:)./X(1,:);

x=0:40:560;

xx=[0:0.1:560];%你希望得到的曲线的精度，0.1是x方向步长

zz1=spline(x,X(1,:),xx);

%subplot(2,2,i1)

plot(xx,zz1);

%title(Tit(i1,:));

xlabel('轴长/mm');ylabel('振型')

hold on;

end

程序有点长，然后出的图我上传在图片附近，我想问问能不能修改一下把这个图的纵坐标改改，让这个图的纵坐标变成正的。有办法吗？谢谢各位了