9.6 達成可能性

定理 9.5.4 は,期待歪みを以下に抑える符号のレートを下から押さえた.残るのは逆向きで,より大きいレートなら実際に符号が作れるか,という問いである.第2章 定理 2.3.2第6章 6.3 節で答えたのと同じ形の問いで,答えも同じく肯定である.本節はその達成可能性を借用として置く.借りる形には歪みの側に小さな余裕が付いていて,レートをより大きくとれば,どのについても,長さを十分大きくとって期待歪みを以下に抑えられる,という形をしている.そのうえで,期待歪みがより小さい再現の作り方があるときには,この余裕を落として期待歪みを以下に抑えられることを見る.

記号を一つ断っておく.本節は符号のレートの値をと書く.本章は定義 9.1.5 のレート歪み関数をつねにのように引数を添えて書くので,引数の有無で二つを見分けられる.

本節はまず,なぜという数がここに現れるのかを,再現の列で情報源の典型系列を覆うという勘定で見る.そのうえで達成可能性を借用として置き,逆定理と突き合わせて本章の主定理に至る.次に長さの符号を二つ作って数を見る.最後に,第6章と本章で得た二つの符号化定理を並べて,向きの違いを見る.

覆う勘定

第2章の圧縮は,典型集合(定義 2.2.1)に番号を振ることで働いた.情報源のブロックはほぼ確実に典型集合に落ち(定理 2.2.3),その元の個数はをわずかに上回るところまでで抑えられる(定理 2.2.5)から,番号を書き下すのに文字あたりビットあれば足りた.歪みを許すと,番号の数をもっと減らせるはずである.一つの番号に一つの系列を対応させるのではなく,一本の再現の列に,そこから歪みの小さい情報源の列をまとめて引き受けさせればよい.

ここから先は見当である.条件付き分布を一つ固定し,対をその同時分布に従うものとして,その独立な複製を組並べたものをと書く.再現の列を一本決めたとき,それが引き受けられる情報源の列は,を知ったあとの世界で典型な列だから,その個数はおよそである,と読む.肩にが来るのは,定理 2.2.5 の勘定をこの条件つきの世界に当てるからである.あちらで肩に乗っていたのは情報源の文字あたりの不確かさで,いまはを知ったぶんだけそれが減り,平均するとになる.覆いたい相手はおよそ個あるから,要る本数はその比

2𝑛𝐻(𝑋)2𝑛𝐻(𝑋ˆ𝑋)=2𝑛𝐼(𝑋;ˆ𝑋)

ほどになる(指数の差が相互情報量になるのは定理 1.3.4 による).本数を少なくしたければを小さくすればよく,歪みの制約のもとでいちばん小さくした値がにほかならない(定義 9.1.5).本書はこの勘定を主張としては述べない.「を知ったあとの典型な列」にあたる集合を,本書は用意していないからである.一本の再現の列が実際に何本を引き受けるかは,長さの場合を例 9.6.2 で数える.

第8章 8.4 節の勘定と,数の出方が同じであることも見ておく.あちらは,半径の球の中に半径の球をいくつ詰められるかを,体積の比として数えた(命題 8.4.1).こちらは,情報源の典型集合が一本の再現の列でどれだけ覆えるかを,個数の比として数えている.違うのは,詰めるのか覆うのかという向きだけである.

借りる道具

解析の事実を一つ借りる.区間の上で定義された凸関数は,その区間の内点で連続であるという事実である.当てる相手は,実数の区間の上の実数値関数として見ただけで,この借用に依存するのは系 9.6.1 の第の主張の証明だけである.本書はこの事実を証明しないが,形式化されていないわけではない.Mathlib に無条件の機械検証済みの定理として置かれている.

道具をもう一つ借りる.レート歪み理論の達成可能性,すなわち「レートがを上回っているかぎり,期待歪みをにいくらでも近いところまで抑える符号がある」という主張である.使う形を書いておく.

を空でない有限集合,上の分布,を歪み尺度(定義 9.1.1)とし,定義 9.1.5 のとおりとする.実数についてが空でないとし,実数を満たし,とする.このとき番号があって,を満たすどのについても,符号語数を満たす長さのレート歪み符号(定義 9.5.1)で,どの成分も分布に従い互いに独立なに値をとる個の確率変数の組に対する期待歪みが以下であるものが存在する.

当てる対象は,情報源アルファベットと再現アルファベットがどちらも空でない有限集合であるときの,定義 9.5.1 のレート歪み符号だけである.情報源の分布には,上の分布であること以外の条件を置かない.この借用を主張の根拠として当てるのは系 9.6.1 の証明だけである.第14章 14.4 節がこの名前を二度引くが,どちらも当てているのではない.一度は 定義 14.4.2 が歪みにの余裕を付ける形をここと同じだと言うところ,もう一度は同節が借りる「副情報つきレート歪みの達成可能性」の筋書きが,覆いの段をここと同じ作り方だと言うところである.

本書がこの主張を証明しないのは,第6章 6.3 節のランダム符号化に使う部品のうち,二つを用意していないからである.一つは,結合典型集合(定義 6.2.3)に落ちる確率を下から抑える評価である.第6章が要ったのは,無関係な対がたまたま結合典型になる確率を上から抑えること(定理 6.2.7)で,そちらは誤りを小さくするために使う.覆う議論で要るのはその逆向き,すなわち独立に引いた再現の列が情報源の列と結合典型になる確率を下から抑えることで,一本も当たらない確率をの形で押さえるのに使う.本書はこの下界を持っていない.もう一つは,結合典型であることからブロックの歪みが期待歪みの近くにあることを引き出す段である.これには第2章 2.4 節の強典型性(定義 2.4.1)を,一つの情報源についてではなく対について立てる必要があり,本書はその形を用意していない.「本書で証明しない」ことと「形式化されていない」ことは別である.本節が借りた主張は,本書が証明を載せないだけで,無条件の機械検証済みの定理として形式化されている.

借りたままにするので,中で何が起きているかの筋書きだけ書いておく.符号帳は設計するのではなく,くじで引く.命題 9.1.8 によりの上で最小にする条件付き分布をとり,それが定める再現の側の分布をと書く.一本の再現の列は,各成分を独立にから引いて作り,それを本ぶん独立に繰り返す.符号化は,情報源の列に対して,対歪み典型 である符号語を探し,見つかればその番号を出し,一本も見つからなければ番号を出す,というものである.ここで歪み典型とは,結合典型(定義 6.2.3の定める同時分布に当てたもの)であって,かつブロックの歪み(定義 9.1.2)が期待歪みの近くにあることをいう.第6章 6.3 節の復号器が結合典型な符号語を探したのと,探す規則は同じで,探す側が復号から符号化へ移っている.

見当のとおり,一本の符号語が当たる確率はおよそで,本のうち一本も当たらない確率はおよそである.ほどあってなら,これはとともにに向かう(のときであり,𝑥 =2𝑛𝑅(𝐷)とおくと肩のが発散するからである).当たったときのブロックの歪みは,歪み典型の定め方から期待歪みの近く,すなわちの近くにある.当たらなかったときの歪みは,アルファベットが有限なのでで抑えられ,その場合の確率がに向かうので,期待歪みへの寄与もに向かう.二つを合わせると,を十分大きくとれば期待歪みは以下に収まる.最後に,引いた符号帳について平均した期待歪みがこの評価を満たすので,少なくとも一つの符号帳が同じ評価を満たす.第6章 6.3 節も最後に同じ形の議論を通った.

形式化: rate_distortion_achievability_operational_general (ソース)

形式化上の注記. 形式化の宣言の仮定は,借りた主張と一つずつ対応している.情報源の分布が確率分布であること,制約集合が空でないこと(本文の「が空でない」にあたる),レートがより真に大きいこと,が正であることの四つで,これがすべてである.結論も同じ形で,ある番号から先のすべての長さについて,符号語数が上の不等式を満たすレート歪み符号があり,その期待歪みが以下である,という形をしている.ただし宣言は符号語数について,を切り上げた整数以上であるという下からの評価も与える.本文が使うのは符号語数が正であることだけで,それは定義 9.5.1 から出るので,借りる形からはこの評価を落としてある.筋書きに書いた歪み典型は,あちらでは distortionTypicalSet (InformationTheory/Shannon/RateDistortion/AchievabilityJointTypicalEncoder.lean) として,結合典型集合と,ブロックの歪みについての条件との共通部分に定めてある.ただし探索の規則は本文の筋書きと違う.形式化の符号化器 jointStronglyTypicalLossyEncoder (InformationTheory/Shannon/RateDistortion/AchievabilityJointStrongTypicality.lean) が探すのは,強典型性(定義 2.4.1)の意味で結合典型な符号語である.ブロックの歪みが期待歪みの近くにあることは,探索の規則からではなく,そこから distortionTypicalSet への包含 jts_subset_dts_of_dist_slack (InformationTheory/Shannon/RateDistortion/AchievabilityUnconditional.lean) を経て出る.

レート歪み定理

これから置く系のうち,達成可能性の側の証明は二段の余裕とりで動くので,筋を先に言っておく.借りた達成可能性が返してくるのは期待歪み以下の符号なので,その余裕をあらかじめ差し引いておく.すなわちより少し小さい歪みに借用を当て,余裕のほうをにとれば,合わせてちょうどになる.そう縮めても借用が使える,つまりがなお成り立つ,というところに借りた凸関数の連続性が要る.レートの側にも同じ形の差し引きがあって,借用のレートをより少し小さいにとることで,符号語数を整数にするぶんの増加を吸収する.

系 9.6.1(レート歪み定理). を空でない有限集合,上の分布,を歪み尺度(定義 9.1.1)とし,定義 9.1.5 のとおりとする.実数についてが空でないとし,各についてを,どの成分も分布に従い互いに独立な個のに値をとる確率変数の組とする.実数について次の二つが成り立つ.

  1. ならば,どのについても,レートが以下でに対する期待歪みが以下である長さのレート歪み符号(定義 9.5.1)は存在しない.
  2. で,かつが空でない実数が存在するならば,番号があって,を満たすどのについても,レートが以下でに対する期待歪みが以下である長さのレート歪み符号が存在する.

証明.

  1. そのような符号があったとし,その符号語数をとする.期待歪みが以下だから定理 9.5.4 が当たりである.いっぽうレートが以下という仮定はだからとなり,に反する.

  2. まず,が空でないどの実数についてもであることを見る.定義 9.1.5は値の集合の下限であり,命題 1.3.2 よりどの値も非負だからはこの集合の下界である.下限は下界のうち最大のものだからである.

歪みの側に余裕を作る.仮定よりが空でない実数がある.命題 9.2.1 より以上のどの実数についてもは空でないから,は区間の上で定まり,命題 9.2.3 よりその上で凸である.だからはこの区間の内点であり,区間の上で定義された凸関数はその区間の内点で連続だから,で連続である.にこの連続性を当てると,があって,の点を満たすかぎりである.そこでとおくとかつであり,とおくとである.

借りたレート歪み理論の達成可能性を,歪みの上限,レート,余裕に当てて番号を得る.を満たす番号をとり,とおく.とすると,符号語数を満たし,に対する期待歪みが以下である長さのレート歪み符号がとれる.

レートを評価する.よりだからである.また定義 9.5.1 よりである.補題 8.2.5 よりは単調だから,の対数をとってで割ると

1𝑛log𝑀𝑅+1𝑛

である.よりだから,レートは以下である.

形式化上の注記. 系 9.6.1 に対応する単独の宣言は無い.しかも片側ずつの機械検証は,同じについてのものになっていない.第の主張に使った達成可能性の宣言が下限をとる相手は,定義 9.1.5 に紐付けた rateDistortionFunctionPmf (InformationTheory/Shannon/RateDistortion/Achievability.lean) である.いっぽう第の主張が経由する定理 9.5.4 の宣言が使うのは,を測度の言葉で書いた rateDistortionFunction (InformationTheory/Shannon/RateDistortion/Converse.lean) で,値を拡張非負実数にとり,相互情報量を相対エントロピーの形で書く.二つの値はどちらも定義 9.1.5にあたるが,単独の宣言としては別のもので,両者を結ぶ宣言は無い.したがって,本文が一つのについて述べている両側の主張は,形式化では別々の量についての二つの主張になっている.

系 9.6.1 の二つはを境として向かい合っている.第の主張は,レートがより小さければ期待歪みを以下に抑える符号がどの長さにも無いと言う.第の主張は,より大きく,しかも期待歪みがより小さい再現の作り方があるなら,長さを十分大きくとってそういう符号がとれると言う.一行で言えば,歪みを許した圧縮に要るレートの限界はで,境目のすぐ上とすぐ下とで符号の有無が入れ替わる.第6章 定理 6.4.9第8章 系 8.4.2 で達成可能性と逆定理が噛み合ったのと,同じ形である.

の主張だけが余分に置いた条件は,より小さい歪みで済む再現の作り方がある,ということである.借りた達成可能性がそのまま与えるのは期待歪みが以下までで,そこからを落とすのに使ったのが,この条件と,の非増加性(命題 9.2.1)・凸性(命題 9.2.3)と,借りた凸関数の連続性である.歪みをいっさい許さないは,どの情報源でもこの条件から外れる.歪み尺度の値は非負(定義 9.1.1)だから,ではが空だからである.外したのは書きぶりの都合ではない.例 9.5.5 の情報源で,レートをのあいだにとってみる(例 9.3.5 よりは約ビットだから,そういうレートはある).期待歪みをにするには,例 9.5.5 の最後の段と同じ理由で符号化写像が単射でなければならず,符号語数は以上,すなわちレートは以上になる.長さを何にとってもそうだから,どれだけ長くとってもそういう符号は存在しない.

数で見る

例 9.6.2(の二値情報源の長さの符号). X =ˆX ={0,1}𝑝(1) =1/4𝑝(0) =3/4𝑑 =𝑑𝐻例 9.1.3)とし,を,どちらの成分も分布に従い互いに独立な確率変数の組とする.長さのレート歪み符号(定義 9.5.1)を二つ考える.

  1. 符号語数で,復号写像が番号を返すもの.そのレートはで,に対する期待歪みはである.
  2. 符号語数で,復号写像が番号,番号を返し,符号化写像がに番号を,残りの三つの列に番号を割り当てるもの.そのレートはビットで,に対する期待歪みはである.

さらに,定理 9.5.4 が期待歪みに置く下界は約ビットであり,第の符号のレートはその倍を超える.

証明. 成分が互いに独立だから,の分布はPr[𝑋2 =(0,1)] =Pr[𝑋2 =(1,0)] =3/16である.また例 9.1.3 より,長さのブロックの歪みは成分が食い違う位置の割合,すなわち食い違った成分の個数の半分である.

  1. レートはである.再現の列はの一本だけだから,ふたたび例 9.1.3 より期待歪みはである.

  2. レートはビットである.四つの列それぞれについてブロックの歪みを数える.は自分と同じ再現の列に移るから歪みはは番号を経てに戻り,食い違う成分が一つだから歪みはである.よって期待歪みはである.

最後の主張を見る.より,定理 9.3.4 の第の場合がに当たるから,は空でなく,その値も同じ場合の式で直接求まる.log23 =1.5849だからであり

𝐻𝑏(316)=316×2.4150+1316×0.2995=0.4528+0.2433=0.6962

である.例 9.3.5そのもので約ビットだから,は約ビットである.レートをこの値で割るとだから,レートは下界の倍を超えている.

形式化上の注記. 例 9.6.2 に対応する宣言は形式化されていない.形式化には歪み尺度を具体的にとった実例が一つも置かれていないので,符号を一つ具体的に作ってその期待歪みを述べた宣言も無い.

二つを見比べると,覆うという言葉に実体が入る.第の符号は再現の列を一本しか持たず,その一本がの四つの列すべてを引き受けている.レートはで番号を書き下す手間がいっさい要らない代わりに,が出るたびにそのぶんの歪みを払うので,期待歪みはの出る割合そのものになる.第の符号はをもう一本足し,四つのうちの一つだけをそちらへ移した.一本が引き受ける本数が減ったぶん期待歪みはまで下がり,代わりにレートがビットに上がる.レートと歪みの交換は,こういう形で起きている.

の符号のレートは,定理 9.5.4 が期待歪みに置く下界にちょうど等しい(例 9.3.5 より).歪みをまで許すなら,長さでも符号語一本で足りるということである.第の符号のレートのほうは,例 9.6.2 の最後の主張のとおり,下界倍を超えたところにある.長さでは再現の列の置き場所が四か所しかないので,この粗さが残る.この情報源では,より定理 9.3.4 の第の場合がにも当たってが空でなく,しかもだから,系 9.6.1 の第の主張がに当たる.長さを大きくとれば,レートをのすぐ上まで下げても期待歪みを以下に抑えられる,ということである.

覆うことと詰めること

第6章と本章で符号化定理を二つ得た.二つを並べると,同じ道具立てが向きを変えて使われていることが見える.以下は前半の勘定と同じく見当の言葉で書く.記号を一つ置く.第6章の通信路の入力を,出力をと書く.本章のはどちらも情報源の側の記号なので,通信路の側には別の字を当てる.例 7.3.5 が一様分布の確率変数に使ったとは別である.

通信路符号化(第6章 定理 6.4.9)が作るのは,入力アルファベットの列に置いた符号語の集まりである.一本の符号語を送ると,出てくる列はその符号語と結合典型な列のどれかで,およそ個の広がりをもつ.二本の符号語の広がりが重なると復号を誤るから,設計の目標は符号語をなるべく重ならないように 散らす ことである.全体の入れ物である出力側の典型集合はおよそ個だから,散らせる本数はその比ほどで,本数を多くしたければを大きくすればよい.入力分布を動かしてを最大にした値が容量である(定義 6.1.4).

情報源符号化(系 9.6.1)が作るのは,再現アルファベットの列に置いた再現の列の集まりである.一本の再現の列は,それと歪み典型な情報源の列を引き受け,およそ個の広がりをもつ.引き受けそこねた列があるとそこで歪みが大きくなるから,設計の目標は情報源の典型集合をなるべく少ない本数で 覆う ことである.覆う相手はおよそ個だから,要る本数はその比ほどで,本数を少なくしたければを小さくすればよい.条件付き分布を動かしてを最小にした値がレート歪み関数である(定義 9.1.5).

三つの向きがそろって裏返っている.散らすことと覆うこと,最大化と最小化,そして符号のレートに対する不等号の向きである.通信路では,達成レートが容量を超えられない(系 6.4.8)ので容量が上からの限界になり,情報源では,レートがを下回れない(定理 9.5.4)のでが下からの限界になる.相互情報量を肩に乗せた同じ形の数が,一方では詰められる本数の上限として,他方では覆うのに要る本数の下限として現れている.

本書は,この対比が定理として成り立つことを示さない.示したのは二つの符号化定理であって,一方から他方を導く主張ではない.並べて見えるのは,二つの証明が同じ部品(典型集合と相互情報量)を使っていて,最適化の向きだけが違う,という事実までである.

形式化上の注記. 形式化の側でも二つの定義は鏡になっていて,容量 capacity (InformationTheory/Shannon/ChannelCoding/ShannonTheorem.lean) が相互情報量の上限,レート歪み関数 rateDistortionFunctionPmf (InformationTheory/Shannon/RateDistortion/Achievability.lean) が相互情報量の下限として書かれている.ただし上限と下限をとる相手は別の宣言である.あちらは通信路についての相互情報量 mutualInfoOfChannel (InformationTheory/Shannon/ChannelCoding/Basic.lean),こちらは同時分布についての相互情報量 mutualInfoPmf (InformationTheory/Shannon/RateDistortion/Achievability.lean) で,二つを同じ主張の中で結ぶ宣言は無い.

第1章で分布から決まる数として定義した相互情報量が,これで二度,符号についての意味を得たことになる.第6章では入力分布を動かして最大にした値が,回の使用あたり運べるビット数の限界として現れた.本章では条件付き分布を動かして最小にした値が,文字あたり記録に要するビット数の限界として現れた.運ぶことと記録することは別の操作だが,限界を決めている量は同じ一つである.

InformationTheory — 形式化検証つき情報理論教科書(レビュー版).数式は MathJax + AMS Euler で事前レンダリング.