24小时热门版块排行榜

返回列表

吴秀彬

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 30.5
帖子: 43
在线: 3.8小时
虫号: 2386670
注册: 2013-03-29
专业: 传动机械学

[求助] matlab求助

w = 1;s0 = 100;s = 150;e = 50;
x11 = 0 : (pi/180) : (pi/4);
z1  = (s + s0)*cos(x11)-e*sin(x11);
y1 = (s0 + s)*sin(x11) - e*cos(x11);
x12 = (pi/4) : (pi/180) : (pi/2);
z2  = (s + s0)*cos(x12)-e*sin(x12);
y2 = (s0 + s)*sin(x12) - e*cos(x12);
x2 = (pi/2) : (pi/180) : (13*pi/18);
z3  = (s + s0)*cos(x2)-e*sin(x2);
y3 = (s0 + s)*sin(x2) - e*cos(x2);
x3 = (13*pi/18) : (pi/180) : (7*pi/6);
z4  = (s + s0)*cos(x3)-e*sin(x3);
y4 = (s0 + s)*sin(x3) - e*cos(x3);
x4 = (7*pi/6) : (pi/180) : (2*pi);
z5 = (s + s0)*cos(x4)-e*sin(x4);
y5 = (s0 + s)*sin(x4) - e*cos(x4); plot(z1,y1,'b',z2,y2,'b',z3,y3,'b',z4,y4,'b',z5,y5,'b');
xlabel('x/mm')
ylabel('y/mm')
title('理轮轮曲线')
>>  v=[];
syms z1 z2 z3 z4 z5
s0 = 100;
e = 50;
s11 = 600*x11.^2/pi.^2;
t11  = (s11 + s0)*cos(z1)-e*sin(z1);
y11 = (s0 + s11)*sin(z1) - e*cos(z1);
tx11=diff(t11,z1);
txx11=diff(t11,z1,2);
yx11=diff(y11,z1);
yxx11=diff(y11,z1,2);
for xx11= 0

pi/100)

pi/4);
k11=subs(abs((tx11^2+yx11^2)^1.5/(tx11*yxx11-txx11*yx11)),{z1},{xx11});
v=[v,1/k11];
end
s12=150 - 1200/pi.^2*(pi/2 - x12);
t12=(s12 + s0)*cos(z2)-e*sin(z2);
y12=(s0 + s12)*sin(z2) - e*cos(z2);
tx12=diff(t12,z2);
txx12=diff(t12,z2,2);
yx12=diff(y12,z2);
yxx12=diff(y12,z2,2);
for xx12= (pi/4) : (pi/180) : (pi/2);
k12=subs(abs((tx12^2+yx12^2)^1.5/(tx12*yxx12-txx12*yx12)),{z2},{xx12});
v=[v,1/k12];
end
s2=150;
t2 = (s2 + s0)*cos(z3)-e*sin(z3);
y2=(s0 + s2)*sin(z3) - e*cos(z3);
tx2=diff(t2,z3);
txx2=diff(t2,z3,2);
yx2=diff(y2,z3);
yxx2=diff(y2,z3,2);
for xx2=(pi/2)

pi/100)

13*pi/18);
k2=subs(abs((tx2^2+yx2^2)^1.5/(tx2*yxx2-txx2*yx2)),{z3},{xx2});
v=[v,1/k2];
end
s3 = 75 + 75*cos(4.5*x3 - 2.75*pi);
t3=(s3 + s0)*cos(z4)-e*sin(z4);
y3 = (s0 + s3)*sin(z4) - e*cos(z4);
tx3=diff(t3,z4);
txx3=diff(t3,z4,2);
yx3=diff(y3,z4);
yxx3=diff(y3,z4,2);
for xx3 = (13*pi/18) : (pi/180) : (7*pi/6);
k3=subs(abs((tx3^2+yx3^2)^1.5/(tx3*yxx3-txx3*yx3)),{z4},{xx3});
v=[v,1/k3];
end
s4 = 0;
t4=(s4 + s0)*cos(z5)-e*sin(z5);
y4 = (s0 + s4)*sin(z5) - e*cos(z5);
tx4=diff(t4,z5);
txx4=diff(t4,z5,2);
yx4=diff(y4,z5);
yxx4=diff(y4,z5,2);
for xx4=(42*pi/36)

pi/100)

2*pi);
k4=subs(abs((tx4^2+yx4^2)^1.5/(tx4*yxx4-txx4*yx4)),{z5},{xx4});
v=[v,1/k4];
end
min(v)
??? Error using ==> mupadmex
Error in MuPAD command: not a square matrix [(Dom::Matrix(Dom::ExpressionField()))::_power]

Error in ==> sym.sym>sym.mpower at 207
B = mupadmex('symobj::mpower',A.s,p.s);

>>

回复此楼