x = fzero(fun,x0,options)
换初始点,换算法,设置更严格的收敛条件,我看了下，好像也就第27,28,91,215,277,462的函数值的绝对值大于1,其他都在0附近吧,自己手动处理下呗

我任意试了几组，好多在0附近获取的都是奇值点，此时的exitflag一般都是-5，通过explot观察后，好多都需要重新设置初始点，主要是数据量太大了，棘手的很

奇值点是什么意思

options = optimset('fzero');
options.Display = 'off';
rng('default');
for i = 1:N
    f=@(y)(y-y.*x(i).^2+x(i).*(y.^2.*x(i).^2-y.^2+b(i).^2.*m(i).^2).^0.5)./(b(i).*(y-y.*x(i).^2+x(i).*(y.^2.*x(i).^2-y.^2+n(i).^2).^0.5))-a(i);
    [xx(i),fx(i),exitflag,output] = fzero(f,0,options);     
    aa = 0;
    while exitflag~=1
        aa = aa + 1;
        x0 = 100*randn(1);
            if ~isinf(f(x0)) && isreal(f(x0))
                [xx(i),fx(i),exitflag,output] = fzero(f,x0,options);                 
            end
        if aa>1000
            break
        end
    end
end
find(isnan(fx)|abs(fx)>1)

    58
    91
    96
   220
   251
   277
   406
   437
   463
   468
   495
   499
   522
   524

太感谢了，解决一大难题对我来说