Question

用matlab解含参数方程组

syms x1 x2
>> [x1,x2]=solve('s1^2=(x+r*sin(x1)*sin(x2))^2+(r*sin(x1)*cos(x2)+d/2)^2+(h-r*cos(x1))^2','s2^2=(x+r*sin(x1)*sin(x2)-l)^2+(r*sin(x1)*cos(x2)-d/2)^2+(h-r*cos(x1))^2','x1','x2')
matlab一直在运行啊 算不出来 求大神指导啊
最后 运行出来的答案有1000多页 太恐怖了

Answer 1

建议放弃，solve一般解不了包含形式sin(x)*cos(y)的符号函数。

追问

那用什么来解呢

追答

你运行出来结果了？
我的结果怎么是这样
x1 =

[ empty sym ]

x2 =

[]

追问

出来一个结果 但是复制到word 里有1000页 我的是x1=atan(1/4*(-(-64*d^2*h^4*x^2*l^2+80*h^4*l^4*d^2-128*s1^.........

追答

你试试simplify(A)函数，如果不能化简，就没办法了。
你也可以根据实际问题尝试数值解法，代入其他变量的数值，用fsolve函数。
数值解法：
%%%%%%%%%%%%%%%%%函数
F = myfun(z)
r=rand;
x=rand;
s1=rand;
d=rand;
h=rand;
s2=rand;
l=rand;%%%%各无关变量赋值
F = [s1^2-(x+r*sin(z(1))*sin(z(2)))^2+(r*sin(z(1))*cos(z(2))+d/2)^2+(h-r*cos(z(1)))^2;
s2^2-(x+r*sin(z(1))*sin(z(2))-l)^2+(r*sin(z(1))*cos(z(2))-d/2)^2+(h-r*cos(z(1)))^2];

%%%%%%%%%%%%%调用
x0 = [-5; -5]; % Make a starting guess at the solution
options=optimset('Display','iter'); % Option to display output
[x,fval] = fsolve(@myfun,x0,options) % Call solve
%%%%注：有时候可能找不到解
祝你好运！！

Answer 2

三角函数是有周期的，当含有一个未知数时，解有多个。所以，解三角函数方程往往得不出一个明确的解析解，当含有多个未知数时，更是如此。