Explicit Ensemble Mean Synchronization for Time Scale Generation with Mixed Atomic Clock Ensembles

技術的概要：混合原子時計アンサンブルを用いた時間スケール生成のための明示的アンサンブル平均同期

問題定義
本論文は、セシウム型（Cs 型）と水素メーザー型（Hm 型）の時計を組み合わせる混合アンサンブルから信頼性の高い時間スケールを生成するという課題に取り組む。従来のカルマンフィルタリング（CKF）は時計偏差の推定に広く用いられているが、この文脈では決定的な限界に直面する。すなわち、測定可能なものが時計の読み値の差のみであるため、混合アンサンブルモデルは検出不可能である。この検出不可能性が、CKF 内の誤差共分散行列の発散を引き起こし、数値的不安定性と不正確な推定をもたらす。さらに、既存の時間スケール生成手法は、しばしば暗黙的なアンサンブル平均に依存するか、均質なアンサンブル（Cs 型のみ、あるいは Hm 型のみのいずれか）に限定されており、混合アンサンブルの相補的な安定性特性を十分に活用できていない。著者は以下の 2 つの主要な課題を特定している。(1) 発散項を含まない混合アンサンブル向けのカルマンフィルタリングアルゴリズムの開発、および (2) 定義された重み付き平均に時計偏差を同期させつつ、短期間または長期間の区間における周波数安定性（ハダマード分散で測定）を最適化する時間スケールの生成。

手法
提案される解決策は、状態空間変換と明示的制御設計を組み合わせた 2 段階のアプローチである。

1. 観測可能標準形分解と変換カルマンフィルタ（TKF）：
   著者はまず、混合アンサンブルモデルに対して観測可能標準形分解を適用する。状態空間を観測可能成分と非観測可能成分に変換することで、新たなカルマンフィルタリングアルゴリズムを導出する。この変換された領域において、発散を引き起こす非観測状態に関連する誤差共分散行列は、ゲイン計算から分離される。その結果得られる変換カルマンフィルタ（TKF）およびその定常状態版は、観測可能状態の収束する誤差共分散行列のみを利用し、発散項を排除して数値的安定性を確保する。

2. 明示的アンサンブル平均同期アルゴリズム：
   TKF からの状態推定値を用いて、著者は明示的アンサンブル平均同期アルゴリズムを提案する。制御入力は、個々の時計の時間偏差を「非観測可能な重み付き平均位相状態」に同期するように設計される。

   • 同期： 状態フィードバック制御則を適用して同期誤差をゼロに駆動し、アンサンブルがユーザー指定の重みベクトル $q$ で定義される共通の時間スケールに収束することを保証する。
   • 周波数安定性の最適化： 短期安定性と長期安定性の間のトレードオフに対処するため、著者は非観測状態の自由運転ダイナミクスを解析する。生成された時間スケールのハダマード分散（HVAR）を最小化する最適な重みベクトル $q$ を導出する。
     • 短期安定性の場合、サンプリング間隔 $\tau \to 0$ において HVAR を最小化するように重みベクトルを導出する。これにより、白色周波数ノイズの低い時計（通常は Hm 型）を優先的に扱う。
     • 長期安定性の場合、$\tau \to \infty$ において HVAR を最小化するように重みベクトルを導出する。これにより、ランダムウォーク周波数ノイズの低い時計（通常は Cs 型）を優先し、周波数ドリフトを有する Hm 型時計にはゼロの重みを割り当てる。

主要な貢献

• 混合アンサンブルのための安定フィルタリング： 本論文は、検出不可能な混合時計アンサンブルにおける従来のフィルタに内在する発散問題を解決する、観測可能標準形分解に基づくカルマンフィルタリングアルゴリズムの厳密な導出を提供する。
• 明示的同期制御： 従来の「暗黙的」アンサンブル平均アプローチとは異なり、この研究はユーザーが重みベクトルを定義して時計を特定の目標状態に同期させる明示的アルゴリズムを導入する。
• 最適化された時間スケール生成： 著者は、特定の時間区間に対して生成された時間スケール（GTS）の周波数安定性を最適化するための重みベクトルの選択方法を実証する。適切な重みを選択することで、個々の時計や単純な平均と比較して、GTS が優れた短期安定性（Hm 型時計を活用）または優れた長期安定性（Cs 型時計を活用）を達成し得ることを示す。

結果
提案アルゴリズムの有効性は、10 個の時計（Cs 型 7 個、Hm 型 3 個）からなるアンサンブルを含む数値例によって検証された。

• 短期最適化： 短期安定性向けに最適化された重みベクトル（$q^0_H$）を適用した際、シミュレーションは制御された時計の時間偏差が、最適化された自由運転ダイナミクスの仮想軌跡に密接に従うことを示した。その結果、1,000 秒の区間において、自由運転時計と比較して理想時間からの発散が著しく減少した GTS が得られた。
• 長期最適化： 長期安定性向けに最適化された重みベクトル（$q^\infty_H$）を適用した際、GTS は $10^5$ 秒の区間において優れた周波数安定性を示した。HVAR 曲線は、生成された時間スケールが長期間にわたり個々の自由運転時計よりも優れた安定性を維持したことを確認した。
• ハダマード分散： HVAR 解析は、提案手法が対象とする時間区間に対して分散を成功裏に最小化し、重みベクトルの理論的最適化を検証したことを確認した。

意義と主張
本論文は、提案する「明示的アンサンブル平均同期アルゴリズム」が、混合原子時計アンサンブルにおける時間スケール生成のための堅牢な解決策を提供すると主張する。誤差共分散行列の発散を排除することで、この手法は従来のフィルタが失敗する場面で数値的安定性を保証する。さらに、ユーザーが重みベクトルを明示的に定義することを可能にすることで、このアルゴリズムは生成された時間スケールの周波数安定性を特定の運用ニーズ（短期の精度を優先するか、長期の精度を優先するか）に合わせて調整する柔軟なメカニズムを提供する。著者は、現在の研究が短期または長期の安定性のいずれかを個別に最適化するものであるが、将来の研究では両方の区間に対して同時にバランスの取れた安定性プロファイルを実現する戦略に焦点を当てることを付記している。