问题标题:
【mathematica怎么求解方程组Reduce[{s1^2==(x+rSin[[Alpha]]Sin[[Beta]])^2+(rSin[[Alpha]]Cos[[Beta]]+d/2)^2+(h-rCos[[Alpha]])^2,s2^2==(x+rSin[[Alpha]]Sin[[Beta]]-l)^2+(rSin[[Alpha]]Cos[[Beta]]-d/2)^2+(h-】
问题描述:
mathematica怎么求解方程组
Reduce[{s1^2==(x+
rSin[[Alpha]]Sin[[Beta]])^2+(rSin[[Alpha]]
Cos[[Beta]]+d/2)^2+(h-rCos[[Alpha]])^2,
s2^2==(x+rSin[[Alpha]]Sin[[Beta]]-
l)^2+(rSin[[Alpha]]Cos[[Beta]]-d/2)^2+(h-
rCos[[Alpha]])^2,0
廖恒成回答:
解析解太复杂,无法求出啊.只能将未知数的实际值代入方程,求得数值解.
林成鲁回答:
我把x,r,d,h,l都用数值带入,然后只求Alpha和Beta关于s1和s2的关系,可是还是求不出来
难道只能也带入s1和s2,最后求出Alpha和Beta?
廖恒成回答:
你给出一组x,r,d,h,l的值,我试着帮你求解一下。
林成鲁回答:
x=4r=0.5d=1h=2l=8真心感谢大神啊
廖恒成回答:
Clear["Global`*"]
x=4;r=0.5;d=1;h=2;l=8;
plot:=FindRoot[{s1^2==(x+
rSin[[Alpha]]Sin[[Beta]])^2+(rSin[[Alpha]]Cos[
[Beta]]+d/2)^2+(h-rCos[[Alpha]])^2,
s2^2==(x+rSin[[Alpha]]Sin[[Beta]]-
l)^2+(rSin[[Alpha]]Cos[[Beta]]-d/2)^2+(h-
rCos[[Alpha]])^2},{{[Alpha],1,Pi/2},{[Beta],1,
Pi/2}}];
Table[plot,{s1,1,3,0.2},{s2,1,3,0.2}]//MatrixForm
修改最后一条语句s1和s2的范围,可以得到符合要求的一系列值。然后对结果拟合即可得到Alpha和Beta关于s1和s2的关系。数值拟合可以参考其它资料。
查看更多