跳到主要内容

東京大学 情報理工学系研究科 電子情報学専攻 2012年8月実施 専門 第5問

Author

祭音Myyura (co-authored with GPT 5.6 SOL)

Description

一定の時間間隔で、赤(R)、緑(G)、青(B)のいずれかのランプの光る装置がある。その光り方には、以下のような特徴がある。

  • R の直後には、R と G のどちらかが、それぞれ 0.80.80.20.2 の確率で光る。
  • G の直後には、G と B のどちらかが、それぞれ 0.90.90.10.1 の確率で光る。
  • B の直後には、B と R のどちらかが、それぞれ 0.90.90.10.1 の確率で光る。

この装置の光の点灯を情報源 SS として、以下の問いに答えよ。

(1) 情報源 SS は、一次のマルコフ情報源である。その状態遷移図を示せ。

(2) R, G, B それぞれが点灯する定常状態確率を求めよ。

(3) 情報源 SS からの点灯系列を効率よく符号化して伝送したい。点灯を 22 つずつまとめて 0,10,1 からなる 22 元符号を割り当てる場合に、最も効率の良い符号を設計せよ。また、その時の 11 点灯あたりの平均符号長を求めよ。

(4) 同じ色の点灯が継続して続く傾向があることを考慮し、以下に示す 99 つの点灯のまとまりに対して 0,10,1 からなる 22 元符号の割り当てを行う場合、最も効率の良い符号を設計せよ。また、その時の 11 点灯あたりの平均符号長を求めよ。

RG, RRG, RRR, GB, GGB, GGG, BR, BBR, BBB\mathrm{RG,\ RRG,\ RRR,\ GB,\ GGB,\ GGG,\ BR,\ BBR,\ BBB}

(5) この情報源 SS の平均符号長はどこまで短くすることができるか。その最小値を求めよ。

(なお、計算に当たっては、log23=1.58, log25=2.32\log_2 3=1.58,\ \log_2 5=2.32 を必要に応じて用いよ。)

Kai

(1)

(2)

定常確率を (πR,πG,πB)(\pi_R,\pi_G,\pi_B) とすると、

0.2πR=0.1πG=0.1πB,πR+πG+πB=1.0.2\pi_R=0.1\pi_G=0.1\pi_B,\qquad \pi_R+\pi_G+\pi_B=1.

したがって、

(πR,πG,πB)=(15,25,25).\boxed{(\pi_R,\pi_G,\pi_B)=\left(\frac15,\frac25,\frac25\right)}.

(3)

P(Xn=i,Xn+1=j)=πiPijP(X_n=i,X_{n+1}=j)=\pi_iP_{ij} に対してハフマン符号を構成する。

点灯の組確率符号
GG0.360.360
BB0.360.3610
RR0.160.16110
RG0.040.041110
GB0.040.0411110
BR0.040.0411111

他の組の確率は 00 である。11 点灯あたりの平均符号長は

L2=0.36+2(0.36)+3(0.16)+(4+5+5)(0.04)2=1.06 bit/点灯.\boxed{\overline L_2= \frac{0.36+2(0.36)+3(0.16)+(4+5+5)(0.04)}2 =1.06\ \text{bit/点灯}}.

(4)

系列を指定されたまとまりに順次、重なりなく区切る。まとまりの先頭色の定常分布を μ\mu とする。先頭色から末尾色への遷移行列 DD と、通常の一段遷移行列 PP

D=(.64.3600.81.19.190.81),P=(.8.200.9.1.10.9).D=\begin{pmatrix}.64&.36&0\\0&.81&.19\\.19&0&.81\end{pmatrix},\qquad P=\begin{pmatrix}.8&.2&0\\0&.9&.1\\.1&0&.9\end{pmatrix}.

次のまとまりの先頭色への遷移は DPDP なので、μ=μDP\mu=\mu DPiμi=1\sum_i\mu_i=1 より

μ=(μR,μG,μB)=162277(12773,24772,24732).\boxed{\mu=(\mu_R,\mu_G,\mu_B) =\frac1{62277}(12773,24772,24732)}.

この分布に対するハフマン符号の一例は次のとおり。

まとまり確率符号
RG.2μR.2\mu_R0101
RRG.16μR.16\mu_R0000
RRR.64μR.64\mu_R011
GB.1μG.1\mu_G0100
GGB.09μG.09\mu_G0010
GGG.81μG.81\mu_G11
BR.1μB.1\mu_B0011
BBR.09μB.09\mu_B0001
BBB.81μB.81\mu_B10

まとまりあたりの平均符号長と平均点灯数はそれぞれ

E[]=3.36μR+2.38μG+2.38μB=803684311385,E[N]=2.8μR+2.9μG+2.9μB=17932662277.\begin{aligned} \mathrm E[\ell]&=3.36\mu_R+2.38\mu_G+2.38\mu_B =\frac{803684}{311385},\\ \mathrm E[N]&=2.8\mu_R+2.9\mu_G+2.9\mu_B =\frac{179326}{62277}. \end{aligned}

したがって、

Lv=E[]E[N]=57406640450.8963 bit/点灯.\boxed{\overline L_{\mathrm{v}}=\frac{\mathrm E[\ell]}{\mathrm E[N]} =\frac{57406}{64045}\simeq0.8963\ \text{bit/点灯}}.

(5)

下限は情報源のエントロピー率である。h2(p)=plog2p(1p)log2(1p)h_2(p)=-p\log_2p-(1-p)\log_2(1-p) とおき、指定の対数近似を用いると

H(S)=15h2(0.2)+45h2(0.1)=log25+12253625log230.5248 bit/点灯.\begin{aligned} H(S)&=\frac15h_2(0.2)+\frac45h_2(0.1)\\ &=\log_2 5+\frac{12}{25}-\frac{36}{25}\log_2 3\\ &\simeq\boxed{0.5248\ \text{bit/点灯}}. \end{aligned}

十分長いブロックを符号化すれば、この値に任意に近づけられる。