Question

如何用MATLAB计算方程组

f(1)=(x(1)+x(2)+x(3)+x(4)+x(5)+22.55*x(1)*x(4)+7.63*x(3)*x(4)+148.222*x(1)^3*x(4)+5.86*10^13*x(1)^12*x(4)^7+67.03*x(1)*x(4)^2+212.17*x(1)*x(4)^6+338.08*x(1)*x(5)+4948932.908*x(1)^3*x(5)^2+3522485477*x(1)^4*x(5)^3+3.46*x(4)*x(5)+4.109*x(3)*x(5)+6.03*x(3)*x(5)^2+4.847*x(3)^2*x(5)+135498.74*x(1)^3*x(2)*x(4)^3+7.55*10^14*x(1)^11*x(2)*x(4)^7)-1;
f(2)=0.5653*(1/3*x(2)+135498.74*x(1)^3*x(2)*x(4)^3+7.55*10^14*x(1)^11*x(2)*x(4)^7)-0.2946*(0.5*x(1)+22.55*x(1)*x(4)+444.666*x(1)^3*x(4)+7.03*10^14*x(1)^12*x(4)^7+67.03*x(1)*x(4)^2+212.17*x(1)*x(4)^6+338.08*x(1)*x(5)+14846798.72*x(1)^3*x(5)^2+1.41*10^10*x(1)^4*x(5)^3+406496.24*x(1)^3*x(2)*x(4)^3+8.305*10^15*x(1)^11*x(2)*x(4)^7);
f(3)=0.2946*(0.5*x(3)+7.63*x(3)*x(4)+4.109*x(3)*x(5)+6.03*x(3)*x(5)^2+9.694*x(3)^2*x(5))-0.05*(1/3*x(2)+135498.74*x(1)^3*x(2)*x(4)^3+7.55*10^14*x(1)^11*x(2)*x(4)^7);
f(4)=0.05*(x(4)+22.55*x(1)*x(4)+7.63*x(3)*x(4)+148.222*x(1)^3*x(4)+4.102*10^14*x(1)^12*x(4)^7+134.06*x(1)*x(4)^2+1273.02*x(1)*x(4)^6+3.46*x(4)*x(5)+406496.24*x(1)^3*x(2)*x(4)^3+5.285*10^15*x(1)^11*x(2)*x(4)^7)-0.3456*(0.5*x(3)+7.63*x(3)*x(4)+4.109*x(3)*x(5)+6.03*x(3)*x(5)^2+9.694*x(3)^2*x(5));
f(5)=0.3456*(x(5)+338.08*x(1)*x(5)+9897865.816*x(1)^3*x(5)^2+1.056*10^10*x(1)^4*x(5)^3+3.46*x(4)*x(5)+4.109*x(3)*x(5)+12.06*x(3)*x(5)^2+4.847*x(3)^2*x(5))-0.05*(x(4)+22.55*x(1)*x(4)+7.63*x(3)*x(4)+148.222*x(1)^3*x(4)+4.102*10^14*x(1)^12*x(4)^7+134.06*x(1)*x(4)^2+1273.02*x(1)*x(4)^6+3.46*x(4)*x(5)+406496.24*x(1)^3*x(2)*x(4)^3+5.285*10^15*x(1)^11*x(2)*x(4)^7);

Answer 1

syms x1 x2 x3 x4 x5 f1 f2 f3 f4 f5
f1=(x1+x2+x3+x4+x5+22.55*x1*x4+7.63*x3*x4+148.222*x1^3*x4+5.86*10^13*x1^12*x4^7+67.03*x1*x4^2+212.17*x1*x4^6+338.08*x1*x5+4948932.908*x1^3*x5^2+3522485477*x1^4*x5^3+3.46*x4*x5+4.109*x3*x5+6.03*x3*x5^2+4.847*x3^2*x5+135498.74*x1^3*x2*x4^3+7.55*10^14*x1^11*x2*x4^7)-1;
f2=0.5653*(1/3*x2+135498.74*x1^3*x2*x4^3+7.55*10^14*x1^11*x2*x4^7)-0.2946*(0.5*x1+22.55*x1*x4+444.666*x1^3*x4+7.03*10^14*x1^12*x4^7+67.03*x1*x4^2+212.17*x1*x4^6+338.08*x1*x5+14846798.72*x1^3*x5^2+1.41*10^10*x1^4*x5^3+406496.24*x1^3*x2*x4^3+8.305*10^15*x1^11*x2*x4^7);
f3=0.2946*(0.5*x3+7.63*x3*x4+4.109*x3*x5+6.03*x3*x5^2+9.694*x3^2*x5)-0.05*(1/3*x2+135498.74*x1^3*x2*x4^3+7.55*10^14*x1^11*x2*x4^7);
f4=0.05*(x4+22.55*x1*x4+7.63*x3*x4+148.222*x1^3*x4+4.102*10^14*x1^12*x4^7+134.06*x1*x4^2+1273.02*x1*x4^6+3.46*x4*x5+406496.24*x1^3*x2*x4^3+5.285*10^15*x1^11*x2*x4^7)-0.3456*(0.5*x3+7.63*x3*x4+4.109*x3*x5+6.03*x3*x5^2+9.694*x3^2*x5);
f5=0.3456*(x5+338.08*x1*x5+9897865.816*x1^3*x5^2+1.056*10^10*x1^4*x5^3+3.46*x4*x5+4.109*x3*x5+12.06*x3*x5^2+4.847*x3^2*x5)-0.05*(x4+22.55*x1*x4+7.63*x3*x4+148.222*x1^3*x4+4.102*10^14*x1^12*x4^7+134.06*x1*x4^2+1273.02*x1*x4^6+3.46*x4*x5+406496.24*x1^3*x2*x4^3+5.285*10^15*x1^11*x2*x4^7);
[solx1 solx2 solx3 solx4 solx5]=solve(f1==0,f2==0,f3==0,f4==0,f5==0,x1,x2,x3,x4,x5);