| 查看: 450 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有12人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有18人回复
-
2027博士申请-临床检验诊断学、病原生物学、微生物方向
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有17人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有10人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有22人回复
-
越来越担心硕士毕业后没有竞争力了
已经有14人回复
-
论文的挂名怎么评判贡献？基金申请中，排名第一的挂名文章
已经有3人回复
-
申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请
已经有3人回复
-
闲聊
已经有6人回复
2楼2016-03-29 10:24:55