変調・符号化#28

シンボル同期 — Gardner検出器でシンボル境界のタイミングを追跡する

整合フィルタの出力は『正しい時刻でサンプリングすれば』2E/N₀を達成できる、が肝心のその時刻はどうやって知るのか。PLLが搬送波の位相を追尾したのに対し、今回はシンボル境界のタイミングを追尾するGardner検出器を導出し、なぜそれがキャリア位相に依存せず動作できるのかを数式で示す。

前提知識: matched-filter

シンボル同期タイミングリカバリGardner検出器早期-遅延サンプリング補間フィルタ

この回で学ぶこと

整合フィルタの回で、私たちは受信SNRを最大化する最適フィルタの形 h(t)=p(Tt)h(t)=p(T-t) を導出し、その出力を時刻 TT でサンプリングすれば、判定変数が

y(T)=0Tr(τ)p(τ)dτ=±Eb+ny(T) = \int_0^T r(\tau)\, p(\tau)\, d\tau = \pm\sqrt{E_b} + n

という理想的な信号空間モデルに帰着し、Pb=Q(2Eb/N0)P_b = Q(\sqrt{2E_b/N_0}) が達成できることを厳密に示しました。しかしこの議論には、あからさまに天下りな仮定が1つ隠れていました。「時刻 TT でサンプリングすれば」というくだりです。実際の受信機は、ある1つのシンボルがいつ始まりいつ終わるのか、あらかじめ知りません。探査機からの信号は伝搬遅延・ドップラーシフト・局部発振器の周波数誤差を経て届くため、受信機のサンプリングクロックとシンボル境界の間には、未知でしかも時間とともにゆっくり変動するずれが存在します。整合フィルタが理論上の性能限界を「達成可能」と示したことと、実際にその限界を「達成する」ことの間には、このタイミングのずれを推定し補正するという、もう1つの仕事が挟まっているのです。

この仕事を担うのがシンボル同期(symbol synchronization)、別名タイミングリカバリ(timing recovery)と呼ばれる受信機のサブシステムです。PLLの回で、私たちは地上局が探査機の残留搬送波の位相 θ(t)\theta(t) を、雑音の中からフィードバックループで追尾し続ける仕組みを学びました。今回のシンボル同期は、これと構造的によく似た問題です。ただし追いかける対象が搬送波の位相ではなく、シンボル境界のタイミング τ(t)\tau(t)(あるいは正規化したタイミングオフセット ε(t)=τ(t)/T\varepsilon(t) = \tau(t)/T)である点が違います。PLLが「位相」の同期ループだったのに対し、今回は「時間」の同期ループを作ります。この回では、そのタイミング誤差を検出する代表的な回路である Gardner検出器 (Gardner Timing Error Detector) を数式で導出し、それがなぜ「決定(判定結果)を使わずに」「搬送波の残留位相に影響されずに」動作できるのかを厳密に示し、最後にPLLとの数学的な類似性・相違点、そして実装上の要である補間フィルタの仕組みまでを見ていきます。

直感的な導入: サンプリング時刻がずれると何が起きるか

ISIとアイダイアグラムの回で、判定回路がサンプリングする時刻 t0t_0 が最適点からずれると、ピーク歪み Dpeak(t0)D_{\text{peak}}(t_0) が増大し、アイダイアグラムの「目」が横方向に閉じていくことを見ました。整合フィルタ理論が保証する 2Eb/N02E_b/N_0 という上限値は、あくまでパルスのピーク(整合フィルタ出力が最大になる瞬間)でサンプリングできたときにのみ達成されます。サンプリング時刻が最適点から τ\tau だけずれると、整合フィルタ出力はそのピークから離れた値しか取れず、実効的な信号振幅が小さくなり、隣接シンボルからの漏れ込み(ISI)も増えるため、実効SNRは劣化します。

つまり受信機は、整合フィルタを正しく実装しただけでは足りず、その出力のどの瞬間をサンプリングすべきかを自分で見つけ、しかも探査機との相対運動によって刻々と変化するそのタイミングを追い続けなければならないわけです。これはPLLが「見失った搬送波の位相を推定し追い続ける」問題と構造的に同じで、パラメータが位相からタイミングに変わっただけです。したがって解決策の骨格もPLLに似せて作ることができます。

  1. タイミング誤差検出器 (Timing Error Detector, TED): 受信サンプルから「今のサンプリングタイミングがどちらの方向にどれだけずれているか」を表す誤差信号を作る。
  2. ループフィルタ: その誤差信号を平滑化・積分する。
  3. 補間器/可変クロック: フィルタ後の値に応じて、実際にサンプリングする時刻(あるいはサンプル列から作り出す仮想的なサンプリング時刻)を調整する。

PLLの位相検波器・ループフィルタ・VCOという3点セットと見比べれば、対応関係は明らかです。以下ではまず、この対応関係の心臓部である「タイミング誤差検出器」をどう作るかを、具体的な数式で詰めていきます。

問題の定式化

整合フィルタ出力(あるいは同じことですが、まだ整合フィルタをかける前でもよい、後述)を、シンボルレート 1/T1/T に対して十分細かくサンプリングした連続時間波形 y(t)y(t) とします。伝送された(未知の)シンボル列を {ak}\{a_k\}、combined pulse(送信パルス整形と受信フィルタを合成した波形、ISIの回で使ったのと同じ記法)を g(t)g(t)、未知のタイミングオフセットを τ\tau とすると、

y(t)=m=amg(tmTτ)+n(t)y(t) = \sum_{m=-\infty}^{\infty} a_m\, g(t - mT - \tau) + n(t)

正規化したタイミングオフセット ετ/T\varepsilon \equiv \tau/T を使うと、受信機の仕事は「ε\varepsilon を雑音の中から推定し、g(tmTτ)g(t-mT-\tau) が最大値(ピーク)を取る瞬間、すなわち t=mT+τt = mT+\tau を狙ってサンプリングすること」だと言い換えられます。

ここで整合フィルタ理論からの重要な前提を思い出しておきます。整合フィルタの回で証明した通り、送信パルスの自己相関関数がそのまま整合フィルタの出力波形になります。NRZ矩形パルス p(t)=Ap(t)=A(0tT0\le t\le T)の自己相関は、よく知られているように三角形になります。以下ではこの具体的なケースを使い、Gardner検出器の誤差信号を最後まで厳密に計算します。合成パルスをピーク位置を原点に取り直して

g(t)={1tTt<T0tTg(t) = \begin{cases} 1 - \dfrac{|t|}{T} & |t| < T \\ 0 & |t| \ge T \end{cases}

と正規化します。この g(t)g(t) は隣接シンボルとしか重ならない(2シンボル以上離れると完全にゼロ)という扱いやすい性質を持っており、以下の計算をきれいに閉じた形にしてくれます。

Gardner検出器の導出

Gardner検出器の特徴は、1シンボルあたり2つのサンプル(シンボル境界上の「オンタイムサンプル」と、その中間点の「中間サンプル」)だけを使い、しかも判定結果(推定シンボル a^k\hat a_k)を一切使わずに誤差信号を作れる点にあります。受信機のサンプリングクロックが kTkT にあると仮定して(実際のずれ ε\varepsilon はこれから推定する)、次の2種類のサンプルを定義します。

xky(kT+τ),xk12y ⁣(kTT2+τ)x_k \equiv y(kT + \tau), \qquad x_{k-\frac12} \equiv y\!\left(kT - \tfrac{T}{2} + \tau\right)

Gardner検出器の誤差信号は、次の単純な積として定義されます。

ek    xk12(xkxk1)e_k \;\equiv\; x_{k-\frac12}\,\big(x_k - x_{k-1}\big)

直感的には、xkx_kxk1x_{k-1}(隣接する2つのオンタイムサンプル)の差分が、その間にある中間サンプル xk12x_{k-\frac12} でどれだけ重み付けされるかを見ています。もしサンプリングタイミングが完全に正しければ(ε=0\varepsilon=0)、中間サンプルはちょうど2つのシンボル遷移の中点にあり、対称性から誤差信号の期待値はゼロになるはずです。タイミングがずれると、この対称性が崩れて誤差信号にゼロでない値が現れる——というのが直感的な仕組みです。これを厳密に確認しましょう。

三角パルスによる厳密な計算

g(t)g(t)t<T|t|<T でしか非ゼロにならないことを使うと、0ε<120 \le \varepsilon < \tfrac12 の範囲で(タイミングオフセットが小さく、隣接する高々2シンボルしか重ならない状況を考えます)、次の値が具体的に計算できます。

xk=ak(1ε)+ak1ε,xk1=ak1(1ε)+ak2εx_k = a_k(1-\varepsilon) + a_{k-1}\,\varepsilon, \qquad x_{k-1} = a_{k-1}(1-\varepsilon) + a_{k-2}\,\varepsilon xk12=ak ⁣(12ε)+ak1 ⁣(12+ε)x_{k-\frac12} = a_k\!\left(\tfrac12-\varepsilon\right) + a_{k-1}\!\left(\tfrac12+\varepsilon\right)

(それぞれ、g(εT)=1εg(-\varepsilon T)=1-\varepsilong(TεT)=εg(T-\varepsilon T)=\varepsilong(T/2εT)=12εg(-T/2-\varepsilon T)=\tfrac12-\varepsilong(T/2εT)=12+εg(T/2-\varepsilon T)=\tfrac12+\varepsilon を代入し、支持区間 t<T|t|<T の外側にあたる項をゼロとして落としたものです。)

これらを ek=xk12(xkxk1)e_k = x_{k-\frac12}(x_k-x_{k-1}) に代入し展開すると、

xkxk1=ak(1ε)+ak1(2ε1)ak2εx_k - x_{k-1} = a_k(1-\varepsilon) + a_{k-1}(2\varepsilon-1) - a_{k-2}\,\varepsilon ek=[ak(12ε)+ak1(12+ε)][ak(1ε)+ak1(2ε1)ak2ε]e_k = \Big[a_k\big(\tfrac12-\varepsilon\big)+a_{k-1}\big(\tfrac12+\varepsilon\big)\Big]\Big[a_k(1-\varepsilon) + a_{k-1}(2\varepsilon-1) - a_{k-2}\varepsilon\Big]

シンボル {ak}\{a_k\} が独立同分布で ak{1,+1}a_k \in \{-1,+1\}E[aiaj]=δijE[a_i a_j] = \delta_{ij}(異なる時刻のシンボルは無相関)であることを使って期待値を取ると、異なる添字の積(akak1a_k a_{k-1}akak2a_k a_{k-2}ak1ak2a_{k-1}a_{k-2})の期待値はすべてゼロになり、ak2=ak12=1a_k^2=a_{k-1}^2=1 の項だけが残ります。

E[ek]=(12ε)(1ε)+(12+ε)(2ε1)E[e_k] = \left(\tfrac12-\varepsilon\right)(1-\varepsilon) + \left(\tfrac12+\varepsilon\right)(2\varepsilon-1)

これを丁寧に展開すると、

E[ek]=(1232ε+ε2)+(2ε212)=32ε+3ε2=32ε(12ε)E[e_k] = \left(\tfrac12 - \tfrac32\varepsilon + \varepsilon^2\right) + \left(2\varepsilon^2 - \tfrac12\right) = -\tfrac32\,\varepsilon + 3\varepsilon^2 = -\frac{3}{2}\,\varepsilon\,(1-2\varepsilon)

これがGardner検出器のS字カーブです。 ε=0\varepsilon=0 でちょうどゼロを通り、ε\varepsilon が小さいうちは

E[ek]32εE[e_k] \approx -\frac{3}{2}\,\varepsilon

と、タイミングオフセットに比例した(符号付きの)誤差信号が得られます。ループフィルタとNCO/補間器がこの誤差信号を使って ε0\varepsilon \to 0 の方向にサンプリングタイミングを修正し続ければ、ループはタイミングオフセットゼロの状態に収束します。PLLの回の位相検波器の特性 vd(t)=KdsinϕKdϕv_d(t)=K_d\sin\phi \approx K_d\phi(小さい ϕ\phi で線形近似)と、まったく同じ構造の「線形化されたS字カーブ」が、ここでもタイミング誤差について再現されていることに注目してください。

なお興味深いことに、ε=0\varepsilon=0 ちょうどのときは、確率的な期待値としてではなく、どんなデータパターンでも厳密に ek=0e_k=0 になります(実際に代入すると xk12(xkxk1)=12(ak+ak1)(akak1)=12(ak2ak12)=0x_{k-\frac12}(x_k-x_{k-1}) = \tfrac12(a_k+a_{k-1})(a_k-a_{k-1}) = \tfrac12(a_k^2-a_{k-1}^2)=0ak2=ak12=1a_k^2=a_{k-1}^2=1 より恒等的にゼロ)。しかし ε0\varepsilon \neq 0 では、eke_k の値そのものはデータパターン(ak,ak1,ak2a_k,a_{k-1},a_{k-2} の組み合わせ)によってばらつきます。この期待値まわりのばらつきは、加法性雑音とは別に生じるセルフノイズ (self-noise) と呼ばれ、Gardner検出器のような「判定結果を使わない (non-data-aided, NDA)」タイプの誤差検出器に特有の、実務上無視できない誤差源になります(後述)。

なぜこの検出器は「非データ判定型」かつ「キャリア位相非依存」で動作できるのか

Gardner検出器の重要な利点は2つあり、両方とも上の導出の構造そのものから直接見て取れます。

(1) 非データ判定型(non-data-aided)であること。 上の E[ek]E[e_k] の計算では、シンボル {ak}\{a_k\}実際の値が何であったかを一切使っていません。使ったのは「aka_k が独立同分布で E[aiaj]=δijE[a_ia_j]=\delta_{ij}」という統計的な性質だけです。これは、受信機が誤り訂正やシンボル判定を行う前の、生のサンプル xk,xk12x_k, x_{k-\frac12} だけから誤差信号を計算できることを意味します。対照的に、いわゆる決定指向型 (decision-directed) の検出器(たとえばMueller-Müller検出器)は、暫定的な判定結果 a^k\hat a_k を使って誤差を計算するため、判定回路より後段に位置する必要があり、判定誤りがあると誤差検出そのものが乱れるという弱点を持ちます。Gardner検出器はこの依存関係を持たないため、受信チェーンの早い段階、極端には判定回路や誤り訂正復号より前でタイミングを確立できます。

(2) キャリア位相に非依存であること。 ここまでは簡単のため実数(ベースバンドで位相回転のない)信号として扱ってきましたが、実際のQPSK/BPSK受信機では、Costasループの回で見たように、搬送波位相の再生が完全に収束しきる前の段階でも、あるいは残留した位相誤差 ϕ0\phi_0 が乗ったままの状態でも、シンボル同期を並行して走らせたいことがよくあります。複素ベースバンド表現で、受信サンプルに未知の一定位相回転 ejϕ0e^{j\phi_0} が乗っている状況を考えます。

y(t)=ejϕ0mamg(tmTτ)+n(t)y(t) = e^{j\phi_0}\sum_m a_m\, g(t-mT-\tau) + n(t)

このとき、Gardner検出器を複素数に自然拡張した形

ek  =  Re ⁣{xk12(xkxk1)}e_k \;=\; \operatorname{Re}\!\left\{\, x_{k-\frac12}^{\,*}\,\big(x_k - x_{k-1}\big) \,\right\}

を使います(実数の場合の xk12(xkxk1)x_{k-\frac12}(x_k-x_{k-1}) を、複素共役を使った内積の形に一般化したものです)。ここで、すべてのサンプルに共通の位相回転 ejϕ0e^{j\phi_0} が乗っているとすると、xkejϕ0xkx_k \to e^{j\phi_0}x_kxk1ejϕ0xk1x_{k-1}\to e^{j\phi_0}x_{k-1}xk12ejϕ0xk12x_{k-\frac12}\to e^{j\phi_0}x_{k-\frac12} となるので、

xk12(ejϕ0xk12)=ejϕ0xk12x_{k-\frac12}^{\,*} \to \big(e^{j\phi_0}x_{k-\frac12}\big)^{*} = e^{-j\phi_0}\,x_{k-\frac12}^{\,*}

したがって、

xk12(xkxk1)    ejϕ0xk12ejϕ0(xkxk1)=ejϕ0ejϕ0xk12(xkxk1)x_{k-\frac12}^{\,*}\big(x_k-x_{k-1}\big) \;\longrightarrow\; e^{-j\phi_0}\,x_{k-\frac12}^{\,*} \cdot e^{j\phi_0}\big(x_k-x_{k-1}\big) = e^{-j\phi_0}e^{j\phi_0}\, x_{k-\frac12}^{\,*}\big(x_k-x_{k-1}\big)

ejϕ0ejϕ0=1e^{-j\phi_0}e^{j\phi_0}=1 ですから、未知の位相回転 ϕ0\phi_0 は完全に打ち消し合い、誤差信号 eke_k にはまったく現れません。 これは複素共役を挟んだ積が本質的に位相回転に対して不変な演算になっている(内積の絶対値・実部が回転で保存される)という一般的な事実の直接の帰結です。この性質のおかげで、Gardner検出器は搬送波再生ループ(Costasループなど)がまだロックし切っていない、あるいはそもそも搬送波再生と独立に動かしたい局面でも、正しくタイミング誤差を検出できます。実務上、多くのディジタル受信機ではシンボル同期ループを搬送波再生ループと並行に(あるいはタイミング同期を先に確立させてから搬送波再生に入る順序で)動かすことができ、この柔軟性の根拠がまさにこの位相不変性です。

ループ全体の構造とPLLとの対比

以上を踏まえ、シンボル同期ループの全体をブロック図として整理すると、PLLの回で見た構造と1対1に対応します。

要素PLL(搬送波位相追尾)シンボル同期ループ(タイミング追尾)
追尾する量搬送波の位相 θ(t)\theta(t)シンボルタイミング τ(t)\tau(t)(正規化 ε(t)\varepsilon(t))
誤差検出器位相検波器: vdKdϕv_d\approx K_d\,\phiGardner検出器(TED): E[ek]KTEDεE[e_k]\approx K_{TED}\,\varepsilon
ループフィルタPIフィルタ F(s)=(1+sτ2)/(sτ1)F(s)=(1+s\tau_2)/(s\tau_1)同様のPIフィルタ(離散時間版)
制御対象VCO(電圧制御発振器、dθo/dt=Kovcd\theta_o/dt = K_o v_c)補間器 + NCO(数値制御発振器、可変サンプリング位相)
ループの積分器数2個(VCO 1個 + フィルタ 1個)→ Type-II2個(NCO 1個 + フィルタ 1個)→ 同じくType-II構造が標準

誤差検出器がロック近傍で線形近似できる(E[ek]KTEDεE[e_k]\approx K_{TED}\varepsilon)という点さえ認めてしまえば、PLLの回で導いた開ループ伝達関数 G(s)=KF(s)/sG(s)=K F(s)/s、閉ループ伝達関数 H(s)H(s)、そしてループ雑音帯域幅

BL=ωn2(ζ+14ζ)B_L = \frac{\omega_n}{2}\left(\zeta+\frac{1}{4\zeta}\right)

の枠組みが(ゲイン KdK_dKTEDK_{TED} に読み替えるだけで)そのまま流用できます。Costasループの回でも同じ手法を使い、Costas検出器の等価線形ゲイン Kd=A2/4K_d=A^2/4 をPLLの式にそのまま代入しました。シンボル同期ループも本質的には**「制御対象がVCOではなく、離散時間で動くタイミング補間器である」PLLの一種**とみなせるわけです。

ただし相違点も明確にしておく必要があります。PLLのVCOは連続時間の物理的な発振器で、電圧に応じて瞬時角周波数を直接偏移させます。これに対しシンボル同期ループの「制御対象」は、実装によって大きく2通りに分かれます。

方式1: アナログ的解法(可変サンプリングクロック)。 VCO(あるいはVCXO, 電圧制御水晶発振器)そのものでA/D変換器のサンプリングクロックの位相を直接調整する方式です。PLLとほぼ同じ回路構成がそのまま流用できますが、A/D変換器のクロックを動的に動かす必要があり、現代の高集積度のディジタル受信機ではあまり好まれません。

方式2: デジタル的解法(補間器 + NCO)。 現代のSDR(Software-Defined Radio)受信機で標準的なのがこちらです。A/D変換器は探査機のシンボルレートとは無関係な、固定の(かつ十分高い)クロックレートで動作し続けます。受信機はこの固定レートのサンプル列 x[n]x[n] から、補間フィルタ (interpolator) を使って、任意の(整数サンプル間隔の狭間にある)非整数時刻のサンプル値を計算で作り出します。

補間器の仕組みの概要

固定レートのサンプル列 x[n]x[n](nn は整数、サンプル周期 TsT_s)から、非整数時刻 t=(n+μ)Tst=(n+\mu)T_s(μ[0,1)\mu \in [0,1) が求めたい小数遅延)におけるサンプル値 x((n+μ)Ts)x\big((n+\mu)T_s\big) を求める問題は、ディジタル信号処理では帯域制限された信号の補間として一般的に扱われます。理論的にはサンプリング定理に基づく sinc 補間が厳密解ですが、無限のタップ数が必要で実装できないため、実際には低次の多項式補間(線形補間、あるいは3次多項式によるキュービック補間)で近似します。

実装上よく使われるのが Farrow構造 と呼ばれる効率的な形です。これは、補間フィルタの係数を μ\mu の多項式として

x((n+μ)Ts)i=0Nμi(kci,kx[nk])x\big((n+\mu)T_s\big) \approx \sum_{i=0}^{N} \mu^{i}\left(\sum_k c_{i,k}\, x[n-k]\right)

という形に分解し、内側の(固定係数の)FIRフィルタ群をあらかじめ設計しておけば、μ\mu が変わるたびにフィルタ係数全体を再計算する必要がなく、μ\mu のべき乗を外側でかけ合わせるだけで済む、という構造です。ループフィルタの出力は、この μ\mu(小数タイミング、シンボル周期に対する位相)を刻々と更新する信号として使われます。

このループの中で、PLLのVCOに相当する役割を果たすのがNCO(Numerically Controlled Oscillator, 数値制御発振器)、より正確には「モジュロ1のカウンタ(アキュムレータ)」です。固定サンプルレートのクロックごとに、ループフィルタの出力(推定されたシンボルレートのわずかなずれに相当する量)をアキュムレータに積算していき、アキュムレータが1を超えて桁上がりする(オーバーフローする)瞬間が「次のシンボル境界に到達した」ことを意味します。オーバーフロー直前のアキュムレータの値がその瞬間の小数タイミング μ\mu を与え、それを使って補間器がその時刻のサンプル値を計算し、これがシンボルレートで出力される「タイミング補正済みサンプル」になります。この一連の流れ(誤差検出 → ループフィルタ → NCO(アキュムレータ) → 補間器)が、PLLの(位相検波器 → ループフィルタ → VCO)の構造を、連続時間の物理発振器を使わずに完全にディジタル演算だけで模倣したものになっているわけです。

実務での使われ方

現代のディジタル受信機・SDR(GNU Radioなどのオープンソースツールチェーンを含む)では、上記の「固定レートADC + Gardner検出器 + Farrow型補間器 + NCO」という構成が、シンボル同期を実装する際の事実上の標準アーキテクチャになっています。典型的な実装では、A/D変換器はシンボルレートの2倍(2サンプル/シンボル)以上でオーバーサンプリングし、Gardner検出器はその中間サンプルを使って誤差を計算し、ループフィルタとNCOがシンボルレートに同期したサンプリング位相を維持し続けます。

DSN(Deep Space Network)や商用衛星地上局のベースバンドプロセッサでは、シンボル同期ブロックはCostasループなどの搬送波再生ブロックと、Viterbi/LDPC復号などの誤り訂正ブロックとの間、整合フィルタの直後に位置づけられ、CCSDS 401.0-B や DSN Telecommunications Link Design Handbook (810-005) にも、この受信チェーン全体の構成が規定されています。

タイミングジッタとBERの関係については、ISIとアイダイアグラムの回で導出したピーク歪みの枠組みがそのまま使えます。シンボル同期ループが残す残留タイミング誤差の標準偏差を σε\sigma_\varepsilon(シンボル周期 TT に対する正規化値)とすると、これはアイダイアグラムの横方向の開口をランダムに侵食する量に対応し、実効的なピーク歪み Dpeak(ε)D_{\text{peak}}(\varepsilon) を増大させ、結果としてビット誤り率を劣化させます。実務上の経験則としては、ループがある程度整定した後の残留ジッタ σε\sigma_\varepsilon を、シンボル周期のおおむね数パーセント以下(たとえば σε/T0.02\sigma_\varepsilon/T \lesssim 0.020.050.05程度)に抑えることが、実用上無視できるBER劣化(1dB未満)の目安として扱われます。この設計では、PLLと同様にループ雑音帯域幅 BLB_L を狭くするほど雑音の取り込みが減り σε\sigma_\varepsilon が下がる一方、探査機側のシンボルレート変動(ドップラーに伴うシンボルレートの伸縮)への追従が遅れるというトレードオフが働き、BLTB_L T(ループ帯域幅とシンボル周期の積)の値がミッションのダイナミクスと信号強度に応じて選定されます。

また前節で触れたセルフノイズは、Gardner検出器のような非データ判定型のTEDに共通する特性で、雑音のない理想的なチャネルでも、データパターンのランダム性そのものに起因する残留ジッタのフロア(下限)を生み出します。これは決定指向型の検出器(受信機が高SNRで安定に動作している局面でより低いセルフノイズを持つことが知られています)との性能上のトレードオフの一因であり、Mengali & D’Andreaなどの標準的な文献で詳しく解析されています。実際のミッション設計では、想定される動作SNRレンジ全体でセルフノイズを含めたジッタ性能をシミュレーションで評価し、Gardner検出器を採用するか、決定指向型や他のNDA検出器(たとえばMueller-Müller検出器)を採用するかが判断されます。

演習問題

  1. 本文では 0ε<120 \le \varepsilon < \tfrac12 の場合について Gardner 検出器の期待値 E[ek]=32ε(12ε)E[e_k] = -\tfrac32\varepsilon(1-2\varepsilon) を導出しました。同様の手順を 12<ε<0-\tfrac12 < \varepsilon < 0 の場合について行い(対称性を使ってもよいですし、g(t)g(t) への代入からやり直してもかまいません)、ε=0\varepsilon=0 の前後で E[ek]E[e_k] が符号を変える(タイミングのずれの方向を正しく検出できる)ことを確認してください。

  2. 複素ベースバンドのGardner検出器 ek=Re{xk12(xkxk1)}e_k = \operatorname{Re}\{x_{k-\frac12}^{*}(x_k-x_{k-1})\} が、共通の位相回転 ejϕ0e^{j\phi_0} に対して不変であることを、本文の導出をなぞって自分の手で再現してください。また、もし複素共役を取らずに ek=xk12(xkxk1)e_k = x_{k-\frac12}(x_k-x_{k-1})(実数の式をそのまま複素数に流用したもの)と誤って定義してしまった場合、位相回転に対してどのような不都合が生じるかを議論してください。

  3. あるシンボル同期ループが ζ=0.707\zeta=0.707、正規化ループ帯域幅 BLT=0.01B_LT = 0.01 で設計されているとします。PLLの回の公式 BL=ωn2(ζ+14ζ)B_L=\dfrac{\omega_n}{2}\left(\zeta+\dfrac{1}{4\zeta}\right) を踏まえ、ωnT\omega_n T(シンボル周期で正規化した自然角周波数)の値を求めてください。またこの値が意味する「ループが整定するのに要するおおよそのシンボル数」を、1/ωn1/\omega_n のオーダーで概算してください。

  4. Gardner検出器(非データ判定型、1シンボルあたり2サンプル必要)と、決定指向型の検出器(判定結果 a^k\hat a_k を使い、原理的には1シンボルあたり1サンプルで動作可能)を比較し、それぞれの利点・欠点を、本文で述べた「非データ判定型であること」「セルフノイズ」「判定回路との依存関係」というキーワードを使って自分の言葉で整理してください。

まとめと次回予告

シンボル同期は、整合フィルタの回で証明した「理論上達成可能な最大SNR」を、実際の受信機で本当に達成するために欠かせない、もう1つのフィードバック制御ループです。PLLが搬送波の位相を追尾するループだったのに対し、今回のGardner検出器を使ったループはシンボルタイミングを追尾するループであり、誤差検出器・ループフィルタ・制御対象(VCOではなく補間器+NCO)という同じ3点セットの構造、そして同じType-II線形ループの解析フレームワークを共有していることを見ました。Gardner検出器そのものは、非データ判定型でありながら三角パルスの具体例で厳密にS字カーブ E[ek]32εE[e_k]\approx-\tfrac32\varepsilon を導出でき、しかも複素共役を使った構成によって未知の搬送波位相回転に対して数学的に不変であることも確認しました。

これで、PCM/PSK/PMによる変調から始まり、PLLCostasループによる搬送波再生、整合フィルタによるSNR最大化、そして今回のシンボル同期まで、地上局が探査機からの微弱な信号を「見つけ、位相を合わせ、タイミングを合わせ、最適に判定する」という受信チェーンの主要な要素をひと通り数式付きで揃えたことになります。次回以降は、測距・追跡カテゴリに戻り、複数の地上局アンテナを使って探査機の到来方向を精密に測る同一ビーム干渉法のような技術に話を進めていきます。

参考文献

  • F. M. Gardner, “A BPSK/QPSK Timing-Error Detector for Sampled Receivers,” IEEE Transactions on Communications, vol. 34, no. 5, 1986
  • U. Mengali, A. N. D’Andrea, Synchronization Techniques for Digital Receivers, Plenum Press
  • H. Meyr, M. Moeneclaey, S. A. Fechtel, Digital Communication Receivers: Synchronization, Channel Estimation, and Signal Processing, Wiley
  • 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
  • DSN Telecommunications Link Design Handbook, DSN No. 810-005