位相同期ループベースの復調回路
コスタス ループは 位相同期回路 (PLL) に基づく回路 で、 抑圧搬送波 変調 信号 (例: 両側 波帯 抑圧搬送波信号) および位相変調信号 (例: BPSK 、 QPSK ) からの 搬送波 周波数 回復 に使用されます。コスタスループは1950 年代に ゼネラルエレクトリック の ジョン P. コスタス によって発明されました。 [1] [2]この発明は [3] 「現代のデジタル通信に大きな影響を与えた」と評されました 。コスタスループの主な用途は無線受信機です。他の PLL ベースの検出器に対するコスタスループの利点は、小さな偏差ではコスタスループの誤差電圧が と比較して であることです 。これは感度が 2 倍になることを意味し、またコスタスループは特に OFDM および GPS 受信機 における ドップラーシフト 搬送波の追跡に特に適しています。 [3]
罪
(
2
(
θ
私
−
θ
ふ
)
)
{\displaystyle \sin(2(\theta _{i}-\theta _{f}))}
罪
(
θ
私
−
θ
ふ
)
{\displaystyle \sin(\theta _{i}-\theta _{f})}
古典的な実装
ロック状態で動作するコスタスループ。
コスタスループの古典的な実装では、 [4] ローカル 電圧制御発振器 (VCO) が 直交出力を 2 つの 位相検出器 ( たとえば 積 検出器) にそれぞれ 1 つずつ提供します。入力 信号 と同じ位相が 両方の位相検出器にも適用され、各 位相検出器の出力は ローパスフィルタ を通過します 。これらのローパスフィルタの出力は別の位相検出器への入力となり、その出力は電圧制御発振器の制御に使用される前にノイズ低減フィルタを通過します。全体的なループ応答は 3 番目の位相検出器の前にある 2 つの個別のローパスフィルタによって制御されますが、3 番目のローパスフィルタはゲインと位相マージンに関して些細な役割を果たします。
上のコスタス ループの図は、コスタス ループ プロセスによって VCO 周波数と入力キャリア周波数が同じになった「ロック」状態で描かれています。この図は「ロック解除」状態を表すものではありません。
数学モデル
時間領域では
BPSKコスタスループの時間領域モデル
最も単純なケースでは 、 ノイズ低減フィルタの入力には影響しません。搬送波と 電圧制御発振器 (VCO) 信号は、 高周波数の周期的な振動です 。ブロックは アナログ乗算器 です 。
メートル
2
(
t
)
=
1
{\displaystyle m^{2}(t)=1}
メートル
2
(
t
)
=
1
{\displaystyle m^{2}(t)=1}
ふ
r
e
ふ
、
ヴ
c
o
(
θ
r
e
ふ
、
ヴ
c
o
(
t
)
)
{\displaystyle f_{ref,vco}(\theta _{ref,vco}(t))}
θ
˙
r
e
ふ
、
ヴ
c
o
(
t
)
{\displaystyle {\dot {\theta }}_{ref,vco}(t)}
⨂
{\displaystyle \bigotimes}
線形フィルタは 、 線形微分方程式のシステムによって数学的に記述できます。
x
˙
=
あ
x
+
b
あなた
d
(
t
)
、
あなた
ら
ふ
=
c
∗
x
。
{\displaystyle {\begin{array}{ll}{\dot {x}}=Ax+bu_{d}(t),&u_{LF}=c^{*}x.\end{array}}}
ここで、 は定数行列、 はフィルタの状態ベクトル、 は 定数ベクトルです。
あ
{\displaystyle A}
x
(
t
)
{\displaystyle x(t)}
b
{\displaystyle b}
c
{\displaystyle c}
VCO のモデルは通常、線形であると想定されます。
θ
˙
ヴ
c
o
(
t
)
=
ω
ヴ
c
o
ふ
r
e
e
+
け
ヴ
c
o
あなた
ら
ふ
(
t
)
、
t
∈
[
0
、
T
]
、
{\displaystyle {\begin{array}{ll}{\dot {\theta }}_{vco}(t)=\omega _{vco}^{free}+K_{vco}u_{LF}(t),&t\in [0,T],\end{array}}}
ここで、 は VCO の自走周波数、は VCO ゲイン係数です。同様に、VCO のさまざまな非線形モデルを検討することも可能です。
ω
ヴ
c
o
ふ
r
e
e
{\displaystyle \omega _{vco}^{free}}
け
ヴ
c
o
{\displaystyle K_{vco}}
マスタージェネレータの周波数が一定であると仮定する
VCOの方程式とフィルタ収量の方程式
θ
˙
r
e
ふ
(
t
)
≡
ω
r
e
ふ
。
{\displaystyle {\dot {\theta }}_{ref}(t)\equiv \omega _{ref}.}
x
˙
=
あ
x
+
b
ふ
r
e
ふ
(
θ
r
e
ふ
(
t
)
)
ふ
ヴ
c
o
(
θ
ヴ
c
o
(
t
)
)
、
θ
˙
ヴ
c
o
=
ω
ヴ
c
o
ふ
r
e
e
+
け
ヴ
c
o
c
∗
x
。
{\displaystyle {\begin{array}{ll}{\dot {x}}=Ax+bf_{ref}(\theta _{ref}(t))f_{vco}(\theta _{vco}(t)),&{\dot {\theta }}_{vco}=\omega _{vco}^{free}+K_{vco}c^{*}x.\end{array}}}
このシステムは 非自律的 であり、調査するのはかなり困難です。
位相周波数領域では
コスタスループの等価位相周波数領域モデル
コスタスループの位相周波数領域モデルのVCO入力
最も単純なケースでは、
ふ
r
e
ふ
(
θ
r
e
ふ
(
t
)
)
=
コス
(
ω
r
e
ふ
t
)
、
ふ
ヴ
c
o
(
θ
ヴ
c
o
(
t
)
)
=
罪
(
θ
ヴ
c
o
(
t
)
)
ふ
r
e
ふ
(
θ
r
e
ふ
(
t
)
)
2
ふ
ヴ
c
o
(
θ
ヴ
c
o
(
t
)
)
ふ
ヴ
c
o
(
θ
ヴ
c
o
(
t
)
−
π
2
)
=
−
1
8
(
2
罪
(
2
θ
ヴ
c
o
(
t
)
)
+
罪
(
2
θ
ヴ
c
o
(
t
)
−
2
ω
r
e
ふ
t
)
+
罪
(
2
θ
ヴ
c
o
(
t
)
+
2
ω
r
e
ふ
t
)
)
{\displaystyle {\begin{aligned}f_{ref}{\big (}\theta _{ref}(t){\big )}=\cos {\big (}\omega _{ref}t{\big )},\ f_{vco}{\big (}\theta _{vco}(t){\big )}&=\sin {\big (}\theta _{vco}(t){\big )}\\f_{ref}{\big (}\theta _{ref}(t){\big )}^{2}f_{vco}\left(\theta _{vco}(t)\right)f_{vco}\left(\theta _{vco}(t)-{\frac {\pi }{2}}\right)&=-{\frac {1}{8}}{\Big (}2\sin(2\theta _{vco}(t))+\sin(2\theta _{vco}(t)-2\omega _{ref}t)+\sin(2\theta _{vco}(t)+2\omega _{ref}t){\Big )}\end{aligned}}}
標準的なエンジニアリングの仮定は、フィルタが入力から上側サイドバンド周波数を除去し、下側サイドバンドは変更しないというものです。したがって、VCO入力は、入力信号とVCO信号の 特定の波形 と 特性に対応する位相検出器を備え た位相ロックループ と 同等のコスタスループになります 。時間領域と位相周波数領域でのフィルタ出力はほぼ等しいことが証明されています。 [5] [6] [7]
φ
(
θ
r
e
ふ
(
t
)
−
θ
ヴ
c
o
(
t
)
)
=
1
8
罪
(
2
ω
r
e
ふ
t
−
2
θ
ヴ
c
o
(
t
)
)
。
{\displaystyle \varphi (\theta _{ref}(t)-\theta _{vco}(t))={\frac {1}{8}}\sin(2\omega _{ref}t-2\theta _{vco}(t)).}
φ
(
θ
)
{\displaystyle \varphi (\theta )}
ふ
r
e
ふ
(
θ
)
{\displaystyle f_{ref}(\theta )}
ふ
ヴ
c
o
(
θ
)
{\displaystyle f_{vco}(\theta )}
したがって、より単純な 微分方程式の
自律系 を研究することが可能である [8]。
x
˙
=
あ
x
+
b
φ
(
Δ
θ
)
、
Δ
θ
˙
=
ω
ヴ
c
o
ふ
r
e
e
−
ω
r
e
ふ
+
け
ヴ
c
o
c
∗
x
、
Δ
θ
=
θ
ヴ
c
o
−
θ
r
e
ふ
。
{\displaystyle {\begin{aligned}{\dot {x}}&=Ax+b\varphi (\Delta \theta ),\\\Delta {\dot {\theta }}&=\omega _{vco}^{free}-\omega _{ref}+K_{vco}c^{*}x,\\\Delta \theta &=\theta _{vco}-\theta _{ref}.\end{aligned}}}
。
クリロフ ・ボゴリュボフ平均化法 では、非自律方程式と自律方程式の解がいくつかの仮定の下で近いことを証明できます。したがって、時間領域のコスタスループブロック図は、位相周波数関係のレベルのブロック図に漸近的に変更できます。
コスタスループの自律動的モデル(非自律モデルの代わりに)の解析への移行により、入力信号の非常に速い時間スケールと信号位相の遅い時間スケールを同時に観察する必要がある、時間領域でのコスタスループのモデリングに関連する困難を克服できます。このアイデアにより、コアパフォーマンス特性( ホールドイン、プルイン、ロックイン範囲) を計算することが可能になります [9] 。
周波数取得
古典的なコスタスループは、搬送波とVCO間の位相差を小さく、理想的にはゼロの値にするように機能します。 [10] [11] [12] 位相差が小さいということは、周波数ロックが達成されたことを意味します。
QPSK コスタスループ
古典的なコスタスループは、より高いデータレートのために QPSK 変調に適応させることができる。 [13]
古典的なQPSKコスタスループ
入力 QPSK 信号は次の通りである
メートル
1
(
t
)
コス
(
ω
参照
t
)
+
メートル
2
(
t
)
罪
(
ω
参照
t
)
、
メートル
1
(
t
)
=
±
1
、
メートル
2
(
t
)
=
±
1.
{\displaystyle m_{1}(t)\cos \left(\omega _{\text{ref}}t\right)+m_{2}(t)\sin \left(\omega _{\text{ref}}t\right),m_{1}(t)=\pm 1,m_{2}(t)=\pm 1.}
ローパスフィルタLPF1とLPF2の入力は
φ
1
(
t
)
=
コス
(
θ
VC
)
(
メートル
1
(
t
)
コス
(
ω
参照
t
)
+
メートル
2
(
t
)
罪
(
ω
参照
t
)
)
、
φ
2
(
t
)
=
罪
(
θ
VC
)
(
メートル
1
(
t
)
コス
(
ω
参照
t
)
+
メートル
2
(
t
)
罪
(
ω
参照
t
)
)
。
{\displaystyle {\begin{aligned}\varphi _{1}(t)&=\cos \left(\theta _{\text{vco}}\right)\left(m_{1}(t)\cos \left(\omega _{\text{ref}}t\right)+m_{2}(t)\sin \left(\omega _{\text{ref}}t\right)\right),\\\varphi _{2}(t)&=\sin \left(\theta _{\text{vco}}\right)\left(m_{1}(t)\cos \left(\omega _{\text{ref}}t\right)+m_{2}(t)\sin \left(\omega _{\text{ref}}t\right)\right).\end{aligned}}}
同期後、LPF1と LPF2 の出力 を使用して復調データ ( および ) を取得します。VCO の周波数を基準周波数に調整するために、信号 とが 制限され、相互乗算されます。
質問
(
t
)
{\displaystyle Q(t)}
私
(
t
)
{\displaystyle I(t)}
メートル
1
(
t
)
{\displaystyle m_{1}(t)}
メートル
2
(
t
)
{\displaystyle m_{2}(t)}
質問
(
t
)
{\displaystyle Q(t)}
私
(
t
)
{\displaystyle I(t)}
あなた
d
(
t
)
=
私
(
t
)
記号
(
質問
(
t
)
)
−
質問
(
t
)
記号
(
私
(
t
)
)
。
{\displaystyle u_{d}(t)=I(t)\operatorname {sgn}(Q(t))-Q(t)\operatorname {sgn}(I(t)).}
次に、信号は ループフィルタによってフィルタリングされ、BPSKコスタスループと同様にVCOのチューニング信号を形成します。したがって、QPSKコスタスは 常微分方程式 のシステムによって 記述できます [14] 。
あなた
d
(
t
)
{\displaystyle u_{d}(t)}
あなた
LF
(
t
)
{\displaystyle u_{\text{LF}}(t)}
x
˙
1
=
あ
LPF
x
1
+
b
LPF
φ
1
(
t
)
、
x
˙
2
=
あ
LPF
x
2
+
b
LPF
φ
2
(
t
)
、
x
˙
=
あ
LF
x
+
b
LF
(
c
LPF
∗
x
1
記号
(
c
LPF
∗
x
2
)
−
c
LPF
∗
x
2
記号
(
c
LPF
∗
x
1
)
)
、
θ
˙
VC
=
ω
VC
無料
+
け
VCO
(
c
LF
∗
x
+
h
LF
(
c
LPF
∗
x
1
記号
(
c
LPF
∗
x
2
)
−
c
LPF
∗
x
2
記号
(
c
LPF
∗
x
1
)
)
)
。
{\displaystyle {\begin{aligned}{\dot {x}}_{1}&=A_{\text{LPF}}x_{1}+b_{\text{LPF}}\varphi _{1}(t),\\{\dot {x}}_{2}&=A_{\text{LPF}}x_{2}+b_{\text{LPF}}\varphi _{2}(t),\\{\dot {x}}&=A_{\text{LF}}x+b_{\text{LF}}{\big (}c_{\text{LPF}}^{*}x_{1}\operatorname {sgn}(c_{\text{LPF}}^{*}x_{2})-c_{\text{LPF}}^{*}x_{2}\operatorname {sgn}(c_{\text{LPF}}^{*}x_{1}){\big )},\\{\dot {\theta }}_{\text{vco}}&=\omega _{\text{vco}}^{\text{free}}+K_{\text{VCO}}{\Big (}c_{\text{LF}}^{*}x+h_{\text{LF}}{\big (}c_{\text{LPF}}^{*}x_{1}\operatorname {sgn}(c_{\text{LPF}}^{*}x_{2})-c_{\text{LPF}}^{*}x_{2}\operatorname {sgn}(c_{\text{LPF}}^{*}x_{1}){\big )}{\Big )}.\\\end{aligned}}}
以下は LPF1 と LPF2 のパラメータであり、以下は ループ フィルタのパラメータです。
A
LPF
,
b
LPF
,
c
LPF
{\displaystyle A_{\text{LPF}},b_{\text{LPF}},c_{\text{LPF}}}
A
LF
,
b
LF
,
c
LF
,
h
LF
{\displaystyle A_{\text{LF}},b_{\text{LF}},c_{\text{LF}},h_{\text{LF}}}
参考文献
^ コスタス 、ジョン P. (1956)。「同期通信」 。IRE 議事録 。44 (12): 1713–1718。doi :10.1109/jrproc.1956.275063 。
^ Costas, John P. (2002年8月) [1956]. 「同期通信」. Proceedings of the IEEE . 90 (8): 1461–1466. doi :10.1109/JPROC.2002.800876.
^ ab Taylor, D. (2002 年 8 月)。「 『 同期通信』入門、John P. Costas の古典論文」。IEEE の議事録 。90 (8): 1459–1460。doi :10.1109/jpcroc.2002.800719 。
^ Feigin, Jeff (2002 年 1 月 1 日). 「実用的な Costas ループ設計」 (PDF) . RF 設計 : 20–36. 2012 年 2 月 11 日時点の オリジナル (PDF)からアーカイブ 。2010 年 2 月 17 日 閲覧 。
^ Leonov, GA; Kuznetsov, NV; Yuldashev, MV; Yuldashev, RV (2012 年 8 月). 「コスタスループの微分方程式」 (PDF) . Doklady Mathematics . 86 (2): 723–728. doi :10.1134/s1064562412050080. S2CID 255276607.
^ Leonov, GA; Kuznetsov, NV; Yuldashev, MV; Yuldashev, RV (2012). 「位相検出器特性の計算のための解析的手法」 (PDF) . IEEE Transactions on Circuits and Systems Part II . 59 (10): 633–637. doi :10.1109/tcsii.2012.2213362. S2CID 2405056. [ 永久リンク切れ ]
^ Leonov, GA; Kuznetsov, NV; Yuldashev, MV; Yuldashev, RV (2015). 「コスタスループの非線形動的モデルと大規模におけるその安定性の分析へのアプローチ」 (PDF) . 信号処理 . 108 . Elsevier: 124–135. doi : 10.1016/j.sigpro.2014.08.033 . 2022-10-09 にオリジナルからアーカイブ (PDF)されました。
^ Kuznetsov, NV; Leonov, GA; Neittaanmaki, P.; Seledzhi, SM; Yuldashev, MV; Yuldashev, RV (2012). 「入力信号の一般波形に対するコスタスループの非線形数学モデル」。2012 IEEE 4th International Conference on Nonlinear Science and Complexity (NSC) 。IEEE Press。pp. 75–80。doi :10.1109 / NSC.2012.6304729。ISBN 978-1-4673-2703-9 . S2CID 5812970。
^ Kuznetsov, NV; Leonov, GA; Seledzhi, SM; Yuldashev, MV; Yuldashev, RV (2017). 「光コスタスループの非線形モデル:プルイン範囲の推定と隠れた振動」. IFAC-PapersOnLine . 50. ELSEVIER: 3325–3330. doi : 10.1016/j.ifacol.2017.08.514 . ISSN 2405-8963.
^ Costas 1956 は、「局部発振器は、上側と下側のサイドバンドのオーディオ出力寄与が互いに強化されるように、適切な位相に維持される必要があります。発振器の位相が最適値から 90° 離れると、オーディオ出力がゼロになります。これは、このタイプの検出器に典型的です。位相制御の実際の方法についてはすぐに説明しますが、この説明では、正しい発振器位相が維持されていると想定します。」と述べています。
^ 積分器を備えたループフィルタを使用すると、定常位相誤差をゼロにすることができます。PID コントローラ§積分項 を参照してください。
^ ベスト、ローランド E. (1997)。 フェーズロックループ (第 3 版)。ニューヨーク:マグロウヒル。pp. 44–45。ISBN 0-07-006051-7 。
^ 米国特許 4,085,378、Carl R. Ryan および James H. Stilwell、「QPSK 復調器」、1976 年 11 月 26 日公開、Motorola Solutions Inc. に譲渡。
^ Best, RE; Kuznetsov, NV; Leonov, GA; Yuldashev, MV; Yuldashev, RV (2016). 「コスタスループの動的解析に関するチュートリアル」. Annual Reviews in Control . 42. ELSEVIER: 27–49. arXiv : 1511.04435 . doi :10.1016/j.arcontrol.2016.08.003. S2CID 10703739.