matlab fsolve 与fzero,fsolve / fzero：找不到解决方案，看似常规

您的系统设置方式，绘制它并观察其行为实际上很方便 . 我向你的函数进行了矢量化并绘制了 f(x) = MLNEfun(x)-x ，其中 MLNE(x) 的输出是 newA . 实际上，您对系统的固定点感兴趣 .

我观察到的是：

在A~3800处有一个奇点和一个根交叉 . 我们可以使用 fzero ，因为它是一个括号中的根求解器，并且在解决方案 fzero(@(x)MLNEfun(x)-x, [3824,3825]) 上给它非常紧的界限，它产生了 3.8243e+03 . 从您的开始猜测开始，这是几个数量级 . 您的系统在~3e5附近没有解决方案 .

Update 在我的匆忙中，我没有放大图，它显示了另一个(表现良好的)根在1.3294e 04.由你决定哪个是物理上有意义的根 . 我在下面说的一切仍然适用 . 在您感兴趣的解决方案附近开始猜测 .

In response to comment 由于你想为 K 的不同值执行此操作，那么你最好的办法就是坚持 fzero ，只要你要求一个变量而不是 fsolve . 这背后的原因是 fsolve 使用牛顿方法的变体，这些方法没有括号，并且很难在像这样的奇点找到解决方案 . 另一方面， fzero 使用布伦特的方法，该方法保证在括号内找到根(如果存在) . 它在奇点附近表现得更好 .

在执行布伦特方法之前，MATLAB的 fzero 的实现也会搜索包围间隔 . 因此，如果您提供足够接近根的单个起始猜测，它应该为您找到它 . 以下 fzero 的样本输出：

fzero(@(x)MLNEfun(x)-x, 3000, optimset('display', 'iter'))

Search for an interval around 3000 containing a sign change:

Func-count a f(a) b f(b) Procedure

1 3000 -616789 3000 -616789 initial interval

3 2915.15 -433170 3084.85 -905801 search

5 2880 -377057 3120 -1.07362e+06 search

7 2830.29 -311972 3169.71 -1.38274e+06 search

9 2760 -241524 3240 -2.03722e+06 search

11 2660.59 -171701 3339.41 -3.80346e+06 search

13 2520 -109658 3480 -1.16164e+07 search

15 2321.18 -61340.4 3678.82 -1.7387e+08 search

17 2040 -29142.6 3960 2.52373e+08 search

Search for a zero in the interval [2040, 3960]:

Func-count x f(x) Procedure

17 2040 -29142.6 initial

18 2040.22 -29158.9 interpolation

19 3000.11 -617085 bisection

20 3480.06 -1.16224e+07 bisection

21 3960 2.52373e+08 bisection

22 3720.03 -4.83826e+08 interpolation

....

87 3824.32 -5.46204e+48 bisection

88 3824.32 1.03576e+50 bisection

89 3824.32 1.03576e+50 interpolation

Current point x may be near a singular point. The interval [2040, 3960]

reduced to the requested tolerance and the function changes sign in the interval,

but f(x) increased in magnitude as the interval reduced.

ans =

3.8243e+03