HDL QPSK 送信機および受信機

この例では次を使用します。

この例では、Simulink® で HDL コード生成とハードウェア実装用に最適化された QPSK 送受信機を実装する方法を示します。

この例で示すモデルは、直交位相偏移変調 (QPSK) に基づいてデータを変調します。この例の目標は、リアルタイムシステム用に情報を送信して復元できる HDL QPSK 通信システムをモデル化することです。受信機で、単一搬送波通信システムに欠かせないシンボルタイミング同期および搬送波の周波数と位相の同期を実装します。

システム仕様

このセクションでは、この例で使用する送信機と受信機の仕様について説明します。フレームフォーマットはパケットベースです。各パケットに長さ 26 ビットのプリアンブルが含まれます。QPSK 変調器で同一のビットが同相と直交位相で変調されるように、13 ビットのバーカーシーケンスの各ビットを 2 回繰り返してプリアンブルシーケンスが生成されます。プリアンブルシーケンスの後に、2240 ビットのペイロードデータが続きます。送信機はルートレイズドコサイン (RRC) パルス整形フィルターを使用して実行されます。ロールオフ係数は 0.5 であるため、帯域幅はシンボルレートの 1.5 倍、シンボルあたりのサンプル数は 4 (サンプルレートはシンボルレートの 4 倍) になります。RRC インパルス応答は隣接する 4 つのシンボルにまたがります。ビットレートは、シンボルレートの 2 倍です。実効平均ビットレートは、ビットレートにフレーム効率を掛けて求められます。フレーム効率は (2240/(2240+26)) = 0.9885 です。

既定のシンボルレートは 1.92 Mbaud に設定されています。したがって、帯域幅は 1.92e6 の 1.5 倍で 2.88 MHz、サンプルレートは 1.92e6 の 4 倍で 7.68 Msps、ビットレートは 1.92e6 の 2 倍で 3.84 Mbps になります。このシステムでサポートされる実効平均ビットレートは、3.84e6 に 0.9885 を掛けて 3.7959 Mbps となります。これらの仕様はシンボルレートを変更すると変わります。

モデルアーキテクチャ

このセクションでは、ブロック線図に示す QPSK 送受信機の大まかなアーキテクチャについて説明します。QPSK 送信機は、シンボルレートの 2 倍のビットレートで入力をサンプリングします。Data Generator & Packetizer は、データビットを収集し、プリアンブルビットを生成し、パケットビットを形成します。HDL Data Scrambler は、各パケットのデータビットをスクランブルしてビット遷移を増やし、同じビットの実行シーケンスが長くならないようにします。QPSK Modulator は、パケットビットを変調して QPSK シンボルを生成します。RRC Transmit Filter は、QPSK シンボルをアップサンプリングしてパルス整形し、シンボルレートの 4 倍のサンプルレートの Tx 波形を生成します。QPSK 受信機は、伝送速度で入力をサンプリングします。Digital AGC は、受信波形が目的の振幅レベルになるようにゲイン制御を実行します。RRC Receive Filter は、AGC の出力に対して整合フィルター処理を実行します。Frequency and Time Synchronizer は、同期処理を実行し、各パケットの QPSK シンボルを生成します。QPSK Demodulator は、QPSK シンボルを復調してパケットビットを生成します。HDL Data Descrambler は、受信機からストリーミングされるパケットのデータビットをデスクランブルします。

ファイル構造

この例は、1 つの Simulink モデルと 3 つの MATLAB® ファイルで構成されます。

commhdlQPSKTxRx.slx — 最上位の Simulink モデル
commhdlQPSKTxRxParameters.m — 初期化に必要な QPSK Tx と QPSK Rx のパラメーターを生成
commhdlQPSKTxRxModelInit.m — モデル commhdlQPSKTxRx.slx を初期化
generateHelloworldMsgBits.m — "Hello world xxx " メッセージビットを生成。"xxx" は 000 から 100 までの値を参照

システムインターフェイス

次の図は、QPSK 送受信機システムの最上位モデルを示しています。

送信機の入力

dataIn — 入力データ。boolean スカラーとして指定します。
validIn — "dataIn" を有効にするための制御信号。boolean スカラーとして指定します。

送信機の出力

dataOut — 出力送信波形。シンボルレートの 4 倍のサンプルレートの 16 ビット複素数データとして返されます。
validOut — "dataOut" を有効にするための制御信号。boolean スカラーとして返されます。
txDiagBus — 診断出力を含むステータス信号。バス信号として返されます。
dataReady — 入力信号に対する準備状況を示す信号。boolean スカラーとして返されます。

送信機は、入力ビットを受け入れる準備ができると、"dataReady" 信号を有効にしてそれを示します。送信機は、そのパケットに対応するすべてのデータビットを受け入れてからパケットを構成します。そのパケットに対応するすべてのデータビットが受信されないと、送信機はダミーパケットを生成します。ダミーパケットのバーカーシーケンスはプリアンブルに使用されず、データビットが内部でランダムに生成されます。入力ビットレートが実効ビットレート以下である間は、入力が "dataReady" からのバックプレッシャーを受け取らないように、"dataReady" 信号は High のままになります。

受信機の入力

dataIn — 入力データ。送信機出力と同じサンプルレートの 16 ビット複素数データとして指定します。
validIn — "dataIn" を有効にするための制御信号。boolean スカラーとして指定します。

受信機の出力

dataOut — 復号化された出力データビット。boolean スカラーとして返されます。
ctrlOut — start、end、valid の信号を含むバス信号。バス信号として返されます。
rxDiagBus — 診断出力を含むステータス信号。バス信号として返されます。

送信機の構造

次の図は、QPSK Tx Subsystem の最上位モデルを示しています。

Bit Packetizer

Bit Packetizer Subsystem は、Packet Controller MATLAB Function、Bits Store Subsystem、および Multiplexer Subsystem で構成されます。プリアンブルシーケンスは、Preamble Bits Store Subsystem 内のルックアップテーブル (LUT) に格納されます。データビットは、Bits Store Subsystem にストリーミングされ、Data Bits Store Subsystem 内の RAM に格納されます。Packet Controller MATLAB Function は、各パケットのプリアンブルシーケンスと RAM に格納された後続のデータビットを読み取ります。Multiplexer Subsystem は、プリアンブルビットとデータビットを効率化します。

Data Bits Store Subsystem は、2 つのパケットを格納できる RAM で構成されます。この RAM により、離散の有効な入力で送信機が動作する柔軟性が提供されます。Packet Controller MATLAB Function は、RAM に少なくとも 1 つのパケットのデータビットが格納されている場合にのみ RAM からデータを読み取ります。RAM がオーバーフローしないように読み取りと書き込みのロジックが設計されています。RAM に少なくとも 1 つのパケットのデータビットが格納されていない場合、送信機はダミーパケットを生成します。プリアンブル検出で検出されないように、ダミーパケットのバーカーシーケンスはプリアンブルに使用されません。

HDL Data Scrambler

HDL Data Scrambler Subsystem は、Bits Generator Subsystem によって生成される制御信号を使用して各パケットのデータビットをスクランブルします。

QPSK Modulator

QPSK Modulator Subsystem は、QPSK Modulator Basebandブロックを使用して、プリアンブルビットとデータビットを変調することで QPSK シンボルを生成します。次の表に示すグレイマッピングが使用されます。

    Bits         Mapping     
    ____    _________________

     00     0.70711+0.70711i 
     01     -0.70711+0.70711i
     11     -0.70711-0.70711i
     10     0.70711-0.70711i

RRC Transmit Filter

RRC Transmit Filter Subsystem は、入力を 4 倍にアップサンプリングし、Discrete FIR Filter (DSP HDL Toolbox)ブロックと RRC インパルス応答を使用して送信機の波形をパルス整形します。QPSK 受信機の受信フィルターで、この送信フィルターに対する整合フィルターが形成されます。

受信機の構造

次の図は、QPSK Rx Subsystem の最上位モデルを示しています。

Automatic Gain Control

シンボル同期装置とキャリア同期装置の位相同期回路 (PLL) の性能には入力信号の振幅が影響するため、それらの前に Automatic Gain Control Subsystem を配置しています。振幅の二乗の出力を AGC の基準と比較して振幅の誤差が生成されます。この誤差にループゲインを乗算して積分器を通過させることで、必要なゲインが計算されます。結果のゲインに AGC の入力を乗算して AGC の出力が生成されます。詳細については、[ 1 ] の Chapter 9.5 を参照してください。

RRC Receive Filter

RRC Receive Filter は、送信機でパルス整形に使用されたフィルターの整合フィルター係数をもつDiscrete FIR Filter (DSP HDL Toolbox)ブロックです。RRC 整合フィルター処理により、波形のアイダイアグラムのアイ開口が最大になる時点においてゼロ ISI 特性となる RC パルス整形波形が生成されます。また、整合フィルター処理の過程で、フィルター出力の信号のノイズパワーに対する割合 (SNR) が最大化されます。

Frequency and Time Synchronizer

Frequency and Time Synchronizer Subsystem は、シンボル同期、大まかな周波数補正、搬送波同期、およびパケット同期用のプリアンブル検出を実行します。また、搬送波同期で補正されずに残る位相不確定性を推定して解決します。

Symbol Synchronizer Subsystem は PLL ベースの実装です。[ 1 ] の Chapter 8.5 で説明されているように、最適な時点 (アイ開口が最大になる時点) でサンプルを生成します。このサブシステムは 4 つの入力サンプルごとに 1 つの出力サンプルを生成します。Interpolation Filter Subsystem は、[ 1 ] の Chapter 8.4.2 で説明されているハードウェアリソース効率に優れた Farrow 構造をもつ区分的な放物線内挿器を実装します。Farrow の係数は Table 8.4.1 (係数の自由パラメーター $\alpha$ は 0.5) に示されています。このフィルターは入力波形に非整数遅延を導入します。Gardner TED Subsystem は、ガードナー法によるタイミング誤差検出器を実装します。タイミング誤差検出器については、[ 1 ] の Chapter 8.4.1 で説明されています。ループフィルターでタイミング誤差がフィルター処理され、そのタイミング誤差が Interpolation Control MATLAB Function ブロックに渡されます。このブロックは、[ 1 ] の Chapter 8.4.3 で説明されている剰余 1 の減算カウンターを実装することで、ループフィルターでフィルター処理されたタイミング誤差に基づいて非整数遅延を計算し、最適なサンプリング時点で内挿を生成します。Rate Handle Subsystem は、ストローブで示される必要な内挿を選択します。このサンプルは、シンボル同期前のアイダイアグラムのアイ開口が最大になる時点に対応します。

Coarse Frequency Compensator Subsystem は、Raise Power to 4 Subsystem で入力シーケンスを 4 のべき乗にします。これによって、入力シーケンスの QPSK 位相マッピングの依存関係は排除されますが、推定範囲が 4 分の 1 になります。このシーケンスが Coarse Frequency Estimator Subsystem にストリーミングされます。Coarse Frequency Estimator Subsystem から取得される推定は、4 のべき乗になっているため、係数 4 を削除して正規化された大まかな周波数推定を得るために 4 で除算されます。この推定で駆動されるNCO (DSP HDL Toolbox)ブロックにより、周波数オフセットを補正するための入力シーケンスで共役と乗算を行った複素指数位相が生成されます。

Coarse Frequency Estimator Subsystem は、差分で入力シーケンスを検出し、入力の複素周波数オフセット推定を抽出します。この推定が Integrator Subsystem で 2^15 の連続する推定で平均化され、最終的な複素推定が得られます。Extract Frequency Subsystem で、Complex to Magnitude-Angle (DSP HDL Toolbox)ブロックによって複素推定から周波数が抽出されます。

周波数推定器で推定されるのは、正規化周波数 (シンボルレート基準) の範囲で -0.125 ～ 0.125 になります。これは、1.92 Mbaud のシンボルレートでは、周波数オフセット範囲 -240 KHz ～ 240 KHz に対応します。推定の精度については、大まかな周波数オフセット補正後の残留周波数オフセットが正規化周波数の範囲 -0.0016 ～ 0.0016 に収まるようになります。これは、1.92 Mbaud のシンボルレートでは、Carrier Synchronizer の PLL が周波数オフセット範囲 -3 KHz ～ 3 KHz に収束することに対応します。

Carrier Synchronizer Subsystem は、操作点が 45 度の正弦波位相誤差検出器をもつ TYPE II PLL です。位相誤差検出器については [ 1 ] の Chapter 7.2.2、設計方程式については Appendix C で説明されています。操作点がゼロの正弦波位相誤差検出器をもつ TYPE II PLL の詳しい解析については、[ 2 ] の Chapter 4 で説明されています。位相検出器の実数部と虚数部の符号関数は、4 象限のすべての角度を第 1 象限の角度 (0 ～ 90 度) に変換します。これにより、第 2 象限 (90 ～ 180 度) の 90 度、第 3 象限 (-180 ～ -90 度) の 180 度、および第 4 象限 (-90 ～ 0 度) の 270 度の角度について不確定性が生じます。位相誤差は、位相検出器の操作点 (45 度) からの偏差として計算されます。Loop Filter Subsystem の比例積分器フィルターで位相誤差がフィルター処理されます。ループフィルターで正規化ループ帯域幅 (サンプルレートで正規化) とループ減衰係数が設定されます。既定の正規化ループ帯域幅は 0.005、既定の減衰係数は 0.7071 に設定されています。フィルター処理された誤差が位相インクリメントのソースとして Direct Digital Synthesis Subsystem に与えられます。このサブシステムでは、複素指数位相の生成にNCO (DSP HDL Toolbox)ブロックを使用します。この複素指数位相を使用して入力の周波数と位相が補正されます。直接デジタル合成の詳しい解析については、[ 1 ] の Chapter 9.2.2 で説明されています。

Preamble Detector Subsystem は、バーカーシーケンスを使用して入力の連続相関を実行します。この相関は、反転したバーカーシーケンスをDiscrete FIR Filter (DSP HDL Toolbox)ブロックの係数とする畳み込みとして実装されます。相関後の出力の振幅は、Correlator Subsystem 内のComplex to Magnitude-Angle (DSP HDL Toolbox)ブロックを使用して求められます。相関の振幅がしきい値と比較されます。Peak Search Subsystem は、しきい値を超えた最大相関ピークの検索を 1 フレーム時間ごとに開始し、タイミングオフセットを記録します。Timing Adjust Subsystem は、パケットのタイミングをタイミングオフセットに基づいて同期し、パケット同期済みのサンプルであることを後続のサブシステムに示す "syncPulse" 信号を生成します。

Phase Ambiguity Estimation and Correction Subsystem は、[ 1 ] の Chapter 7.7.1 で説明されているユニークワード方式に基づいて位相不確定性の解決を行います。この方式では、プリアンブルシーケンスを参照シーケンスとして使用します。入力のプリアンブルシーケンスで参照シーケンスの共役と乗算を行い、残留位相を位相不確定性の推定として抽出します。この推定を使用して、コンスタレーションを不確定性と反対方向に回転させることで不確定性を補正します。

Packet Controller Subsystem は、パケット境界の制御信号を生成します。

QPSK Demodulator

QPSK Demodulator Subsystem は、QPSK Demodulator Basebandブロックを使用して、パケット同期済みのシンボルを復調することでビットを生成します。

HDL Data Descrambler

HDL Data Descrambler Subsystem は、復調されたビットをデスクランブルしてユーザービットを生成します。このサブシステムは送信機側で使用されるスクランブラと同じです。

モデルの実行

commhdlQPSKTxRx.slx モデルを実行してシミュレートします。Input Data サブシステムマスクでカスタムデータを設定し、Channel サブシステムマスクでチャネル構成を設定できます。QPSKTxRxVerification.m スクリプトに、commhdlQPSKTxRx.slx モデルを検証する手続きが記述されています。この検証スクリプトは、スクリプト内で基準波形を生成してから、基準波形と送信機出力を比較し、送信されたビットと受信機の復号化されたユーザービットを比較します。

検証と結果

sim commhdlQPSKTxRx を実行してモデルを実行します。

>> sim commhdlQPSKTxRx;

Simulation completed
Running the verification script

QPSK Tx:
Maximum absolute symbol error: Real:1.4496e-05 Imaginary:1.4496e-05

Maximum absolute RRC output error: Real:7.8708e-05 Imaginary:7.8708e-05

QPSK Rx:
Initial frames not compared : 29

Number of packets missed = 0 out of 30

Number of packets false detected = 0 out of 30

Number of bits errored = 0 out of 67200

HDL コード生成

モデル全体にパイプラインレジスタ (シアンで表示) が追加されており、送信機と受信機のサブシステムに長いクリティカルパスがないことを確認できます。

この例で参照している HDL コードをチェックして生成するには、HDL Coder™ 製品が必要です。

送信機と受信機のサブシステムの HDL コードを生成するには、モデルを更新し、次のコマンドを使用します。

  makehdl('commhdlQPSKTxRx/QPSK Tx')  and   makehdl('commhdlQPSKTxRx/QPSK Rx')

テストベンチを生成するには、次のコマンドを使用します。

  makehdltb('commhdlQPSKTxRx/QPSK Tx')  and   makehdltb('commhdlQPSKTxRx/QPSK Rx')

テストベンチの生成時間はシミュレーション時間に依存します。

結果の HDL コードは、Xilinx® Zynq®-7000 ZC706 評価ボード用に合成されたものです。配置配線後のリソース使用状況を次の表に示します。最大動作周波数は、送信機が 276 MHz、受信機が 194 MHz です。

       Resources       Tx Usage    Rx Usage
    _______________    ________    ________

    Slice Registers      321        10441  
    Slice LUT            160        7490   
    RAMB36               0          10     
    RAMB18               1          0      
    DSP48                18         110

その他の調査

次の表に示す変数を commhdlQPSKTxRx.slx モデルの Channel サブシステムマスクで調整することで、チャネル条件を変更してモデルを実行できます。

        Variable Name                                               Description                                      
    ______________________    _______________________________________________________________________________________

    dataBits                  Data bits to the transmitter                                                           
    Rsym                      Symbol rate specified in symbols per second                                            
    fractionalTimingOffset    Normalized timing phase offset specified in the range >= 0 and < 1                     
    timingFrequencyOffset     Timing frequency offset specified in PPM                                               
    EbN0dB                    Energy per information bit to single sided noise power spectral density specified in dB
    CFO                       Carrier frequency offset specified in Hz                                               
    CPO                       Carrier phase offset specified in degrees