跳到主要内容

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

Author

Yu

Description

【問 1】

情報源アルファベットが X={a,b,c,d}\mathcal{X}=\{a,b,c,d\} の離散無記憶情報源があり, 記号 xXx \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} 情報源 X1X2X_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 について, v1||v||_1 でその要素の絶対値の総和を表す。行ベクトル vv について, vP1v1||vP||_1 \le ||v||_1 を示せ。

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

题目描述

【问题 1】某离散无记忆信源的信源字母表为 X={a,b,c,d}\mathcal{X}=\{a,b,c,d\},符号 xXx\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 信源 X1X2X_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=jXt=i)p_{ij}=P(X_{t+1}=j\mid X_t=i)。回答:

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

考点

  • 离散信源熵:由含参数的符号分布计算熵,并分析其随 α\alpha 的变化。
  • 二元哈夫曼编码:针对概率次序随参数区间改变的情况分别建树、求平均码长并画分段函数。
  • 马尔可夫信源稳态分布:理解随机矩阵特征值 11 的左右特征向量及稳态概率的意义。
  • 随机矩阵的谱性质:借助 1\ell_1 范数的不增性证明转移矩阵全部特征值的模不超过 11

Kai

【問 1】

(1)

H(S)=1α2log2α21α2log21α2H(S) = 1 - \frac{\alpha}{2}\log_2\frac{\alpha}{2} - \frac{1 - \alpha}{2} \log_2\frac{1-\alpha}{2}

(2)

(3)

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

(4)

(5)

【問 2】

(1)

T=λEP=λ0.250.250.50.4λ0.6000.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.750.250.50.40.4000.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 はマルコフ情報源の定常分布

(3)

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

(4)

vP1=v1p11+v2p21+v3p31+v1p12+v2p22+v3p32+v1p13+v2p23+v3p33(v1p11+v2p21+v3p31)+(v1p12+v2p22+v3p32)+(v1p13+v2p23+v3p33)=v1j=13p1j+v2j=13p2j+v3j=13p3j=v1+v2+v3=v1\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)

Pv1=p11v1+p12v1+p13v1+p21v2+p22v2+p23v2+p31v3+p32v3+p33v3(p11v1+p12v1+p13v1)+(p21v2+p22v2+p23v2)+(p31v3+p32v3+p33v3)=v1j=13p1j+v2j=13p2j+v3j=13p3j=v1+v2+v3=v1Pv=λvとすると,Pv1=λv1=λv1+λv2+λv3=λv1Pv1v1λv1v1λ1\begin{aligned} ||Pv||_1 &= |p_{11}v_1 + p_{12}v_1 + p_{13}v_1| + |p_{21}v_2 + p_{22}v_2 + p_{23}v_2| + |p_{31}v_3 + p_{32}v_3 + p_{33}v_3| \\ &\le (|p_{11}v_1 | + |p_{12}v_1| + |p_{13}v_1|) + (|p_{21}v_2| + |p_{22}v_2| + |p_{23}v_2|) + (|p_{31}v_3| + |p_{32}v_3| + |p_{33}v_3|) \\ &= |v_1|\sum_{j = 1}^3 p_{1j} + |v_2|\sum_{j = 1}^3 p_{2j} + |v_3|\sum_{j = 1}^3 p_{3j} \\ &= |v_1| + |v_2| + |v_3| \\ &= ||v||_1 \\ Pv &= \lambda v \text{とすると,}||Pv||_1 = ||\lambda v||_1 = |\lambda v_1| + |\lambda v_2| + |\lambda v_3| = |\lambda| \cdot ||v||_1 \\ ||Pv||_1 &\le ||v||_1 \Rightarrow |\lambda| \cdot ||v||_1 \le ||v||_1 \Rightarrow |\lambda| \le 1 \end{aligned}