| 查看: 570 | 回复: 1 | |||
[求助]
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); >> |
回复此楼» 猜你喜欢
UltView：一个比 ParaView 更简单的 Abaqus ODB 查看工具（也支持 VTU/VTK）
已经有0人回复
-
求一个新加坡的机械领域访问交流邀请函
已经有1人回复
-
金属材料论文润色/翻译怎么收费?
已经有172人回复
-
定向工程博士如何申请
已经有2人回复
-
紧急招收2026年秋季入学博士生1名（湘潭大学 固体废弃物低碳利用湖南省工程研究中心）
已经有20人回复
-
【急招】“双一流”高校-新能源材料课题组招收2026年秋季入学博士生1名(湘潭大学)
已经有22人回复
-
祈祷面上基金能中
已经有129人回复
-
自用分享｜桌面深层文件夹分组映射工具 HyperDesk 免费体验版
已经有2人回复
jc1984
金虫 (正式写手)
- 应助: 35 (小学生)
- 金币: 3002.6
- 散金: 290
- 帖子: 958
- 在线: 281.8小时
- 虫号: 562940
- 注册: 2008-05-24
- 性别: GG
- 专业: 传动机械学
2楼2013-05-11 11:02:21
