请各位大神帮忙看下下边的程序,一直求不到正确的结果。matlab提示no solution found
function x = ceshi
global lamda xs xm
lamda = 10;
xs = 1e3*[0.108885324773886;
0.574079999999998;
1.807059352577014;];
xm = [-0.625083880696937, -0.497901804628957 , -0.630229765520985;
3.106230625461430 , 3.194555748931001, 3.273121443089857;];
x0 = 1e3*[1.805059352577014,1.807059352577014,1.809059352577014];
options = optimoptions('fsolve','Algorithm','trust-region-reflective');
[x,fval] = fsolve(@myfun,x0,options);
end
function F = myfun(x)
global xs lamda xm
F(1) = ( (lamda-x(1))/lamda * xm(1, 1) - xs(1) )^2 + ( (lamda-x(1))/lamda * xm(2, 1) - xs(2) )^2 + ( x(1) - xs(3) )^2 -400;
F(2) = ( (lamda-x(2))/lamda * xm(1, 2) - xs(1) )^2 + ( (lamda-x(2))/lamda * xm(2, 2) -