| 查看: 426 | 回复: 1 | |||
[求助]
这个高次方程组数学模型怎么求解呢？
|
|
我的数学模型是利用一个方程组求出θ1，θ2，θ3的表达式。数学模型如上传的图片。其中H1,H2,H3,l1,l2,l3,r,R均为已知。θ1，θ2，θ3取值范围为0-60°。 我编写的求解函数如下，总是报错，还请各位高手帮助 syms H1 H2 H3 l1 l2 l3 R r theta1 theta2 theta3 eq1=(3/2*R-l1.*sin(theta1)-1/2*l2.*sin(theta2)).^2+(sqrt(3)/2*(l2.*sin(theta2)-R))^2+(H1-l1.*cos(theta1)-H2+l2.*cos(theta2)).^2-3*(r.^2); eq2=(1/2*(l2.*sin(theta2)-l3.*sin(theta3))).^2+(sqrt(3)/2*(2*R-l2.*sin(theta2)-l3.*sin(theta3))).^2+(H2-l2.*cos(theta2)-H3+l3.*cos(theta3)).^2-3*(r.^2); eq3=(3/2*R-l1.*sin(theta1)-1/2*l3.*sin(theta3)).^2+(sqrt(3)/2*(R-l3.*sin(theta3))).^2+(H1-l1.*cos(theta1)-H3+l3.*cos(theta3)).^2-3*(r.^2); theta0=[rand[1,3]*60] f = matlabFunction([eq1;eq2;eq3],'vars',{[theta1 theta2 theta3]}); options=optimset('Display','iter','MaxFunEvals',1000,'TolFun',1e-10); [x,fval,exitflag,output]=fsolve(f,theta0,options)  111.jpg |
回复此楼» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
2楼2016-03-29 10:24:55
20