这样解:
syms s theta;
eqv1=248*1500*65*(s+67.5*theta)^2+248*1500*90*(s+230*theta)^2+1500*(0.9388*s^2+2.6132*s)-16*600000;
eqv2=47.5*248*1500*65*(s+67.5*theta)^2-115*248*1500*90*(s+230*theta)^2+115*1500*(0.9388*s^2+2.6132*s)-16*(500000*100-100000*35);
S=vpasolve([eqv1==0,eqv2==0],[s,theta])
S.s
S.theta
----------------------
解各为两对共轭复数。
>> syms s theta
eqv1=248*1500*65*(s+67.5*theta)^2+248*1500*90*(s+230*theta)^2+1500*(0.9388*s^2+2.6132*s)-16*600000;
eqv2=47.5*248*1500*65*(s+67.5*theta)^2-115*248*1500*90*(s+230*theta)^2+115*1500*(0.9388*s^2+2.6132*s)-16*(500000*100-100000*35);
[eqv1,eqv2]=vpasolve(eqv1,eqv2,[s,theta])
eqv1 =
0.97036579567261583982564910600157 + 0.095475542191596131154371622941616i
0.97036579567261583982564910600157 - 0.095475542191596131154371622941616i
- 0.97081699073775224189564920423288 + 0.095605986164757021108557270903887i
- 0.97081699073775224189564920423288 - 0.095605986164757021108557270903887i
eqv2 =
- 0.0042190562249738357301920133576129 - 0.0014148541929321459107515710370507i
- 0.0042190562249738357301920133576129 + 0.0014148541929321459107515710370507i
0.0042209301541532865472000341497865 - 0.0014163133261377556791725464228764i
0.0042209301541532865472000341497865 + 0.0014163133261377556791725464228764i
有结果的。