mathematica怎么求解方程组
mathematica怎么求解方程组
Reduce[{s1^2 == (x +
r Sin[\[Alpha]] Sin[\[Beta]])^2 + (r Sin[\[Alpha]] \
Cos[\[Beta]] + d/2)^2 + (h - r Cos[\[Alpha]])^2,
s2^2 == (x + r Sin[\[Alpha]] Sin[\[Beta]] -
l)^2 + (r Sin[\[Alpha]] Cos[\[Beta]] - d/2)^2 + (h -
r Cos[\[Alpha]])^2,0 0 用 mathematica 始终在运行啊
解析解太复杂,无法求出啊.只能将未知数的实际值代入方程,求得数值解.我把x,r,d,h,l都用数值带入,然后只求Alpha和Beta关于s1和s2的关系,可是还是求不出来
难道 只能也带入s1和s2,最后求出Alpha和Beta?你给出一组x,r,d,h,l的值,我试着帮你求解一下。x=4r=0.5 d=1 h=2 l=8 真心感谢大神啊Clear["Global`*"]
x = 4 ; r = 0.5 ; d = 1 ; h = 2 ; l = 8;
plot := FindRoot[{s1^2 == (x +
r Sin[\[Alpha]] Sin[\[Beta]])^2 + (r Sin[\[Alpha]] Cos[\
\[Beta]] + d/2)^2 + (h - r Cos[\[Alpha]])^2,
s2^2 == (x + r Sin[\[Alpha]] Sin[\[Beta]] -
l)^2 + (r Sin[\[Alpha]] Cos[\[Beta]] - d/2)^2 + (h -
r Cos[\[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的关系。数值拟合可以参考其它资料。