Optimal Calibration of the Endpoint-corrected Hilbert Transform

技術要約：エンドポイント補正型ヒルベルト変換の最適キャリブレーション

問題提起
振動の瞬時位相を正確かつ低遅延で推定することは、クローズドループ型のセンシングおよびアクチュエーション、特に位相ロック型神経刺激やリアルタイム制御システムにおいて極めて重要である。標準的な手法では、スライディングウィンドウを用いたヒルベルト変換（HT）により解析信号を構成する。しかし、離散時間実装における離散フーリエ変換（DFT）は、ウィンドウ内における周期性を仮定している。ウィンドウ境界において信号が完全に周期的でない場合、この仮定によって不連続性が生じ、ギブス現象（Gibbs ringing）を引き起こし、ウィンドウ端における振幅および位相の著しい歪みを招く。リアルタイムシステムは最新のサンプル（エンドポイント）に依存するため、これらの境界アーティファクトは、極めて重要な位相推定値を直接的に損なうことになる。

これらの問題を軽減するために、因果的な狭帯域フィルタを解析スペクトルに適用するエンドポイント補正型ヒルベルト変換（ecHT）が導入されたが、その系統的なエンドポイント歪みについては、原理的な閉じた形式による解析が行われていなかった。従来の評価は主に経験的なものに留まっており、残留バイアスの性質や最適な補正戦略は未定義のままであった。

手法
本研究は、DFTに基づく解析的手法を用いてecHTの特性評価を行う。著者らは、単一正弦波信号に対してエンドポイント演算子を解析的に導出し、ecHTの出力を以下の2つの明確な成分に分解した：

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ベースのecHTパイプラインに対する単一の複素数乗算のみであり、無視できる程度である。これにより、リアルタイム互換性が維持される。

結果
cecHTの有効性は、シミュレーションおよび実世界のデータを通じて検証された：

• シミュレーション: 単一正弦波シミュレーションにおいて、cecHTは平均位相誤差を約8.53°（未補正のecHT）から0.40°へと一貫して減少させ、最大誤差削減幅を11.70°から1.23°へと減少させた。このキャリブレーションにより、アルゴリズムはパラメータ（帯域幅、フィルタ次数、周波数ドリフト）の選択に対して実質的に独立し、位相精度を犠牲にすることなく、他の制約（例：ノイズ耐性）のためにパラメータを調整することが可能となった。
• EEGアルファ波: EEGアルファ波データセットに適用した結果、cecHTは平均位相誤差を4.9°から0.0°（$p < 10^{-6}$）へと減少させると同時に、位相ロック値（PLV）を維持し、位相ラグ指数（PLI）を大幅に減少させた。これは位相の非対称性が除去されたことを示している。
• 本態性振戦: 振戦データ解析において、キャリブレーションは平均位相誤差を-15.9°から5.2°へと減少させた。スペクトルドリフトを追跡するために中心周波数（$f_0$）を周期的に更新する手法と組み合わせたところ、平均誤差はさらに0.0°へと減少した。これは、残留誤差が主にドリフトに起因するものであり、キャリブレーション手法自体の固有の性質ではないことを示している。

意義および主張
本論文は、cecHTが既存のecHTパイプラインに対する「実用的なドロップイン・リプレースメント（そのまま置き換え可能な代替手段）」を提供し、明確な設計規則を持つ統計的に堅牢な推定器であることを主張している。補正可能なバイアスと固有のリークを分離することで、本フレームワークは、実用的なツールとしての側面と、位相推定の限界を理解するための診断的な言語の両方を提供する。

著者らは、キャリブレーションはフィルタの位相応答にエンコードされた系統的バイアスを排除できる一方で、スペクトルリークや非定常性に起因する補正不可能な誤差フロアを除去することはできない点を強調している。したがって、本手法は安定した条件下で偏りのない推定を実現するものであり、残るジッターは推定器の限界または信号自体の不安定性を反映している。本研究は、cecHTを、単純な因果的フィルタと複雑な適応モデルの中間に位置する実用的なソリューションとして位置づけており、周波数の安定性とアーティファクトへの堅牢性が管理されている限り、リアルタイムの神経信号処理に適しているとしている。