Optimal Calibration of the Endpoint-corrected Hilbert Transform

技术摘要：端点修正希尔伯特变换的最优校准

问题陈述
精确、低延迟地估计振荡的瞬时相位，对于闭环感知与驱动（特别是相位锁定神经刺激和实时控制系统）至章重要。标准方法是通过在滑动窗口上构建希尔伯特变换（HT）来构建解析信号。然而，使用离散傅里叶变换（DFT）的离散时间实现假设窗口内具有周期性。当信号在窗口边界处并非完美周期时，该假设会产生不连续性，从而导致吉布斯振铃现象（Gibbs ringing），这会在窗口边缘引起显著的幅度与相位失真。由于实时系统依赖于最新的样本（即端点），这些边界伪影会直接破坏关键的相位估计。

虽然端点修正希尔伯特变换（ecHT）已被引入，通过应用因果窄带滤波器来缓解这些问题，但其系统的端点失真一直缺乏原则性的闭式分析。以往的评估主要基于经验，导致残余偏差的性质以及最优校准策略仍未定义。

方法论
本研究通过基于 DFT 的分析方法对 ecHT 进行了表征。作者针对窄带单音信号推导了端点算子，将 ecHT 输出分解为两个截然不同的部分：

1. 确定性复增益 ($G_+$)： 一个正频率项，会诱发可校准的幅度和相位偏差。
2. 残余泄漏项 ($G_-$)： 一个负频率项，它逃逸了解析掩模和滤波器，产生随瞬时相位变化的相位相关纹波。

基于这种分解，本文推导出了均方误差（MSE）最优标量校准。作者证明，对于任何 ecHT 实现，都存在一个唯一的复标量 $C_{opt}$，能够最小化校准输出与理想解析信号之间的误差。该标量通过下式解析得出：
$$C_{opt} = \frac{E[\hat{Z}^* Z]}{E[|\hat{Z}|^2]}$$
其中 $Z$ 是真实的端点，$\hat{Z}$ 是 ecHT 估计值。从几何角度看，这代表了一个正交投影（标量维纳滤波器），用于提取最大线性可恢复的信息。作者证明，$C_{opt}$ 仅取决于正频率增益和负频率泄漏，可以在无需数据驱动拟合的情况下，针对任何窗口长度、滤波器带宽、阶数和中心频率进行预计算。

核心贡献

1. 解析分解： 本文确立了 ecHT 误差由可校准的系统偏差（由 $G_+$ 决定）和不可消除的方差底限（由 $G_-$ 决定）组成。
2. 最优校准 (cecHT)： 作者引入了“校准后的 ecHT”（cecHT），通过应用推导出的标量 $C_{opt}$ 来消除系统相位偏差。这使得在稳定条件下能获得近乎为零的平均相位误差。
3. 设计规则： 研究提供了端点相位和幅度误差的显式界限，并将实际设计参数（窗口长度、滤波器带宽/阶数以及中心频率失配）通过端点群延迟与残余偏差联系起来。
4. 计算效率： 校准过程增加的开销极小（仅为一次复数乘法），保持了与现有的 $O(L \log L)$ FFT-based ecHT 流程的实时兼容性。

结果
通过模拟和真实数据验证了 cecHT 的有效性：

• 模拟： 在单音模拟中，cecHT 始终将平均相位误差从约 8.53°（未校准的 ecHT）降低至 0.40°，并将最大误差从 11.70° 降低至 1.23°。校准使算法几乎不受参数选择（带宽、滤波器阶数、频率漂移）的影响，允许在不牺牲相位精度的情况下针对其他约束（如噪声鲁棒性）进行参数调节。
• EEG Alpha 波： 在应用到 EEG Alpha 波数据集时，cecHT 将平均相位误差从 4.9° 降至 0.0° ($p < 10^{-6}$)，同时保留了相位锁定值（PLV）并显著降低了相位滞后指数（PLI），表明消除了相位不对称性。
• 特发性震颤： 在震颤数据分析中，校准将平均相位误差从 -15.9° 降至 5.2°。当结合定期更新中心频率 ($f_0$) 以追踪频谱漂移时，平均误差进一步降至 0.0°，这证明残余误差主要是由漂移引起的，而非校准方法本身的固有问题。

意义与主张
论文声称，cecHT 为现有的 ecHT 流程提供了一个“实用的掉入式替代方案（drop-in replacement）”，提供了一个具有明确设计规则的统计稳健估计器。通过将可校减偏差与本质泄漏分离，该框架既提供了一个用于改进的实用工具，也提供了一种理解相位估计极限的诊断语言。

作者强调，虽然校准可以消除由滤波器相位响应编码的系统偏差，但无法消除由频谱泄漏和非平稳性引起的不可消除的误差底限。因此，该方法在稳定条件下实现无偏估计，剩余的抖动反映的是估计器的局限性或真实的信号不稳定性。该研究将 cecHT 定位为简单因果滤波器与复杂自适应模型之间的中间解决方案，只要管理好频率稳定性和伪影鲁棒性，它便适用于实时神经信号处理。