跳到主要内容

金沢大学 自然科学研究科 電子情報通信学専攻 2022年8月実施 専門科目 情報理論

Author

金沢大学

Description

無記憶情報源 S={s1,s2}S=\{s_1,s_2\} において s1s_1 の発生確率 P(s1)P(s_1)p, 0p1p,\ 0\leq p\leq 1 であった。この無記憶情報源 SS の出力を用いて a1=s1, a2=s2a_1=s_1,\ a_2=s_2 とおき,送信記号集合 A={a1,a2}A=\{a_1,a_2\} を構成する。

この送信記号集合 AA に属する記号 aia_i を,通信路行列

T=(1qqq1q),0q1T= \begin{pmatrix} 1-q & q\\ q & 1-q \end{pmatrix}, \qquad 0\leq q\leq 1

である通信路 CC を介して送信したとき,受信記号集合 B={b1,b2}B=\{b_1,b_2\} に属する記号 bjb_j が受信されるとする。通信路 CC の通信路行列 TTiijj 列要素は,aia_i を送信したときに bjb_j が受信される事象の発生確率を表す。

ただし,エントロピー関数 H(x), 0x1H(x),\ 0\leq x\leq 1

H(x)=xlog2x(1x)log2(1x)H(x)=-x\log_2 x-(1-x)\log_2(1-x)

であり,x=1/2x=1/2 のとき最大値 11 となる。

問1

SS の 2 次拡大情報源 S2S^2 のエントロピー H(S2)H(S^2) を求め,エントロピー関数 H(p)H(p) を用いて表しなさい。

問2

符号 CC では S2S^2 の情報源記号が符号化され,情報源記号 s1s1, s1s2s_1s_1,\ s_1s_2 は下表に示す c1, c2c_1,\ c_2 にそれぞれ符号化される。

情報源記号符号語
s1s1s_1s_1c1=000000c_1=000000
s1s2s_1s_2c2=001110c_2=001110
s2s1s_2s_1c3=c_3=
s2s2s_2s_2c4=c_4=

符号語 ck,cCc_k,c_\ell\in C のハミング距離 h(ck,c)h(c_k,c_\ell) を用いて

dmin(C)=minckc, ck,cCh(ck,c)d_{\min}(C)=\min_{c_k\neq c_\ell,\ c_k,c_\ell\in C}h(c_k,c_\ell)

と定義される。CC の最小ハミング距離 dmin(C)d_{\min}(C)33 であり,s2s1, s2s2s_2s_1,\ s_2s_2 それぞれの符号語 c3, c4c_3,\ c_4 のハミング重み w(c3),w(c4)w(c_3),w(c_4) がいずれも 33 であるとする。

このとき,符号語 c3, c4c_3,\ c_4 となりえる符号語の組み合わせの数を求めなさい。ただし,ある符号語 ca,cbc_a,c_bCC として採用できるとき,c3=ca, c4=cbc_3=c_a,\ c_4=c_b という割り当てと c3=cb, c4=cac_3=c_b,\ c_4=c_a という割り当ては,符号語 c3,c4c_3,c_4 となりえる符号語の組み合わせの数としては 1 つと数えることとする。

問3

H(B)H(B) を求め,エントロピー関数 H(α)H(\alpha) を用いて表しなさい。ただし,

α=p+q2pq\alpha=p+q-2pq

とする。

問4

事後確率 P(a1b1)P(a_1|b_1)P(a1b2)P(a_1|b_2) を求めなさい。

問5

H(AB)H(A|B) を求め,エントロピー関数を用いて表しなさい。そして,通信路 CC の通信路容量を求めなさい。

Kai

問1

S2S^2 の各情報源記号の発生確率は

P(s1s1)=p2P(s_1s_1)=p^2
P(s1s2)=P(s2s1)=p(1p)P(s_1s_2)=P(s_2s_1)=p(1-p)
P(s2s2)=(1p)2P(s_2s_2)=(1-p)^2

である。

したがって,

H(S2)=p2log2p22p(1p)log2{p(1p)}(1p)2log2(1p)2=2p2log2p2p(1p)log2p2p(1p)log2(1p)2(1p)2log2(1p)=2plog2p2(1p)log2(1p)=2H(p)\begin{aligned} H(S^2) &=-p^2\log_2 p^2 -2p(1-p)\log_2 \{p(1-p)\} -(1-p)^2\log_2 (1-p)^2 \\ &=-2p^2\log_2 p -2p(1-p)\log_2 p -2p(1-p)\log_2(1-p) -2(1-p)^2\log_2(1-p) \\ &=-2p\log_2 p-2(1-p)\log_2(1-p) \\ &=2H(p) \end{aligned}

よって,

H(S2)=2H(p)H(S^2)=2H(p)

である。

問2

c1=000000c_1=000000 であり,w(c3)=w(c4)=3w(c_3)=w(c_4)=3 であるため,

h(c1,c3)=h(c1,c4)=3h(c_1,c_3)=h(c_1,c_4)=3

となり,これは dmin(C)=3d_{\min}(C)=3 を満たす。

次に,c2=001110c_2=001110 とのハミング距離を考える。c2c_2 の 1 が立っている位置は第 3, 4, 5 ビットである。c3,c4c_3,c_4 はハミング重み 3 であるから,c2c_2 と 1 の位置が 2 個以上重なると,c2c_2 とのハミング距離が 2 以下となり,dmin(C)=3d_{\min}(C)=3 を満たさない。

したがって,c3,c4c_3,c_4 の候補は,c2c_2 の 1 の位置と高々 1 個だけ重なる符号語である。

この条件を満たす符号語は次の 10 個である。

c3,c4c_3,c_4 の符号語の候補符号語
ccand1,1c_{\mathrm{cand}1,1}111000111000
ccand1,2c_{\mathrm{cand}1,2}101001101001
ccand1,3c_{\mathrm{cand}1,3}011001011001
ccand2,1c_{\mathrm{cand}2,1}110100110100
ccand2,2c_{\mathrm{cand}2,2}100101100101
ccand2,3c_{\mathrm{cand}2,3}010101010101
ccand3,1c_{\mathrm{cand}3,1}110010110010
ccand3,2c_{\mathrm{cand}3,2}100011100011
ccand3,3c_{\mathrm{cand}3,3}010011010011
ccand4c_{\mathrm{cand}4}110001110001

ただし,ccand4=110001c_{\mathrm{cand}4}=110001 は,他の 9 個の候補のいずれとも 1 の位置が 2 個重なる。したがって,他の候補とのハミング距離は

h=2h=2

となり,c3,c4c_3,c_4 の組として同時に採用することはできない。

一方,残りの 9 個の候補については,1 つの符号語に対して,ハミング距離が 3 以上となる相手が 4 個存在する。

したがって,順序を区別しない組み合わせの数は

9×42=18\frac{9\times 4}{2}=18

である。

よって,求める符号語 c3,c4c_3,c_4 となりえる符号語の組み合わせの数は

1818

である。

問3

受信記号 b1b_1 の発生確率は

P(b1)=P(a1)P(b1a1)+P(a2)P(b1a2)=p(1q)+(1p)q=p+q2pq\begin{aligned} P(b_1) &=P(a_1)P(b_1|a_1)+P(a_2)P(b_1|a_2)\\ &=p(1-q)+(1-p)q\\ &=p+q-2pq \end{aligned}

である。

ここで

α=p+q2pq\alpha=p+q-2pq

より,

P(b1)=αP(b_1)=\alpha

また,

P(b2)=1αP(b_2)=1-\alpha

である。

したがって,

H(B)=αlog2α(1α)log2(1α)=H(α)\begin{aligned} H(B) &=-\alpha\log_2\alpha-(1-\alpha)\log_2(1-\alpha)\\ &=H(\alpha) \end{aligned}

よって,

H(B)=H(α)H(B)=H(\alpha)

である。

問4

ベイズの定理より,

P(a1b1)=P(a1)P(b1a1)P(b1)P(a_1|b_1) = \frac{P(a_1)P(b_1|a_1)}{P(b_1)}

である。ここで,

P(a1)=p,P(b1a1)=1q,P(b1)=αP(a_1)=p,\qquad P(b_1|a_1)=1-q,\qquad P(b_1)=\alpha

だから,

P(a1b1)=p(1q)αP(a_1|b_1)=\frac{p(1-q)}{\alpha}

である。

同様に,

P(a1b2)=P(a1)P(b2a1)P(b2)P(a_1|b_2) = \frac{P(a_1)P(b_2|a_1)}{P(b_2)}

であり,

P(b2a1)=q,P(b2)=1αP(b_2|a_1)=q,\qquad P(b_2)=1-\alpha

より,

P(a1b2)=pq1αP(a_1|b_2)=\frac{pq}{1-\alpha}

である。

問5

問4より,

P(a1b1)=p(1q)α,P(a2b1)=(1p)qαP(a_1|b_1)=\frac{p(1-q)}{\alpha}, \qquad P(a_2|b_1)=\frac{(1-p)q}{\alpha}

であり,

P(a1b2)=pq1α,P(a2b2)=(1p)(1q)1αP(a_1|b_2)=\frac{pq}{1-\alpha}, \qquad P(a_2|b_2)=\frac{(1-p)(1-q)}{1-\alpha}

である。また,

P(b1)=α,P(b2)=1αP(b_1)=\alpha,\qquad P(b_2)=1-\alpha

である。

したがって,条件付きエントロピー H(AB)H(A|B)

H(AB)= P(b1)H(Ab1)+P(b2)H(Ab2)= α(p(1q)αlog2p(1q)α(1p)qαlog2(1p)qα)+(1α)(pq1αlog2pq1α(1p)(1q)1αlog2(1p)(1q)1α)\begin{aligned} H(A|B) =&\ P(b_1)H(A|b_1)+P(b_2)H(A|b_2)\\ =&\ \alpha \left( -\frac{p(1-q)}{\alpha}\log_2\frac{p(1-q)}{\alpha} -\frac{(1-p)q}{\alpha}\log_2\frac{(1-p)q}{\alpha} \right)\\ &+(1-\alpha) \left( -\frac{pq}{1-\alpha}\log_2\frac{pq}{1-\alpha} -\frac{(1-p)(1-q)}{1-\alpha} \log_2\frac{(1-p)(1-q)}{1-\alpha} \right) \end{aligned}

となる。これを整理すると,

H(AB)=p(1q)log2p(1q)α(1p)qlog2(1p)qαpqlog2pq1α(1p)(1q)log2(1p)(1q)1α\begin{aligned} H(A|B) =&-p(1-q)\log_2\frac{p(1-q)}{\alpha} -(1-p)q\log_2\frac{(1-p)q}{\alpha}\\ &-pq\log_2\frac{pq}{1-\alpha} -(1-p)(1-q)\log_2 \frac{(1-p)(1-q)}{1-\alpha} \end{aligned}

である。

さらに対数を展開すると,

H(AB)=p(1q){log2p+log2(1q)log2α}(1p)q{log2(1p)+log2qlog2α}pq{log2p+log2qlog2(1α)}(1p)(1q){log2(1p)+log2(1q)log2(1α)}\begin{aligned} H(A|B) =&-p(1-q)\{\log_2 p+\log_2(1-q)-\log_2\alpha\}\\ &-(1-p)q\{\log_2(1-p)+\log_2 q-\log_2\alpha\}\\ &-pq\{\log_2 p+\log_2 q-\log_2(1-\alpha)\}\\ &-(1-p)(1-q) \{\log_2(1-p)+\log_2(1-q)-\log_2(1-\alpha)\} \end{aligned}

となる。各項をまとめると,

H(AB)=plog2p(1p)log2(1p)qlog2q(1q)log2(1q)+αlog2α+(1α)log2(1α)\begin{aligned} H(A|B) =&-p\log_2 p-(1-p)\log_2(1-p)\\ &-q\log_2 q-(1-q)\log_2(1-q)\\ &+\alpha\log_2\alpha+(1-\alpha)\log_2(1-\alpha) \end{aligned}

である。

したがって,

H(AB)=H(p)+H(q)H(α)H(A|B)=H(p)+H(q)-H(\alpha)

となる。

次に,相互情報量は

I(A;B)=H(A)H(AB)I(A;B)=H(A)-H(A|B)

であるから,

I(A;B)=H(p){H(p)+H(q)H(α)}=H(α)H(q)\begin{aligned} I(A;B) &=H(p)-\{H(p)+H(q)-H(\alpha)\}\\ &=H(\alpha)-H(q) \end{aligned}

である。

よって,通信路 CC の通信路容量は

Ccap=maxpI(A;B)=maxp{H(α)H(q)}C_{\mathrm{cap}} = \max_p I(A;B) = \max_p \{H(\alpha)-H(q)\}

である。

ここで,H(q)H(q)pp に依存しないので,H(α)H(\alpha) が最大となるように pp を選べばよい。エントロピー関数 H(α)H(\alpha)

α=12\alpha=\frac{1}{2}

のとき最大値 11 をとる。

また,

α=p+q2pq\alpha=p+q-2pq

であるから,q12q\neq \frac{1}{2} のとき,

p=12p=\frac{1}{2}

とすれば

α=12\alpha=\frac{1}{2}

となる。したがって,

Ccap=1H(q)C_{\mathrm{cap}}=1-H(q)

である。