Origin中HLN公式拟合求助：参数未识别且拟合不收敛

问题概述 #

在进行HLN拟合时遇到两类问题：

• 公式评估时提示 x, y, A, L not used
• 强行拟合显示 Fit did not converge, and reason unknown

需拟合参数为前置因子A（预期范围-1到1）和相位相干长度L，仅拥有7个磁导率数据点。

脚本公式问题修正 #

原Labtalk脚本存在冗余代码，导致变量关联逻辑异常：

psi_deriv=real_polygamma(x,k)  // 定义后未使用，属于冗余代码
k=0  // 定义后未使用，冗余
y=(A*1.23314709*10^-5)*(real_polygamma(0.5+1.64557041*10^(-16)/(x*L^2),1)-ln(1.64557041*10^(-16)/(x*L^2)))

修正后的脚本 #

删除冗余代码，优化格式提升兼容性：

y=(A*1.23314709e-5)*(real_polygamma(0.5+1.64557041e-16/(x*L^2),1)-ln(1.64557041e-16/(x*L^2)))

• 用科学计数法简写10^-5为e-5，符合Labtalk的常规写法
• 移除未使用的psi_deriv和k定义，消除“变量未使用”的提示

拟合收敛问题解决建议 #

1. 参数初始值设置
   • 给A设置初始值（如0），给L设置物理合理的初始值（比如1e-8 m，根据样品特性调整），避免拟合器从随机值开始迭代
2. 参数范围约束
   • 严格锁定A的范围为-1到1，在拟合设置中开启参数范围限制并填入对应数值
   • 给L设置合理的物理范围（如1e-10到1e-6 m），避免拟合器进入无意义的参数空间
3. 小样本数据优化
   • 检查7个数据点是否覆盖HLN公式的敏感区间（如低场区域），确保数据能反映公式的变化趋势
   • 尝试加权拟合，给误差较小的数据点更高权重，提升拟合稳定性
4. 拟合算法切换
   • 更换Labtalk拟合器的算法（比如从Levenberg-Marquardt改为Nelder-Mead），不同算法对小样本和复杂公式的适应性存在差异

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