| 查看: 601 | 回复: 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); >> |
回复此楼» 猜你喜欢
分享一个可以用于高精度模拟液态碳氢燃料相变的LAMMPS程序
已经有0人回复
-
有幸中面上，特来还愿
已经有114人回复
-
工程热物理与能源利用论文润色/翻译怎么收费?
已经有114人回复
-
英国拉夫堡大学-CSC 合作渠道 |电池方向招收 2027 博士生（免学费+生活费, 30个名额）
已经有18人回复
-
分享力学性能测试标准
已经有0人回复
-
意大利都灵理工大学AFM 课题组招收2027年CSC公派博士研究生
已经有0人回复
jc1984
金虫 (正式写手)
- 应助: 35 (小学生)
- 金币: 3002.6
- 散金: 290
- 帖子: 958
- 在线: 281.8小时
- 虫号: 562940
- 注册: 2008-05-24
- 性别: GG
- 专业: 传动机械学
2楼2013-05-11 11:02:21
