欠定方程求解 Matlab,500论坛币悬赏求解欠定方程的最小二乘解

function y = myfunc(w)

y= [w(1)+w(3)-0.6;

5.85*w(1)+6.12*w(2)+3.42*w(3)+3.88*w(4)+4.34*w(5)+0.99*w(6)+1.89*w(7)+0.5*w(8)-5;

w(3)+w(4)+w(5)+w(6)+w(7)+w(8)-0.4;

0.72*w(1)^2+1.25*w(2)^2+2.5*w(3)^2+0.21*w(4)^2+0.44*w(5)^2+16.1*w(6)^2+9.69*w(7)^2+20.27*w(8)^2+...

0.55*w(1)*w(2)+0.7*w(1)*w(3)+0.21*w(1)*w(4)+0.23*w(1)*w(5)-1.4*w(1)*w(6)-0.04*w(1)*w(7)-...

6.46*w(1)*w(8)+1.65*w(2)*w(3)+0.53*w(2)*w(4)+0.67*w(2)*w(5)-1.56*w(2)*w(6)-0.04*w(2)*w(7)+...

0.4*w(2)*w(8)+0.87*w(3)*w(4)+0.77*w(3)*w(5)+2.36*w(3)*w(6)+1.77*w(3)*w(7)+4.12*w(3)*w(8)+...

0.58*w(4)*w(5)+0.15*w(4)*w(6)+0.31*w(4)*w(7)+0.77*w(4)*w(8)-0.16*w(5)*w(6)+...

0.06*w(5)*w(7)+0.61*w(5)*w(8)+19*w(6)*w(7)+29.88*w(6)*w(8)+22.67*w(7)*w(8)-1.3776^2];

end

命令行: fsolve(@myfunc,zeros(8,1))

结果:

0.2667

0.1872

0.3333

0.1488

0.1607

-0.0855

0.0308

-0.1882