跳到主要内容

東京工業大学 工学院 情報通信系 2017年8月実施 S2 エントロピーとJensenの不等式

Author​

思齐塾, 祭音Myyura (co-authored with GPT 5.6 SOL)

Description​

S2.

XX と YY をそれぞれ 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_X と PYP_Y によって表す。ただし、 M≥2,N≥2M \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(Y∣X)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(Y∣X)H(X, Y) = H(X) + H(Y|X)
I(X;Y)=H(Y)−H(Y∣X)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 つの不等式
0≤H(X)0 \le H(X)
H(X)≤log⁡2MH(X) \le \log_2 M
H(Y∣X)≤H(Y)H(Y|X) \le H(Y)

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

题目描述​

随机变量 XX、YY 分别取 MM 个实数值 {x1,…,xM}\{x_1,\ldots,x_M\} 和 NN 个实数值 {y1,…,yN}\{y_1,\ldots,y_N\},其中 M,N≥2M,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(Y∣X)H(Y\mid X)、联合熵 H(X,Y)H(X,Y) 和互信息 I(X;Y)I(X;Y);熵的单位取 bit。
  2. 证明恒等式
H(X,Y)=H(X)+H(Y∣X),H(X,Y)=H(X)+H(Y\mid X),
I(X;Y)=H(Y)−H(Y∣X).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. 注意到 log⁡2x\log_2x 是上凸函数,证明
0≤H(X),H(X)≤log⁡2M,H(Y∣X)≤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=1Mpilog⁡2pi,\boxed{H(X)=-\sum_{i=1}^{M}p_i\log_2p_i},
H(Y∣X)=−∑i=1M∑j=1Npijlog⁡2pijpi,\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=1M∑j=1Npijlog⁡2pij,\boxed{H(X,Y)=-\sum_{i=1}^{M}\sum_{j=1}^{N}p_{ij}\log_2p_{ij}},
I(X;Y)=∑i=1M∑j=1Npijlog⁡2pijpiqj.\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)=log⁡pij−log⁡pi\log(p_{ij}/p_i)=\log p_{ij}-\log p_i と周辺化 ∑jpij=pi\sum_jp_{ij}=p_i を使えば

H(X)+H(Y∣X)=−∑ipilog⁡2pi−∑i,jpij(log⁡2pij−log⁡2pi)=−∑i,jpijlog⁡2pij=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(Y∣X)=−∑jqjlog⁡2qj+∑i,jpijlog⁡2pijpi=∑i,jpijlog⁡2pijpiqj=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=1−pM+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)=q∑i=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] なので −log⁡2pi≥0-\log_2p_i\ge0 であり,直ちに

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

を得る。

次に 3) を上に凸な log⁡2\log_2 と zi=1/piz_i=1/p_i に適用すると

H(X)=∑ipilog⁡21pi≤log⁡2(∑ipi1pi)=log⁡2M.H(X)=\sum_ip_i\log_2\frac1{p_i} \le\log_2\left(\sum_ip_i\frac1{p_i}\right) =\boxed{\log_2M}.

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

−I(X;Y)=∑i,jpijlog⁡2piqjpij≤log⁡2(∑i,jpijpiqjpij)=log⁡2(∑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(Y∣X)I(X;Y)=H(Y)-H(Y\mid X) から

H(Y∣X)≤H(Y)\boxed{H(Y\mid X)\le H(Y)}

が従う。