跳到主要内容

京都大学 情報学研究科 知能情報学専攻 2024年8月実施 専門科目 S-3

Author

itsuitsuki

Description

大学公表の原題

Q.1

A discrete memoryless channel C consists of two discrete memoryless channels D and E, which are connected serially as shown in the following figure.

XDYEZX\xrightarrow{D}Y\xrightarrow{E}Z

The input alphabet of D is Σa={a1,a2,a3,a4}\Sigma_a = \{a_1, a_2, a_3, a_4\}. Both of the output alphabet of D and the input alphabet of E are Σb={b1,b2,b3,b4}\Sigma_b = \{b_1, b_2, b_3, b_4\}. The output alphabet of E is Σc={c1,c2}\Sigma_c = \{c_1, c_2\}. Let random variables X,YX, Y, and ZZ be respectively on Σa,Σb\Sigma_a, \Sigma_b, and Σc\Sigma_c. The channel transition matrix p(YX)p(Y|X) for D and the channel transition matrix q(ZY)q(Z|Y) for E are given as

p(YX)=(121200012120001212120012),andq(ZY)=(10100101).p(Y|X) = \begin{pmatrix} \frac{1}{2} & \frac{1}{2} & 0 & 0 \\ 0 & \frac{1}{2} & \frac{1}{2} & 0 \\ 0 & 0 & \frac{1}{2} & \frac{1}{2} \\ \frac{1}{2} & 0 & 0 & \frac{1}{2} \end{pmatrix}, \quad \text{and} \quad q(Z|Y) = \begin{pmatrix} 1 & 0 \\ 1 & 0 \\ 0 & 1 \\ 0 & 1 \end{pmatrix}.

Answer the following questions.

(1) Compute the channel capacity of D.

(2) Assume that XX follows the probability distribution r(ai)r(a_i) given below. Compute the mutual information I(X;Z)I(X; Z). You must show its derivation.

r(a1)=r(a3)=12,andr(a2)=r(a4)=0.r(a_1) = r(a_3) = \frac{1}{2}, \quad \text{and} \quad r(a_2) = r(a_4) = 0.

Q.2

We consider only the AND operation and the XOR (exclusive or) operation \oplus for the elements in {0,1}\{0, 1\}. We define a word as an element in {0,1}k (k1)\{0, 1\}^k \ (k \ge 1), each of which is represented as a row vector. Consider the liner codes generated with a matrix GG of kk rows and nn columns (k,n1)(k, n \ge 1) as x=wG\boldsymbol{x} = \boldsymbol{w}G, where w{0,1}k\boldsymbol{w} \in \{0, 1\}^k is a word and x{0,1}n\boldsymbol{x} \in \{0, 1\}^n is a codeword. Let C(G)={wGw{0,1}k}C(G) = \{\boldsymbol{w}G \mid \boldsymbol{w} \in \{0, 1\}^k\}. Answer the following questions.

(1) Show that xyC(G)\boldsymbol{x} \oplus \boldsymbol{y} \in C(G) holds for all x,yC(G)\boldsymbol{x}, \boldsymbol{y} \in C(G), where xy\boldsymbol{x} \oplus \boldsymbol{y} is the row vector obtained by element-wise XOR of two row vectors x\boldsymbol{x} and y\boldsymbol{y}.

(2) For the set C(G)C(G), prove that

minx,yC(G) and xyd(x,y)=minxC(G) and x0d(x,0),\min_{\boldsymbol{x}, \boldsymbol{y} \in C(G) \text{ and } \boldsymbol{x} \neq \boldsymbol{y}} d(\boldsymbol{x}, \boldsymbol{y}) = \min_{\boldsymbol{x} \in C(G) \text{ and } \boldsymbol{x} \neq \boldsymbol{0}} d(\boldsymbol{x}, \boldsymbol{0}),

where dd is the Hamming distance and 0=(0,0,,0)\boldsymbol{0} = (0, 0, \dots, 0).

(3) For the case that GG is given below, find a matrix FF of nn rows and nn columns such that G=GFG' = GF generates a systematic code. Moreover, by using the matrices FF and GG', compute the parity check matrix HH for C(G)C(G).

G=(0001101001101001101001101000).G = \begin{pmatrix} 0 & 0 & 0 & 1 & 1 & 0 & 1 \\ 0 & 0 & 1 & 1 & 0 & 1 & 0 \\ 0 & 1 & 1 & 0 & 1 & 0 & 0 \\ 1 & 1 & 0 & 1 & 0 & 0 & 0 \end{pmatrix}.

题目描述

  1. 考虑图示的串联系统,其中离散无记忆信道 D 的输入、输出字母表均为四元集合,信道 E 把四元输入映射为二元输出:

    XDYEZX\xrightarrow{D}Y\xrightarrow{E}Z

    D 与 E 的转移概率矩阵分别为

    p(YX)=(121200012120001212120012),q(ZY)=(10100101).p(Y\mid X)= \begin{pmatrix} \frac12&\frac12&0&0\\ 0&\frac12&\frac12&0\\ 0&0&\frac12&\frac12\\ \frac12&0&0&\frac12 \end{pmatrix}, \qquad q(Z\mid Y)= \begin{pmatrix} 1&0\\ 1&0\\ 0&1\\ 0&1 \end{pmatrix}.

    (1)求信道 D 的信道容量;(2)若输入符号 a1,a3a_1,a_3 的概率各为 12\frac12,而 a2,a4a_2,a_4 的概率为 0,请写出推导过程并计算串联后输入 XX 与最终输出 ZZ 的互信息 I(X;Z)I(X;Z)

  2. 在二元域上考虑由 k×nk\times n 生成矩阵 GG 定义的线性码。对 kk 位行向量信息字 w\boldsymbol{w},码字为 x=wG\boldsymbol{x}=\boldsymbol{w}G,所有码字组成 C(G)C(G)

    (1)证明任意 x,yC(G)\boldsymbol{x},\boldsymbol{y}\in C(G) 的按位异或 xy\boldsymbol{x}\oplus\boldsymbol{y} 仍属于 C(G)C(G)

    (2)证明

    minx,yC(G)xyd(x,y)=minxC(G)x0d(x,0),\min_{\substack{\boldsymbol{x},\boldsymbol{y}\in C(G)\\\boldsymbol{x}\ne\boldsymbol{y}}} d(\boldsymbol{x},\boldsymbol{y}) = \min_{\substack{\boldsymbol{x}\in C(G)\\\boldsymbol{x}\ne\boldsymbol{0}}} d(\boldsymbol{x},\boldsymbol{0}),

    其中 dd 为 Hamming 距离,0=(0,0,,0)\boldsymbol{0}=(0,0,\ldots,0)

    (3)当

    G=(0001101001101001101001101000)G= \begin{pmatrix} 0&0&0&1&1&0&1\\ 0&0&1&1&0&1&0\\ 0&1&1&0&1&0&0\\ 1&1&0&1&0&0&0 \end{pmatrix}

    时,求一个 n×nn\times n 矩阵 FF,使 G=GFG'=GF 生成系统码;再利用 FFGG' 计算原码 C(G)C(G) 的校验矩阵 HH

Kai

All logarithms in Q.1 have base 22, and all matrix operations in Q.2 are over F2\mathbb F_2.

Q.1

(1)

Every row of p(YX)p(Y\mid X) has entropy 11, so H(YX)=1H(Y\mid X)=1 for any input distribution. Since YY has four possible values, H(Y)2H(Y)\le2. Uniform input makes the output uniform because each column sums to one, attaining the bound. Therefore

CD=maxI(X;Y)=21=1 bit.C_D=\max I(X;Y)=2-1=1\text{ bit}.

(2)

Multiplying the two channel matrices gives

p(ZX)=(101/21/2011/21/2).p(Z\mid X)= \begin{pmatrix}1&0\\1/2&1/2\\0&1\\1/2&1/2\end{pmatrix}.

Only a1,a3a_1,a_3 occur. Their outputs are deterministically c1,c2c_1,c_2, respectively, each with probability 1/21/2. Thus H(Z)=1H(Z)=1, H(ZX)=0H(Z\mid X)=0, and I(X;Z)=1I(X;Z)=1 bit.

Q.2

(1)

Write x=uG\boldsymbol x=\boldsymbol uG and y=vG\boldsymbol y=\boldsymbol vG. Distributivity over F2\mathbb F_2 gives

xy=(uv)GC(G).\boldsymbol x\oplus\boldsymbol y=(\boldsymbol u\oplus\boldsymbol v)G\in C(G).

(2)

For distinct codewords, xy\boldsymbol x\oplus\boldsymbol y is a nonzero codeword and

d(x,y)=d(xy,0).d(\boldsymbol x,\boldsymbol y) =d(\boldsymbol x\oplus\boldsymbol y,\boldsymbol0).

Conversely, every nonzero codeword occurs in the pair (x,0)(\boldsymbol x,\boldsymbol0). The two sets of distances therefore coincide and have equal minima. This statement assumes C(G){0}C(G)\ne\{\boldsymbol0\}; for the zero code both sets are empty (or both minima are ++\infty under that convention).

(3)

Partition G=[A B]G=[A\ B], where AA consists of the first four columns. It is invertible over F2\mathbb F_2. Choose

F=diag(A1,I3)=(0111000111000011000001000000000010000000100000001).F=\operatorname{diag}(A^{-1},I_3)= \begin{pmatrix} 0&1&1&1&0&0&0\\ 1&1&1&0&0&0&0\\ 1&1&0&0&0&0&0\\ 1&0&0&0&0&0&0\\ 0&0&0&0&1&0&0\\ 0&0&0&0&0&1&0\\ 0&0&0&0&0&0&1 \end{pmatrix}.

Then

G=GF=[I4 B]=(1000101010001000101000001000).G'=GF=[I_4\ B]= \begin{pmatrix} 1&0&0&0&1&0&1\\ 0&1&0&0&0&1&0\\ 0&0&1&0&1&0&0\\ 0&0&0&1&0&0&0 \end{pmatrix}.

A parity-check matrix for C(G)C(G') is H=[BT I3]H'=[B^T\ I_3]. Since a row codeword transforms as x=xF\boldsymbol x'=\boldsymbol xF, the original code has parity-check matrix H=HFTH=H'F^T:

H=(101110011100100111001).\boxed{H= \begin{pmatrix} 1&0&1&1&1&0&0\\ 1&1&1&0&0&1&0\\ 0&1&1&1&0&0&1 \end{pmatrix}}.

Indeed GHT=0GH^T=0 and rankH=3\operatorname{rank}H=3, so the kernel consists of exactly the 242^4 original codewords.