一、 为什么实验数据需要降噪与频域分析？
在理想状态下，我们期望获得的仪器数据是一条光滑、流畅且规律明显的曲线，能够完美地反映变量之间的物理或化学关系。然而现实中的实验环境往往充满着不可避免的干扰，在采集数据时会受到各种因素的影响。这些干扰在数据图表上表现为叠加在主要趋势线上的杂乱毛刺或剧烈波动，如果直接使用这些带有大量噪声的原始数据进行作图或分析，不仅会导致图表显得粗糙、不专业，难以达到SCI期刊发表的美观标准，更严重的是，噪声会掩盖数据的真实特征，导致峰值定位偏差、积分面积计算错误或导数分析失真，从而得出错误的科学结论。因此，对数据进行降噪处理，还原信号的本来面目，是数据处理中至关重要的一步。
为了有效地去除噪声，我们首先需要理解噪声的本质。对于大多数实验人员来说，我们习惯于在时域观察数据，即观察信号随时间变化的波形。在时域图中，噪声和有效信号往往纠缠在一起，很难通过简单的观察将它们剥离。此时我们需要转换视角，利用数学工具将信号从时域转换到频域。通过频域分析，我们可以清晰地看到信号的内部结构。实验中的有效信号（如化学反应的浓度变化、温度的缓慢升降）变化相对平缓，主要集中在低频区域；而噪声（如电子元件的抖动）通常变化极快，主要分布在高频区域。还有一些特定的干扰，如50Hz的市电干扰，会在频域图上表现为特定的尖峰。一旦我们在频域中识别出了代表噪声的频率范围，就可以有针对性地设计滤波器，保留低频的有效信号颗粒，筛除高频的噪声杂质。
通过Origin等软件进行FFT变换和滤波，实际上就是给实验数据做一次精细的外科手术，切除病灶并保留健康的组织（有效信号），从而确保后续的数据分析结果准确可靠，最终呈现出一张既美观又严谨的SCI科研图表。
二、 导入带有噪声的实验数据与初步预览
打开 Origin 软件，保持一个空白的工作界面，直接在文件夹中选中附件的示例数据文件，将其拖入 Origin 的工作区窗口。软件会自动识别分隔符并将数据填入表格中。此时，我们能看到一个包含三列数据的工作簿，分别是代表时间的"Time (s)"、代表含噪实验数据的"Amplitude (a.u.)"以及作为参考的纯净信号"Clean Signal"：
图1 导入后的工作表
数据导入后需要确认列属性是否设置正确，这是绘图的基础。观察表格的列标题，第一列"Time"应当被标记为 X（表示横轴/自变量），第二列"Amplitude"应当被标记为 Y（表示纵轴/因变量）。为了模拟真实的实验场景，我们暂时忽略第三列"Clean Signal"，仅关注前两列，因为在实际科研中，我们往往只能测得带有噪声的混合信号，而无法直接预知完美的纯净曲线。确认数据无误后，我们进行数据的初步可视化预览，以便直观感受噪声的强度。选中 A 列（Time）和 B 列（Amplitude），点击顶部菜单栏选择绘图------基础2D图------折线图。Origin 会弹出一个图形窗口，展示出一条剧烈波动的曲线。
图2 选中数据列后打开菜单绘制折线图
观察这张初步生成的预览图，可以发现曲线虽然整体呈现出某种周期性的波动趋势，但线条本身非常粗糙，充满了密集的毛刺和锯齿。在时域视角下，噪声的幅度已经严重干扰了对信号细节的观察，如果直接对这条曲线求导或寻找峰值，结果将充满误差。因此，接下来需要进行 FFT 变换与滤波处理：
图3 绘制出的效果
三、 快速傅里叶变换（FFT）的操作流程
首先需要确保当前的活动窗口是刚刚绘制好的含噪信号波形图，然后在 Origin 顶部的菜单栏中，依次点击分析（Analysis）------信号处理（Signal Processing）------ FFT------FFT...。软件会自动锁定当前图层中的曲线作为输入信号，无需手动寻找数据列：
图4 打开FFT分析
在弹出的 FFT 对话框中，选项参数保持默认即可，点击确定开始分析：
图5 点击确定开始分析
此时，Origin 会生成一个新的工作表。在这个新工作表中，按住 Ctrl 选择E 列和 A 列，然后再次选择绘制折线图，即可获得频谱图：
图6 选择新工作表的E和A列并绘制折线图
图7 生成的频谱图
在这个新图中，横坐标不再是时间（Time），而是变成了 Frequency（频率）；纵坐标则是 Amplitude（幅度）。可以发现，原本在时域上杂乱无章的波形，现在变成了一系列竖立的"柱子"或峰。这个转换过程，就像是把一杯混合了各种配料的果汁（时域信号），瞬间分离成了原本的配方表（频域信号），让我们一眼就能看清里面到底加了什么料。
四、 如何通过频谱图识别有效信号与高频噪声
刚刚生成的 FFT 频谱图中，横轴 Frequency (Hz) 代表频率，即信号变化的快慢；纵轴 Amplitude 代表幅度，即该频率成分在原始信号中的强度。为了精准定位信号，我们需要使用 Origin 的 Data Reader（数据读取） 工具。点击工具栏上的"十字准星"图标，鼠标指针会变成一个十字架。将鼠标移动到频谱图中最高的那根柱子顶端，屏幕上会浮动显示该点的坐标。可以发现，最高的峰值精确对应在 X = 5 (Hz) 附近，而第二高的峰值对应在 X = 20 (Hz) 附近。这两个尖锐、突出的峰值就是有效信号，它们代表了实验中真实的物理变化规律：
图8 有效信号
接下来，观察 20Hz 右侧的区域，可以看到一条紧贴着底部、密密麻麻、起伏不定的"草丛"。这些幅度很低、频率分布极广（一直延伸到高频）的杂乱波形，就是要剔除的"高斯白噪声"。在频域分析中，噪声往往表现为这种全频段的低幅"底噪"。识别的关键在于鹤立鸡群的对比：有效信号通常表现为特定频率下的高耸尖峰，而噪声则是弥漫在底部的杂草。
图9 噪声识别
通过这种目视观察，就可以确定滤波的策略了。既然有效信号截止在 20Hz，而 20Hz 之后全是噪声，我们只需要在两者之间划一道界线。比如可以选择 30Hz 或 40Hz 作为截止频率，保留 30Hz 以下的所有内容，并切除 30Hz 以上的所有内容。这就是接下来将要进行的低通滤波（Low-Pass Filtering）的核心逻辑。
五、 使用FFT滤波器（Low-Pass）剔除噪声
首先，需要切回到最初展示原始含噪波形的图形窗口。在 Origin 的顶部菜单栏中，依次点击分析------信号处理------FFT滤波器。这个功能模块专门用于根据频率特性来清洗数据：
图10 确保含噪波形在前并选择FFT滤波器
在弹出的对话框中，最重要的设置项是滤波器类型。点击下拉菜单，可以看到高通、带通等选项。根据之前的分析，有效信号位于低频区，因此需要让低频通过，阻挡高频。确保选项选择了低通。接下来是整个操作的核心------设置截止频率。这个数值就是刚刚设定的门禁标准，回顾频谱图，有效信号最高频率为 20Hz。为了保险起见，不要将截止频率死死地卡在 20Hz，以免误伤有效信号的边缘；同时也不能设得太高，否则会放进太多噪声。根据经验，留出一点余量是明智的，这里输入 40。这意味着，0 到 40Hz 的信号会被完整保留，而大于 40Hz 的剧烈震荡将被强行抹平：
图11 修改参数并确定
观察结果可以发现，灰色的原始曲线上覆盖了一条红色的新曲线。这条红线既完美地贴合了原始数据的波动趋势，又消除了锯齿状毛刺：
图12 过滤后的效果
六、 Savitzky-Golay平滑处理与FFT滤波的效果对比
虽然 FFT 低通滤波在处理周期性波动信号时表现出色，但在科研绘图中还有一个方案可供选择------Savitzky-Golay (S-G) 平滑。对于光谱学研究者而言，S-G 平滑往往是首选方案。为了对比两者的差异，我们需要撤销刚刚的处理，然后重新回到波形图：
图13 选择撤销处理结果
图14 返回该界面
这一次我们直接针对时域进行修整，点击菜单栏的分析------信号处理------平滑：
图15 选择平滑处理
在弹出的平滑设置对话框中，方法（Method） 下拉菜单里默认的通常就是 Savitzky-Golay，如果不是，请手动选择它。S-G 平滑的原理与 FFT 完全不同：FFT 是直接去除高频成分，而 S-G 平滑则像是一个在曲线上滑动的"熨斗"。它利用多项式在滑动窗口内对数据进行最小二乘法拟合。这里最关键的参数是窗口点数（Points of Window），这个数值决定了"熨斗"的大小：数值太小，降噪不明显；数值太大，虽然曲线变滑了，但可能会把原本尖锐的信号峰也给熨平了，导致峰高降低、峰宽变大。针对 1000Hz 的采样率，为了获得与 FFT 类似的平滑度，可以将窗口点数尝试设置为 50 左右，多项式阶数（Polynomial Order）保持默认的 2 阶即可：
图16 平滑参数设置
点击确定后，可以发现结果已经生成：
图17 平滑处理结果
对比之前的结果可以发现，FFT 低通滤波的结果呈现出一种丝滑感，因为它从物理本质上彻底消灭了高频振动，生成的曲线像流体一样流畅，非常适合处理这种由正弦波构成的周期性物理信号。S-G 平滑的结果则更像是在描红，它非常忠实地贴合原始数据的局部走势。虽然它也去除了大部分毛刺，但在某些局部可能不如 FFT 那么圆润。
总结两者的适用情况可以明确，如果你的实验数据是周期性的波形（比如交流电信号、震动波、光波），或者你明确知道噪声是特定频率的高频干扰，FFT 滤波便能够还原出最纯净的波形。如果你的实验数据包含尖锐的峰（例如色谱图、X射线衍射图、红外光谱），需要极力避免峰值高度被削减或峰位置发生偏移，那么 Savitzky-Golay 平滑是更好的选择。因为它在降噪的同时，能最大程度地保留信号的形状特征和峰面积，这对于后续的定量计算至关重要。
七、 原始信号与降噪后信号的对比图绘制
经过前面的处理，现在工作表中已经包含了两种数据：原始的含噪数据，以及经过滤波或平滑处理后的数据。现在我们需要将它们整合在一张图表中，用直观的方式向审稿人展示降噪的效果。
首先回到包含所有数据列的工作表，按住键盘上的 Ctrl 键，依次点击选中三列数据：X轴（Time）、原始数据列（Amplitude） 以及处理后的数据列。选中后，再次选择绘图------基础2D图------折线图：
图18 选择三列数据并再次画图
生成的图表中，两条曲线会重叠在一起。为避免更改一条曲线导致两条曲线一同变化，需要先选择一条曲线，然后将其移出组：
图19 移出组设置
接下来双击灰色的原始含噪数据曲线，在打开的绘图属性------线条选项卡中，将透明度调整为 50% 左右。这一步非常关键，因为我们希望原始噪声作为背景存在，既要能看见，又不能喧宾夺主，抢了处理后数据的风头：
图20 绘图细节
接着选中 FFT 滤波后的红色线条，在绘图细节中将宽度适当加粗，设为 3.0。这样一条光滑、清晰的主线就会浮现在灰色的噪声背景之上，达到了突出显示的效果：
图21 修改线条宽度
最后，直接双击图例，标明"Raw Data"和"FFT Filtered (40Hz Cut-off)"。这样的对比图，一眼就能让人看出结果的意义：
图22 修改图例
至此，我们完成了一次从混乱到清晰的数据进化之旅。从理解频域的意义，到通过 FFT 诊断噪声，再到利用滤波器精准切除，最后通过对比图完美呈现，实现了科研制图的目标。