[x1 x2 x3 y1 y2 y3 x y z] = solve(...
'(x2 - x1) * (x2 - x3) + (y2 - y1) * (y2 - y3) = 0',...
'(x2 - x1)^2 + (y2 - y1)^2 = d1^2',...
'(x2 - x3)^2 + (y2 - y3)^2 = d2^2',...
'k * z * a * (1 - 3 * l0 / ((x3 - a)^2 + (y3 - b)^2 + z^2)^(1/2)) = m * g * x',...
'k * z * b * (1 - 1 * l0 / (x1^2 + y1^2 + z^2)^(1/2)) = m * g * (b - y)',...
'x1 * b * (1 - l0 / (x1^2 + y1^2 + z^2)^(1/2)) + a * (y3 - b) * (1 - 3 * l0 / ((x3 - a)^2 + (y3 - b)^2 + z^2)^(1/2)) = 0',...
'm * g = k * z * (3 - l0 / (x1^2 + y1^2 + z^2)^(1/2) - 2 * l0 / (x2^2 + (y2 - b)^2 + z^2)^(1/2) - 3 * l0 / ((x3 - a)^2 + (y3 - b)^2 + z^2)^(1/2))',...
'(1 - l0/(x1^2 + y1^2 + z^2)^(1/2)) * x1 + (1 - 2 * l0 / (x2^2 + (y2 - b)^2 + z^2)^(1/2)) * x2 + (1 - 3 * l0 / ((x3 - a)^2 + (y3 - b)^2 + z^2)^(1/2)) * (x3 - a) = 0',...
'(1 - l0/(x1^2 + y1^2 + z^2)^(1/2)) * y1 + (1 - 2 * l0 / (x2^2 + (y2 - b)^2 + z^2)^(1/2)) * (y2 - b) + (1 - 3 * l0 / ((x3 - a)^2 + (y3 - b)^2 + z^2)^(1/2)) * (y3 - b) = 0',...
'x1','x2','x3','y1','y2','y3','x','y','z');
为什么运行一小时都不出结果 也没有报错啊 不是有问题吧
哪位大哥CPU强劲 帮忙运行一下呗~
matlab 符号计算本来就橡旁银启和不强
尽量不要求复杂的符号方程，如果参数知道，可以转化成求数值梁宴解的问题
好晕，这是神马
我得去开电脑。
