跳到主要内容

東京工業大学 工学院 情報通信系 2018年8月実施 概率统计

Author

思齐塾, 祭音Myyura

Description

S2.

XXYY をそれぞれ MM 個の実数値 {x1,x2,,xM}\{x_1, x_2, \dots, x_M\}NN 個の実数値 {y1,y2,,yN}\{y_1, y_2, \dots, y_N\} のいずれかの値を取る確率変数とし、その同時確率分布を PXYP_{XY} によって、またその周辺分布を PXP_XPYP_Y によって表す。ただし、 M2,N2M \ge 2, N \ge 2 とし、 PXY(xi,yj)>0(i=1,2,,M,j=1,2,,N)P_{XY}(x_i, y_j) > 0 (i = 1, 2, \dots, M, j = 1, 2, \dots, N) とする。

  1. エントロピー H(X)H(X) 、条件付きエントロピー H(YX)H(Y|X) 、同時エントロピー H(X,Y)H(X, Y) 、相互情報量 I(X;Y)I(X; Y) を確率分布 PXY,PXP_{XY}, P_X ならびに PYP_Y を用いて書け。ただし、エントロピーの単位はビットとする。

  2. 2 つの恒等式

H(X,Y)=H(X)+H(YX)H(X, Y) = H(X) + H(Y|X)
I(X;Y)=H(Y)H(YX)I(X; Y) = H(Y) - H(Y|X)

が成り立つことを示せ。

  1. 実数値関数 ff は、任意の λ[0,1]\lambda \in [0, 1] と任意の正の実数 x,yx, y について
λf(x)+(1λ)f(y)f(λx+(1λ)y)\lambda f(x) + (1 - \lambda)f(y) \le f(\lambda x + (1 - \lambda)y)

が成り立つとき、上に凸の関数であるという。関数 ff が上に凸の関数であるとき、任意の MM に対して

i=1MPX(xi)f(zi)f(i=1MPX(xi)zi)\sum_{i=1}^M P_X(x_i)f(z_i) \le f\left( \sum_{i=1}^M P_X(x_i)z_i \right)

が成り立つことを数学的帰納法を用いて証明せよ。ただし、 z1,z2,,zMz_1, z_2, \dots, z_M は任意の正の実数とする。

  1. 3 つの不等式
0H(X)0 \le H(X)
H(X)log2MH(X) \le \log_2 M
H(YX)H(Y)H(Y|X) \le H(Y)

が成り立つことを証明せよ。 log2x\log_2 x が上に凸の関数であることに注意せよ。

题目描述

随机变量 XXYY 分别取 MM 个实数值 {x1,,xM}\{x_1,\ldots,x_M\}NN 个实数值 {y1,,yN}\{y_1,\ldots,y_N\},其中 M,N2M,N\geq2。用 PXYP_{XY} 表示联合分布,PX,PYP_X,P_Y 表示边缘分布,并假设

PXY(xi,yj)>0(i=1,,M; j=1,,N).P_{XY}(x_i,y_j)>0 \quad (i=1,\ldots,M;\ j=1,\ldots,N).
  1. PXY,PX,PYP_{XY},P_X,P_Y 写出熵 H(X)H(X)、条件熵 H(YX)H(Y\mid X)、联合熵 H(X,Y)H(X,Y) 和互信息 I(X;Y)I(X;Y);熵的单位取 bit。
  2. 证明恒等式
H(X,Y)=H(X)+H(YX),H(X,Y)=H(X)+H(Y\mid X),
I(X;Y)=H(Y)H(YX).I(X;Y)=H(Y)-H(Y\mid X).
  1. 若实函数 ff 对每个 λ[0,1]\lambda\in[0,1] 及任意正实数 x,yx,y 都满足
λf(x)+(1λ)f(y)f(λx+(1λ)y),\lambda f(x)+(1-\lambda)f(y) \leq f\bigl(\lambda x+(1-\lambda)y\bigr),

则称 ff 为上凸函数。假设 ff 上凸,用数学归纳法证明对任意 MM 及任意正实数 z1,,zMz_1,\ldots,z_M

i=1MPX(xi)f(zi)f ⁣(i=1MPX(xi)zi).\sum_{i=1}^MP_X(x_i)f(z_i) \leq f\!\left(\sum_{i=1}^MP_X(x_i)z_i\right).
  1. 注意到 log2x\log_2x 是上凸函数,证明
0H(X),H(X)log2M,H(YX)H(Y).0\leq H(X),\qquad H(X)\leq\log_2M,\qquad H(Y\mid X)\leq H(Y).

Kai

S2

以下, pij=PXY(xi,yj)p_{ij}=P_{XY}(x_i,y_j) , pi=PX(xi)p_i=P_X(x_i) , qj=PY(yj)q_j=P_Y(y_j) と略記する。仮定よりこれらは正である。

1)

ビットを単位とする各量は

H(X)=i=1Mpilog2pi,\boxed{H(X)=-\sum_{i=1}^{M}p_i\log_2p_i},
H(YX)=i=1Mj=1Npijlog2pijpi,\boxed{H(Y\mid X)=-\sum_{i=1}^{M}\sum_{j=1}^{N}p_{ij} \log_2\frac{p_{ij}}{p_i}},
H(X,Y)=i=1Mj=1Npijlog2pij,\boxed{H(X,Y)=-\sum_{i=1}^{M}\sum_{j=1}^{N}p_{ij}\log_2p_{ij}},
I(X;Y)=i=1Mj=1Npijlog2pijpiqj.\boxed{I(X;Y)=\sum_{i=1}^{M}\sum_{j=1}^{N}p_{ij} \log_2\frac{p_{ij}}{p_iq_j}}.

2)

log(pij/pi)=logpijlogpi\log(p_{ij}/p_i)=\log p_{ij}-\log p_i と周辺化 jpij=pi\sum_jp_{ij}=p_i を使えば

H(X)+H(YX)=ipilog2pii,jpij(log2pijlog2pi)=i,jpijlog2pij=H(X,Y).\begin{aligned} H(X)+H(Y\mid X) &=-\sum_i p_i\log_2p_i -\sum_{i,j}p_{ij}(\log_2p_{ij}-\log_2p_i)\\ &=-\sum_{i,j}p_{ij}\log_2p_{ij}=H(X,Y). \end{aligned}

また

H(Y)H(YX)=jqjlog2qj+i,jpijlog2pijpi=i,jpijlog2pijpiqj=I(X;Y).\begin{aligned} H(Y)-H(Y\mid X) &=-\sum_jq_j\log_2q_j +\sum_{i,j}p_{ij}\log_2\frac{p_{ij}}{p_i}\\ &=\sum_{i,j}p_{ij}\log_2\frac{p_{ij}}{p_iq_j} =I(X;Y). \end{aligned}

よって二つの恒等式が示された。

3)

重みを pi=PX(xi)p_i=P_X(x_i) と書く。 M=1M=1 は等号であり, M=2M=2 は上に凸の定義そのものである。ある MM で成立すると仮定し, M+1M+1 個の正の重みについて

q=i=1Mpi=1pM+1(0<q<1),zˉ=i=1Mpiqziq=\sum_{i=1}^{M}p_i=1-p_{M+1}\quad(0<q<1),\qquad \bar z=\sum_{i=1}^{M}\frac{p_i}{q}z_i

とおく。帰納法の仮定と二点に対する凸性から

i=1M+1pif(zi)=qi=1Mpiqf(zi)+pM+1f(zM+1)qf(zˉ)+pM+1f(zM+1)f(qzˉ+pM+1zM+1)=f(i=1M+1pizi).\begin{aligned} \sum_{i=1}^{M+1}p_if(z_i) &=q\sum_{i=1}^{M}\frac{p_i}{q}f(z_i)+p_{M+1}f(z_{M+1})\\ &\le qf(\bar z)+p_{M+1}f(z_{M+1})\\ &\le f(q\bar z+p_{M+1}z_{M+1})\\ &=f\left(\sum_{i=1}^{M+1}p_iz_i\right). \end{aligned}

したがって数学的帰納法により任意の MM で主張が成り立つ。

4)

pi(0,1]p_i\in(0,1] なので log2pi0-\log_2p_i\ge0 であり,直ちに

0H(X)\boxed{0\le H(X)}

を得る。

次に 3) を上に凸な log2\log_2zi=1/piz_i=1/p_i に適用すると

H(X)=ipilog21pilog2(ipi1pi)=log2M.H(X)=\sum_ip_i\log_2\frac1{p_i} \le\log_2\left(\sum_ip_i\frac1{p_i}\right) =\boxed{\log_2M}.

最後に相互情報量の非負性を同じ不等式で示す。 log2\log_2 の上への凸性より

I(X;Y)=i,jpijlog2piqjpijlog2(i,jpijpiqjpij)=log2(i,jpiqj)=0.\begin{aligned} -I(X;Y) &=\sum_{i,j}p_{ij}\log_2\frac{p_iq_j}{p_{ij}}\\ &\le\log_2\left(\sum_{i,j}p_{ij}\frac{p_iq_j}{p_{ij}}\right) =\log_2\left(\sum_{i,j}p_iq_j\right)=0. \end{aligned}

ゆえに I(X;Y)0I(X;Y)\ge0 。2) の恒等式 I(X;Y)=H(Y)H(YX)I(X;Y)=H(Y)-H(Y\mid X) から

H(YX)H(Y)\boxed{H(Y\mid X)\le H(Y)}

が従う。