fsolve is not returning the right solutions

1 Ansicht (letzte 30 Tage)
Eliyah Hawila
Eliyah Hawila am 31 Dez. 2020
Bearbeitet: Matt J am 31 Dez. 2020
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

Antworten (2)

Alex Sha
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;

Matt J
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 = 1×3
0 0 2
F=eqns([1,0,1])
F = 1×3
0 0 2
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

Kategorien

Mehr zu Systems of Nonlinear Equations finden Sie in Help Center und File Exchange

Produkte


Version

R2020b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by