Extended Kalman Smoothing of Free Spin Precession Signals for Accurate Magnetic Field Determination

技術要約：高精度な磁場決定のための自由スピン歳差運動信号に対する拡張カルマン平滑化

問題提起
静的な磁束密度の正確な決定は、関係式 $f_L = \frac{\gamma}{2\pi}B_0$ を通じて、偏極核（$^3$Heなど）の自由スspin歳差運動周波数（$f_L$）を測定することに依存している。$^3$Heの磁気回転比 $\gamma$ は極めて高い精度（相対不確かさ $<10^{-9}$）で既知であるが、磁場測定の全体的な精度は、現在、核磁気共鳴（NMR）信号からの周波数推定の精度によって制限されている。従来の主要な手法である非線形最小二乗法（具体的には正弦・余弦適合、またはSCF）は、自由歳差運動減衰（FPD）信号に適用する場合、重大な課題に直面する。これらの信号は、指数関数的な振幅減衰と、磁場のドリフトによる時間依存の周波数シフトを示す。従来のアプローチでは、一定の振幅と周波数を仮定するために、データを固定長のブロックに分割する必要がある場合が多い。これはトレードオフを生み出す。すなわち、短いブロックはドリフトに適応できるが、ノイズによる分散が高くなる。一方で、長いブロックはノイズを低減できるが、急速な周波数変化や振幅減衰の追従に失敗し、モデルの不一致やバイアスを引き起こす。

手法
著者らは、FPD信号から時変の周波数および振幅を推定するために、拡張カルマン平滑化（EKS）の新規な応用を提案している。リアルタイムで動作する標準的なカルマンフィルタとは異なり、EKSは利用可能なすべての測定値（過去および未来）を利用して状態推定を洗練させるため、センサデータの後処理に最適である。

本手法は、以下の主要なコンポーネントで構成される：

1. 状態空間モデリング： システムは時間領域ではなく周波数領域でモデル化される。信号は重複しないブロックに分割される。各ブロックについて、状態ベクトル $x_k$ は振幅 $A_k$、周波数オフセット $\delta f_k$、および位相 $\phi_k$ を含む。これらの状態の進化は、振幅と周波数については積分ランダムウォークとして、位相については数値積分としてモデル化される。
2. 観測方程式： 時領域の信号を直接フィッティングするのではなく、EKSはモデルの予測値と信号ブロックの離散フーリエ変換（DFT）係数を比較する。これにより、時領域の信号における加法的ガウスノイズが、周波数領域における独立したガウスノイズ成分へと変換され、ノイズ分散が大幅に減少する。
3. パラメータ最適化： モデルの静的パラメータ（プロセスノイズ共分散 $Q$、測定ノイズ共分散 $R$、および初期値）は、期待値最大化（EM）アルゴリズムを用いて自動的に最適化される。これにより、ブロック長やノイズパラメータの手動チューニングなしに、信号特性に適応することが可能となる。
4. 参照手法： EKSは、線形パラメータと非線形パラメータを分離するために可変投影を用いる標準的なSine-Cos-Fit（SCF）法と比較検証される。SCFアプローチでは、シミュレーションにおいて、平均二乗誤差（MSE）を最小化するように最適なブロック長が事後的に決定され、これが従来手法の「ベストケース」となるシナリオを代表する。

主な貢献

• アルゴリズムの革新： 本論文は、フーリエ係数を観測モデルとして利用することで、緩やかに変化する周波数と減衰する振幅を持つ調和信号の追跡に特化した、周波数領域EKSフレームワークを導入している。
• 自動適応： 本手法は、最小二乗法において重要かつ困難なチューニングステップである、ブロック長やノイズパラメータの手動選択を排除する。EMアルゴリズムがプロセスノイズおよび測定ノイズ共分散の最適化を処理する。
• 包括的なベンチマーク： 本研究は、広範なシミュレーション（信号対雑音比（SNR）および周波数ドリフト拡散定数の変化を含む）と、$^3$He磁力計からの実際の実験データを用いて、EKSとSCFの厳格な比較を行っている。

結果

• シミュレーション研究： 幅広いSNRおよび周波数ドリフト率の範囲において、EKSはSCFと比較して系統的に推定誤差を低減した。中程度から高いSNRおよび進化する周波数の領域において、EKSは周波数推定を約2倍改善した（$\rho_{RMSE,f} \approx -1$ で示される）。SCFは、非常に低いSNRかつ周波数ドリフトが無視できる領域においてのみ、クラメール・ラオの下限に近づくことでEKSを上回った。
• 振幅推定： EKSは、信号が指数関数的に減衰する場合に固定ブロック内での一定振幅の仮定が劣化するため、振幅推定においてSCFを大幅に上回った。
• 実験的検証： $^3$Heスピン歳差運動測定（継続時間 約4時間）の実験データに適用したところ、EKSは周波数変動の追跡において優れた性能を示した。大きなブロック長を用いたSCFの推定は、急速な周波数変化（例：約300秒時点での上方への急上昇）を解像できず、信頼性の低い不確かさ区間を示したが、EKSは一貫した追跡を維持した。特筆すべきは、EKSが、ブロック長をほぼ2桁小さくしてもなお、SCFと同等の不確かさ区間を達成したことであり、これは精度を犠牲にすることなく高い時間分解能を実現したことを示している。

意義
本論文は、EKSベースの解析が、NMRベースの磁力計における周波数推定の精度を大幅に向上させ、これは絶対磁場決定の基礎となるステップであることを主張している。ノイズを効果的にフィルタリングし、信号の減衰と周波数ドリフトに自動的に適応することで、本手法は高精度測定の安定性と正確性を高める。著者らは、この多用途なアプローチが、磁場変動の追跡において特に価値が高いことを強調している。これは現代の磁力計アプリケーションにおける重要な要件である。本研究は、システム依存のパラメータに依存しない、堅牢で使いやすい実装を提供しており、超高精度NMR測定における現在の最先端の最小二乗法に代わる有望な候補となっている。