这里是16个方程解16个未知数的代码
syms x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 x11 x12 x13 x14 x15 x16
eqns=[31541.3*x1-9800-x2*sin(x3)==0,
x2*cos(x3)-180*(2-x1)==0,
19.7133+x2*sin(x3)-98-x4*sin(x5)==0,
x4*cos(x5)-x2*cos(x3)==0,
19.7133+x4*sin(x3)-98-x6*sin(x7)==0,
x4*cos(x3)-x6*cos(x7)==0,
19.7133+x6*sin(x5)-98-x8*sin(x9)==0,
x4*cos(x5)-x8*cos(x9)==0,
19.7133+x8*sin(x7)-98-x10*sin(x11)==0,
x8*cos(x7)-x10*cos(x11)==0,
2838.717+x10*sin(x11)-980-x12*sin(x13)-1200*9.8==0,
x10*cos(x11)-x12*cos(x13)==0,
x12*sin(x13)-691.488-x14*sin(x15)==0,
x12*cos(x13)-x14*cos(x15)==0,
x1+sin(x3)+sin(x5)+sin(x7)+sin(x9)+sin(x11)+x16-18==0,
x16-x14*cos(x15)/(7*9.8)*(-1+sqrt(1+(7*9.8*22.05/(x14*cos(x15)))^2))==0
%180*(2-x1)-x14*cos(x15)==0
];
vars=[x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 x11 x12 x13 x14 x15 x16];
[x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 x11 x12 x13 x14 x15 x16] = vpasolve(eqns,vars);
x3=x3*180/3.14159;
x5=x5*180/3.14159;
x7=x7*180/3.14159;
x9=x9*180/3.14159;
x11=x11*180/3.14159;
x13=x13*180/3.14159;
x15=x15*180/3.14159;
x3=mod(-28349.73579141389,360);
x5=mod(1350.2893382549,360);
x7=mod(-6029.716981624863,360);
x9=mod(27450.3113840014925,360);
x11=mod(-2969.71435382043738,360);
x13=mod(35910.3240477991,360);
x15=mod(-809.706571880199685235,360);
x=double([x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 x11 x12 x13 x14 x15 x16])