9.3 二値情報源
9.2 節は𝑅( ⋅)のグラフの両端を決めた.左の端では,情報源アルファベットが再現アルファベットに含まれ,しかも歪み尺度が「0になるのは一致するときだけ」を満たすとき,値は情報源のエントロピーである(命題 9.2.4).右の端では,𝐷max以上の歪みを許すと値が0になる(命題 9.2.5).残っているのは両端のあいだの形で,本節はいちばん小さい情報源についてそれを決める.情報源も再現も0と1の二文字にとり,歪みは当たったか外れたかだけを見る Hamming 歪み(例 9.1.3)にとる.答えは,歪みを許しすぎない範囲では,情報源のエントロピーから,許した歪みのぶんのエントロピーを引いた形になる.
本節のアルファベットは有限だから,定義 9.1.5 と 9.2 節の主張はそのまま当てられる.記号を一つ断っておく.本節はPr[𝑋 =1]を𝜋と書く.補題 9.2.2 の証明が同時分布を𝜋𝑗(𝑥,ˆ𝑥)と書いていて字が重なるが,あちらはつねに二つの引数をとるので,引数のない𝜋はつねに二値情報源のパラメータである.9.2 節が両端を確かめたときに使った𝜋と同じものである.
二値エントロピー関数の形
答えを書き下すには,二値エントロピー関数𝐻𝑏(例 1.1.2)について三つのことが要る.左右対称であること,凹であること,そして左半分で増えることである.𝐻𝑏が𝑝 =1/2で最大になるのだから三つめもそうなっているはずで,実際に前の二つから出る.出発点は 1.1 節で認めた𝜑(𝑡) = −𝑡log𝑡の凹性ひとつである.
補題 9.3.1(二値エントロピー関数の対称性・凹性・単調性). 𝐻𝑏を例 1.1.2 の二値エントロピー関数とする.次の三つが成り立つ.
- どの𝑝 ∈[0,1]についても𝐻𝑏(1 −𝑝) =𝐻𝑏(𝑝)である.
- 𝑎,𝑏 ∈[0,1]と𝜆 ∈[0,1]について𝐻𝑏(𝜆𝑎 +(1 −𝜆)𝑏) ≥𝜆𝐻𝑏(𝑎) +(1 −𝜆)𝐻𝑏(𝑏)である.すなわち𝐻𝑏は[0,1]の上で凹である.
- 0 ≤𝑠 ≤𝑡 ≤1/2を満たす実数𝑠,𝑡について𝐻𝑏(𝑠) ≤𝐻𝑏(𝑡)である.
証明.
-
例 1.1.2 の定め方から𝐻𝑏(𝑝) =𝜑(𝑝) +𝜑(1 −𝑝)である(𝜑(𝑡) = −𝑡log𝑡は 1.1 節の記号である).右辺は𝑝と1 −𝑝を入れ替えても変わらない.
-
𝑐 :=𝜆𝑎 +(1 −𝜆)𝑏とおく.1.1 節で認めた𝜑の凹性を2点の形で𝑎,𝑏に当てると𝜑(𝑐) ≥𝜆𝜑(𝑎) +(1 −𝜆)𝜑(𝑏)である.同じ凹性を1 −𝑎,1 −𝑏に当てると,𝜆(1 −𝑎) +(1 −𝜆)(1 −𝑏) =1 −𝑐だから𝜑(1 −𝑐) ≥𝜆𝜑(1 −𝑎) +(1 −𝜆)𝜑(1 −𝑏)である.二つを足すと第1の等式で使った形から𝐻𝑏(𝑐) ≥𝜆𝐻𝑏(𝑎) +(1 −𝜆)𝐻𝑏(𝑏)を得る.
-
𝑠 =1/2ならば𝑡 =1/2で両辺が等しいから,𝑠 <1/2としてよい.𝑡を𝑠と1 −𝑠の凸結合として書けるように重みを選びたい.そうできれば,第2の主張の凹性より𝐻𝑏(𝑡)は𝐻𝑏(𝑠)と𝐻𝑏(1 −𝑠)の重み付き平均以上であり,第1の主張より後者は𝐻𝑏(𝑠)に等しいから,𝐻𝑏(𝑡) ≥𝐻𝑏(𝑠)が出る.𝜆𝑠 +(1 −𝜆)(1 −𝑠) =𝑡を𝜆について解くと
𝜆:=1−𝑠−𝑡1−2𝑠とおく.𝑠 <1/2より分母は正である.𝑠 <1/2と𝑡 ≤1/2より𝑠 +𝑡 <1だから分子も正であり,𝑠 ≤𝑡より1 −𝑠 −𝑡 ≤1 −2𝑠だから𝜆 ∈(0,1]である.また
𝜆𝑠+(1−𝜆)(1−𝑠)=(1−𝑠)−𝜆(1−2𝑠)=(1−𝑠)−(1−𝑠−𝑡)=𝑡である.𝑠と1 −𝑠はどちらも[0,1]に入るから,第2の主張を𝑎 :=𝑠,𝑏 :=1 −𝑠ととって当てると𝐻𝑏(𝑡) ≥𝜆𝐻𝑏(𝑠) +(1 −𝜆)𝐻𝑏(1 −𝑠)であり,第1の主張より𝐻𝑏(1 −𝑠) =𝐻𝑏(𝑠)だから,右辺は𝐻𝑏(𝑠)に等しい.◻
下界
定義 9.1.5 が最小にするのは𝐼(𝑝;𝑞)である.定理 1.3.4 でこれを𝐻(𝑋) −𝐻(𝑋 ∣ˆ𝑋)と読むと,第1項は𝑞を含まない定数だから,問いは「𝑋とˆ𝑋の結びつきをいちばん弱くする」から「再現を見たあとに残る𝑋の不確かさをいちばん大きくする」に移る.そして歪みの制約が言っているのは,二つの文字が食い違う確率が𝐷以下だ,ということだけである.食い違ったかどうかを表す二値の確率変数を一つ置けば,そこに二値エントロピー関数が出てくる.
命題 9.3.2(相互情報量の下界). 𝜋 ∈[0,1],0 ≤𝐷 ≤1/2とする.X =ˆX ={0,1}とし,𝑝を𝑝(1) =𝜋,𝑝(0) =1 −𝜋で定まるX上の分布,𝑑𝐻を例 9.1.3 の Hamming 歪みとし,Q( ⋅)と𝐼(𝑝; ⋅)を定義 9.1.5 のとおりとする.このとき,どの𝑞 ∈Q(𝐷)についても
𝐼(𝑝;𝑞)≥𝐻𝑏(𝜋)−𝐻𝑏(𝐷)である(𝐻𝑏は例 1.1.2 の二値エントロピー関数).
証明. 𝑞 ∈Q(𝐷)をとり,(𝑋,ˆ𝑋)を定義 9.1.5 のとおり同時分布𝑝(𝑥)𝑞(ˆ𝑥 ∣𝑥)をもつ対とする.定義 9.1.5 より𝐼(𝑝;𝑞) =𝐼(𝑋;ˆ𝑋)であり,定理 1.3.4 より𝐼(𝑋;ˆ𝑋) =𝐻(𝑋) −𝐻(𝑋 ∣ˆ𝑋)である.𝑋は確率𝜋で1,確率1 −𝜋で0をとるから,例 1.1.2 より𝐻(𝑋) =𝐻𝑏(𝜋)である.
食い違いを表す確率変数を置く.𝑋 ≠ˆ𝑋のとき1,𝑋 =ˆ𝑋のとき0をとる確率変数を𝐸とする.例 9.1.3 を𝑛 =1に当てると,𝑞の期待歪みは𝔼[𝑑𝐻(𝑋,ˆ𝑋)] =Pr[ˆ𝑋 ≠𝑋] =Pr[𝐸 =1]であり,𝑞 ∈Q(𝐷)よりPr[𝐸 =1] ≤𝐷である.
𝐻(𝑋 ∣ˆ𝑋) =𝐻(𝐸 ∣ˆ𝑋)を見る.Pr[ˆ𝑋 =ˆ𝑥] >0となるˆ𝑥 ∈{0,1}を固定し,𝑒ˆ𝑥 :=Pr[𝑋 ≠ˆ𝑥 ∣ˆ𝑋 =ˆ𝑥]とおく(Pr[ˆ𝑋 =ˆ𝑥] =0の項は,定義 1.2.2 のとおりどちらの量にも寄与しない).𝑋は{0,1}に値をとるから,ˆ𝑋 =ˆ𝑥のもとでの𝑋の条件付き分布は,ˆ𝑥でないほうの文字に𝑒ˆ𝑥,ˆ𝑥に1 −𝑒ˆ𝑥を置いたものである.定義 1.2.2 の𝐻(𝑋 ∣ˆ𝑋 =ˆ𝑥)はこの分布に定義 1.1.1 を当てたもので,項は二つだけだから−𝑒ˆ𝑥log𝑒ˆ𝑥 −(1 −𝑒ˆ𝑥)log(1 −𝑒ˆ𝑥)に等しく,例 1.1.2 よりこれは𝐻𝑏(𝑒ˆ𝑥)である.一方ˆ𝑋 =ˆ𝑥のもとで𝐸 =1となるのは𝑋 ≠ˆ𝑥のときだから,ˆ𝑋 =ˆ𝑥のもとでの𝐸の条件付き分布は1に𝑒ˆ𝑥,0に1 −𝑒ˆ𝑥を置いたもので,同じ理由で𝐻(𝐸 ∣ˆ𝑋 =ˆ𝑥) =𝐻𝑏(𝑒ˆ𝑥)である.定義 1.2.2 はこれらをPr[ˆ𝑋 =ˆ𝑥]で平均したものだから,二つの条件付きエントロピーは等しい.
条件を落として抑える.定理 1.2.4 の基本形(条件の一方を自明な定数にとった形)より𝐻(𝐸 ∣ˆ𝑋) ≤𝐻(𝐸)である.𝐸は確率Pr[𝐸 =1]で1,確率1 −Pr[𝐸 =1]で0をとるから,例 1.1.2 より𝐻(𝐸) =𝐻𝑏(Pr[𝐸 =1])である.Pr[𝐸 =1] ≤𝐷 ≤1/2だから補題 9.3.1 の第3の主張より𝐻𝑏(Pr[𝐸 =1]) ≤𝐻𝑏(𝐷)である.
以上を合わせると𝐻(𝑋 ∣ˆ𝑋) ≤𝐻𝑏(𝐷)であり,𝐼(𝑝;𝑞) =𝐻𝑏(𝜋) −𝐻(𝑋 ∣ˆ𝑋) ≥𝐻𝑏(𝜋) −𝐻𝑏(𝐷)を得る.◼
逆向きに作る
命題 9.3.2 の証明で不等号が入ったのは二か所である.定理 1.2.4 で条件ˆ𝑋を落としたところと,補題 9.3.1 の第3の主張でPr[𝐸 =1]を𝐷に置き換えたところである.下界を達成する条件付き分布を作る筋は,抑えたこの二つを抑えた値ちょうどにすることである.食い違いの確率をちょうど𝐷にし,しかも再現ˆ𝑋を知っても食い違いについて何も分からないようにする,ということである.
後者が逆向きの構成の理由になる.素直に思いつくのは,情報源𝑋を先に置いて,そこに確率𝐷の反転を独立に加えて再現ˆ𝑋を作る向きである.この作り方では,食い違いは作った時点で𝑋と独立になる.欲しいのはそちらではなく,ˆ𝑋と独立な食い違いのほうである.順序を入れ替えて,再現ˆ𝑋を先に置き,そこに確率𝐷の反転を独立に加えて情報源𝑋を作れば,食い違いは作った時点でˆ𝑋と独立である.反転確率𝐷の二元対称通信路(例 6.1.8)を,情報源から再現へではなく,再現から情報源へ向けて使う,と言ってもよい.残るのは,こうして作った同時分布の第1周辺が情報源の分布𝑝に一致するように,ˆ𝑋の側のパラメータを選ぶことだけである.
命題 9.3.3(下界を達成する条件付き分布). 𝜋を実数,𝐷を実数とし,0 ≤𝐷 <min(𝜋,1 −𝜋)とする.X =ˆX ={0,1}とし,𝑝を𝑝(1) =𝜋,𝑝(0) =1 −𝜋で定まるX上の分布,𝑑𝐻を例 9.1.3 の Hamming 歪みとし,Q( ⋅)と𝐼(𝑝; ⋅)を定義 9.1.5 のとおりとする.このとき,期待歪みがちょうど𝐷に等しい条件付き分布𝑞∗ ∈Q(𝐷)であって
𝐼(𝑝;𝑞∗)=𝐻𝑏(𝜋)−𝐻𝑏(𝐷)を満たすものが存在する(𝐻𝑏は例 1.1.2 の二値エントロピー関数).
証明. 仮定からmin(𝜋,1 −𝜋) >0,すなわち0 <𝜋 <1である.またmin(𝜋,1 −𝜋) ≤1/2だから𝐷 <1/2で,1 −2𝐷 >0である.ˆ𝑋が確率𝑟で1をとるとして,第1周辺が𝜋になるという条件𝑟(1 −𝐷) +(1 −𝑟)𝐷 =𝜋を𝑟について解くと
𝑟:=𝜋−𝐷1−2𝐷とおく.分母は正であり,𝐷 <𝜋より分子も正だから𝑟 >0である.𝐷 <1 −𝜋より𝜋 −𝐷 <1 −2𝐷だから𝑟 <1である.
同時分布を書き下す.X ׈X上の数の組𝐽を
𝐽(1,1):=𝑟(1−𝐷),𝐽(0,1):=𝑟𝐷,𝐽(1,0):=(1−𝑟)𝐷,𝐽(0,0):=(1−𝑟)(1−𝐷)で定める.どれも非負であり,総和は𝑟 +(1 −𝑟) =1である.これは「ˆ𝑋を確率𝑟で1にとり,そこに確率𝐷の反転を独立に加えて𝑋を作る」という手順が与える同時分布にほかならない.
第1周辺が𝑝であることを見る.𝐽(1,1) +𝐽(1,0) =𝑟(1 −𝐷) +(1 −𝑟)𝐷 =𝐷 +𝑟(1 −2𝐷) =𝐷 +(𝜋 −𝐷) =𝜋であり,𝐽(0,1) +𝐽(0,0) =𝑟𝐷 +(1 −𝑟)(1 −𝐷) =(1 −𝐷) −𝑟(1 −2𝐷) =(1 −𝐷) −(𝜋 −𝐷) =1 −𝜋である.すなわち𝑥について足すと𝑝(𝑥)になる.
条件付き分布を取り出す.𝑝(1) =𝜋 >0,𝑝(0) =1 −𝜋 >0だから𝑞∗(ˆ𝑥 ∣𝑥) :=𝐽(𝑥,ˆ𝑥)/𝑝(𝑥)と定めることができ,値は非負で,各𝑥について∑ˆ𝑥𝑞∗(ˆ𝑥 ∣𝑥) =𝑝(𝑥)/𝑝(𝑥) =1である.よって𝑞∗は定義 9.1.5 の条件付き分布であり,𝑝(𝑥)𝑞∗(ˆ𝑥 ∣𝑥) =𝐽(𝑥,ˆ𝑥)だから,𝑞∗に対応する対(𝑋,ˆ𝑋)の同時分布は𝐽である.
期待歪みを求める.𝑑𝐻は𝑥 ≠ˆ𝑥の対でだけ値1をとるから,定義 9.1.5 の期待歪みは𝐽(0,1) +𝐽(1,0) =𝑟𝐷 +(1 −𝑟)𝐷 =𝐷である.よって𝑞∗ ∈Q(𝐷)であり,期待歪みはちょうど𝐷である.
相互情報量を求める.定理 1.3.4 より𝐼(𝑝;𝑞∗) =𝐻(𝑋) −𝐻(𝑋 ∣ˆ𝑋)である.第1周辺が𝑝だから例 1.1.2 より𝐻(𝑋) =𝐻𝑏(𝜋)である.第2周辺はPr[ˆ𝑋 =1] =𝐽(1,1) +𝐽(0,1) =𝑟(1 −𝐷) +𝑟𝐷 =𝑟,Pr[ˆ𝑋 =0] =1 −𝑟で,0 <𝑟 <1よりどちらも正である.Pr[𝑋 ≠1 ∣ˆ𝑋 =1] =𝐽(0,1)/𝑟 =𝐷であり,Pr[𝑋 ≠0 ∣ˆ𝑋 =0] =𝐽(1,0)/(1 −𝑟) =𝐷である.したがって,どちらのˆ𝑥についても,ˆ𝑋 =ˆ𝑥のもとでの𝑋の条件付き分布はˆ𝑥でないほうの文字に𝐷,ˆ𝑥に1 −𝐷を置いたものであり,定義 1.2.2 と定義 1.1.1 により𝐻(𝑋 ∣ˆ𝑋 =ˆ𝑥)は二項の和−𝐷log𝐷 −(1 −𝐷)log(1 −𝐷),すなわち例 1.1.2 の𝐻𝑏(𝐷)に等しい.定義 1.2.2 はこれをPr[ˆ𝑋 =ˆ𝑥]で平均したものだから𝐻(𝑋 ∣ˆ𝑋) =𝐻𝑏(𝐷)であり,𝐼(𝑝;𝑞∗) =𝐻𝑏(𝜋) −𝐻𝑏(𝐷)を得る.◼
レート歪み関数
定理 9.3.4(二値情報源のレート歪み関数). 𝜋 ∈[0,1]とする.X =ˆX ={0,1}とし,𝑝を𝑝(1) =𝜋,𝑝(0) =1 −𝜋で定まるX上の分布,𝑑𝐻を例 9.1.3 の Hamming 歪みとし,Q( ⋅)と𝑅( ⋅)を定義 9.1.5 のとおりとする.実数𝐷について次が成り立つ(𝐻𝑏は例 1.1.2 の二値エントロピー関数).
- 0 ≤𝐷 <min(𝜋,1 −𝜋)ならばQ(𝐷)は空でなく,𝑅(𝐷) =𝐻𝑏(𝜋) −𝐻𝑏(𝐷)である.
- 𝐷 ≥min(𝜋,1 −𝜋)ならばQ(𝐷)は空でなく,𝑅(𝐷) =0である.
証明.
-
命題 9.3.3 より,𝐼(𝑝;𝑞∗) =𝐻𝑏(𝜋) −𝐻𝑏(𝐷)を満たす𝑞∗ ∈Q(𝐷)が存在する.とくにQ(𝐷)は空でなく,𝑅(𝐷)は値の集合{𝐼(𝑝;𝑞) :𝑞 ∈Q(𝐷)}の下限だから𝑅(𝐷) ≤𝐻𝑏(𝜋) −𝐻𝑏(𝐷)である.逆向きを見る.min(𝜋,1 −𝜋) ≤1/2だから0 ≤𝐷 ≤1/2であり,命題 9.3.2 よりこの値の集合のどの元も𝐻𝑏(𝜋) −𝐻𝑏(𝐷)以上である.下限は下界のうち最大のものだから𝑅(𝐷) ≥𝐻𝑏(𝜋) −𝐻𝑏(𝐷)である.二つを合わせて等号を得る.
-
まず命題 9.2.5 の𝐷maxを計算する.ˆ𝑥 =0に対しては∑𝑥𝑝(𝑥) 𝑑𝐻(𝑥,0) =𝑝(1) =𝜋,ˆ𝑥 =1に対しては∑𝑥𝑝(𝑥) 𝑑𝐻(𝑥,1) =𝑝(0) =1 −𝜋だから,𝐷max =min(𝜋,1 −𝜋)である.命題 9.2.5 よりQ(𝐷max)は空でないから,𝐷max ≤𝐷に命題 9.2.1 を当てるとQ(𝐷)も空でない.そのうえで命題 9.2.5 の後半を𝐷に当てると,𝐷 ≥𝐷maxだから𝑅(𝐷) =0である.
◼
両端で確かめる. 左の端は𝐷 =0で,定理 9.3.4 の第1の場合の値は𝐻𝑏(𝜋) −𝐻𝑏(0) =𝐻𝑏(𝜋)である.X =ˆXであり,𝑑𝐻(𝑥,ˆ𝑥) =0がˆ𝑥 =𝑥と同値だから,命題 9.2.4 の仮定が満たされていて,これはその𝑅(0) =𝐻(𝑋)と一致している.右の端は𝐷max =min(𝜋,1 −𝜋)で,第1の場合の式にこの値を入れると,𝜋 ≤1/2なら𝐻𝑏(𝜋) −𝐻𝑏(𝜋) =0,𝜋 >1/2なら𝐻𝑏(𝜋) −𝐻𝑏(1 −𝜋) =0である(後者は補題 9.3.1 の第1の主張による).すなわち第1の場合の式と第2の場合の値は,境目で段差なくつながっている.
通信路容量と同じ式になる. 公平なコイン(𝜋 =1/2)では𝐷max =1/2で,0 ≤𝐷 <1/2のとき𝑅(𝐷) =1 −𝐻𝑏(𝐷)ビットである.第6章 例 6.1.8 の二元対称通信路の容量と同じ式が,ここでは最小にした相互情報量として現れている.偶然ではない.命題 9.3.3 の作り方は反転確率𝐷の二元対称通信路を再現から情報源へ向けて使うもので,𝜋 =1/2のときは再現の側が一様分布になる(命題 9.3.3 の証明の𝑟が1/2になる).その通信路の容量を達成する入力がまさに一様入力なので,下界を達成する対の相互情報量が,例 6.1.8 の容量の値そのものになる.
得られた式が 9.2 節の一般論と食い違わないことも見ておく.𝐷について非増加であることは,補題 9.3.1 の第3の主張が[0,1/2]の上で𝐻𝑏が増えると言っているので,−𝐻𝑏(𝐷)が増えないことから読める(命題 9.2.1).凸であることは,補題 9.3.1 の第2の主張により𝐻𝑏が凹だから,−𝐻𝑏が凸であることから読める(命題 9.2.3).一般論が形だけを言っていたところに,二値の情報源では閉じた式が入る.
数で見る
例 9.3.5(𝜋 =1/4の二値情報源). X =ˆX ={0,1}とし,𝑝(1) =1/4,𝑝(0) =3/4,𝑑 =𝑑𝐻(例 9.1.3)として𝑅( ⋅)を定義 9.1.5 のとおりとする.このとき𝑅(0)は約0.8113ビット,𝑅(0.05)は約0.5249ビット,𝑅(0.1)は約0.3423ビットであり,𝐷 ≥1/4では𝑅(𝐷) =0である.
証明. min(1/4,3/4) =1/4だから,定理 9.3.4 より0 ≤𝐷 <1/4では𝑅(𝐷) =𝐻𝑏(1/4) −𝐻𝑏(𝐷)であり,𝐷 ≥1/4では𝑅(𝐷) =0である.
𝐻𝑏(1/4)を求める.log23 =1.5849…だから
𝐻𝑏(1/4)=14log24+34log243=12+34(2−log23)=0.5+0.75×0.4150…=0.8112…である.𝐷 =0では𝐻𝑏(0) =0だから𝑅(0) =𝐻𝑏(1/4)で,約0.8113ビットである.
𝐷 =0.05ではlog20.05 = −4.3219…,log20.95 = −0.0740…だから𝐻𝑏(0.05) =0.05 ×4.3219… +0.95 ×0.0740… =0.2160… +0.0703… =0.2863…であり,𝑅(0.05) =0.8112… −0.2863… =0.5248…で約0.5249ビットである.
𝐷 =0.1ではlog20.1 = −3.3219…,log20.9 = −0.1520…だから𝐻𝑏(0.1) =0.1 ×3.3219… +0.9 ×0.1520… =0.3321… +0.1368… =0.4689…であり,𝑅(0.1) =0.8112… −0.4689… =0.3422…で約0.3423ビットである.◼
規模感. 例 9.3.5 の情報源では,歪みをいっさい許さないときの値が約0.8113ビットである.文字の5パーセントの食い違いを許すだけで約0.5249ビットまで下がり,減り方はおよそ35パーセントである.10パーセントの食い違いまで許すと約0.3423ビットで,およそ58パーセント減る.最初のうちは,わずかな歪みを許すだけで値が大きく落ちる.1/4まで許すと0になる.情報源を見ずにいつも0を出す再現でも,食い違うのは文字の1/4だけだからである(命題 9.2.5 の𝐷maxの読み方).
ここまでの主張は,どれも最適化問題の値についてのものである.9.5 節と 9.6 節へ進む前に,9.4 節でもう一つの具体例,すなわち実数値をとる情報源を扱う.
InformationTheory — 形式化検証つき情報理論教科書(レビュー版).数式は MathJax + AMS Euler で事前レンダリング.