九州大学 システム情報科学府 情報理工学専攻 2019年8月実施 情報理論
Author
Yu
Description
【問 1】
下図は, 定常 2 重マルコフ情報源 S の状態遷移図である. 下記の設問に答えよ.
(1) 上図の状態遷移図を元に, このマルコフ情報源 S の遷移確率行列 A を求めよ. ただし, 状態 "00","10", "01" の順に, 行を記せ.
(2) 時点 0 で状態 "00" にいたとして, 時点 2 で状態 "10" にいる確率はいくらか答えよ.
(3) 上記のマルコフ情報源 S の定常分布 w=(w1,w2,w3) を求めよ. ただし, wi の添え字 i=1,2,3 は状態 "00", "10", "01" にそれぞれ対応するものとする.
(4) 定常 2 重マルコフ情報源 X1,X2,… のエントロピーレート limn→∞n1H(X1,X2,…,Xn) が H(X3∣X1,X2) に一致することを示せ. ただし, H(X3∣X1,X2) は条件付きエントロピーである.
(5) 上記のマルコフ情報源 S に対するエントロピーレート H(S) を求めよ.
【問 2】
X,Y を {0,1} に値をとる確率変数とする. パラメータ α,β,γ∈[0,1] に対し,
P(X=0)=αP(Y=0∣X=0)=βP(Y=0∣X=1)=γ,P(X=1)=1−α,,P(Y=1∣X=0)=1−β,,P(Y=1∣X=1)=1−γ,
とする. 2 値エントロピー関数を
h(p)={−plogp−(1−p)log(1−p),0,for0<p<1,forp=0,1
とする. 以下の問いに答えよ.
(1) 条件付きエントロピー H(Y∣X) を 2 値エントロピー関数を用いて表現せよ.
(2) β=1−γ とする. このとき, 相互情報量 I(X;Y)を最大化するα を求めよ. また I(X;Y) の最大値を 2 値エントロピー関数と α,β を用いて表現せよ.
(3) α,β をある値に固定する. ただし, 0<α<1 とする. 相互情報量 I(X;Y) を最小化する γ の値を α,β を用いて表せ. また, その最小値を示せ.
(4) α,β をある値に固定する. ただし, 0<α<1,β>21 とする. 相互情報量 I(X;Y) を最大化する γ の値を求めよ.
题目描述
【问题 1】原题状态迁移图给出了平稳二阶 Markov 信源 S,其状态为 00、10、01。回答:
-
根据状态图求信源 S 的转移概率矩阵 A,矩阵行按 00、10、01 的顺序排列。
-
若时刻 0 处于状态 00,求时刻 2 处于状态 10 的概率。
-
求信源 S 的平稳分布 w=(w1,w2,w3),其中 w1,w2,w3 依次对应状态 00、10、01。
-
对平稳二阶 Markov 信源 X1,X2,…,证明其熵率
n→∞limn1H(X1,X2,…,Xn)
等于条件熵 H(X3∣X1,X2)。
-
求上图 Markov 信源 S 的熵率 H(S)。
【问题 2】令 X,Y 为取值于 {0,1} 的随机变量。对参数
α,β,γ∈[0,1],给定
P(X=0)P(Y=0∣X=0)P(Y=0∣X=1)=α,=β,=γ,P(X=1)P(Y=1∣X=0)P(Y=1∣X=1)=1−α,=1−β,=1−γ.
定义二元熵函数
h(p)={−plogp−(1−p)log(1−p),0,0<p<1,p=0,1.
回答:
- 用二元熵函数表示条件熵 H(Y∣X)。
- 设 β=1−γ,求使互信息 I(X;Y) 最大的 α;并用二元熵函数及参数 α,β 表示 I(X;Y) 的最大值。
- 固定 α,β 且 0<α<1,用 α,β 表示使 I(X;Y) 最小的 γ,并给出该最小值。
- 固定 α,β 且 0<α<1、β>21,求使 I(X;Y) 最大的 γ。
- 二阶马尔可夫信源:由状态图构造转移矩阵,计算多步转移概率和平稳分布。
- 熵率与条件熵:利用二阶马尔可夫性和熵的链式法则证明熵率公式,并据稳态权重求具体熵率。
- 条件熵与互信息:用二元熵函数表达给定二元信道的条件熵和输出熵。
- 信道容量与参数优化:在对称或固定转移参数下对输入分布或信道参数优化互信息,求最大值或最小值及其达到条件。
Kai
【問 1】
(1)
4321000141210
(2)
A2=1698321412101638121⇒P=41
(3)
{w1+w2+w3=1wA=w⇒w=(21,41,41)
(4)
n→∞limn1H(X1,X2,…,Xn)=n→∞limn1[H(X1)+H(X2∣X1)+H(X3∣X1,X2)+⋯+H(Xn∣X1,…,Xn−1)]=n→∞limn1[H(X1)+H(X2∣X1)+H(X3∣X1,X2)+⋯+H(Xn∣Xn−2,Xn−1)]=n→∞limn1[H(X3∣X1,X2)+⋯+H(Xn∣Xn−2,Xn−1)]=n→∞limnn−2H(X3∣X1,X2)=H(X3∣X1,X2)
(5)
H(S)=w1H(41)+w2H(21)+w3H(1)=45−83log23
【問 2】
(1)
H(Y∣X)=αh(β)+(1−α)h(γ)
(2)
α=21
αmaxI(X;Y)=1−h(β)
(3)
γminI(X;Y)=0
(4)
I(X;Y)f(γ)=H(Y)−H(Y∣X)=h(αβ+(1−α)γ)−αh(β)−(1−α)h(γ)=h(αβ+(1−α)γ)−(1−α)h(γ)の最大値問題を解くことになる