w=0:0.1:6;
M=11.36;
a=1;
L0=0.25;
r1=0.75;
r2=1.375;
kv=800;
wn=(kv/M)^0.5;
o=w./wn;
d=a*kv*(r2-1)/(wn^2*M*L0^2*r1^3);
n=(kv/(wn^2*M))*(1-2*a*(1-r2+r1)/r1);
y=0.02;
A=U^2*(3*d*U^2+4*(n-2*y^2));
B=(3*d.*U.^2+4*n.*U).^2-64.*U^4*y^2.*(n-y^2)-48*d.*U.^6*y^2;
o.^2==(A+B.^0.5)./4./(U.^2-1);
Ta=(U.^2+1+2.*U./o.^2.*(3/4*d.*U.^3+n.*U-U.*o.^2)).^0.5;
fplot(w,Ta);
代码贴上啊,像这种
w=;
M=;
不然谁给你打代码,debug