每次处理都保存输入与参数
处理步骤会改变线形和面积，结果要能回到原始数据重算。
读图时保留这项条件：自动相位和基线结果仍需检查，不能只凭谱图整齐判断正确。
转动相位，看看面积为何会失真
切换正确相位、90° 相位错误和斜基线，再分别拖动两个参数。
4.6.1 从复数 FID 到可检查的谱图
本节接续4.5的采样概念。读者已能区分谱宽、采集时间和谱线宽度；现在要判断一个处理操作改变了什么，并保存能够重做的步骤。常规一维氢谱的处理顺序通常为读取原始数据与参数、检查数字滤波延迟、加窗、零填充、傅里叶变换、相位校正、基线校正和化学位移定标。软件导入时可能自动完成其中数步，因此打开一个谱图后应先查看处理历史。Oxford 的培训资料明确指出，自动导入包含加窗、变换及相位操作，随后仍可能需要手动调整（Mycroft，2024）。
展开说明与推导
复数 FID（complex FID）保留接收器的两路正交信号，可写为实部加虚部。虚数单位满足 \(i^2=-1\)，它用于记录旋转相位，不表示另一类氢。以下采用一个教学模型：单条谱线的信号从 \(t=0\) 开始，以正号复指数旋转，并按单指数衰减。
其中 \(a\) 为初始振幅，\(R\) 是有效横向衰减率，单位为 \(\mathrm{s^{-1}}\)；\(f_0\) 是相对接收载频的频率偏置，单位为 Hz；\(\phi\) 是初相位，单位为弧度。实际 FID 是多条谱线与噪声的叠加。保存时不能只留实部，因为后面的相位校正需要在复平面旋转。原始 FID、采集参数和已有处理参数应分别保留，处理结果另存为新版本。
处理参数描述从原始数组到结果数组的变换；采集参数描述测量如何获得原始数组。修改显示纵轴倍率只改变画法；修改窗函数、相位或基线会改变供分析的数组。记录时应将二者分开。
4.6.2 数字滤波与群延迟校正
接收信号经数字滤波和降采样后，输出的时间起点可能滞后于理想 FID。这个群延迟（group delay）可包含非整数个采样间隔。对于纯延迟模型，在本节负号傅里叶约定下，延迟 \(\tau\) 会引入频率因子 \(\exp(-i2\pi f\tau)\)，即随频率变化的相位。滤波器在有限记录起点的行为还可能产生瞬态，不能仅靠调整一个零阶相位处理。
展开说明与推导
校正应先读厂商参数，判断导入程序是否已处理延迟，再执行与该数据格式相符的算法。例如 nmrglue 官方文档中的 Bruker 数据示例提供数字滤波校正，接口同时说明其截断选择可能移除有用数据（nmrglue 项目，在线文档）。这项功能有明确的数据格式和参数条件，不能对所有 FID 统一删去若干点。已校正的数据若再执行一次，时间起点会再次移动，之后补相位也未必恢复原始信息。
实际检查至少包含原始点数、校正后点数、使用的延迟参数以及首段 FID 的变化。厂商的“数据点数”可能把实部、虚部分别计数，必须先换成复数采样点数，再计算时间间隔和补零长度。若只读到一份处理后的实数谱，通常无法完整重做这一步。应回到原始文件核对，避免在频域不断增加一阶相位的绝对值。
4.6.3 窗函数怎样改变线形
窗函数（window function）是乘在 FID 上的权重。后段 FID 的信号已经衰减，噪声仍然存在；减小后段权重可改善某些信号的可见性，也会改变线形。常见的指数加窗（exponential apodization）用正的线宽增量 \(L\) 表示权重。对式（4.6-1）的理想无限长单指数信号，有
展开说明与推导
这里 \(L\) 和吸收线的半高全宽 \(\Gamma\) 均以 Hz 表示。因为两指数相乘使衰减率变为 \(R+\pi L\)，线宽便增加 \(L\)；此推导要求单指数衰减、充分长的记录和正确相位（Keeler，2004）。设原线宽为 1.2 Hz，选择 \(L=0.5\ \mathrm{Hz}\)，理想线宽变为 1.7 Hz。加窗削弱了长时间的信息，峰的底部与邻峰更容易重叠，不能将较平滑的外观直接解释为更高的分辨能力。
窗函数的参数单位必须核查。nmrglue 的 em 接口要求将 Hz 表示的增量除以 Hz 表示的谱宽；谱宽 4000 Hz、增量 0.5 Hz 对应参数 0.000125，直接传入 0.5 会产生完全不同的权重。高斯窗、正弦窗和含上升指数的分辨增强处理各有线形与噪声代价。比较小耦合间隔时，应同时保留轻度加窗或未加窗的结果，并记录所用函数，避免把窗函数产生的形状变化当成新的裂分证据。
4.6.4 傅里叶变换、数组排序与零填充
离散傅里叶变换（discrete Fourier transform，DFT）在离散频率位置计算各旋转分量的贡献；快速傅里叶变换（fast Fourier transform，FFT）是计算同一变换的算法。令加窗后的 \(N\) 个复数样本为 \(x_n\)，采样间隔为 \(\Delta t\) 秒；补零至 \(M\) 点后，按本节约定计算
展开说明与推导
\(k=0,\ldots,M-1\) 是频率数组索引，式中的 \(\Delta t\) 用于与连续变换的量纲衔接；NumPy 默认正变换没有该系数。正号旋转的输入在负号指数变换下出现在正频率侧。软件可能采用另一个指数符号、共轭输入或不同归一化，必须把相应频率轴一并定义（NumPy 开发者，在线文档）。改变这些约定而沿用旧坐标，可能把谱图镜像或改变振幅尺度。
标准 FFT 输出先放零频率，再放正频率和负频率；将零频率移到中央只是数组重排。NMR 程序还常按习惯将高化学位移画在左侧。重排或反转数据时，坐标数组必须同时变换；只反转曲线会制造错误的位移指认。由复数 FID 计算氢谱应保留复数变换，不能仅因最终显示实部就改用只接收实数的变换接口。
零填充（zero filling）只在已有 FID 后补上零。频域取样间隔从 \(1/(N\Delta t)\) 变为 \(1/(M\Delta t)\)，可以更密地描述同一条线，却没有新增采集信息，也没有延长真实的观测。教学设定 \(N=8192\)、谱宽 \(1/\Delta t=4000\ \mathrm{Hz}\)，采集时间按 \(N\Delta t\) 计为 2.048 s；补到 32768 点使频域网格由约 0.4883 Hz 变为约 0.1221 Hz。两个本来无法分开的宽峰，不会因此自动变为两个独立信号。
傅里叶变换保留实部、虚部与相位关系。取模 \(|X|\) 可使线形不再随整体相位旋转，但所得模谱具有不同的线形与噪声分布；常规吸收谱的相位、积分和裂分判读不能直接照搬到模谱。
习题 4.6-1
沿用谱宽 4000 Hz、8192 个复数点、补至 32768 点的教学设定。计算采样间隔、采集时间和补零前后的频域网格；对原半高全宽 1.2 Hz 的理想单指数线施加 0.5 Hz 指数加窗，求新线宽和 \(t=1\ \mathrm{s}\) 时的窗权重。
参考解答采样间隔为 \(1/4000=0.00025\ \mathrm{s}\)，即 0.25 ms。按约定采集时间为 \(8192\times0.00025=2.048\ \mathrm{s}\)。两网格分别为 \(4000/8192=0.48828125\ \mathrm{Hz}\) 与 \(4000/32768=0.1220703125\ \mathrm{Hz}\)，相差四倍。零填充没有改变原来的 8192 个测量值。式（4.6-2）给新线宽 \(1.2+0.5=1.7\ \mathrm{Hz}\)，窗权重为 \(\exp(-\pi\times0.5\times1)=0.20788\)。它只把后段信号乘小；增加的 0.5 Hz 线宽不能由补零抵消。实际谱若存在截断、非单指数衰减或相位误差，应另外检查这些条件。
4.6.5 相位校正的符号与枢轴
吸收线形（absorption lineshape）在中心两侧呈同号对称分布；色散线形（dispersion lineshape）在中心两侧异号。按式（4.6-1）与负号连续变换，令角频率差 \(q=2\pi(f-f_0)\)，单位为 \(\mathrm{rad\,s^{-1}}\)。零初相位时有 \(S(f)=a/(R+iq)=A+iD\)，其中 \(A=aR/(R^2+q^2)\)、\(D=-aq/(R^2+q^2)\)。对这个零初相位谱乘 \(\exp(i\theta)\) 后，实部成为 \(A\cos\theta-D\sin\theta\)；若原谱已有初相位 \(\phi\)，总旋转角为 \(\phi+\theta\)。一侧抬高、另一侧下凹的尾部因而可能来自相位，而非杂质峰。
展开说明与推导
零阶相位（zero-order phase）让所有频率旋转相同角度；一阶相位（first-order phase）让旋转角度随频率线性变化。本页为可核验算例定义
\(f_{\mathrm{p}}\) 是相位枢轴（phase pivot），\(W\) 为完整谱宽，二者以 Hz 计；\(\theta_0\) 是枢轴上的角度，\(\theta_1\) 是从低频边到高频边的总相位变化，均以弧度代入指数。练习以度显示时先乘 \(\pi/180\)。这一定义便于推理，并不等同于所有软件的参数：枢轴、索引方向和一阶相位的跨度约定都可能不同。某一软件的两个角度值不能脱离其操作顺序转用于另一个软件。
手动相位时，可先选一条不重叠、信噪比较高的窄线放置枢轴，调零阶使其底部与尾部接近吸收形，再在远离枢轴的清楚谱线检查一阶相位，往返小幅调整。应覆盖整个有效位移区，不能只使最高峰好看。Oxford 的教程将零阶和一阶分别描述为频率无关和频率相关（Mycroft，2024）。自动相位的目标函数可能受强溶剂峰、重叠、宽信号及基线牵引；算法结束后仍应检查弱峰和谱图两端。若不同区域始终需要互相矛盾的校正，回查滤波、导入约定与采集问题，避免用基线拟合掩盖色散尾部。
习题 4.6-2
以式（4.6-4）定义 \(u=(f-f_{\mathrm{p}})/W\)。设两条清楚窄线在 \(u=-0.25\) 和 \(u=0.25\) 的残余相位分别为 \(+5^\circ\) 与 \(+35^\circ\)，相位误差可用直线表示。求应施加的 \(\theta_0\)、\(\theta_1\)，并在两条线及枢轴核验。
参考解答误差的斜率为 \((35-5)/(0.25-(-0.25))=60^\circ\)，枢轴误差为两值的平均 \(20^\circ\)，所以误差为 \(20^\circ+60^\circ u\)。本节校正因子采用正号指数，应施加 \(\theta_0=-20^\circ\)、\(\theta_1=-60^\circ\)。在 \(u=-0.25\) 处校正角为 \(-20+15=-5^\circ\)，与原误差相加为零；在 \(u=0.25\) 处为 \(-20-15=-35^\circ\)，也相加为零；枢轴处为 \(-20^\circ\)。实现指数运算时这些角度必须换成弧度。若软件使用负号旋转或相反索引，应根据其定义换算，不能照填这两个值。
4.6.6 基线校正先确认无峰区域
基线校正（baseline correction）估计并扣除非信号背景。常见方法包括常数扣除、无峰节点之间的插值和低阶多项式拟合。先完成相位校正，再放大观察峰外区域；如果尾部呈异号并随峰出现，首先怀疑相位。Oxford 教程的节点法要求把节点放在真实基线上（Mycroft，2024）。节点落在峰底或宽信号中，拟合就会把一部分真实面积扣掉。
展开说明与推导
基线对积分的影响可直接计算。若以 ppm 为横坐标 \(x\)，区域 \([x_1,x_2]\) 内有恒定背景 \(c\)，该区域会额外获得面积
这里 \(c\) 的单位为谱强度，\(I\) 为谱强度乘 ppm。相同背景在更宽的积分窗口内产生更大的偏差，弱而宽的信号尤其敏感。若拟合曲线追随了宽峰，扣除的面积还可能为负方向的系统误差。提高多项式阶数能让曲线更贴合局部变化，却不证明局部变化属于仪器背景；只比较校正后曲线是否平直是不充分的。
OH、NH、交换中的信号或聚合物信号可能很宽。不能以“看起来像隆起”作为删除依据。保留校正前后谱及拟合背景，避开可疑区域，并比较替代节点或模型；若没有可信的无峰区域，记录基线不确定性。Mestrelab 的官方处理指南提醒，强平滑型基线操作可能损失峰尾或扰动积分（Mestrelab Research，2015）。
4.6.7 化学位移定标与强度尺度
化学位移定标（chemical-shift referencing）用可信参考物建立横轴。沿用本书与 IUPAC 的定义，\(\nu\) 和 \(\nu_{\mathrm{ref}}\) 分别为样品核与参考核的共振频率，单位一致：
展开说明与推导
\(\delta\) 是无量纲比值。参考物、溶剂、温度和参考方式应与结果一起报告；TMS 和水相常用的 DSS 都可能受介质条件影响，参考峰不能当作脱离条件的固定自然常数（Harris et al.，2008）。用残留溶剂定标时先确认溶剂及峰身份，避免把附近水峰、杂质或卫星误选为参考。改变横轴零点不改变线间 Hz 距离，因而不能用整体平移修正错误的谱宽。
定标还须与纵轴处理区分。显示放大、将某个积分赋值为 3、将最大峰归一为 1，分别服务于读图、相对氢数和强度展示；这些操作均不能补偿采集时的饱和、交换或选择性抑制。常规处理后的积分怎样划界，在4.7展开；对含量测量还要检查7.1的定量采集条件。
习题 4.6-3
一积分区域宽 0.40 ppm，观测面积为 100 强度单位·ppm，恒定基线为 10 强度单位。另有可信参考峰在旧坐标 0.07 ppm，应设为 0 ppm；待指认峰在旧坐标 3.27 ppm。求扣除背景后的面积、背景相对真实面积的偏差、定标后位移，并以 500 MHz 标称频率近似换算相对参考的 Hz 偏置。
参考解答背景面积为 \(10\times0.40=4\)，真实面积为 \(100-4=96\) 强度单位·ppm。不扣背景会高估 \(4/96\times100\%=4.1667\%\)；分母要用校正后的面积。按固定频率比例尺作通常的零点平移近似，横轴统一减去 0.07 ppm，待指认峰成为 \(3.27-0.07=3.20\ \mathrm{ppm}\)，即 \(\delta=3.20\times10^{-6}\)。近似取 \(\nu_{\mathrm{ref}}=500\times10^6\ \mathrm{Hz}\)，得到 \(\Delta\nu=1600\ \mathrm{Hz}\)。面积扣除与横轴平移是独立操作。若还更换定义式中的精确参考频率，位移须按 \(\delta_{\mathrm{new}}=(\delta_{\mathrm{old}}-\delta_{\mathrm{ref,old}})/(1+\delta_{\mathrm{ref,old}})\) 换算；本例的微小比例修正不影响所列有效位数。实际参考频率与标称频率有差别时，应使用采集元数据计算。
4.6.8 保存可重放的处理记录
处理记录应能让另一个人在同一原始数据上重建数组，随后再核对坐标与图像。以上 8192 点教学设定可逐步记录为：确认每点含两路正交分量；谱宽取 4000 Hz；按格式检查并处理数字滤波；选择 0.5 Hz 指数窗；在确认校正后点数不超过目标长度后补至 32768 点；执行已定义符号的复数 FFT；调整并保存相位；选无峰节点估计背景；用已确认参考峰设置 ppm 轴。每步都保留输入条件和输出检查，不能将一张导出的图片当成完整记录。
| 记录对象 | 必须保存的内容 | 重放时核查 |
|---|---|---|
| 原始输入 | 数据来源、文件标识或哈希、核种、复数点数、谱宽、载频 | 读取的数组与采集条件对应 |
| 滤波与时域 | 是否校正、参数和算法；校正后点数；窗函数、单位、首点缩放 | 没有重复校正或额外截断 |
| 变换与相位 | 补零长度、指数符号、归一化、排序；零阶、一阶、枢轴 | 谱线位置和整个谱区的相位 |
| 背景与定标 | 节点或掩膜、拟合方法、参考物与赋值、纵轴尺度 | 宽信号面积与坐标没有误改 |
| 结果版本 | 软件版本、操作顺序、处理数组、峰表和图像的对应关系 | 版本间差异可追溯到具体操作 |
展开说明与推导
重放并不要求所有处理参数适用于所有样品。一个模板适合窄线小分子，换到宽峰样品后可能需要减弱基线拟合；某次自动相位成功，也不能把同样的角度强加到下一次采集。保留两套合理参数的结果，检查关键位移、弱峰与面积结论是否稳定，比只保存一个“最终”图更有助于判断处理不确定性。
本书已在固定环境中用 nmrglue 0.12 运行 BMRB 乙醇条目 bmse000297 的一个真实 FID，执行滤波校正、1 Hz 加窗、补零与相位处理；数据来源、环境和输出对照见8.2 用 Python 处理真实 FID。它是单个存档采集的处理练习，积分比例未经定量方法验证，也不能代表所有厂商格式均已测试。本节的数值练习和两幅解析线形图仍为教学设定。保存原始输入与处理步骤后，再进入寻峰和积分。
关键术语
- 复数 FID （complex FID）
- 包含实部和虚部、保留旋转相位信息的时域采样信号。
- 群延迟 （group delay）
- 滤波系统引入的有效时间延迟；校正须依赖数据格式与相关参数。
- 窗函数 （window function）
- 乘在 FID 上、改变不同时刻权重的函数。
- 指数加窗 （exponential apodization）
- 以指数衰减权重削弱后段 FID；理想单指数线的半高全宽随之增加。
- 离散傅里叶变换 （DFT）
- 将有限采样数组变换到离散频率网格的复数运算。
- 零填充 （zero filling）
- 在 FID 后补零以增加频域取样密度，不增加测量信息。
- 零阶相位 （zero-order phase）
- 对整个复数谱施加相同旋转角度的校正。
- 一阶相位 （first-order phase）
- 旋转角度随频率线性变化的校正，数值须附带跨度约定。
- 相位枢轴 （phase pivot）
- 定义零阶相位、且一阶相位相对该点展开的频率位置。
- 基线校正 （baseline correction）
- 估计并扣除非信号背景，同时保留真实宽信号的处理。
- 吸收线形 （absorption lineshape）
- 理想单线中心两侧同号对称的谱线成分。
- 化学位移定标 （chemical-shift referencing）
- 按已确认的参考频率或参考峰赋值建立化学位移坐标。
参考文献与延伸阅读
- Keeler J. 2004. Fourier transformation and data processing. Understanding NMR Spectroscopy 公开讲义，第 4 章，University of Cambridge. 官方讲义 PDF。
- nmrglue 项目. Processing 1D Bruker Data；remove_digital_filter. 官方在线文档，核查于 2026-10-11。处理示例；数字滤波接口。
- nmrglue 项目. proc_base.em；proc_base.ps. 官方在线文档，核查于 2026-10-11。指数窗定义与单位；相位接口。
- NumPy 开发者. Discrete Fourier Transform. NumPy 官方文档，核查于 2026-10-11。DFT 符号、归一化与频率排序。
- Mycroft C. 2024. Introduction to MestreNova for NMR Analysis. University of Oxford NMR Facility，2024-03-06 培训资料。官方教程 PDF。
- Mestrelab Research. 2015. Starting Guide for the Automatic Analysis of Ligand Screening Data. 官方指南。本文只援引其对相位、基线、重处理模板及过强平滑的操作条件说明。
- Harris RK, Becker ED, Cabral de Menezes SM, Granger P, Hoffman RE, Zilm KW. 2008. Further conventions for NMR shielding and chemical shifts (IUPAC Recommendations 2008). Pure and Applied Chemistry, 80: 59–84. IUPAC 官方条目；DOI。