6.1 通信路と通信路容量

第2章と第3章は,情報源のもつ冗長さを取り除く話だった.同じ内容をより短く書けるのはどこまでか.その限界がエントロピーだった.本章で扱うのは逆向きの操作である.通信路は送った記号を確率的に取り違えるので,受け取った側が誤りを直せるように,送る側があらかじめ冗長さを足す.足しすぎれば運べる情報が減り,足りなければ誤りが残る.この綱引きの決着がどこにつくかが本章の主題である.

素朴に考えると,雑音のある通信路では誤り確率を 0 に近づけようとするたびに伝送速度が0 に落ちそうに思える.同じ記号を何度も繰り返して多数決をとれば誤りは減らせるが,繰り返した回数だけ速度は落ちる.Shannon の発見は,そうはならないということである.通信路ごとにある数が定まっていて,より小さい速度である限り,誤り確率はいくらでも 0 に近づけられる.本節ではその(通信路容量)を定義し,それが上限ではなく最大値として達成されることを示し,二つの通信路で値を計算する.

本章ではの底を 2 にとる.第2章・第3章と同じくの形の量を扱うので,指数と対数の底をそろえておくと式が読みやすい.

形式化上の注記(本章共通)

単位は本文と形式化で違う.形式化はを自然対数にとるので,相互情報量も容量もナットが単位であり,本文がと書く量はあちらではにあたる.底をそろえれば同じ主張である.

形式化の Channel は条件付き分布の族そのもの(各入力に対して出力の分布を返す核)であって,各入力についての総和が 1 という条件は含まない.その条件は,Channel を使う主張の側に別の前提として付いている.

通信路と,その相互情報量

定義 6.1.1(離散無記憶通信路). 入力アルファベットと出力アルファベットはどちらも有限とし,は空でないとする.通信路 とは,各に対して上の分布を与える対応である.すなわちであって,各についてを満たす.長さの入力に対する出力の分布が

𝑊𝑛(𝑦𝑛∣𝑥𝑛)=𝑛−1∏𝑖=0𝑊(𝑦𝑖∣𝑥𝑖)

で与えられるとき,この通信路を 記憶のない(memoryless)通信路 という.アルファベットが有限なので,これを 離散無記憶通信路 とも呼ぶ.上の分布を入力に与えたときの入力と出力の 結合分布 を,その第 2 成分の周辺分布を出力分布 と呼ぶ.

は「を送ったときにが受け取られる確率」である.通信路が定めるのはこの条件付き分布だけで,何を送るかは通信路の側では決まっていない.入力分布を与えてはじめて入力と出力が対になり,結合分布と出力分布が定まる.積の形の条件が言っているのは,各時刻の雑音がそれ以前の入出力に一切依存しない,ということである.これが「記憶のない」の内容で,第3章の記憶のある情報源とちょうど逆向きの仮定になっている.あちらは情報源に記憶があってよかったが,こちらは通信路に記憶がない.

形式化: Channel,結合分布 jointDistribution,出力分布 outputDistribution (ソース)

いちばん小さい通信路で確かめる. とし,送ったビットが確率で反転する通信路を考える.遷移確率は

𝑊(0∣0)=𝑊(1∣1)=0.9,𝑊(1∣0)=𝑊(0∣1)=0.1

である.同じ入力に対する二つの値の和が 1 になっていることが定義 6.1.1 の条件で,同じ出力に対する和には何の条件もない.通信路が定めるのは入力ごとの分布だからである.ここに入力分布を与えると,結合分布はとに,とにを置き,出力分布はになる.この通信路は,反転する確率を一般の値にした形で例 6.1.8 にふたたび現れる.

定義 6.1.2(通信路の相互情報量). 上の分布と通信路に対し,出力分布をとおく.入力分布のもとでのの相互情報量 を

𝐼(𝑝;𝑊):=∑𝑥,𝑦𝑝(𝑥)𝑊(𝑦∣𝑥)log⁡𝑊(𝑦∣𝑥)𝑞(𝑦)

で定める.和はの項についてとる.

対数の中身は,「が送られたと知ることで,が受け取られるという見込みが何倍になったか」である.何を送っても出力の分布が変わらない通信路(雑音が入力を完全に潰してしまう通信路)では,この比が全点で 1 になってになる.逆にごとに出力の分布がはっきり違えば比は 1 から離れ,は大きくなる.1.3 節の言葉でいえば,は入力と出力が独立からどれだけ隔たっているかを一つの数にまとめた量である.

形式化: mutualInfoOfChannel (ソース)

定義 6.1.2 は通信路の側から書いた式で,第1章の相互情報量とは見た目が違う.両者が同じものであることを先に確かめておく.

命題 6.1.3. 上の分布と通信路をとり,を結合分布に従う対とする.このとき

𝐼(𝑝;𝑊)=𝐼(𝑋;𝑌)=𝐻(𝑋)+𝐻(𝑌)−𝐻(𝑋,𝑌)=𝐻(𝑌)−𝐻(𝑌∣𝑋)=𝐻(𝑋)−𝐻(𝑋∣𝑌).

証明. の同時分布はであり,その周辺分布はとである.和をとる範囲ではなので約分してよく,定義 1.3.1 の対数の中身は

𝑝(𝑥,𝑦)𝑝(𝑥)𝑞(𝑦)=𝑝(𝑥)𝑊(𝑦∣𝑥)𝑝(𝑥)𝑞(𝑦)=𝑊(𝑦∣𝑥)𝑞(𝑦)

になる.したがって定義 6.1.2 の右辺は定義 1.3.1 の右辺そのもので,である.残る三つの等式は,この対に定理 1.3.4 を当てたものである.◼

形式化: 結合分布の相互情報量との一致 mutualInfoOfChannel_eq_mutualInfo_prod,三つのエントロピーの差の形 mutualInfoOfChannel_eq_HX_add_HY_sub_HZ (ソース),出力の不確かさから雑音の分を引いた形 mutualInfoOfChannel_toReal_eq_outputEntropy_sub (ソース)

注記

最後の等式に対応する単独の宣言はない.定理 1.3.4 の形式化との合成で得られる.

雑音の分だけ減る,と読む. 以下でいちばん使うのはの形である.は受け取った側から見た出力の不確かさ,は何を送ったかを知ってもなお残る不確かさ,すなわち雑音そのものの不確かさである.運べた情報とは,出力の不確かさから雑音の分を差し引いた残りだ,と読める.送る側に選べるのはだけで,を動かせばもも動く.その中で差を最大にするを探す,というのが次の定義である.

通信路容量

定義 6.1.4(通信路容量). 通信路の 通信路容量 を,入力分布全体にわたる上限

𝐶(𝑊):=sup𝑝𝐼(𝑝;𝑊)

で定める.は上の分布全体,すなわちかつを満たす関数の全体を動く.

形式化: capacity (ソース)

は通信路だけで決まる数で,入力分布にも符号にも依存しない.定義の時点では「1 回の使用あたりどれだけ運べるか」という操作的な意味はまだ持っておらず,単に相互情報量の最大化問題の値でしかない.それが誤りなく送れるレートの限界と一致する,というのが 6.3 節以降の内容である.

上限と書いたが,実際には最大値である.それを示すには,が入力分布のつくる集合の上で連続であることと,その集合が有界閉集合であることの二つが要る.前者は,一次結合・合成・有限和が連続性を保つという微積分の計算規則と,補題 1.1.8 が示したの上の連続性から下で組み立てる(計算規則のほうは既知とする).後者から最大値の存在を引き出すために,次の一つを証明せずに借りる.

Weierstrass の最大値定理を借りる. 借りるのは「有限次元の実ベクトル空間の空でない有界閉集合の上の実数値連続関数は最大値をとる」という形である.当てる相手は,実ベクトル空間の部分集合である確率単体

Δ:={𝑝:X→ℝ∣𝑝(𝑥)≥0 (∀𝑥∈X), ∑𝑥𝑝(𝑥)=1}

の上の実数値関数である.この借用に依存するのは定理 6.1.5 と,同じ最大値定理を 9.1 節で借り直している第9章 命題 9.1.8,11.5 節で借り直している第11章 命題 11.5.7・定理 11.6.2,および 15.2 節で借り直している第15章 命題 15.2.2・命題 15.8.5 である.本章の以降の節と,第12章 定理 12.3.4 が事前分布について最大値をとるところは,定理 6.1.5 の結論だけを使う.借り直した先が当てる相手は,確率単体そのものではない.第9章は実ベクトル空間の部分集合である制約集合に,第11章はの閉区間と,確率単体の閉部分集合とに,第15章は銘柄の番号の上の確率単体,すなわちの部分集合であるポートフォリオ全体に当てる.本書はこの最大値定理を証明しないが,形式化されていないわけではない.定理 6.1.5 の形式化は,Mathlib にある無条件の機械検証済みのこの定理をそのまま呼び出しているからである.「本書で証明しない」ことと「形式化されていない」ことは別である.

定理 6.1.5. を通信路(定義 6.1.1)とし,上の分布全体のなす集合をと書く.このときはの上で連続であり,を最大にする入力分布が存在してである.

証明. 連続性から示す.を結合分布に従う対とすると,命題 6.1.3 よりである.三つの項をの関数として書き下すと

𝐻(𝑋)=∑𝑥𝜑(𝑝(𝑥)),𝐻(𝑌)=∑𝑦𝜑(∑𝑥𝑝(𝑥)𝑊(𝑦∣𝑥)),𝐻(𝑋,𝑌)=∑𝑥,𝑦𝜑(𝑝(𝑥)𝑊(𝑦∣𝑥))

となる.どの項も,の成分の一次結合にを合成したものの有限和である.一次結合は連続で,は上で連続だから,はの上で連続である.

最大値の存在に移る.はの有界閉集合である(有界性はから,閉であることは定義が等式と不等式で書けていることから従う).よって借用した Weierstrass の最大値定理により,はの上で最大値をとる.それを与える点をとする.

最後にを見る.は値の集合の上界であり,かつその集合の元でもある.上限は上界のうち最小のものだから,この二つからである.◼

形式化: 連続性 continuous_mutualInfoOfChannel_left,最大化する入力分布の存在exists_capacity_achiever (ソース)

注記

等号に対応する単独の宣言はない.容量が上限として定義されていることcapacity と,その値の集合が上に有界であること capacity_bddAbove (InformationTheory/Shannon/ChannelCoding/ShannonTheorem.lean) との合成で得られる.

命題 6.1.6. 任意の通信路に対し.

証明. 定理 6.1.5 のをとるとであり,命題 6.1.3 よりこれは対の相互情報量に等しい.命題 1.3.2 よりそれは非負である.◼

形式化: capacity_nonneg (ソース)

命題 6.1.7. 任意の通信路に対し.

証明. 任意の入力分布をとり,を結合分布に従う対とする.命題 6.1.3 よりであり,条件付きエントロピーは非負量の平均なので非負(定義 1.2.2 と命題 1.1.4).よってで,定理 1.1.5 よりである.同じことを命題 6.1.3 のもう一つの形で行えばを得る.は任意だったので,は値の集合の上界であり,上限はそれ以下である.◼

読み方は素直である.1 回の使用で運べる情報は,入れ物の大きさを超えられない.送れる記号が通りしかなければ,雑音がまったくなくても 1 回あたりより多くは運べないし,受け取る側の区別が通りしかないなら同じことが出口の側で起きる.命題 6.1.7 は,容量が雑音の少なさだけでなく,入口と出口の狭いほうにも押さえられていることを言っている.

形式化上の注記

命題 6.1.7 に対応する宣言は形式化されていない.命題 6.1.7 に付した証明が,この主張の保証のすべてである.

二つの通信路で計算する

例 6.1.8(二元対称通信路). を固定する.とし,

𝑊(1∣0)=𝑊(0∣1)=𝜌,𝑊(0∣0)=𝑊(1∣1)=1−𝜌

で定まる通信路を,反転確率 の 二元対称通信路 という.送ったビットが確率で反転し,確率でそのまま届く通信路である(図 6.1.1).その容量は

𝐶(𝑊)=1−𝐻𝑏(𝜌)

であり,一様入力で達成される(は例 1.1.2 の二値エントロピー).さらに,どの入力分布についても,とおくと

𝐼(𝑝;𝑊)=𝐻𝑏((1−𝜌)𝜋+𝜌(1−𝜋))−𝐻𝑏(𝜌)

である.

図 6.1.1. 反転確率の二元対称通信路.左が入力,右が出力で,矢印に添えた数がその遷移の確率である.

証明. まずが入力分布によらないことを見る.を固定すると,は確率で,確率での反対をとるので,その分布のエントロピーはによらずである.定義 1.2.2 はこれをで平均したものだから,どんな入力分布に対してもである.

次にを入力分布の式で書く.とおくと

𝑞(1)=𝑝(0)𝑊(1∣0)+𝑝(1)𝑊(1∣1)=(1−𝜋)𝜌+𝜋(1−𝜌)

であり,である.は二つの値しかとらないから例 1.1.2 よりであり,命題 6.1.3 よりとなって,主張の一般の式を得る.

上界はここから出る.に定理 1.1.5 を当てるとだからである.最後にこれが達成されることを見る.とするとでであり,出力は一様だから例 1.1.2 よりで,に等しい.上界が達成されたので,定義 6.1.4 の上限はこの値である.◼

形式化上の注記

例 6.1.8 に対応する宣言は形式化されていない.例 6.1.8 に付した証明が,この主張の保証のすべてである.

例 6.1.9(二元消失通信路). を固定する.X ={0,1},とし,

𝑊(𝑥∣𝑥)=1−𝜂,𝑊(∗∣𝑥)=𝜂,𝑊(1−𝑥∣𝑥)=0(𝑥∈{0,1})

で定まる通信路を,消失確率 の 二元消失通信路 という.送ったビットは確率でそのまま届き,確率で「読めなかった」ことを示す記号に化ける.反転して届くことはない(図 6.1.2).その容量は

𝐶(𝑊)=1−𝜂

であり(二元対称通信路の容量との比較は 図 6.1.3),一様入力で達成される.さらに,どの入力分布についても,とおくとである(は例 1.1.2 の二値エントロピー).

図 6.1.2. 消失確率の二元消失通信路.出力のは読めなかったことを示す記号で,とが入れ替わる遷移はない.

証明. 入力分布を任意にとり,とおく.二元対称通信路のときと同じく,を固定したときのの分布は「確率で,確率で」で,そのエントロピーはによらずである.よって.

出力分布は,𝑞(1) =𝜋(1 −𝜂),である.定義 1.1.1 に代入し,の積を和に分けると

𝐻(𝑌)=−(1−𝜂)(𝜋log⁡𝜋+(1−𝜋)log⁡(1−𝜋))−(1−𝜂)log⁡(1−𝜂)−𝜂log⁡𝜂=(1−𝜂)𝐻𝑏(𝜋)+𝐻𝑏(𝜂)

となる(確率 0 の項はの約束により両辺で消えるので,分け方に注意は要らない).命題 6.1.3 より

𝐼(𝑝;𝑊)=𝐻(𝑌)−𝐻(𝑌∣𝑋)=(1−𝜂)𝐻𝑏(𝜋)

である.例 1.1.2 よりだからであり,一様入力ではこの値が達成される.上界が達成されたので,定義 6.1.4 の上限はこの値である.◼

図 6.1.3. 二元対称通信路の容量と,二元消失通信路の容量.横軸は反転確率または消失確率である.

形式化上の注記

例 6.1.9 に対応する宣言は形式化されていない.例 6.1.9 に付した証明が,この主張の保証のすべてである.

二つを数で比べる. 反転確率の二元対称通信路の容量はビットである.いっぽう消失確率の二元消失通信路はビットで,こちらのほうがずっと大きい.どちらも10 回に 1 回は届いた記号が信用できないのに,差がこれだけ開くのは,どこが壊れたかを受け取る側が知っているかどうかの違いによる.消失通信路では壊れた位置がとして見えるので,残りの割はそのまま信用してよい.対称通信路では壊れた位置が見えないので,届いたビットのどれもが疑わしくなる.反転確率をまで上げると容量は約ビットになり,雑音がないときの 1 ビットのちょうど半分に落ちる.

両端で確かめる. 二元対称通信路でならで容量は 1 ビットである.雑音がないので 1 回に 1 ビット運べる.でもなので容量はやはり 1 ビットで,これは「必ず反転する」通信路が受け取る側で読み替えれば無雑音と同じだからである.いちばん悪いのはで,より容量は 0 になる.このとき出力は入力に依らず一様で,何も運べない.二元消失通信路のほうはで1 ビット,で 0 ビットと,消失確率にそのまま比例して減る.

容量を達成する入力分布は定理 6.1.5 で存在だけが分かっている.それが満たさなければならない条件は 6.6 節で与える(逆向き,つまりその条件を満たせば最大化子であるという主張は扱わない).二つの例で一様入力が最大化したのは通信路の対称性のおかげで,一般にはそうならない.

情報理論入門(Lean による形式証明つき)・GitHub