このページの翻訳は最新ではありません。ここをクリックして、英語の最新版を参照してください。

ldpcDecode

バイナリ LDPC 符号の復号化

R2021b 以降

構文

[y,actualnumiter,finalparitychecks] = ldpcDecode(llr,decodercfg,maxnumiter)

[y,actualnumiter,finalparitychecks] = ldpcDecode(llr,decodercfg,maxnumiter,Name,Value)

説明

関数 ldpcDecode は、4 つのアルゴリズムのうちのいずれかを使用して入力コードワードを復号化します。詳細については、アルゴリズムを参照してください。LDPC 符号は、スパース パリティ チェック行列およびシャノン限界に近いパフォーマンスを達成できる長いブロック長をもつ、線形の誤り制御符号です。

[y,actualnumiter,finalparitychecks] = ldpcDecode(llr,decodercfg,maxnumiter) は、入力 ldpcDecoderConfig 構成オブジェクト decodercfg で指定された LDPC 行列を使用して、入力対数尤度比 (LLR) llr を復号化します。正の LLR は、対応するビットがゼロである可能性が高いことを示します。入力 maxnumiter で指定された反復の最大数以内ですべてのパリティ チェックが満たされると、復号化は終了します。LDPC 符号は、スパース パリティ チェック行列およびシャノン限界に近いパフォーマンスを達成できる長いブロック長をもつ、線形の誤り制御符号です。

[y,actualnumiter,finalparitychecks] = ldpcDecode(llr,decodercfg,maxnumiter,Name,Value) は、名前と値の引数を 1 つ以上使用してオプションを指定します。たとえば、DecisionType='soft' は、軟判定復号化を指定し、LLR を出力します。

例

符号化率 3/4 の LDPC コードワードの復号化

IEEE® 802.11 で指定された符号化率 3/4 の LDPC 符号を構成するために、プロトタイプ行列とブロック サイズのパラメーターを初期化します。関数 ldpcQuasiCyclicMatrix を使用して、パリティチェック行列を作成します。

P = [
    16 17 22 24  9  3 14 -1  4  2  7 -1 26 -1  2 -1 21 -1  1  0 -1 -1 -1 -1
    25 12 12  3  3 26  6 21 -1 15 22 -1 15 -1  4 -1 -1 16 -1  0  0 -1 -1 -1
    25 18 26 16 22 23  9 -1  0 -1  4 -1  4 -1  8 23 11 -1 -1 -1  0  0 -1 -1
     9  7  0  1 17 -1 -1  7  3 -1  3 23 -1 16 -1 -1 21 -1  0 -1 -1  0  0 -1
    24  5 26  7  1 -1 -1 15 24 15 -1  8 -1 13 -1 13 -1 11 -1 -1 -1 -1  0  0
     2  2 19 14 24  1 15 19 -1 21 -1  2 -1 24 -1  3 -1  2  1 -1 -1 -1 -1  0
    ];
blockSize = 27;
pcmatrix = ldpcQuasiCyclicMatrix(blockSize,P);

LDPC 符号化器と LDPC 復号化器の構成オブジェクトを作成し、それらのプロパティを表示します。

cfgLDPCEnc = ldpcEncoderConfig(pcmatrix)

cfgLDPCEnc = 
  ldpcEncoderConfig with properties:

     ParityCheckMatrix: [162x648 logical]

   Read-only properties:
           BlockLength: 648
    NumInformationBits: 486
    NumParityCheckBits: 162
              CodeRate: 0.7500

cfgLDPCDec = ldpcDecoderConfig(pcmatrix)

cfgLDPCDec = 
  ldpcDecoderConfig with properties:

     ParityCheckMatrix: [162x648 logical]
             Algorithm: 'bp'

   Read-only properties:
           BlockLength: 648
    NumInformationBits: 486
    NumParityCheckBits: 162
              CodeRate: 0.7500

LDPC 符号化され、QPSK 変調されたビット ストリームを AWGN チャネル経由で送信します。信号を復調し、受信コードワードを復号化してから、ビット エラーをカウントします。入れ子の for ループを使用し、複数の SNR 設定、および送信データの LDPC 前方誤り訂正 (FEC) 符号化がある場合とない場合のフレームについて処理します。

M = 4;
maxnumiter = 10;
snr = [3 6 20];
numframes = 10;
qpskmod = comm.PSKModulator(M,'BitInput',true);
qpskmod2 = comm.PSKModulator(M);

ber = comm.ErrorRate;
ber2 = comm.ErrorRate;

for ii = 1:length(snr)
    qpskdemod = comm.PSKDemodulator(M,'BitOutput',true, ...
        'DecisionMethod','Approximate log-likelihood ratio', ...
        'Variance',1/10^(snr(ii)/10));
    qpskdemod2 = comm.PSKDemodulator(M);
    for counter = 1:numframes
        data = randi([0 1],cfgLDPCEnc.NumInformationBits,1,'int8');
        % Transmit and receive with LDPC coding
        encodedData = ldpcEncode(data,cfgLDPCEnc);
        modSignal = qpskmod(encodedData);
        receivedSignal = awgn(modSignal,snr(ii));
        demodSignal = qpskdemod(receivedSignal);
        receivedBits = ldpcDecode(demodSignal,cfgLDPCDec,maxnumiter);
        errStats = ber(data,receivedBits);
        % Transmit and receive with no LDPC coding
        noCoding = qpskmod2(data);
        rxNoCoding = awgn(noCoding,snr(ii));
        rxBitsNoCoding = qpskdemod2(rxNoCoding);
        errStatsNoCoding = ber2(data,int8(rxBitsNoCoding));
    end
    fprintf(['SNR = %2d\n   Coded: Error rate = %1.2f, ' ...
        'Number of errors = %d\n'], ...
        snr(ii),errStats(1),errStats(2))
    fprintf(['Noncoded: Error rate = %1.2f, ' ...
        'Number of errors = %d\n'], ...
        errStatsNoCoding(1),errStatsNoCoding(2))
    reset(ber);
    reset(ber2);
end

SNR =  3
   Coded: Error rate = 0.07, Number of errors = 335

Noncoded: Error rate = 0.15, Number of errors = 714

SNR =  6
   Coded: Error rate = 0.00, Number of errors = 0

Noncoded: Error rate = 0.04, Number of errors = 196

SNR = 20
   Coded: Error rate = 0.00, Number of errors = 0

Noncoded: Error rate = 0.00, Number of errors = 0

入力引数

llr — 対数尤度比

行列

対数尤度比。入力 decodercfg の BlockLength プロパティと等しい行数をもつ行列として指定します。llr の各列は、1 つのコードワードに対応します。関数は各列を個別に復号化します。正の LLR は、対応するビットがゼロである可能性が高いことを示します。

データ型: single | double

decodercfg — LDPC 復号化器の構成

ldpcDecoderConfig オブジェクト

LDPC 復号化器の構成。ldpcDecoderConfig オブジェクトとして指定します。

maxnumiter — 復号化反復の最大数

正のスカラー

復号化反復の最大数。正のスカラーとして指定します。

データ型: double

名前と値の引数

オプションの引数のペアを Name1=Value1,...,NameN=ValueN として指定します。ここで、Name は引数名で、Value は対応する値です。名前と値の引数は他の引数の後に指定しなければなりませんが、ペアの順序は重要ではありません。

例: Termination='max'

OutputFormat — 出力形式

'info' (既定値) | 'whole'

出力形式。以下のいずれかの値として指定します。

• 'info' — 復号化された情報ビットのみを出力します。関数が出力する行数は、入力 decodercfg の NumInformationBits プロパティと同じになります。

• 'whole' — 情報ビットとパリティチェック ビットを含む、復号化された LDPC コードワード ビットすべてを出力します。関数が出力する行数は、入力 decodercfg の BlockLength プロパティと同じになります。

データ型: char | string

DecisionType — Decision type

'hard' (既定値) | 'soft'

LDPC 復号化の判定タイプ。次のいずれかの値として指定します。

• 'hard' — 硬判定復号化を実行し、復号化されたビットを int8 データ型の値として出力します。

• 'soft' — 軟判定復号化を実行し、入力と同じデータ型で LLR を出力します。

データ型: char | string

MinSumScalingFactor — 正規化された min-sum 復号化アルゴリズムのスケーリング係数

0.75 (既定値) | 範囲 (0, 1] のスカラー

正規化された min-sum 復号化アルゴリズムのスケーリング係数。範囲 (0, 1] のスカラーとして指定します。詳細については、正規化された min-sum 復号化を参照してください。

依存関係

このプロパティを有効にするには、入力 decodercfg の Algorithm プロパティを 'norm-min-sum' に設定します。

データ型: double

MinSumOffset — min-sum 復号化アルゴリズムのオフセット

0.5 (既定値) | スカラー

min-sum 復号化アルゴリズムのオフセット。スカラーとして指定します。詳細については、オフセット min-sum 復号化を参照してください。

依存関係

このプロパティを有効にするには、入力 decodercfg の Algorithm プロパティを 'offset-min-sum' に設定します。

データ型: double

Termination — 復号化終了基準

'early' (既定値) | 'max'

復号化終了基準。次のいずれかの値として指定します。

• 'early' — 入力 maxnumiter で指定された反復の最大数以内ですべてのパリティ チェックが満たされると、復号化の反復を終了します。

• 'max' — 最大数 maxnumiter の反復が完了したときに復号化を終了します。

データ型: char | string

Multithreaded — マルチスレッド実行の有効化

true または 1 (既定値) | false または 0

マルチスレッド実行の有効化。logical 1 (true) または logical 0 (false) として指定します。インタープリター モードで MATLAB^® を実行し、この引数を true に設定すると、関数は複数のスレッドで復号化アルゴリズムを実行します。

ヒント

大きなパリティチェック行列では、マルチスレッド実行により、LDPC 復号化の処理時間が大幅に短縮されます。

依存関係

このプロパティを有効にするには、MATLAB をインタープリター モードで実行します。

出力引数

y — 復号化されたコードワード

行列

復号化されたコードワード。llr(1:K,:) の復号化されたビットを表す K 行の行列として返されます。K は、入力 decodercfg の NumInformationBits プロパティと同じです。復号化操作の場合、llr の各列は 1 つのコードワードに対応します。関数は各列を個別に復号化します。名前と値の引数 'OutputFormat' は、出力に復号化された情報ビットを含める (既定) か、LDPC コードワード ビット全体を含めるかを指定します。名前と値の引数 'DecisionType' は、復号化判定タイプとこの出力のデータ型を指定して判定します。

詳細については、アルゴリズムを参照してください。

データ型: int8 | double | single

actualnumiter — 復号化の実際の反復回数

行ベクトル

復号化の実際の反復回数。行ベクトルとして返されます。コードワードのすべてのパリティ チェックが満たされると、反復の最大数 maxnumiter に達する前でも復号化が停止することがあります。この出力は、関数がコードワードに対して実行した実際の反復回数から成る行ベクトルです。

データ型: double

finalparitychecks — 各コードワードに対する最終パリティ チェック

行列

各コードワードに対する最終パリティ チェック。入力 decodercfg の ParityCheckBits プロパティと等しい行数をもつ行列として返されます。復号化操作の場合、この出力の各列は、対応するコードワードに対する最終パリティ チェックになります。

データ型: double

アルゴリズム

以下のメッセージ伝達アルゴリズムのいずれかを使用する LDPC 復号化。

確率伝播復号化

確率伝播アルゴリズムの実装は、Gallager によって提示された復号化アルゴリズムに基づいています [2]。

LDPC 符号化されて送信されたコードワード c = c₀、c₁、...、c_n-1 に対して、LDPC 復号化器への入力は、対数尤度比 (LLR) の値 $$L(c_i)=\log \left(\frac{\mathrm{Pr}(c_i=0|\text{channel output for }{c}_i)}{\mathrm{Pr}(c_i=1|\text{channel output for }{c}_i)}\right)$$ になります。

各反復で、アルゴリズムの主なコンポーネントが以下の式に基づいて更新されます。

$$L(r_{ji})=2\text{atanh}\left({\displaystyle \prod _{i^{\prime }\in V_j\backslash i}\tanh \left(\frac{1}{2}L(q_{i^{\prime }j})\right)}\right)$$,

$$L(q_{ij})=L(c_i)+{\displaystyle \sum _{j^{\prime }\in C_i\backslash j}L(r_{j^{\prime }i})}$$ (最初の反復の前に $$L(q_{ij})=L(c_i)$$ として初期化される)、および

$$L(Q_i)=L(c_i)+{\displaystyle \sum _{j^{\prime }\in C_i}L(r_{j^{\prime }i})}$$.

それぞれの反復の最後で、L(Q_i) には、送信ビット c_i に対する LLR 値の更新された推定値が含まれます。値 L(Q_i) は c_i の軟判定出力です。L(Q_i) < 0 の場合、c_i の硬判定出力は 1 です。それ以外の場合、c_i の硬判定出力は 0 です。

すべてのパリティ チェックが満たされた場合に停止するように復号化が設定されている場合、アルゴリズムは、各反復の終わりにパリティ チェック式 (H c' = 0) を検証します。すべてのパリティ チェックが満たされた場合、または最大反復回数に到達した場合、復号化は停止します。

インデックス集合 $$C_i\backslash j$$ と $$V_j\backslash i$$ は、パリティ チェック行列 (PCM) に基づいています。インデックス集合 C_i と V_j はそれぞれ、PCM の列 i と行 j のすべての非ゼロ要素に対応します。

次の図は、与えられた PCM の i = 5 と j = 3 に対するインデックス集合の計算を示しています。

無限数がアルゴリズム式で使用されないように、atanh(1) は 19.07 に、atanh(–1) は –19.07 に設定されます。有限値で演算を行うために、MATLAB は tanh(19.07) に対して 1 を、tanh(-19.07) に対して –1 を返します。

階層的確率伝播復号化

階層的確率伝播アルゴリズムの実装は、Hocevar [3] の II.A 節にある復号化アルゴリズムに基づいています。復号化ループは、PCM の行 (層) のサブセットを反復します。レイヤー内の各行 m と各ビット インデックス j について、実装ではアルゴリズムの主なコンポーネントが以下の式に基づいて更新されます。

(1) $$L(q_{mj})=L(q_j)-R_{mj}$$,

(2) $$A_{mj}={\displaystyle \sum _{\begin{array}{c}n \in N\left(m\right)\\ n\ne j\end{array}}\text{ψ}(L(q_{mn}))}$$,

(3) $$s_{mj}={\displaystyle \prod _{\begin{array}{c}n \in N\left(m\right)\\ n\ne j\end{array}}\text{sign}(L(q_{mn}))}$$,

(4) $$R_{mj}=-s_{mj}\text{ψ}(A_{mj})$$、および

(5) $$L(q_j)=L(q_{mj})+R_{mj}$$.

各レイヤーについて、復号化の式 (5) が、現在の LLR 入力 $$L(q_{mj})$$ と前のレイヤーの更新 $$R_{mj}$$ から得られる結合された入力に対して適用されます。

ノードのサブセットのみがレイヤーで更新されるため、階層的確率伝播アルゴリズムは、確率伝播アルゴリズムと比較して高速です。階層的確率伝播アルゴリズムを使用すれば、確率伝播復号化が達成するのと同じエラー レートを達成するために使用する復号化反復回数が半分で済みます。

正規化された min-sum 復号化

正規化された min-sum 復号化アルゴリズムの実装は、式 (2) を次のように置き換えた階層的確率伝播アルゴリズムに従います。

$$A_{mj}=\underset{\begin{array}{c}n \in N\left(m\right)\\ n\ne j\end{array}}{\min }(\left|L(q_{mn}) \right|\cdot \alpha )$$,

ここで、α は、関数 ldpcDecode への入力引数 MinSumScalingFactor で指定される、範囲 (0, 1] のスケーリング係数です。この式は、Chen [4] にある式 (4) を応用したものです。

オフセット min-sum 復号化

オフセット min-sum 復号化アルゴリズムの実装は、式 (2) を次のように置き換えた階層的確率伝播アルゴリズムに従います。

$$A_{mj} = \max (\underset{\begin{array}{c}n \in N\left(m\right)\\ n\ne j\end{array}}{\min } (\left|L\left(q_{mn}\right)\right|- \beta ), 0)$$,

ここで、β は、0 以上のオフセットで、関数 ldpcDecode への入力引数 MinSumOffset で指定されます。この式は、Chen [4] にある式 (5) を応用したものです。

参照

[1] IEEE Std 802.11™-2020 (Revision of IEEE Std 802.11-2016). "Part 11: Wireless LAN Medium Access Control (MAC) and Physical Layer (PHY) Specifications." IEEE Standard for Information technology — Telecommunications and information exchange between systems. Local and metropolitan area networks — Specific requirements.

[2] Gallager, Robert G. Low-Density Parity-Check Codes. Cambridge, MA: MIT Press, 1963.

[3] Hocevar, D.E. "A reduced complexity decoder architecture via layered decoding of LDPC codes." In IEEE Workshop on Signal Processing Systems, 2004. SIPS 2004. doi: 10.1109/SIPS.2004.1363033

[4] Chen, Jinghu, R.M. Tanner, C. Jones, and Yan Li. "Improved min-sum decoding algorithms for irregular LDPC codes." In Proceedings. International Symposium on Information Theory, 2005. ISIT 2005. doi: 10.1109/ISIT.2005.1523374

拡張機能

C/C++ コード生成
MATLAB® Coder™ を使用して C および C++ コードを生成します。

バージョン履歴

R2021b で導入

参考

関数

• ldpcEncode | ldpcQuasiCyclicMatrix | dvbs2ldpc

オブジェクト

• ldpcDecoderConfig | ldpcEncoderConfig | comm.gpu.LDPCDecoder

×

例を開く

この例の変更されたバージョンがあります。編集された方の例を開きますか?

いいえ、変更されたバージョンを上書きします。 はい