【情報理論】生成行列の求め方と誤り訂正方法、問題

問題

下記の検査行列Hを考える。次の問いに答えよ。

(1)H=(101110011100100111001)

(1)生成行列Gを求めよ。
(2)Gにより生成する符号語を全て求めよ。
(3)[1 1 1 1 1 1 0]を受信したときの誤り訂正を実施せよ。

生成行列とは

情報ビットから符号を生成するための行列です。

大前提として、本記事では情報理論の内、符号の誤りを検出する方法の一つを紹介します。

符号とは、情報を表す情報ビットと検査部分を表す検査ビットで分かれています。

検査行列を活用し、誤りを訂正する方法もありますが、本記事では生成行列から誤りを確認する方法を紹介します。

入力符号に生成行列を積算すると、取り得る出力を網羅的に入手できます。取り得る出力どれにも該当しない出力を検出する=誤りであることが分かる大方針の元、説明していきます。

生成行列の求め方

手順1:検査行列Hを成分分解する

下記のように分解できます。

(2)H=[PtI]

P情報・検査ビット関連行列と言います。検査行列から生成行列を求める。またはその逆を行うために必ず出てくる行列です。キーになりますので覚えておきましょう。

検査行列は正方行列ではなく、列の数が行の数より多い行列です。

そのため、行の個数分を単位行列Iで表しても、単位行列に属さない列が出てきます。その分を情報・検査ビット関連行列で表せば良いです。

手順2:Ptを転置する

生成行列Gは、下記のように表すことができます。

(3)G=[IP]

手順1で求めた情報・検査ビット関連行列(転置)をもう一度転置し、行の成分の数だけ単位行列を追加します。

生成する符号の組み合わせを生成行列を用いて求める

入力する符号ベクトルをcとすれば、下記の計算で生成する符号の候補を求めることができます。

(4)cG

求め方の例

例)3行1列ベクトル c=(c0,c1,c2)

生成行列G=(100010101100111)を考える。

このとき、(4)式を計算すると

(5)(c0,c1,c2,c1c2,c0c1c2)になる。

c0 c2は0,1しか取らないのだから、23=8通りの入力に対し(5)式は下記の出力しか取りえないです。

c0c1c2c1c2c0c1c2
00000
00111
01011
01100
10000
10110
11010
11101

もし、上記の出力以外が出てきたら、それは誤りであることが分かります。

注意点

本問は線形符号であることが前提です。そのため、ciの和は排他的論理和で表すことが必要です。

線形関数についての詳しい記事は、こちらで説明しています。

解答例

(1)生成行列の求め方

行の個数は3つのため、単位行列は3*3成分である。

Hは(7,3)成分のため、後半の5列目~7列目が単位行列成分

一方で、前半の1~4列目の部分は情報・検査ビット関連行列成分。

(2)式より

(6)Pt=(101111100111)

が得られる。(3)式により、Gを求める。

(7)P=(110011111101)

行の個数が4つのため、単位行列は4*4成分。

(8)G=[I,P]=(1000110010001100101110001101)

(2)符号語の生成

(1)より、行の成分は4つ。よって、入力符号c

(9)c=(c0,c1,c2,c3)で4成分定義する。

この時の生成符号は

(10)cG=(c0,c1,c2,c3,c0c2c3,c0c1c2,c1c2c3)

ci(0,1)のため、考えられる生成符号は下記のパターンになる。

c0c1c2c3c0c2c3c0c1c2c1c2c3
0000000
0001101
0010111
0011010
0100011
0101110
0110100
0111001
1000110
1001011
1010001
1011100
1100101
1101000
1110010
1111111

(3)誤り訂正

[1 1 1 1 1 1 0]は、(2)の符号語の一覧表に無いため、誤りである。

最もハミング距離が小さい符号は、[1 1 1 1 1 1 1]のため、これに訂正する。

タイトルとURLをコピーしました