用MATLAB的fsolve做运动分析时，初始值不同，解出来的答案也不同，怎么选择呢？

多解问题，不奇怪。例如，直线与圆的交点。

function hah
x0=[-1 -1];
%x0=[1 1];
x=fsolve(@myfun,x0)
ezplot('x^2+y^2-1'),hold
ezplot('x-y')
axis equal
function f=myfun(X)
x=X(1);y=X(2);
f(1)=x-y;
f(2)=x^2+y^2-1;

结果：

Optimization terminated: first-order optimality is less than options.TolFun.

x =

-0.7071 -0.7071

function hah
%x0=[-1 -1];
x0=[1 1];
x=fsolve(@myfun,x0)
function f=myfun(X)
x=X(1);y=X(2);
f(1)=x-y;
f(2)=x^2+y^2-1;

结果：
Optimization terminated: first-order optimality is less than options.TolFun.

x =

0.7071 0.7071

都是答案。