変調・符号化#171

搬送波周波数捕捉 — ロックする前に、まず見つけ出す

PLLは位相誤差が小さい前提の追尾装置だが、その前提はどこから来るのか。打ち上げ直後など周波数が数十kHzのオーダーで不明な局面で、地上局受信機が広い不確かさ範囲から搬送波をまず粗く見つけ出す「捕捉」問題を、周波数掃引法とFFTベース高速捕捉の数式で理解する。

前提知識: PLL — 弱い搬送波を追いかけ続けるフィードバック制御

周波数捕捉アクイジション周波数掃引FFTドップラー不確かさ

この回で学ぶこと

PLLの回では、位相検波器の非線形特性 $v_d(t) = K_d\sin\phi(t)$ を $\sin\phi(t)\approx\phi(t)$ と線形近似するところから、すべての解析が始まりました。この近似が成り立つのは、位相誤差 $\phi(t)$ がすでに小さい——つまり VCOの周波数がすでに受信搬送波の近くにあり、ループがすでにロックしている(あるいはロックに近い)状態にある——ときだけです。ダンピング係数、ループ帯域幅、終値定理による定常誤差の議論は、すべてこの「すでに近くにいる」という前提の上に築かれていました。

しかし、この前提はどこから来るのでしょうか。地上局の受信機の局部発振器(VCO)は、電源を入れた瞬間、探査機の送信周波数についていったい何を知っているのでしょうか。実は、何も知りません。あるいは、せいぜい「大まかにこのあたりだろう」という予測値しか持っていません。

その予測がどれくらい不確かなのかを考えてみましょう。探査機の送信周波数 $f_T$ 自体は搭載発振器の仕様で数Hz〜数十Hzオーダーの精度で分かっています。しかし受信機が実際に地上で観測する周波数は、ドップラー効果を受けた

f_R = f_T\left(1 + \frac{v_{\text{rel}}}{c}\right)

であり、 $v_{\text{rel}}$ (探査機と地上局の視線方向相対速度)を事前にどれだけ正確に予測できるかが鍵になります。軌道が精密に確定していれば $v_{\text{rel}}$ の予測誤差は小さく、ドップラーシフトの不確かさも数百Hz程度に収まります。しかし、打ち上げ直後(まだ軌道決定が済んでいない)、軌道投入直後(まさに軌道が変わりつつある)、あるいは通信途絶からの再捕捉(セーフモード復帰など、探査機の姿勢や軌道の履歴が不明な局面)では、 $v_{\text{rel}}$ の予測誤差が大きくなり、受信周波数の不確かさは容易に数kHz〜数十kHzのオーダーに達します。

この不確かさの幅は、PLLの回で見たループの線形動作範囲——位相誤差が $\sin\phi\approx\phi$ とみなせる範囲、周波数で言えばループ帯域幅 $B_L$ (せいぜい数Hz〜数十Hz)——と比べて、桁違いに広いことに注意してください。VCOの周波数がこの広い不確かさ範囲のどこか(受信搬送波から数kHz以上離れた場所)にある状態で通常のPLLの電源を入れても、位相検波器は $\sin\phi(t)$ の非線形領域で高速に振動するだけで、ループフィルタはその振動を平均してほぼゼロの補正しか出さず、VCOはいつまで経っても正しい周波数に近づいてくれません。線形PLLは「近くにいる」ことを仮定した追尾装置であって、「探す」ための装置ではないのです。

そこで必要になるのが、PLLの線形追尾に先立つ捕捉(アクイジション、acquisition)フェーズです。これは「周波数が数kHz〜数十kHzの範囲のどこかにあるらしい」という粗い情報だけを頼りに、まず受信搬送波のおおよその周波数を見つけ出し、PLLがロック(追尾)を始められる線形動作範囲内までVCOを導く問題です。この回では、この捕捉問題を2つの代表的な手法——周波数掃引法とFFTベースの高速捕捉——を通じて数式で理解します。

直感的な導入: 暗闇でラジオの周波数を合わせる

日常の感覚に近い例で考えてみましょう。古いアナログラジオのダイヤルを回して放送局を探す作業を思い浮かべてください。目当ての局がどのあたりの周波数にあるかまったく分からないとき、あなたは2つの方法を取れます。

ダイヤルをゆっくり端から端まで回し、音が聞こえた瞬間に手を止める。 これが周波数掃引(スイープ)法に相当します。速く回しすぎると聞こえるはずの局を一瞬で通り過ぎてしまい、逆に遅すぎると全部の周波数を確認し終えるまでに時間がかかりすぎます。
もし魔法のような装置があって、ダイヤル全域の電波強度を一瞬で全部同時に測定し、いちばん強いところを教えてくれたら。 これがFFTベースの捕捉法に相当します。全周波数を並列に評価できるため、原理的にはスイープよりずっと速く見つけられます。

深宇宙受信機における「捕捉」は、まさにこの探索問題です。以下、それぞれの手法を定式化していきます。

問題の定式化: 周波数不確かさとSNR

受信信号(搬送波成分のみ、簡略化のため位相は無視)を

r(t) = A\cos\big(2\pi (f_0+\Delta f) t + \theta_0\big) + n(t)

と書きます。 $f_0$ は受信機が予測している(公称の)周波数、 $\Delta f$ は未知の周波数オフセット、 $n(t)$ は雑音です。捕捉問題とは、事前情報として

|\Delta f| \le \Delta f_{\max}

(不確かさの半幅 $\Delta f_{\max}$ 、たとえば数十kHz)しか分かっていない状態で、 $\Delta f$ の真値を、PLLが線形にロックできる範囲(おおよそループ帯域幅 $B_L$ の数倍程度)まで絞り込むことです。信号対雑音比の指標として、PLLの回で導入したのと同様に、受信搬送波電力 $P_C$ と片側雑音スペクトル密度 $N_0$ を使います。捕捉アルゴリズムが動作するために必要な観測時間を $T$ とすると、その間に積分できる実効的なSNRは、周波数が正しく合っている(あるいはFFTのビンが正しい周波数を捉えている)と仮定すれば近似的に

\text{SNR} \approx \frac{P_C\, T}{N_0}

のオーダーで決まります。観測時間 $T$ を長く取るほど積分利得が上がり弱い信号でも検出できますが、逆に $T$ が長いほど「1つの周波数(あるいは1つの掃引ステップ)を調べるのに時間がかかる」ため、探索全体が遅くなります。この検出感度と探索速度のトレードオフが、以下で見る2つの手法に共通して現れる基本原理です。

周波数掃引(スイープ)法

もっとも古典的な捕捉方式は、VCOの周波数を時間とともに一定の速度で連続的に変化させながら、各瞬間の周波数でどれだけ強い信号が「見えているか」を評価する方法です。

掃引信号のモデル

VCOの瞬時周波数を時間とともに直線的に変化させます。

f_{\text{VCO}}(t) = f_0 + \dot f_{\text{sw}}\, t

$\dot f_{\text{sw}}$ が**掃引速度(スイープレート)**で、単位はHz/sです。これは、PLLの回で扱った「一定のドップラーレート $\dot f$ に対する定常位相誤差」の議論と数学的に同じ形の入力ですが、ここでの $\dot f_{\text{sw}}$ はドップラーによって受信機に押し付けられるものではなく、受信機自身が能動的にVCOへ加える探索用の掃引である点が異なります。

このVCO出力と受信信号をミキシングし、ローパスフィルタを通すと、ビート周波数(差周波数)

f_{\text{beat}}(t) = \Delta f - \dot f_{\text{sw}}\, t

を持つ信号が得られます(VCOの走査開始点を $\Delta f$ が測られる原点に合わせています)。この $f_{\text{beat}}(t)$ がゼロに近づく瞬間、すなわち $t^* \approx \Delta f/\dot f_{\text{sw}}$ の付近で、ビート周波数がループのごく低い周波数域(直流に近い成分)に入り込み、エネルギー検出器の出力にピークが現れます。このピークを検出した時点でVCOの走査を止め(あるいはピーク近傍の周波数を推定して直接VCOをそこにセットし)、PLLの線形追尾ループへと制御を引き渡します。

掃引速度と検出感度のトレードオフ

掃引法の性能を左右する核心的な問題は、**「信号のそばをどれだけの速さで通り過ぎるか」**です。検出器(典型的には帯域幅 $B_i$ のバンドパスフィルタとその後段のエネルギー検出、あるいは相関器)が信号に反応するには、その検出器の帯域内に信号がとどまる時間がある程度必要です。検出器の実効帯域幅を $B_i$ とすると、信号がその帯域を横切るのにかかる時間(通過時間、dwell timeに相当)は

\tau_{\text{dwell}} \approx \frac{B_i}{\dot f_{\text{sw}}}

で近似できます。この $\tau_{\text{dwell}}$ の間に検出器が積分できるエネルギーがSNRを決めるので、掃引速度 $\dot f_{\text{sw}}$ が大きい(速く走査する)ほど $\tau_{\text{dwell}}$ は短くなり、検出に使える積分時間が減って感度が落ちます。定性的には、検出に必要な最小SNRを確保するための条件は

\dot f_{\text{sw}} \lesssim \frac{P_C}{N_0}\cdot\frac{1}{\text{SNR}_{\min}}

のような形で、受信電力が弱いほど掃引速度を遅くしなければならないという関係になります(具体的な係数は検出器の構成やしきい値設計に依存しますが、この反比例の関係は掃引法に共通する本質です)。

一方で、掃引速度を遅くすると、不確かさ範囲全体 $2\Delta f_{\max}$ を一往復するのにかかる時間

T_{\text{acq}} \approx \frac{2\Delta f_{\max}}{\dot f_{\text{sw}}}

が長くなります。つまり掃引法には、

\dot f_{\text{sw}} \text{ を上げる} \Rightarrow T_{\text{acq}} \text{ は短くなるが検出感度が落ち、信号を見逃す(通り過ぎる)危険が増す}

\dot f_{\text{sw}} \text{ を下げる} \Rightarrow \text{検出感度は上がるが } T_{\text{acq}} \text{ が長くなる}

という直接のトレードオフがあり、掃引速度は「見逃さない最大速度」の範囲内で、捕捉時間をなるべく短くするように選ばれます。さらに、 $\Delta f_{\max}$ (不確かさの半幅)が大きいミッションフェーズほど、同じ掃引速度でも探索範囲が広がる分、捕捉時間は線形に伸びます。これが「捕捉範囲と捕捉時間のトレードオフ」の掃引法における具体的な現れです。

FFTベースの高速捕捉

周波数掃引法の弱点は、ある瞬間には不確かさ範囲のうちのごく一部(検出器の帯域幅 $B_i$ の分)しか調べておらず、残りの周波数については何も情報を得ていない、という点です。順番に一つずつ調べていく逐次探索である以上、原理的に時間がかかります。

これに対して、受信信号を短時間だけサンプリングし、その離散フーリエ変換(DFT、実装上は高速フーリエ変換FFT)を計算すれば、不確かさ範囲に含まれるすべての周波数ビンのエネルギーを一度の演算で同時に評価できます。これがFFTベースの高速捕捉の核心です。

定式化

受信信号を複素ベースバンドに変換した後、時間長 $T_{\text{obs}}$ にわたって(実効サンプリング間隔 $T_s$ 、サンプル数 $N=T_{\text{obs}}/T_s$ )サンプリングした離散信号を $x[n]$ とします。そのDFTは

X[k] = \sum_{n=0}^{N-1} x[n]\, e^{-j2\pi kn/N}, \qquad k = 0, 1, \dots, N-1

で、各ビン $k$ は周波数

f_k = k\cdot \frac{1}{N T_s} = k\cdot \frac{1}{T_{\text{obs}}}

に対応します。すなわち、周波数分解能はビン幅

\delta f = \frac{1}{T_{\text{obs}}}

によって決まり、観測時間 $T_{\text{obs}}$ を長く取るほど周波数分解能は細かくなります。搬送波が実際に存在すれば、 $|X[k]|^2$ (パワースペクトル)は真の周波数オフセット $\Delta f$ に対応するビン付近で鋭いピークを示すので、全ビンにわたって $|X[k]|^2$ を走査し、最大値を与える $k^*$ を見つければ、

\widehat{\Delta f} = \frac{k^*}{T_{\text{obs}}}

がオフセットの推定値になります。

なぜ掃引法より速いのか

掃引法との本質的な違いは、並列性にあります。掃引法は不確かさ範囲 $2\Delta f_{\max}$ を検出器帯域幅 $B_i$ の「窓」で逐次的に(時間をずらしながら)見ていくのに対し、FFTは1回の観測(長さ $T_{\text{obs}}$ )で、幅 $\delta f = 1/T_{\text{obs}}$ の $N$ 個のビンすべてを同時に評価します。これはちょうど、暗い部屋を懐中電灯で一点ずつ照らして探すか(掃引)、部屋全体を一瞬で照らすフラッシュを焚いて写真を撮るか(FFT)、という違いに相当します。

必要な周波数分解能 $\delta f$ (捕捉後にPLLへ引き渡すために十分な精度、あるいはさらに細かい探索を続けるための粗い絞り込みの精度)を得るのに必要な観測時間は $T_{\text{obs}} = 1/\delta f$ で決まっており、この $T_{\text{obs}}$ 一回分の観測だけで不確かさ範囲全体をカバーできます。これに対して掃引法は、同程度の周波数分解能(検出器帯域幅 $B_i\sim \delta f$ )を得ようとすると、不確かさ範囲全体をカバーするのに $2\Delta f_{\max}/\delta f$ 回分の「窓」を順番に見て回る必要があり、総観測時間はおおよそ

T_{\text{acq}}^{\text{sweep}} \sim \frac{2\Delta f_{\max}}{\delta f}\cdot T_{\text{obs,窓}} \sim \frac{2\Delta f_{\max}}{\delta f^2}

のオーダーになります(窓1つあたりの観測時間も概ね $1/\delta f$ 程度必要になるため)。一方FFTベースの方法は

T_{\text{acq}}^{\text{FFT}} \sim T_{\text{obs}} = \frac{1}{\delta f}

で済みます。同じ分解能 $\delta f$ を得るのに、掃引法はおよそ $2\Delta f_{\max}/\delta f$ 倍(不確かさ範囲をビン幅で割った、探索すべきビンの本数)だけFFTより時間がかかる計算になり、これがFFTベースの捕捉が「高速」と呼ばれる理由です。もちろん現実の実装では、FFTの演算量(点数 $N$ に対して $O(N\log N)$ )や、 $N$ が大きくなったときのメモリ・計算資源の制約が別のボトルネックになりうるため、単純に「FFTなら常に高速」とは言い切れませんが、探索原理としての並列性の優位は明確です。

捕捉範囲と捕捉時間のトレードオフ(一般形)

掃引法・FFT法いずれにおいても共通する一般的な関係をまとめておきます。捕捉に使える時間 $T_{\text{acq}}$ 、達成したい周波数分解能(捕捉精度) $\delta f$ 、探索すべき不確かさ範囲 $\Delta f_{\max}$ の間には、おおよそ

T_{\text{acq}} \propto \frac{\Delta f_{\max}}{\delta f} \times (\text{1ステップあたりの必要観測時間})

という関係があり、不確かさ範囲が広いほど、また要求される捕捉精度(分解能)が細かいほど、捕捉には長い時間がかかります。実務上は、この捕捉時間を、要求される検出感度(信号がどれだけ弱くても取りこぼさないか)とのトレードオフの中で許容範囲に収めるように、観測時間・掃引速度・FFTの点数といったパラメータを設計します。捕捉が済んだ後は、粗く絞り込まれた周波数(あるいはその周辺のさらに細かい探索を経た値)にVCOをセットし、そこから先はPLLの回で確立した線形追尾モデルにバトンタッチします。この「引き渡し」がうまくいくためには、捕捉フェーズで達成する周波数精度が、PLLの線形動作範囲(ループ帯域幅 $B_L$ のオーダー)に収まっている必要があり、これが捕捉フェーズに要求される最低限の分解能の目安になります。

実務での使われ方

深宇宙探査機のミッションでは、周波数捕捉は打ち上げ直後の**初期取得(initial acquisition)**フェーズで特に重要な役割を果たします。打ち上げ直後は探査機の軌道がまだ精密に確定しておらず、DSN(Deep Space Network)の運用者が事前に用意する周波数予測モデル(打ち上げ計画と概略軌道から算出される予測ドップラープロファイル)の誤差が、巡航フェーズの精密軌道決定後と比べて大きくなります。この結果、受信すべき搬送波周波数の不確かさが数kHz〜数十kHzのオーダーに達することがあり、DSNの受信機(Block V Receiverおよびそのデジタル後継機)は、この不確かさ範囲にわたる周波数捕捉戦略を実行するように設計されています。具体的には、事前の周波数予測(打ち上げ軌道解析から得られる時刻同期したドップラープロファイル)によって不確かさ範囲をできる限り絞り込んだ上で、残る不確かさに対して掃引法やFFTベースのサーチを適用し、ピークを検出した時点でPLLの線形追尾へ切り替える、という段階的な戦略が一般的です。近年のデジタル受信機では、広帯域A/D変換後にデジタル信号処理でFFTベースの捕捉を行い、掃引法よりも短時間で搬送波の粗い位置を特定してからPLL(あるいはCostasループ)に引き渡す実装が主流になりつつあります。

同様の捕捉問題は、通信途絶からの回復時にも生じます。探査機がセーフモードに入ったり、姿勢異常でアンテナが地球を向かなくなったりした後の再捕捉では、直近の軌道情報自体は分かっていても、探査機側の発振器の温度ドリフトや姿勢変化に伴う予期しないドップラー変動が加わるため、やはりある程度の不確かさ範囲を捜索する必要があります。運用上、こうした緊急時のシナリオに備えて、捕捉戦略のパラメータ(掃引速度、探索範囲、FFTの積分時間)は、ミッションの信号強度(リンクバジェット)とダイナミクスの想定に応じてあらかじめ複数のプロファイルが用意され、状況に応じて切り替えられます。この設計思想はDSN Telecommunications Link Design Handbook (810-005) の受信機捕捉に関するモジュールに規定されています。

演習問題

周波数不確かさの半幅が $\Delta f_{\max} = 20$ kHz、掃引速度 $\dot f_{\text{sw}} = 500$ Hz/s のとき、不確かさ範囲全体( $2\Delta f_{\max}$ )を一往復するのにかかる時間 $T_{\text{acq}}$ を求めてください。
検出器の実効帯域幅が $B_i = 10$ Hz であるとき、問1の掃引速度 $\dot f_{\text{sw}} = 500$ Hz/s のもとで信号が検出器帯域を通過する時間(dwell time) $\tau_{\text{dwell}} \approx B_i/\dot f_{\text{sw}}$ を求めてください。この時間が短すぎると考えられる場合、掃引速度をどう調整すべきか、感度とのトレードオフを踏まえて議論してください。
FFTベースの捕捉で周波数分解能 $\delta f = 10$ Hz を達成したい場合、必要な観測時間 $T_{\text{obs}} = 1/\delta f$ を求めてください。またサンプリング周波数を $f_s = 1$ MHz とするとき、必要なFFT点数 $N = T_{\text{obs}}/T_s$ ( $T_s = 1/f_s$ )のオーダーを見積もってください。
なぜPLLの回で確立した線形PLLモデルは、そのままでは「周波数がまったく分からない」捕捉問題に使えないのか、位相検波器の特性 $v_d(t)=K_d\sin\phi(t)$ が非線形であることに触れながら、この回で学んだ捕捉フェーズとPLLの追尾フェーズの役割分担を自分の言葉で説明してください。

まとめと次回予告

PLLの回で扱った線形追尾モデルは、「位相誤差がすでに小さい」という前提の上に成り立っていましたが、探査機のドップラーシフトの不確かさが大きい局面(打ち上げ直後や再捕捉時など)では、その前提そのものを満たすために、まず広い周波数範囲から搬送波のおおよその位置を見つけ出す捕捉フェーズが必要でした。周波数掃引法は、局部発振器を一定速度で走査しながらエネルギーのピークを探す古典的な方式で、掃引速度と検出感度の間に明確なトレードオフを持ちます。FFTベースの捕捉は、短時間の観測ですべての周波数ビンを並列に評価できるため、原理的に掃引法より高速に周波数を絞り込めます。いずれの手法でも、不確かさ範囲が広いほど、また要求精度が高いほど捕捉には時間がかかるという一般的なトレードオフがあり、捕捉が完了した時点でPLLの線形追尾に制御を引き渡す、という2段構えの設計が深宇宙受信機の標準的な戦略になっています。

次回は、この「周波数領域での探索」というアイデアをさらに一般化し、信号がスパース(疎)である——つまり有用な情報がごく一部の成分に集中している——という性質を積極的に利用して、通常より少ない観測から信号を復元する圧縮センシングの考え方に軽く触れます。

参考文献

F. M. Gardner, Phaselock Techniques, 3rd ed., Wiley
A. J. Viterbi, Principles of Coherent Communication, McGraw-Hill
DSN Telecommunications Link Design Handbook, DSN No. 810-005(受信機捕捉・キャリアトラッキングに関するモジュール)
J. H. Yuen (ed.), Deep Space Telecommunications Systems Engineering, JPL Publication 82-76
J. G. Proakis, M. Salehi, Digital Communications, 5th ed., McGraw-Hill
CCSDS 401.0-B, Radio Frequency and Modulation Systems, Part 1: Earth Stations and Spacecraft