如何让最小二乘（Levenberg-Marquardt）法在工程代码中正常运行

问题表现 #

• 目标为最小化残差向量F使其趋近于0，迭代过程出现结果无变化、输出荒谬值的现象，初步怀疑存在变量类型不匹配、向量输入维度错误类问题，暂未定位根因
• 优化目标为调整各级速度变量c1、w2，使各级动叶通流流量G_perfil_M、喷嘴通流流量G_perfil_F尽可能一致，二者允许误差由代码末尾的Levenberg-Marquardt（LM）算法参数定义，目前调试受阻

已定位的核心隐患 #

• 全局可变数组状态污染：代码中c1、w2、p1、p2、i0、i2等中间计算变量为全局作用域的numpy数组，LM算法每轮迭代调用目标函数时，上一轮的计算结果会残留在数组中，没有做逐次迭代的状态重置，直接导致梯度计算偏离真实值，引发迭代不收敛、结果跳变
• 变量定义冲突：p2_fixo初始定义为长度50的float数组，后续直接被覆盖为标量2.5，虽然当前运算逻辑能正常执行，但存在数组维度广播的潜在风险，后续扩展代码时极易触发维度不匹配报错
• 变量作用域风险：G_estagio在循环内的if x_func == stage分支中赋值，一旦循环范围调整，会直接触发未定义变量报错
• 初始值非物理：所有c1、w2的初始猜测值直接设为600m/s，远超工业汽轮机级蒸汽速度的常规范围（200~400m/s），极易触发IAPWS水属性计算的非物理状态区间，返回无效值或计算报错
• 无异常值拦截机制：目标函数计算过程中没有校验IAPWS模块的输入输出合理性，一旦计算得到非物理的压力、比容、流量值，会直接传入残差计算环节，导致LM算法向荒谬的解空间迭代

调试修正方向 #

• 将目标函数内部用到的所有中间计算变量改为函数内的局部变量，每次迭代进入目标函数时重新初始化计算状态，彻底消除全局数组的状态残留问题
• 统一常量定义逻辑：直接将p2_fixo = 2.5作为标量常量提前声明，删除此前长度50的无意义数组初始化逻辑
• 提前估算合理的初始猜测值：根据通流面积、预估初始流量反推c1、w2的合理取值范围，避免初始值落在非物理区间
• 增加输入输出校验：每次调用IAPWS水属性接口前，判断给定的焓、熵参数是否落在合法的水/蒸汽状态区间内，对非物理计算结果直接返回大残差，避免异常值污染迭代过程
• 增加残差维度校验：在目标函数返回前，确认残差向量长度与优化变量长度满足LM算法的要求（残差维度≥优化变量维度），避免维度不匹配导致的迭代失效

内容的提问来源于stack exchange，提问作者Marcus R Marincek