用matlab求解非线性方程组

来源:百度知道 编辑:UC知道 时间:2024/06/15 14:42:30
y(1)=-32827140*(k3+k2)+487243*(k2)+33208127*(k3)-13571428*(k3)*(w2)-23022959183*(k2)*(w3)
-152015*(k3)*(w1)+2888577*(k2+k3)*(w3)-438976000000*(k2)*(w1)-4899775691*(k2+k3)*(w2)=0
y(2)=-30537926*(k3+k2)+489610*(k2)+31024732*(k3)-13597282*(k3)*(w2)-23110760209*(k2)*(w3)
-152320*(k3)*(w1)+2900173*(k2+k3)*(w3)-439940937937*(k2)*(w1)-4928324169*(k2+k3)*(w2)=0
y(3)=-35288379*(k3+k2)+488148*(k2)+35777442*(k3)-13623200*(k3)*(w2)-23198947280*(k2)*(w3)
-152595*(k3)*(w1)+2910658*(k2+k3)*(w3)-442461458018*(k2)*(w1)-4956063936*(k2+k3)*(w2)=0
y(4)=-31038486*(k3+k2)+489139*(k2)+31524125*(k3)-13710163*(k3)*(w2)-23496071969*(k2)*(w3)
-152883*(k3)*(w1)+2921674*(k2+k3)*(w3)-448975631075*(k2)*(w1)-5006578944*(k2+k3)*(w2)=0
y(5)=61269*(k3)-90250*(k2+k3)-15.535*(k2)-10123802*(k3)-1721938*(k3)*(w2)-13571428*(k3)*(a2)-292*(k3)*(w1)-152018*(k3)*(a1)+10914*(w3)*(k2+k3)+7062*(a3)*(k2+k3)+255478267518*(k2)*(w3)+60590293*(k2)*(a3)-2*(k2)*56402229485*(w1)-2

建立fun.m
function y= myfun(x)
y = (-32827140*(x(2)+x(1))+487243*(x(1))...
+33208127*(x(2))-13571428*(x(2))*(x(4))-23022959183*(x(1))*(x(5))-152015*(x(2))*(x(3))+2888577*(x(1)...
+x(2))*(x(5))-438976000000*(x(1))*(x(3))-4899775691*(x(1)+x(2))*(x(4)))^2 ...
+(-30537926*(x(2)+x(1))+489610*(x(1))+31024732*(x(2))-13597282*(x(2))*(x(4))...
-23110760209*(x(1))*(x(5))-152320*(x(2))*(x(3))+2900173*(x(1)+x(2))*(x(5))...
-439940937937*(x(1))*(x(3))-4928324169*(x(1)+x(2))*(x(4)))^2+(-35288379*(x(2)+x(1))...
+488148*(x(1))+35777442*(x(2))-13623200*(x(2))*(x(4))-23198947280*(x(1))*(x(5))-152595*(x(2))*(x(3))...
+2910658*(x(1)+x(2))*(x(5))-442461458018*(x(1))*(x(3))-4956063936*(x(1)+x(2))*(x(4)))^2 ...
+(-31038486*(x(2)+x(1))+489139*(x(1))+31524125*(x(2))-13710163*(x(2))*(x(4))...
-23496071969*(x(1))*(x(5))-152883*(x(2))*(x(3))+2921674*(x(1)+x(2))*(x(5))-448975631075*(x(1))*(x(3))...
-5006578944*(x(1)+x(2))*(x(4)))^2+(61269*(x(