Question

matlab中求解方程组,用vpasolve解不出来，哪个大神能够帮忙解决一下，不甚感激。

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);

Answer 1

这样解：

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

----------------------

解各为两对共轭复数。

追问

能不能求解最接近的实根？谢谢啦，很急，在线等。。。

追答

确定的复数解已经求出来，要接近的实数解的话就在实部附近找找看吧，找到误差最小的。

Answer 2

>> 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
有结果的。