9.5 逆定理
9.4 節までで,レート歪み関数𝑅(𝐷)の値は二つの具体例について決まった.けれどもそれはまだ最適化問題の値であって,そこに符号についての意味を与えるのが本節と 9.6 節である.本節はまず,情報源を圧縮して再現するとはどういう操作かを符号として定義し,そのうえで,期待歪みを𝐷以下に抑える符号のレートが𝑅(𝐷)を下回れないことを示す.ここで初めて𝑅(𝐷)に,符号についての限界という操作的な意味が付きはじめる.
置く符号は,第2章 定義 2.3.1 のブロック情報源符号を,歪みを許す形に緩めたものになる.第6章 定義 6.2.1 のブロック通信路符号ではなくこちらに近いのは,動かす相手が違うからである.第6章の符号は,こちらが選んだメッセージを符号語に変えて通信路へ送り,出てきたものからメッセージを当てる.第2章の符号は,向こうから来た列を番号に潰し,番号から列を作り直す.本節の符号もこの後者の形で,第2章と違うのは,作り直す先が情報源アルファベットではなく再現アルファベットであることと,要求が「一致する」ではなく「歪みが𝐷以下」であることの二つだけである.
論証の形も第2章 定理 2.3.4 に近いが,一つ落ちる道具がある.ファノの不等式を使わない.第2章の逆定理は,復号に失敗する確率という一つの数を,ファノの不等式で情報量の言葉に翻訳する必要があった.本節では歪みという量が最初から𝑅(𝐷)の定義(定義 9.1.5)の中に書き込まれているので,その翻訳が要らない.そのぶん結論の形も変わり,𝑛を大きくする極限をとらずに,各𝑛の符号一つ一つについて不等式が成り立つ.
歪みを許す符号
定義 9.5.1(レート歪み符号). XとˆXを空でない有限集合,𝑛 ≥1,𝑀 ≥1とする.長さ𝑛,符号語数𝑀の レート歪み符号 とは,符号化写像𝑐 :X𝑛 →{1,…,𝑀}と復号写像𝑔 :{1,…,𝑀} →ˆX𝑛の組(𝑐,𝑔)のことをいう.その レート を
1𝑛log𝑀で定める.さらに𝑑を歪み尺度(定義 9.1.1),𝑋𝑛をX𝑛に値をとる確率変数の組とするとき,この符号の𝑋𝑛に対する 期待歪み を
𝔼[𝑑(𝑋𝑛, 𝑔(𝑐(𝑋𝑛)))]で定める(𝑑(𝑥𝑛,ˆ𝑥𝑛)は定義 9.1.2 のブロックの歪みである).
符号化写像は長さ𝑛のブロックを𝑀個の番号のどれかに潰し,復号写像はその番号から再現の列を作る.レートの読み方は第2章 定義 2.3.1 と同じで,番号を書き下すのに要るlog𝑀ビットを𝑛文字で割った,1文字あたりの取り分である.期待歪みは,情報源のブロックと,それが番号を経て戻ってきた再現のブロックとのあいだの歪みを平均した値である.定義 9.1.2 のブロックの歪みが成分についての平均だったから,この量も1文字あたりの歪みになっている.符号化写像も復号写像も確率を使わない.同じブロックには同じ番号が付き,同じ番号からは同じ再現の列が出る.
ブロックを成分に分ける
定理 9.5.4 の証明は,ブロックについての量を成分についての量の和に落とすところで進む.落とし方を先に二つ切り出しておく.一つは相互情報量についてのもので,情報源の成分が互いに独立であることだけを使う.
補題 9.5.2(相互情報量の優加法性). XとˆXを空でない有限集合,𝑛 ≥1とする.𝑋𝑛 =(𝑋0,…,𝑋𝑛−1)をXに値をとる互いに独立な確率変数の組,ˆ𝑋𝑛 =(ˆ𝑋0,…,ˆ𝑋𝑛−1)をˆXに値をとる確率変数の組とすると
𝑛−1∑𝑖=0𝐼(𝑋𝑖;ˆ𝑋𝑖)≤𝐼(𝑋𝑛;ˆ𝑋𝑛)である.
証明. 𝑋<𝑖 :=(𝑋0,…,𝑋𝑖−1)と書く.以下,𝑖 =0については条件にとる組が空であると読み,条件付きの量を無条件の量と読む(定理 1.5.2 の右辺の第0項も同じ読み方である).
まず全体を成分ごとの増分に分ける.定理 1.5.2 を,変数列𝑋0,…,𝑋𝑛−1と,相手側の変数をˆ𝑋𝑛として当てると
𝐼(𝑋𝑛;ˆ𝑋𝑛)=𝑛−1∑𝑖=0𝐼(𝑋𝑖;ˆ𝑋𝑛∣𝑋<𝑖)である.
各項を下から抑える.𝑖を一つ固定する.定理 1.4.3 を,条件を𝑋<𝑖,相手をˆ𝑋𝑛ととって当てると
𝐼(𝑋𝑖;ˆ𝑋𝑛∣𝑋<𝑖)=𝐻(𝑋𝑖∣𝑋<𝑖)−𝐻(𝑋𝑖∣𝑋<𝑖,ˆ𝑋𝑛)である.第1項は𝐻(𝑋𝑖)に等しい.実際,𝑋0,…,𝑋𝑛−1は互いに独立だから𝑋𝑖と𝑋<𝑖は独立であり,命題 1.3.2 の等号条件より𝐼(𝑋𝑖;𝑋<𝑖) =0,これと定理 1.3.4 のエントロピー表現を合わせると𝐻(𝑋𝑖 ∣𝑋<𝑖) =𝐻(𝑋𝑖)となる.
第2項を上から抑える.条件を減らしてもエントロピーは減らない(定理 1.2.4)から,条件(𝑋<𝑖,ˆ𝑋𝑛)のうちˆ𝑋𝑖だけを残して
𝐻(𝑋𝑖∣𝑋<𝑖,ˆ𝑋𝑛)≤𝐻(𝑋𝑖∣ˆ𝑋𝑖)である.定理 1.2.4 を当てる形にするには,条件をˆ𝑋𝑖とそれ以外に分けておけばよい.条件にとる組の成分を並べ替えても,条件の各値のもとでの𝑋𝑖の条件付き分布は変わらないから,定義 1.2.2 の値も変わらない.そこでˆ𝑋𝑖を先頭に,残りの成分を後ろに並べ替えたうえで,𝑌 ←ˆ𝑋𝑖,𝑍 ←残りの成分を並べた組ととればよい(𝑛 =1で残りの成分が無いときは,𝑍を自明な定数にとれば両辺が一致する).二つを合わせると,ふたたび定理 1.3.4 より
𝐼(𝑋𝑖;ˆ𝑋𝑛∣𝑋<𝑖)≥𝐻(𝑋𝑖)−𝐻(𝑋𝑖∣ˆ𝑋𝑖)=𝐼(𝑋𝑖;ˆ𝑋𝑖)を得る.𝑖について足し合わせれば主張が従う.◻
補題 9.5.2 が言っているのは,ブロックとして見たときの結びつきは,成分ごとに見たときの結びつきの和より弱くならない,ということである.独立性はちょうど一か所,𝐻(𝑋𝑖 ∣𝑋<𝑖) =𝐻(𝑋𝑖)の段で効いている.情報源の成分に相関があれば,先行する成分を知るだけで𝑋𝑖の不確かさが減るので,この等号が崩れて証明は通らなくなる.結論のほうも落ちる.𝑛 =2で𝑋0と𝑋1がつねに等しく,ˆ𝑋0 =ˆ𝑋1 =𝑋0ととると,左辺は2𝐻(𝑋0),右辺は𝐻(𝑋0)になり,𝐻(𝑋0) >0なら不等号の向きが逆になるからである.
もう一つは,成分ごとの歪みの平均にレート歪み関数を当てる段である.命題 9.2.3 の凸性は2点についてのものなので,𝑛点の形にしておく.
補題 9.5.3(レート歪み関数の𝑛点凸性). XとˆXを空でない有限集合,𝑝をX上の分布,𝑑を歪み尺度(定義 9.1.1)とし,Q( ⋅)と𝑅( ⋅)を定義 9.1.5 のとおりとする.𝑛 ≥1とし,実数𝐷0,…,𝐷𝑛−1についてQ(𝐷0),…,Q(𝐷𝑛−1)がどれも空でないとする.このとき¯𝐷 :=1𝑛∑𝑛−1𝑖=0𝐷𝑖についてQ(¯𝐷)は空でなく
𝑅(¯𝐷)≤1𝑛𝑛−1∑𝑖=0𝑅(𝐷𝑖)である.
証明. Q(𝐷)が空でない実数𝐷の全体をIと書く.Iが区間であることを見る.𝑠 ∈Iと𝑠 ≤𝑡ならば,命題 9.2.1 よりQ(𝑡)も空でないから𝑡 ∈Iである.よってIは,二つの元をとればそのあいだの実数をすべて含む集合,すなわち区間である.
¯𝐷 ∈Iであることを見る.𝐷0,…,𝐷𝑛−1の最小値を𝐷minと書くと,平均は最小値以上だから𝐷min ≤¯𝐷である.𝐷min ∈Iだから,いま見たことより¯𝐷 ∈Iであり,とくにQ(¯𝐷)は空でない.
−𝑅( ⋅)がIの上で凹であることを見る.Iの二つの元𝑎,𝑏と𝜆 ∈[0,1]をとると,Q(𝑎)とQ(𝑏)はどちらも空でないから命題 9.2.3 が当たり,𝑅(𝜆𝑎 +(1 −𝜆)𝑏) ≤𝜆𝑅(𝑎) +(1 −𝜆)𝑅(𝑏)である.両辺の符号を変えれば,これは−𝑅( ⋅)についての凹性の2点形にほかならない.
そこで補題 1.1.9 を,区間I上の凹関数−𝑅( ⋅),点𝐷0,…,𝐷𝑛−1 ∈I,重みをどれも1/𝑛ととって当てると
1𝑛𝑛−1∑𝑖=0(−𝑅(𝐷𝑖))≤−𝑅(¯𝐷)である.両辺の符号を変えて主張を得る.◻
逆定理
定理 9.5.4(レート歪み符号の逆定理). XとˆXを空でない有限集合,𝑝をX上の分布,𝑑を歪み尺度(定義 9.1.1)とし,Q( ⋅)と𝑅( ⋅)を定義 9.1.5 のとおりとする.𝑛 ≥1,𝑀 ≥1とし,𝑋𝑛 =(𝑋0,…,𝑋𝑛−1)を,どの成分も分布𝑝に従い互いに独立なXに値をとる確率変数の組とする.長さ𝑛,符号語数𝑀のレート歪み符号(定義 9.5.1)の𝑋𝑛に対する期待歪みが実数𝐷以下ならば,Q(𝐷)は空でなく
1𝑛log𝑀≥𝑅(𝐷)である.
証明. 符号を(𝑐,𝑔)と書き,ˆ𝑋𝑛 :=𝑔(𝑐(𝑋𝑛))とおく.その第𝑖成分をˆ𝑋𝑖とし,𝐷𝑖 :=𝔼[𝑑(𝑋𝑖,ˆ𝑋𝑖)]とおく.
期待歪みは成分ごとの歪みの平均である. 定義 9.1.2 より𝑑(𝑋𝑛,ˆ𝑋𝑛) =1𝑛∑𝑖𝑑(𝑋𝑖,ˆ𝑋𝑖)だから,期待値の線形性により,定義 9.5.1 の期待歪みは1𝑛∑𝑖𝐷𝑖に等しい.仮定よりこれは𝐷以下である.
各成分は制約集合の元を与える. 𝑖を一つ固定する.ˆXは空でないので,その元を一つ選んでˆ𝑥0と書く.𝑝(𝑥) >0である𝑥については𝑞𝑖(ˆ𝑥 ∣𝑥) :=Pr[ˆ𝑋𝑖 =ˆ𝑥 ∣𝑋𝑖 =𝑥]と定め,𝑝(𝑥) =0である𝑥についてはˆ𝑥 =ˆ𝑥0のとき1,そうでないとき0と定める.どちらの場合も値は非負で,各𝑥についてˆ𝑥にわたる総和は1だから,𝑞𝑖は定義 9.1.5 の条件付き分布である.𝑋𝑖の分布は𝑝だから,𝑝(𝑥) >0では𝑝(𝑥)𝑞𝑖(ˆ𝑥 ∣𝑥) =Pr[𝑋𝑖 =𝑥, ˆ𝑋𝑖 =ˆ𝑥]であり,𝑝(𝑥) =0では両辺とも0である.すなわち𝑞𝑖に対応する対の同時分布は(𝑋𝑖,ˆ𝑋𝑖)の同時分布そのものだから,定義 9.1.5 より𝐼(𝑝;𝑞𝑖) =𝐼(𝑋𝑖;ˆ𝑋𝑖)であり,その期待歪みは𝔼[𝑑(𝑋𝑖,ˆ𝑋𝑖)] =𝐷𝑖である.よって𝑞𝑖 ∈Q(𝐷𝑖)であり,とくにQ(𝐷𝑖)は空でない.𝑅(𝐷𝑖)は値の集合{𝐼(𝑝;𝑞) :𝑞 ∈Q(𝐷𝑖)}の下限だから
𝑅(𝐷𝑖)≤𝐼(𝑋𝑖;ˆ𝑋𝑖)である.
符号語数がブロックの相互情報量を抑える. ˆ𝑋𝑛の値は𝑔の値域に属する.その値域の大きさを𝑘と書くと,𝑀 ≥1だから値域は空でなく1 ≤𝑘である.ˆ𝑋𝑛は大きさ𝑘の集合に値をとる確率変数だから,定理 1.1.5 より𝐻(ˆ𝑋𝑛) ≤log𝑘である.𝑔の定義域は{1,…,𝑀}だから𝑘 ≤𝑀であり,補題 8.2.5 の単調性よりlog𝑘 ≤log𝑀である.よって𝐻(ˆ𝑋𝑛) ≤log𝑀である.また定理 1.3.4 を対(𝑋𝑛,ˆ𝑋𝑛)に当てると𝐼(𝑋𝑛;ˆ𝑋𝑛) =𝐻(ˆ𝑋𝑛) −𝐻(ˆ𝑋𝑛 ∣𝑋𝑛)であり,条件付きエントロピーは非負の量を非負の重みで平均したものだから非負である(定義 1.2.2 と命題 1.1.4).よって𝐼(𝑋𝑛;ˆ𝑋𝑛) ≤𝐻(ˆ𝑋𝑛) ≤log𝑀である.
三つの評価が一つの連鎖になる. 補題 9.5.2 を組𝑋𝑛と組ˆ𝑋𝑛に当てると∑𝑖𝐼(𝑋𝑖;ˆ𝑋𝑖) ≤𝐼(𝑋𝑛;ˆ𝑋𝑛)である.ここまでを合わせると
log𝑀≥𝐼(𝑋𝑛;ˆ𝑋𝑛)≥𝑛−1∑𝑖=0𝐼(𝑋𝑖;ˆ𝑋𝑖)≥𝑛−1∑𝑖=0𝑅(𝐷𝑖)である.¯𝐷 :=1𝑛∑𝑖𝐷𝑖とおくと,各Q(𝐷𝑖)が空でないので補題 9.5.3 が当たり,Q(¯𝐷)は空でなく𝑅(¯𝐷) ≤1𝑛∑𝑖𝑅(𝐷𝑖)である.よって1𝑛log𝑀 ≥𝑅(¯𝐷)である.最後に,期待歪みが成分ごとの歪みの平均であることから得た¯𝐷 ≤𝐷に命題 9.2.1 を当てると,Q(𝐷)は空でなく𝑅(𝐷) ≤𝑅(¯𝐷)である.二つを合わせて1𝑛log𝑀 ≥𝑅(𝐷)を得る.◼
左辺が誰のものかを見る. 左辺の1𝑛log𝑀は符号を設計する側が決める量で,右辺の𝑅(𝐷)は情報源の分布と歪み尺度だけで決まる量である.したがって定理 9.5.4 は,設計する側が下げられるところに,情報源の側から下限が置かれている形に読める.裏返して言えば,レートを𝑅(𝐷)より下げた符号は,期待歪みを𝐷以下に抑えられない.レートを下げたければ𝐷を大きくとる,すなわち歪みを多く許すことになり,そのぶん右辺は上がらない(命題 9.2.1).第2章 定理 2.3.4 と違うのは,ここに極限が出てこないことである.あちらは誤り確率が0に向かう族についてしか言えず,結論も下極限の形だった.こちらは各𝑛の符号1本ごとに不等式が成り立ち,𝑛を大きくしても限界は緩まない.違いを生んでいるのは,歪みが「𝐷以下」という形で最初から符号に課された条件になっていることである.
数で見る
例 9.5.5(𝜋 =1/4の二値情報源の符号語数). X =ˆX ={0,1},𝑝(1) =1/4,𝑝(0) =3/4,𝑑 =𝑑𝐻(例 9.1.3)とし,𝑋100を,どの成分も分布𝑝に従い互いに独立な100個の確率変数の組とする.長さ100のレート歪み符号(定義 9.5.1)の符号語数を𝑀とし,𝑋100に対する期待歪みを考える.期待歪みが0ならばlog𝑀 ≥100 𝑅(0)であり,右辺は約81.13である.期待歪みが0.05以下ならばlog𝑀 ≥100 𝑅(0.05)であり,右辺は約52.49である.さらに,期待歪みが0ならばlog𝑀 ≥100である.
証明. min(1/4, 3/4) =1/4であり0 ≤0.05 <1/4だから,定理 9.3.4 の第1の場合が𝐷 =0と𝐷 =0.05のどちらにも当たり,Q(0)とQ(0.05)はどちらも空でない.期待歪みが0のときは定理 9.5.4 を𝑛 =100,𝐷 =0に当て,期待歪みが0.05以下のときは同じ定理を𝑛 =100,𝐷 =0.05に当てて,どちらも両辺を100倍すれば主張の不等式を得る.二つの右辺の値は,例 9.3.5 が求めた𝑅(0)と𝑅(0.05)の値のちょうど100倍である.
最後の主張は定理 9.5.4 を通さずに直接見る.符号を(𝑐,𝑔)と書く.歪み尺度の値は非負(定義 9.1.1)だからブロックの歪み𝑑(𝑋100,𝑔(𝑐(𝑋100)))も非負であり,その期待値は非負の値を非負の重みで足し合わせたものだから,期待歪みが0になるのは,正の確率をもつどの列でもブロックの歪みが0のときである.𝑝(0)と𝑝(1)はどちらも正だから𝑋100はどの列も正の確率でとり,よってどの𝑥100 ∈X100についてもブロックの歪みは0である.𝑑𝐻の値が0になるのは二つの文字が一致するときに限る(例 9.1.3)から,成分ごとに一致して𝑔(𝑐(𝑥100)) =𝑥100である.したがって𝑐は単射で𝑀 ≥2100であり,補題 8.2.5 の単調性よりlog𝑀 ≥100である.◼
規模感. 例 9.5.5 の情報源を素朴に1文字1ビットで書けば,長さ100のブロックに番号は2100通り要り,log𝑀は100である.歪みをいっさい許さないなら,例 9.5.5 よりlog𝑀は約81.13を下回れない.ただし同じ例が示しているとおり,この場合に実際に要るのはlog𝑀 =100で,素朴なやり方から1ビットも減らせない.定理 9.5.4 の下界は,各𝑛で最良の値を与えるものではない,ということである.第2章の下界にも同じ開きがあった.定理 2.3.4 が与えるのは,誤り確率を0に向かわせる族について𝐻(𝑋)という値までで,誤りをいっさい許さない族についても,それより強い下界は与えない.歪みを0ちょうどに固定したときに定理 9.5.4 が𝑅(0)までしか与えないのと,同じ格好である.文字の5パーセントの食い違いまで許すと,限界が約52.49まで下がり,素朴なやり方の半分ほどになる.ただし例 9.5.5 が言っているのは限界だけで,この本数の符号語で実際に歪みを0.05以下に抑えられるかは,まだ何も言っていない.下界と実際に要る量の開きが縮むかどうかも,長さを大きくしたときに何が言えるかによる.
定理 9.5.4 は片側だけの主張である.𝑅(𝐷)より小さいレートでは期待歪みを𝐷以下に抑えられないことは言えたが,𝑅(𝐷)より大きいレートで実際に抑えられるかは,まだ問うていない.それが 9.6 節の内容で,そこで初めて𝑅(𝐷)の両側がそろう.
InformationTheory — 形式化検証つき情報理論教科書(レビュー版).数式は MathJax + AMS Euler で事前レンダリング.