跳到主要内容

九州大学 システム情報科学府 情報理工学専攻 2018年8月実施 情報理論

Author​

Yu, 祭音Myyura

Description​

【問 1】​

情報源アルファベットが X={a,b,c,d}\mathcal{X}=\{a,b,c,d\} の離散無記憶情報源があり, 記号 x∈Xx \in \mathcal{X} の出現確率 p(x)p(x) が以下の表で与えられている。ただし, 0<α≤120 <\alpha \le \frac{1}{2} とする。

xxaabbccdd
p(x)p(x)α2\frac{\alpha}{2}1−α2\frac{1-\alpha}{2}14\frac{1}{4}14\frac{1}{4}

以下の問いに答えよ。

(1) この情報源のエントロピーを α\alpha の関数として求めよ。

(2) 0<α<140 < \alpha < \frac{1}{4} のときのこの情報源に対する 22 元ハフマン符号を 11 つ求めよ。

(3) (2)の符号の平均符号長 (符号長の期待値)を α\alpha の関数として求めよ。

(4) 14≤α≤12\frac{1}{4} \le \alpha \le \frac{1}{2} のときのこの情報源に対する 22 元ハフマン符号の平均符号長を求めよ。

(5) この情報源に対する 22 元ハフマン符号の平均符号長を L(α)L(\alpha) とおく。 α∈(0,12]\alpha \in (0,\frac{1}{2}] に対する L(α)L(\alpha) のガラフをかけ。

【問 2】​

アルファベット {1,2,3}\{1,2,3\} 上の Markov\text{Markov} 情報源 X1X2⋯X_1X_2\cdots の遷移確率行列が

P=(pij)=(0.250.250.50.40.6000.20.8)P = (p_{ij}) = \begin{pmatrix} 0.25 & 0.25 & 0.5 \\ 0.4 & 0.6 & 0 \\ 0 & 0.2 & 0.8 \end{pmatrix}

で与えられているとき, 次の各問に答えよ。ただし, pijp_{ij} は Xt=iX_t = i のもとで Xt+1=jX_{t + 1} = j となる条件付き確率を表す。

(1) PP の固有値の一つが 11 であることを示せ。

(2) 行ベクトル vv が固有値 11 に対応する PP の左固有ベクトルであり, その要素の総和が 11 であるとする。vv の持つ意味を述べよ。ただし, 固有値 λ\lambda に対応する PP の左固有ベクトルとは, vP=λvvP = \lambda v を満たすゼロでない行ベクトル vv のことである。

(3) 固有値 11 に対応する PP の右固有ベクトルを求めよ。ただし, 固有値 λ\lambda に対応する PP の右固有ベクトルとは, Pu=λuPu = \lambda u を満たすゼロでない列ベクトル uu のことである。

(4) 任意のベクトル vv について, ∣∣v∣∣1||v||_1 でその要素の絶対値の総和を表す。行ベクトル vv について, ∣∣vP∣∣1≤∣∣v∣∣1||vP||_1 \le ||v||_1 を示せ。

(5) PP の任意の固有値の絶対値は 11 以下であることを示せ。

题目描述​

【问题 1】某离散无记忆信源的信源字母表为 X={a,b,c,d}\mathcal{X}=\{a,b,c,d\},符号 x∈Xx\in\mathcal{X} 的出现概率如下,其中 0<α≤120<\alpha\le\frac12:

xxaabbccdd
p(x)p(x)α2\frac{\alpha}{2}1−α2\frac{1-\alpha}{2}14\frac1414\frac14

回答:

  1. 求该信源的熵,将结果表示为 α\alpha 的函数。
  2. 当 0<α<140<\alpha<\frac14 时,给出该信源的一种二元 Huffman 编码。
  3. 求 (2) 中编码的平均码长(码长的期望),表示为 α\alpha 的函数。
  4. 当 14≤α≤12\frac14\le\alpha\le\frac12 时,求该信源二元 Huffman 编码的平均码长。
  5. 记该信源二元 Huffman 编码的平均码长为 L(α)L(\alpha),画出 α∈(0,12]\alpha\in(0,\frac12] 上的 L(α)L(\alpha) 图像。

【问题 2】字母表为 {1,2,3}\{1,2,3\} 的 Markov 信源 X1X2⋯X_1X_2\cdots 的转移概率矩阵为

P=(pij)=(0.250.250.50.40.6000.20.8),P=(p_{ij})= \begin{pmatrix} 0.25&0.25&0.5\\ 0.4&0.6&0\\ 0&0.2&0.8 \end{pmatrix},

其中 pij=P(Xt+1=j∣Xt=i)p_{ij}=P(X_{t+1}=j\mid X_t=i)。回答:

  1. 证明 11 是 PP 的一个特征值。
  2. 设行向量 vv 是 PP 对应特征值 11 的左特征向量,且各分量之和为 11,说明 vv 的概率意义。这里左特征向量满足 vP=λvvP=\lambda v 且非零。
  3. 求 PP 对应特征值 11 的右特征向量;右特征向量为满足 Pu=λuPu=\lambda u 的非零列向量 uu。
  4. 对任意向量 vv,以 ∥v∥1\|v\|_1 表示其各分量绝对值之和。证明对行向量 vv 有 ∥vP∥1≤∥v∥1\|vP\|_1\le\|v\|_1。
  5. 证明 PP 的任意特征值的绝对值均不超过 11。

Kai​

【問 1】​

(1)​

H(S)=1−α2log⁡2α2−1−α2log⁡21−α2H(S) = 1 - \frac{\alpha}{2}\log_2\frac{\alpha}{2} - \frac{1 - \alpha}{2} \log_2\frac{1-\alpha}{2}

(2)​

例えば

b↦0,d↦10,a↦110,c↦111b\mapsto0,\qquad d\mapsto10,\qquad a\mapsto110,\qquad c\mapsto111

とすればよい。

(3)​

Lˉ=1−α2+14⋅2+(14+α2)⋅3=74+α(ビット)\bar{L} = \frac{1 - \alpha}{2} + \frac{1}{4} \cdot 2 + (\frac{1}{4} + \frac{\alpha}{2}) \cdot 3 = \frac{7}{4} + \alpha \quad(\text{ビット})

(4)​

全記号の符号長を 22 にできるので、平均符号長は 22 ビットである。

(5)​

従って

L(α)={74+α(0<α<14),2(14≤α≤12).L(\alpha)= \begin{cases} \dfrac74+\alpha & (0<\alpha<\dfrac14),\\ 2 & (\dfrac14\leq\alpha\leq\dfrac12). \end{cases}

【問 2】​

(1)​

∣T∣=∣λE−P∣=∣λ−0.25−0.25−0.5−0.4λ−0.600−0.2λ−0.8∣|T| = |\lambda E - P| = \begin{vmatrix} \lambda - 0.25 & -0.25 & -0.5\\ -0.4 & \lambda - 0.6 & 0 \\ 0 & -0.2 & \lambda - 0.8 \end{vmatrix}
λ=1とすると, ∣T∣=∣0.75−0.25−0.5−0.40.400−0.20.2∣=0\lambda = 1 \text{とすると, } |T| = \begin{vmatrix} 0.75 & -0.25 & -0.5 \\ -0.4 & 0.4 & 0 \\ 0 & -0.2 & 0.2 \end{vmatrix} = 0

(2)​

vv はマルコフ情報源の定常分布である。この行列では正規化した左固有ベクトルは v=(8,15,20)/43v=(8,15,20)/43 であり、vP=vvP=v を満たす。初期分布が vv なら、すべての時刻で同じ分布が保たれる。

(3)​

Tu=0⇒u=[111]Tu = 0 \Rightarrow u = \begin{bmatrix} 1 \\ 1 \\ 1 \\ \end{bmatrix}

(4)​

∣∣vP∣∣1=∣v1p11+v2p21+v3p31∣+∣v1p12+v2p22+v3p32∣+∣v1p13+v2p23+v3p33∣≤(∣v1p11∣+∣v2p21∣+∣v3p31∣)+(∣v1p12∣+∣v2p22∣+∣v3p32∣)+(∣v1p13∣+∣v2p23∣+∣v3p33∣)=∣v1∣∑j=13p1j+∣v2∣∑j=13p2j+∣v3∣∑j=13p3j=∣v1∣+∣v2∣+∣v3∣=∣∣v∣∣1\begin{aligned} ||vP||_1 &= |v_1p_{11} + v_2p_{21} + v_3p_{31}| + |v_1p_{12} + v_2p_{22} + v_3p_{32}| + |v_1p_{13} + v_2p_{23} + v_3p_{33}| \\ &\le (|v_1p_{11}| + |v_2p_{21}| + |v_3p_{31}|) + (|v_1p_{12}| + |v_2p_{22}| + |v_3p_{32}|) + (|v_1p_{13}| + |v_2p_{23}| + |v_3p_{33}|) \\ &= |v_1|\sum_{j = 1}^3p_{1j} + |v_2|\sum_{j = 1}^3p_{2j} +|v_3|\sum_{j = 1}^3p_{3j} \\ &= |v_1| + |v_2| + |v_3| \\ &= ||v||_1 \end{aligned}

(5)​

λ\lambda を PP の任意の固有値とする。PP と PTP^T の固有値は 一致するので、vP=λvvP=\lambda v を満たす非零の左固有ベクトル vv が 存在する。(4) は複素ベクトルに対しても同様に成り立つから、

∣λ∣ ∥v∥1=∥λv∥1=∥vP∥1≤∥v∥1.|\lambda|\,\|v\|_1=\|\lambda v\|_1=\|vP\|_1 \leq\|v\|_1.

∥v∥1>0\|v\|_1>0 より ∣λ∣≤1|\lambda|\leq1 である。