9.2 レート歪み関数の性質
9.1 節は𝑅(𝐷)を,歪みを𝐷以下に抑えるという制約のもとでの相互情報量の下限として定め,それが最小値でもあること(命題 9.1.8)まで確かめた.値の計算はまだ一つもしていない.本節では,具体的な情報源を決める前に,𝐷の関数として𝑅( ⋅)がどんな形をしているかを調べる.示すのは四つで,𝐷について非増加であること,凸であること,𝐷 =0での値,そして𝑅(𝐷)が0になる𝐷のうち最小のものである.どれも 9.3 節と 9.4 節で具体的な情報源のレート歪み関数を計算するときに,答え合わせに使える.記号は定義 9.1.1 と定義 9.1.5 のものをそのまま引き継ぐ.
歪みを多く許すほど下がる
命題 9.2.1. XとˆXを空でない有限集合,𝑝をX上の分布,𝑑を歪み尺度(定義 9.1.1)とし,Q( ⋅)と𝑅( ⋅)を定義 9.1.5 のとおりとする.実数𝐷1 ≤𝐷2についてQ(𝐷1)が空でないならば,Q(𝐷1) ⊆Q(𝐷2)であり,Q(𝐷2)も空でなく
𝑅(𝐷2)≤𝑅(𝐷1)である.
証明. 𝑞 ∈Q(𝐷1)とすると,その期待歪みは𝐷1以下であり,𝐷1 ≤𝐷2だから𝐷2以下でもある.よって𝑞 ∈Q(𝐷2)であり,包含が従う.この段はQ(𝐷1)が空でも通るので,包含だけなら空でないという仮定は要らない.Q(𝐷1)が空でないのでQ(𝐷2)も空でなく,𝑅(𝐷2)が定まる.
𝑅(𝐷2)は値の集合{𝐼(𝑝;𝑞) :𝑞 ∈Q(𝐷2)}の下界である.いま見た包含より{𝐼(𝑝;𝑞) :𝑞 ∈Q(𝐷1)}はその部分集合だから,𝑅(𝐷2)は後者の下界でもある.下限は下界のうち最大のものだから𝑅(𝐷2) ≤𝑅(𝐷1)である.◼
読み方は素直である.歪みを多く許すほど,選べる再現の作り方は増える.増えたぶんだけ小さい値が選べるかもしれず,少なくとも大きくはならない.効いているのは制約集合の包含だけで,相互情報量の性質は一つも使っていない.
凸である
命題 9.2.1 の証明は制約集合の包含だけで通ったが,凸性の証明はそれでは済まない.二つの再現の作り方を混ぜたときに相互情報量がどう動くかが要るので,それを先に切り出しておく.道具は第1章 定理 1.7.1 の対数和不等式ひとつである.
補題 9.2.2(条件付き分布についての相互情報量の凸性). XとˆXを空でない有限集合,𝑝をX上の分布,𝑞1と𝑞2を条件付き分布(定義 9.1.5),𝜆 ∈[0,1]とする.𝑞𝜆(ˆ𝑥 ∣𝑥) :=𝜆 𝑞1(ˆ𝑥 ∣𝑥) +(1 −𝜆) 𝑞2(ˆ𝑥 ∣𝑥)で定まる𝑞𝜆もまた条件付き分布であり
𝐼(𝑝;𝑞𝜆)≤𝜆𝐼(𝑝;𝑞1)+(1−𝜆)𝐼(𝑝;𝑞2)である.
証明. 筋は,対(𝑥,ˆ𝑥)ごとに対数和不等式を当てて,それを足し上げるだけである.先に,確率0の項が現れない対について本線を書き切り,確率0の対の扱いは末尾にまとめる.
𝑞𝜆が条件付き分布であることは,非負の数の非負係数の和が非負であることと,各𝑥について∑ˆ𝑥𝑞𝜆(ˆ𝑥 ∣𝑥) =𝜆 +(1 −𝜆) =1であることによる.𝜆 =0と𝜆 =1では両辺が一致するので,以下では0 <𝜆 <1とする.
記号を用意する.𝑗 =1,2,𝜆について,𝑞𝑗に対応する同時分布を𝜋𝑗(𝑥,ˆ𝑥) :=𝑝(𝑥) 𝑞𝑗(ˆ𝑥 ∣𝑥),その第2周辺分布をˆ𝑝𝑗(ˆ𝑥) :=∑𝑥𝜋𝑗(𝑥,ˆ𝑥)と書く.どの𝑗でも第1周辺分布は𝑝である.𝑞𝜆の定め方から,各(𝑥,ˆ𝑥)について𝜋𝜆(𝑥,ˆ𝑥) =𝜆 𝜋1(𝑥,ˆ𝑥) +(1 −𝜆) 𝜋2(𝑥,ˆ𝑥)であり,これを𝑥について足して,各ˆ𝑥についてˆ𝑝𝜆(ˆ𝑥) =𝜆 ˆ𝑝1(ˆ𝑥) +(1 −𝜆) ˆ𝑝2(ˆ𝑥)である.
対(𝑥,ˆ𝑥)を一つ固定し,
𝑎1:=𝜆𝜋1(𝑥,ˆ𝑥),𝑎2:=(1−𝜆)𝜋2(𝑥,ˆ𝑥),𝑏1:=𝜆𝑝(𝑥)ˆ𝑝1(ˆ𝑥),𝑏2:=(1−𝜆)𝑝(𝑥)ˆ𝑝2(ˆ𝑥)とおく.これらは非負で,𝑎1 +𝑎2 =𝜋𝜆(𝑥,ˆ𝑥),𝑏1 +𝑏2 =𝑝(𝑥) ˆ𝑝𝜆(ˆ𝑥)である.𝑏1と𝑏2がどちらも正である対では,定理 1.7.1 の対数和不等式がそのまま当たって
𝜋𝜆(𝑥,ˆ𝑥)log𝜋𝜆(𝑥,ˆ𝑥)𝑝(𝑥)ˆ𝑝𝜆(ˆ𝑥)≤𝜆𝜋1(𝑥,ˆ𝑥)log𝜋1(𝑥,ˆ𝑥)𝑝(𝑥)ˆ𝑝1(ˆ𝑥)+(1−𝜆)𝜋2(𝑥,ˆ𝑥)log𝜋2(𝑥,ˆ𝑥)𝑝(𝑥)ˆ𝑝2(ˆ𝑥)を得る(対数和不等式の右辺の𝑎𝑗log(𝑎𝑗/𝑏𝑗)では,𝑎𝑗と𝑏𝑗に掛かる𝜆と1 −𝜆が対数の中で約分される).
この不等式がどの対でも成り立てば,あとは足すだけである.定義 1.3.1 の和は同時分布が正である項についてとるもので,同時分布が0である項の寄与は0だから,上の不等式をX ׈Xの全体で足し合わせると,左辺は𝐼(𝑝;𝑞𝜆),右辺は𝜆𝐼(𝑝;𝑞1) +(1 −𝜆)𝐼(𝑝;𝑞2)になる.
残るのは確率0の対で,そこでも同じ不等式が成り立つ.まず,𝑏𝑗 =0となる𝑗があれば,そこでは𝑎𝑗 =0でもある.実際0 <𝜆 <1だから,𝑏𝑗 =0は𝑝(𝑥) =0またはˆ𝑝𝑗(ˆ𝑥) =0を意味する.前者なら𝜋𝑗(𝑥,ˆ𝑥) =𝑝(𝑥)𝑞𝑗(ˆ𝑥 ∣𝑥) =0であり,後者なら𝜋𝑗(𝑥,ˆ𝑥)は非負の数の和ˆ𝑝𝑗(ˆ𝑥)の一項だからやはり0である.したがって,𝑏𝑗 =0である𝑗を落としても𝑎1 +𝑎2と𝑏1 +𝑏2は変わらない.落とした𝑗については𝜋𝑗(𝑥,ˆ𝑥) =0だから,定義 1.3.1 の和でこの対が𝐼(𝑝;𝑞𝑗)に与える寄与も0である.そこで𝑏𝑗 >0である𝑗だけを残して対数和不等式を当て,落とした𝑗の項を右辺に0として書き足せば,上と同じ不等式を得る.残る𝑗が一つも無ければ𝑎1 =𝑎2 =0,すなわち𝜋1(𝑥,ˆ𝑥) =𝜋2(𝑥,ˆ𝑥) =𝜋𝜆(𝑥,ˆ𝑥) =0となり,この対が不等式の両辺に与える寄与はどちらも0である.◻
補題 9.2.2 が言っているのは,二つの再現の作り方を混ぜると,情報源との結びつきは混ぜる前の平均より強くならない,ということである.混ぜているのは𝑞だけで,情報源の分布𝑝は固定していることに注意したい.不等号の向きは,相互情報量を隔たりとして読むと見える.定義 1.3.1 の式をそのまま読めば,相互情報量は同時分布と二つの周辺分布の積との相対エントロピー(1.6 節)であり,情報源と再現が独立からどれだけ離れているかを測っている.証明で対数和不等式を当てた位置を見ると,まとめているのは「どちらの作り方を使ったか」という区別である.対数和不等式が言うのは,区別を捨ててまとめると値が小さくなることだった(定理 1.7.1).つまり,混ぜて区別を捨てた作り方が独立から隔たっている量は,二つの隔たりを同じ重みで混ぜた値を超えない.
命題 9.2.3. XとˆXを空でない有限集合,𝑝をX上の分布,𝑑を歪み尺度(定義 9.1.1)とし,Q( ⋅)と𝑅( ⋅)を定義 9.1.5 のとおりとする.実数𝐷1,𝐷2についてQ(𝐷1)とQ(𝐷2)がどちらも空でないならば,𝜆 ∈[0,1]に対してQ(𝜆𝐷1 +(1 −𝜆)𝐷2)も空でなく
𝑅(𝜆𝐷1+(1−𝜆)𝐷2)≤𝜆𝑅(𝐷1)+(1−𝜆)𝑅(𝐷2)である.
証明. 命題 9.1.8 により,𝐼(𝑝; ⋅)をQ(𝐷1)の上で最小にする𝑞1と,Q(𝐷2)の上で最小にする𝑞2をとることができ,𝑅(𝐷1) =𝐼(𝑝;𝑞1),𝑅(𝐷2) =𝐼(𝑝;𝑞2)である.𝑞𝜆 :=𝜆𝑞1 +(1 −𝜆)𝑞2を補題 9.2.2 のとおりに定める.
𝑞𝜆がQ(𝜆𝐷1 +(1 −𝜆)𝐷2)に属することを見る.補題 9.2.2 より𝑞𝜆は条件付き分布である.定義 9.1.5 の期待歪みは𝑞の成分について一次だから
∑𝑥,ˆ𝑥𝑝(𝑥)𝑞𝜆(ˆ𝑥∣𝑥)𝑑(𝑥,ˆ𝑥)=𝜆∑𝑥,ˆ𝑥𝑝(𝑥)𝑞1(ˆ𝑥∣𝑥)𝑑(𝑥,ˆ𝑥)+(1−𝜆)∑𝑥,ˆ𝑥𝑝(𝑥)𝑞2(ˆ𝑥∣𝑥)𝑑(𝑥,ˆ𝑥)であり,𝜆 ≥0と1 −𝜆 ≥0より右辺は𝜆𝐷1 +(1 −𝜆)𝐷2以下である.よって𝑞𝜆は制約を満たし,とくにQ(𝜆𝐷1 +(1 −𝜆)𝐷2)は空でない.
𝑅( ⋅)は値の集合の下限だから𝑅(𝜆𝐷1 +(1 −𝜆)𝐷2) ≤𝐼(𝑝;𝑞𝜆)であり,補題 9.2.2 より
𝐼(𝑝;𝑞𝜆)≤𝜆𝐼(𝑝;𝑞1)+(1−𝜆)𝐼(𝑝;𝑞2)=𝜆𝑅(𝐷1)+(1−𝜆)𝑅(𝐷2)である.二つを合わせて主張を得る.◼
命題 9.2.3 は,二つの歪みの上限のあいだを線形に補間した点で,𝑅( ⋅)の値が二つの端を結ぶ弦より下にあることを言っている.証明が使ったのは二つだけで,𝐷1と𝐷2で最小を与える作り方を混ぜると,期待歪みは二つの期待歪みを同じ重みで混ぜた値になり(期待歪みが𝑞について一次だから),結びつきは二つの相互情報量の平均より強くならない(補題 9.2.2)ということである.
なぜ弦より下になると期待できるのかは,時間で分けてみると見当がつく.符号を定義するのは 9.5 節なので符号の話としてはまだ言えないが,長さ𝑛のブロックのうち𝜆の割合を歪み𝐷1用の符号で,残りを歪み𝐷2用の符号で処理したとする.歪みは成分についての平均だから全体で𝜆𝐷1 +(1 −𝜆)𝐷2になり,1文字あたりのレートも同じ重みの平均𝜆𝑅(𝐷1) +(1 −𝜆)𝑅(𝐷2)になる.中間の歪みをこれだけのレートで実現する手立てが一つある以上,中間の点で弦より上には行きようがない,というわけである.
両端
残るは両端である.左の端では歪みをいっさい許さず,右の端では𝑅( ⋅)が0になる.情報源アルファベットが再現アルファベットに含まれ,かつ歪み尺度が「0になるのは一致するときだけ」を満たすならば,左の端で𝑅( ⋅)はエントロピーに戻る.どちらも命題 9.2.4 の仮定に置く.歪み尺度の条件が落とせないことは,命題のあとで見る.
命題 9.2.4. XとˆXを空でない有限集合でX ⊆ˆXを満たすものとし,𝑝をX上の分布,𝑋を分布𝑝に従うXに値をとる確率変数とする.歪み尺度𝑑(定義 9.1.1)が「𝑑(𝑥,ˆ𝑥) =0となるのはˆ𝑥 =𝑥のとき,かつそのときに限る」を満たすとし,Q( ⋅)と𝑅( ⋅)を定義 9.1.5 のとおりとする.このときQ(0)は空でなく
𝑅(0)=𝐻(𝑋)である.
証明. Q(0)が空でないことを見る.X ⊆ˆXだから,ˆ𝑥 =𝑥のとき𝑞(ˆ𝑥 ∣𝑥) :=1,そうでないとき𝑞(ˆ𝑥 ∣𝑥) :=0と定めることができ,これは条件付き分布である.仮定より𝑑(𝑥,𝑥) =0だから,その期待歪みは∑𝑥𝑝(𝑥) 𝑑(𝑥,𝑥) =0であり,𝑞 ∈Q(0)である.
次にQ(0)の元がどれも同じ相互情報量をもつことを見る.𝑞 ∈Q(0)とすると,期待歪み∑𝑥,ˆ𝑥𝑝(𝑥)𝑞(ˆ𝑥 ∣𝑥)𝑑(𝑥,ˆ𝑥)は0以下であり,各項は非負だから,すべての項が0である.すなわち𝑝(𝑥)𝑞(ˆ𝑥 ∣𝑥) >0ならば𝑑(𝑥,ˆ𝑥) =0であり,仮定によりそれはˆ𝑥 =𝑥を意味する.よって𝑝(𝑥) >0である𝑥については,𝑞(ˆ𝑥 ∣𝑥) >0となるˆ𝑥が𝑥に限られ,∑ˆ𝑥𝑞(ˆ𝑥 ∣𝑥) =1と合わせて𝑞(𝑥 ∣𝑥) =1である.したがって𝑞に対応する同時分布は,ˆ𝑥 =𝑥のとき𝑝(𝑥),そうでないとき0であり,第2周辺分布はˆ𝑥 ∈Xのとき𝑝(ˆ𝑥),そうでないとき0である.
その相互情報量を定義 1.3.1 で計算する.同時分布が正になるのはˆ𝑥 =𝑥かつ𝑝(𝑥) >0のときだけで,そこでの値は𝑝(𝑥),二つの周辺分布の積は𝑝(𝑥) 𝑝(𝑥)だから
𝐼(𝑝;𝑞)=∑𝑥:𝑝(𝑥)>0𝑝(𝑥)log𝑝(𝑥)𝑝(𝑥)𝑝(𝑥)=−∑𝑥:𝑝(𝑥)>0𝑝(𝑥)log𝑝(𝑥)=𝐻(𝑋)である(最後の等号は,𝑝(𝑥) =0の項を0log0 =0と約束した定義 1.1.1 による).
値の集合が一点{𝐻(𝑋)}なので,その下限は𝐻(𝑋)である.◼
歪み尺度の条件は落とせない.たとえば𝑑が恒等的に0なら,どの条件付き分布もQ(0)に属するので,ˆXの一点に集中する条件付き分布をとれば,情報源と再現は独立になって相互情報量は0になる.相互情報量は非負(命題 1.3.2)だから,このとき𝑅(0) =0である.𝐻(𝑋) >0である情報源では,これは命題 9.2.4 の結論と違う値である.落とせないのは条件そのもので,条件を課す範囲のほうには余地がある.証明が歪み尺度の条件を使うのは情報源が正の確率を与える文字についてだけなので,確率0の文字での条件は落としてよい.
命題 9.2.4 の値は,第2章 定理 2.3.6 のinfR =𝐻(𝑋)と同じ𝐻(𝑋)である.ただし命題 9.2.4 が述べているのは最適化問題の値についてであって,符号については何も言っていない.二つを結ぶには𝑅(𝐷)が符号のレートの限界であることが要り,それが 9.5 節と 9.6 節の内容である.
命題 9.2.5(𝑅( ⋅)が0になる歪み). XとˆXを空でない有限集合,𝑝をX上の分布,𝑑を歪み尺度(定義 9.1.1)とし,Q( ⋅)と𝑅( ⋅)を定義 9.1.5 のとおりとする.
𝐷max:=minˆ𝑥∈ˆX ∑𝑥∈X𝑝(𝑥)𝑑(𝑥,ˆ𝑥)とおくと,Q(𝐷max)は空でなく𝑅(𝐷max) =0である.さらに,Q(𝐷)が空でない実数𝐷について,𝑅(𝐷) =0となるのは𝐷 ≥𝐷maxのとき,かつそのときに限る.
証明. ˆXは空でない有限集合だから𝐷maxは定まる.最小を与える点を一つとってˆ𝑥0と書き,ˆ𝑥 =ˆ𝑥0のとき𝑞0(ˆ𝑥 ∣𝑥) :=1,そうでないとき𝑞0(ˆ𝑥 ∣𝑥) :=0と定める.これは条件付き分布であり,その期待歪みは∑𝑥𝑝(𝑥) 𝑑(𝑥,ˆ𝑥0) =𝐷maxである.
𝐷 ≥𝐷maxならば𝑅(𝐷) =0であることを見る(𝐷 =𝐷maxを含む).𝑞0の期待歪みは𝐷max ≤𝐷だから𝑞0 ∈Q(𝐷)であり,とくにQ(𝐷)は空でない.𝑞0に対応する同時分布は,ˆ𝑥 =ˆ𝑥0のとき𝑝(𝑥),そうでないとき0であって,これは第1周辺分布𝑝と,ˆ𝑥0に1を置く第2周辺分布との積にほかならない.よって𝑋とˆ𝑋は独立で,命題 1.3.2 より𝐼(𝑝;𝑞0) =0である.したがって𝑅(𝐷) ≤0である.一方,命題 1.3.2 よりQ(𝐷)のどの元でも𝐼(𝑝;𝑞) ≥0だから,0は値の集合の下界であり,下限はそれ以上である.二つを合わせて𝑅(𝐷) =0を得る.
逆に,Q(𝐷)が空でなく𝑅(𝐷) =0とする.命題 9.1.8 より,𝐼(𝑝;𝑞∗) =𝑅(𝐷) =0を満たす𝑞∗ ∈Q(𝐷)がある.命題 1.3.2 の等号条件より,𝑞∗に対応する対(𝑋,ˆ𝑋)は独立であり,同時分布は𝑝(𝑥) ˆ𝑝(ˆ𝑥)である(ˆ𝑝は第2周辺分布).よってその期待歪みは
∑𝑥,ˆ𝑥𝑝(𝑥)ˆ𝑝(ˆ𝑥)𝑑(𝑥,ˆ𝑥)=∑ˆ𝑥ˆ𝑝(ˆ𝑥)(∑𝑥𝑝(𝑥)𝑑(𝑥,ˆ𝑥))≥∑ˆ𝑥ˆ𝑝(ˆ𝑥)𝐷max=𝐷maxである(不等号は𝐷maxが内側の和の最小値でˆ𝑝が非負なこと,最後の等号は∑ˆ𝑥ˆ𝑝(ˆ𝑥) =1による).𝑞∗ ∈Q(𝐷)よりこの期待歪みは𝐷以下だから,𝐷 ≥𝐷maxである.◼
命題 9.2.5 により,𝑅(𝐷) =0となる𝐷のうち最小のものが𝐷maxである.最小のものなのに𝐷maxという最大を思わせる名前が付くのは,これが歪みを許す意味のある範囲の右端だからである.ここから先はいくら歪みを許しても値が0のままで,見るべきものが残っていない.𝑞0のとり方が読み方を与える.情報源を見ずに,あらかじめ決めた一文字ˆ𝑥0をいつも出す,というのが「何も送らない」再現であり,その中でいちばん歪みの小さいものが払う歪みが𝐷maxである.それより多くの歪みを許してよいなら,情報源について何も知らなくて済む.
両端で確かめる. 二値の情報源と Hamming 歪み(例 9.1.3)で,二つの端を見ておく.X =ˆX ={0,1}とし,Pr[𝑋 =1] =𝜋とする.補題 9.2.2 の証明が同時分布を𝜋𝑗(𝑥,ˆ𝑥)と書いていて字が重なるが,あちらはつねに二つの引数をとるので,引数のない𝜋はつねに二値情報源のパラメータである.X =ˆXであり,𝑑𝐻(𝑥,ˆ𝑥) =0はˆ𝑥 =𝑥と同値だから,命題 9.2.4 の二つの仮定はどちらも満たされ,𝑅(0) =𝐻(𝑋) =𝐻𝑏(𝜋)である(𝐻𝑏は例 1.1.2 の二値エントロピー).右の端は,ˆ𝑥 =0に対して∑𝑥𝑝(𝑥)𝑑𝐻(𝑥,0) =𝜋,ˆ𝑥 =1に対して1 −𝜋だから,命題 9.2.5 より𝐷max =min(𝜋,1 −𝜋)である.公平なコイン(𝜋 =1/2)なら𝑅(0) =1ビット,𝐷max =1/2になる.後者は,情報源を見ずにいつも0を出せば文字の半分が食い違う,という勘定と合っている.
本節の四つを並べると,𝑅( ⋅)のグラフの形はかなり絞られる.命題 9.2.4 の仮定のもとでは,𝐷 =0で𝐻(𝑋)から始まり,𝐷が増えるにつれて非増加で(命題 9.2.1)凸(命題 9.2.3)であり,𝐷maxで0に達して,その先は0のままである(命題 9.2.5).残っているのは両端のあいだの形だけで,9.3 節は二値の情報源についてそれを決める.9.4 節はアルファベットが実数値の場合で,本節の枠の外に出て量を定義し直したうえで,同じ形の答えを求める.
InformationTheory — 形式化検証つき情報理論教科書(レビュー版).数式は MathJax + AMS Euler で事前レンダリング.