Origin非线性拟合与峰值分析操作视频详解

一、视频讲解脉络与前置准备

在打开 Origin 开展拟合前，必须对工作表（Worksheet）中的原始数据进行合规性整理。很多视频教程直接从点击菜单开始，容易忽视原始数据质量对非线性回归算法收敛速度的决定性影响。

1. 视频实操时间轴速查

操作视频时长约 15 分钟，建议读者按照以下时间节点与模块规划对应学习：

• 00:00 - 02:40
  数据导入与异常点处理：散点图绘制、采样点步长核查、使用 Mask 工具屏蔽仪器飞点。
• 02:41 - 06:15
  基线（Baseline）构建与背景剔除：AsLS 非对称最小二乘平滑法的参数调优与底线贴合技巧。
• 06:16 - 10:30
  峰值分析器多峰解卷积：二阶导数识别重叠肩峰、高斯与洛伦兹模型分配、设置参数共享与下界。
• 10:31 - 13:20
  NLFit 迭代收敛与残差诊断：Levenberg-Marquardt 阻尼单步调试、收敛判定与残差图分析。
• 13:21 - 15:00
  参数导出与物理有效性检验：Reduced χ²、Adj. R² 评估与分峰积分面积报告输出。

2. 数据导入与绘图规范

将实验仪器导出的 CSV 或 TXT 文件拖入 Origin 工作表。在绘制图形时，严禁直接使用平滑样条线（Spline）绘图，必须选择 Plot > Symbol > Scatter（散点图）或 Line + Symbol。散点图能真实呈现测量点密度与离散度，便于在拟合过程中直观对比原始点与理论拟合曲线的贴合情况。

二、非线性曲线拟合（NLFit）核心操作分步拆解

非线性曲线拟合的数学本质是通过 Levenberg-Marquardt (L-M) 阻尼最小二乘法，使加权残差平方和（Chi-Square, χ²）达到极小值。按下快捷键 Ctrl + Y 或点击菜单 Analysis > Fitting > Nonlinear Curve Fit 即可唤出 NLFit 主对话框。

1. 拟合函数选择（Function Selection）

在 Function Selection 选项卡中，左侧为函数类别（Category），右侧为具体函数（Function）。选型依据应当来自体系的物理或化学机理，而非单纯比拼数学相关性：

• 动力学与弛豫过程：选择 Exponential 类别下的 ExpDec1（单指数衰减）或 ExpDec2（双指数衰减）。
• 吸附等温线与受体结合：选择 Growth/Sigmoidal 类别下的 Hill 或 LangmuirEXT。
• 连续光谱与色谱单峰：选择 Peak Functions 类别下的 Gauss（高斯分布）或 Lorentz（洛伦兹分布）。

2. 参数初始化与边界硬约束（Bounds）

非线性拟合对参数初值高度敏感。如果初始值偏离真实极值较远，算法可能陷入局部极小值或直接发散导致奇异矩阵错误。在视频实操中，按以下三步进行参数锁定：

1. 初始估算（Initialize Parameters）：切换至 Parameters 标签，点击“Initialize Parameters”图标。Origin 内置的快速傅里叶变换（FFT）或几何估算模块会自动填入振幅（Amplitude）、中心位（Center）和半峰宽（Width）。若自动赋初值偏差明显，可点击工具栏的 Active Adjust 交互按钮，在图形窗口中手动拖动垂直游标对齐峰顶。
2. 施加边界约束（Bounds 标签）：切换至 Bounds 标签。对于光谱与色谱分析，必须强制设置面积参数 A > 0 和宽度参数 w > 0。若不设定下界，L-M 算法在探索多维参量空间时极易将峰宽收敛为负数或使面积反向抵消，造成数学残差很小但物理上完全失真的荒谬拟合。
3. 参数共享与线性约束（Constraints）：如果同一样品包含多个同位素峰或双峰，且仪器展宽保持一致，可在 Share 列勾选宽度共享（如令 w1 = w2），或者在 Constraints 标签输入比例关系式（如 XPS 谱图中 A1 = 0.5 * A2）。约束能大幅降低待求解的未知自由度，使原本发散的迭代迅速收敛。

3. 单步试探与收敛判定

初学者切忌一打开窗口就直接点击“Fit until converged”。规范的视频操作演示流程为：

• 第一步：点击“1 Iteration”（单步迭代）：观察图形窗口中红色的拟合曲线如何移动。如果红线向真实谱峰靠拢且 χ² 下降，说明初值和约束方向正确；若红线瞬间拉平、剧烈震荡或冲出坐标边界，应立即点击“Reset Code/Parameters”恢复初值并检查 Bounds。
• 第二步：点击“Fit until converged”（迭代至收敛）：L-M 算法自适应调整阻尼因子 λ。当连续两次迭代之间的相对卡方改变量满足 Δχ² / χ² < 10⁻⁹ 时，状态栏显示绿色的“Fit converged”（拟合已收敛），此时方可点击“OK”生成报告。

三、峰值分析器（Peak Analyzer）重叠峰解卷积与基线扣除

当谱图由多个相互重叠、伴有倾斜背底的复杂包络峰构成时（如拉曼光谱的 D/G 峰、XRD 衍射重叠峰、色谱共流出组分），单独使用 NLFit 难以精确定位子峰。此时需调用 Analysis > Peaks and Baseline > Peak Analyzer 向导。

1. 四大目标模式（Goal）的选择

2. 基线类型（Baseline Mode）的合理配置

在向导的 Baseline Mode 步骤中，切忌盲目选择高阶多项式（Polynomial），高阶多项式容易在图谱两端发生龙格现象（Runge phenomenon），导致两端翘起并切削边缘峰底。推荐策略如下：

• 非对称最小二乘平滑（Asymmetric Least Squares, AsLS）：对于连续漂移的荧光背底，AsLS 是最为稳健的算法。通过设置平滑参数 λ（通常在 10³ ~ 10⁶）和非对称权重 p（通常在 0.001 ~ 0.01），基线能紧密贴合峰谷下边界而绝不上凸切削真实峰。
• 用户自定义锚点（User Defined）：勾选 Enable Auto Find，设置适当的锚点数量（如 8~12 个点），点击 Modify/Del 按钮在谱图低谷处手动移动锚点，采用 B-Spline（B样条曲线） 插值平滑连接，避免锚点落在谱峰侧翼。

3. 重叠肩峰识别：二阶导数寻峰法

在向导的 Find Peaks 步骤中，如果重叠子峰仅表现为主峰侧翼的微小突起（肩峰），默认的 Local Maximum（局部极大值法）由于一阶导数未过零点而必定漏检。在视频演示中，解决这一难题的标准操作为：

• 将 Peak Finding Method 切换为 Second Derivative (二阶导数法)。二阶导数的负极值点能敏锐捕获谱线曲率发生突变的位置，从而自动锁定隐藏肩峰的中心横坐标。
• 适度调节 Smooth Window Size（平滑窗口宽度）以滤除高频随机白噪声对二阶导数的假峰干扰。
• 若个别弱峰仍未标出，点击 Add 按钮并在图形对应位置双击鼠标左键手动补齐寻峰标记。

4. 分峰函数选型（Gaussian vs Lorentzian vs Voigt）

在 Fit Peaks 页面点击 Fit Control 进入拟合控制器。为不同物理成因的峰选择匹配的理论峰型：

• 高斯函数（Gaussian）：适用于仪器分辨率或热多普勒效应主导展宽的场景，峰顶较圆润、尾翼衰减极快（如固态常温吸收光谱）。
• 洛伦兹函数（Lorentzian）：适用于量子激发态自然寿命主导展宽的场景，峰顶尖锐且具有明显的宽尾翼（如红外光谱、核磁共振谱）。
• 沃伊特函数（Voigt）：高斯与洛伦兹的数学卷积，同时兼顾仪器展宽和本征寿命展宽。在 XRD 晶粒细化与应变分析、XPS 高分辨谱分峰中，Voigt 函数能提供最精准的谱形还原。

四、拟合优度评价与结果物理意义验证

拟合收敛生成报告后，判断分峰结果是否可信，绝不能仅仅观察曲线在视觉上是否贴合。必须从统计学指标和实验物理机制两个维度进行双重严谨复核。

1. 核心统计指标审查

1. 调整后判定系数（Adjusted R²）：在多峰拟合中，单纯增加峰的数量必然使普通 R² 趋向于 1。Adjusted R² 引入了自由度惩罚项，只有新加峰带来的残差下降足以弥补参数增加时，Adj. R² 才会提升。通常要求 Adj. R² > 0.98，高精细谱图要求 > 0.995。
2. 简化卡方（Reduced Chi-Square, χ²ᵥ）：计算公式为 $\chi^2 / DOF$（残差平方和除以自由度）。在比较高斯拟合与 Voigt 拟合时，哪个模型的 χ²ᵥ 更小，说明哪个模型在统计惩罚下更贴近真实物理分布。
3. 残差分布形态（Residual Plot）：在报告中勾选 Residual Plot。合格拟合的残差图必须呈现纯随机白噪声分布，散点在零线上下均匀、致密、无规律地跳动。如果残差图呈现明显的波浪纹、正弦波动或马鞍形 U 弯，绝对表明存在未识别的隐藏子峰，或者峰尾展宽函数严重不匹配。
4. 参数标准误（Standard Error）：打开 Summary Table，审查每个参数的 Error 项。若某个峰面积的标准误与其数值处于相同量级（例如 Area = 120 ± 95），说明两个重叠峰间距过窄，算法产生了严重的强共线性（Collinearity），此分峰结果在物理上不可靠。

2. 物理合理性自洽检验

• 中心峰位漂移检查：解卷积得到的各峰中心坐标 $x_c$ 必须与文献报道的标准能级或振动基团完全吻合，不得为了迁就残差而漂移到毫无理论依据的位置。
• 半高宽（FWHM）物理限值：同一批次测试中，仪器的硬件分辨率决定了最小物理谱宽。若分峰结果中出现远小于仪器分辨率的“超窄尖刺峰”，或者跨越整个谱带的“超宽扁平峰”，均为数学死锁造成的伪峰，必须删除该峰或对 FWHM 施加更严格的 Bounds 约束。

五、视频实操高频故障诊断与关键注意项

在跟随视频进行实际软件操作时，以下三类报错与异常现象出现频率最高，下表总结了成因诊断与标准排查处方：