Matlab中fminsearch停止准则问题：拟合微分方程组参数无法终止

我来帮你拆解下这个问题——你说fminsearch已经找到了可接受的最小值，但就是停不下来，结合你给出的optimset设置，大概率是这几个原因在搞鬼：

1. 目标函数的平坦区域卡壳了 #

fminsearch用的是Nelder-Mead单纯形算法，这货在遇到目标函数变化极小的“平坦区”时容易犯轴：虽然函数值已经低于你设的TolFun=1e-4，但单纯形的顶点移动幅度还没达到TolX=1e-4的要求，结果就一直在这个区域里反复试探，迟迟不收敛。

• 解决办法：
  • 既然你不追求极致拟合，不如把TolX调宽松点，比如改成1e-3，让参数的收敛阈值放宽，这样算法更容易判定“已经到位了”；
  • 给参数加个隐形边界：fminsearch本身不支持直接设参数范围，但你可以在目标函数里加惩罚项——比如某个参数超出合理范围时，直接返回一个超大的函数值，让算法自动避开这些无效区域，减少无意义的迭代。

2. 微分方程求解的数值噪声在捣乱 #

拟合微分方程时，每次计算目标函数都要解DEs，数值求解的误差可能会让目标函数值出现微小的随机波动。这就导致fminsearch误以为“函数值还在变化，还没收敛”，一直在迭代。

• 解决办法：
  • 优化DEs的求解精度：比如用ode45时，通过odeset设置更小的RelTol和AbsTol，减少数值求解带来的噪声；
  • 给目标函数做“平滑处理”：比如对DEs的求解结果重复计算3-5次取平均值，稳定目标函数的输出，让算法更容易判断收敛状态。

3. MaxFunEvals的设置可能不匹配你的场景 #

你设了MaxFunEvals=1000，但如果你的参数维度较高（比如3个及以上），1000次函数评估可能还没让单纯形满足TolX的要求。不过反过来，既然你已经得到了可接受的结果，完全可以主动终止迭代。

• 解决办法：
  • 在目标函数里加个“提前终止开关”：当函数值低于你的预期阈值（比如1e-4）时，直接返回一个极小值，让fminsearch认为已经找到最优解；
  • 直接调小MaxFunEvals，比如设为500，强制在迭代次数达标后停止——反正你本来就不追求完美拟合，够用就行。

4. 参数迭代后期在最优值附近震荡 #

如果初始参数离最优区域较远，fminsearch前期找到可接受最小值后，后期可能会在这个值附近来回微调，始终达不到TolX的要求，导致停不下来。

• 解决办法：
  • 把当前找到的可接受参数作为新的初始值，重新运行fminsearch，算法从更接近最优解的位置开始，大概率能快速收敛；
  • 手动监控迭代过程：当你看到函数值连续几次变化都小于TolFun时，直接中断程序就行，没必要死等算法自动停止。

内容的提问来源于stack exchange，提问作者Esme_