fsolve is not returning the right solutions
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
My code is below. The correct solutions must be (1,0,1) and (-1,0,1) but I keep getting weird answers regardless of my initial guess
function solveeqs()
x0 = [1 1 1];
x = fsolve(@eqns,x0);
x;
eqns(x0)
end
function F = eqns(z)
x1 = z(1);
x2 = z(2);
x3 = z(3);
F(1) = x1 - x1*x2 - x2.^3 + x3*x1.^2 + x3*x2.^2 - x3 - x3*x1 + x3*x1*x2 + x3*x2.^3;
F(2) = x1 - x3*x1 + x3*x2 + 2*x3*x1*x2;
F(3) = 1 + 2*x2.^2 + x3.^2;
end
0 Kommentare
Antworten (2)
Alex Sha
am 31 Dez. 2020
Hi, if you think the solutions should be (1,0,1) or (-1,0,1) , then one equation:
F(3) = 1 + 2*x2.^2 + x3.^2;
should be:
F(3) = 1 + 2*x2.^2 - x3.^2;
0 Kommentare
Matt J
am 31 Dez. 2020
Bearbeitet: Matt J
am 31 Dez. 2020
The solutions you've cited definitely do not solve the equations you've provided,
F=eqns([-1,0,1])
F=eqns([1,0,1])
function F = eqns(z)
x1 = z(1);
x2 = z(2);
x3 = z(3);
F(1) = x1 - x1*x2 - x2.^3 + x3*x1.^2 + x3*x2.^2 - x3 - x3*x1 + x3*x1*x2 + x3*x2.^3;
F(2) = x1 - x3*x1 + x3*x2 + 2*x3*x1*x2;
F(3) = 1 + 2*x2.^2 + x3.^2;
end
0 Kommentare
Siehe auch
Kategorien
Mehr zu Systems of Nonlinear Equations finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!