6.1 通信路と通信路容量
第2章と第3章は,情報源のもつ冗長さを取り除く話だった.同じ内容をより短く書けるのはどこまでか.その限界がエントロピーだった.本章で扱うのは逆向きの操作である.通信路は送った記号を確率的に取り違えるので,受け取った側が誤りを直せるように,送る側があらかじめ冗長さを足す.足しすぎれば運べる情報が減り,足りなければ誤りが残る.この綱引きの決着がどこにつくかが本章の主題である.
素朴に考えると,雑音のある通信路では誤り確率を 0 に近づけようとするたびに伝送速度が0 に落ちそうに思える.同じ記号を何度も繰り返して多数決をとれば誤りは減らせるが,繰り返した回数だけ速度は落ちる.Shannon の発見は,そうはならないということである.通信路ごとにある数𝐶が定まっていて,𝐶より小さい速度である限り,誤り確率はいくらでも 0 に近づけられる.本節ではその𝐶(通信路容量)を定義し,それが上限ではなく最大値として達成されることを示し,二つの通信路で値を計算する.
本章ではlogの底を 2 にとる.第2章・第3章と同じく2𝑛𝑅の形の量を扱うので,指数と対数の底をそろえておくと式が読みやすい.
通信路と,その相互情報量
定義 6.1.1(離散無記憶通信路). 入力アルファベットXと出力アルファベットYはどちらも有限とし,Xは空でないとする.通信路 とは,各𝑥 ∈Xに対してY上の分布𝑊( ⋅ ∣𝑥)を与える対応である.すなわち𝑊(𝑦 ∣𝑥) ≥0であって,各𝑥について∑𝑦∈Y𝑊(𝑦 ∣𝑥) =1を満たす.長さ𝑛の入力𝑥𝑛 =(𝑥0,…,𝑥𝑛−1)に対する出力の分布が
𝑊𝑛(𝑦𝑛∣𝑥𝑛)=𝑛−1∏𝑖=0𝑊(𝑦𝑖∣𝑥𝑖)で与えられるとき,この通信路を 記憶のない(memoryless)通信路 という.アルファベットが有限なので,これを 離散無記憶通信路 とも呼ぶ.X上の分布𝑝を入力に与えたときの入力と出力の 結合分布 を𝑝(𝑥)𝑊(𝑦 ∣𝑥),その第 2 成分の周辺分布𝑞(𝑦) :=∑𝑥𝑝(𝑥)𝑊(𝑦 ∣𝑥)を出力分布 と呼ぶ.
𝑊(𝑦 ∣𝑥)は「𝑥を送ったときに𝑦が受け取られる確率」である.通信路が定めるのはこの条件付き分布だけで,何を送るかは通信路の側では決まっていない.入力分布𝑝を与えてはじめて入力と出力が対になり,結合分布と出力分布が定まる.積の形の条件が言っているのは,各時刻の雑音がそれ以前の入出力に一切依存しない,ということである.これが「記憶のない」の内容で,第3章の記憶のある情報源とちょうど逆向きの仮定になっている.あちらは情報源に記憶があってよかったが,こちらは通信路に記憶がない.
いちばん小さい通信路で確かめる. X =Y ={0,1}とし,送ったビットが確率0.1で裏返る通信路を考える.遷移確率は
𝑊(0∣0)=𝑊(1∣1)=0.9,𝑊(1∣0)=𝑊(0∣1)=0.1である.同じ入力に対する二つの値の和が 1 になっていることが定義 6.1.1 の条件で,同じ出力に対する和には何の条件もない.通信路が定めるのは入力ごとの分布だからである.ここに入力分布𝑝(0) =𝑝(1) =1/2を与えると,結合分布は(0,0)と(1,1)に0.45,(0,1)と(1,0)に0.05を置き,出力分布は𝑞(0) =𝑞(1) =1/2になる.この通信路は,裏返る確率を一般の値にした形で例 6.1.8 にふたたび現れる.
定義 6.1.2(通信路の相互情報量). X上の分布𝑝と通信路𝑊に対し,出力分布を𝑞(𝑦) :=∑𝑥𝑝(𝑥)𝑊(𝑦 ∣𝑥)とおく.入力分布𝑝のもとでの𝑊の相互情報量 を
𝐼(𝑝;𝑊):=∑𝑥,𝑦𝑝(𝑥)𝑊(𝑦∣𝑥)log𝑊(𝑦∣𝑥)𝑞(𝑦)で定める.和は𝑝(𝑥)𝑊(𝑦 ∣𝑥) >0の項についてとる.
対数の中身𝑊(𝑦 ∣𝑥)/𝑞(𝑦)は,「𝑥が送られたと知ることで,𝑦が受け取られるという見込みが何倍になったか」である.何を送っても出力の分布が変わらない通信路(雑音が入力を完全に潰してしまう通信路)では,この比が全点で 1 になって𝐼(𝑝;𝑊) =0になる.逆に𝑥ごとに出力の分布がはっきり違えば比は 1 から離れ,𝐼(𝑝;𝑊)は大きくなる.1.3 節の言葉でいえば,𝐼(𝑝;𝑊)は入力と出力が独立からどれだけ隔たっているかを一つの数にまとめた量である.
定義 6.1.2 は通信路の側から書いた式で,第1章の相互情報量とは見た目が違う.両者が同じものであることを先に確かめておく.
命題 6.1.3. X上の分布𝑝と通信路𝑊をとり,(𝑋,𝑌)を結合分布𝑝(𝑥)𝑊(𝑦 ∣𝑥)に従う対とする.このとき
𝐼(𝑝;𝑊)=𝐼(𝑋;𝑌)=𝐻(𝑋)+𝐻(𝑌)−𝐻(𝑋,𝑌)=𝐻(𝑌)−𝐻(𝑌∣𝑋)=𝐻(𝑋)−𝐻(𝑋∣𝑌).
証明. (𝑋,𝑌)の同時分布は𝑝(𝑥,𝑦) =𝑝(𝑥)𝑊(𝑦 ∣𝑥)であり,その周辺分布は∑𝑦𝑝(𝑥,𝑦) =𝑝(𝑥)と∑𝑥𝑝(𝑥,𝑦) =𝑞(𝑦)である.和をとる範囲𝑝(𝑥,𝑦) >0では𝑝(𝑥) >0なので約分してよく,定義 1.3.1 の対数の中身は
𝑝(𝑥,𝑦)𝑝(𝑥)𝑞(𝑦)=𝑝(𝑥)𝑊(𝑦∣𝑥)𝑝(𝑥)𝑞(𝑦)=𝑊(𝑦∣𝑥)𝑞(𝑦)になる.したがって定義 6.1.2 の右辺は定義 1.3.1 の右辺そのもので,𝐼(𝑝;𝑊) =𝐼(𝑋;𝑌)である.残る三つの等式は,この対(𝑋,𝑌)に定理 1.3.4 を当てたものである.◼
雑音の分だけ減る,と読む. 以下でいちばん使うのは𝐼(𝑝;𝑊) =𝐻(𝑌) −𝐻(𝑌 ∣𝑋)の形である.𝐻(𝑌)は受け取った側から見た出力の不確かさ,𝐻(𝑌 ∣𝑋)は何を送ったかを知ってもなお残る不確かさ,すなわち雑音そのものの不確かさである.運べた情報とは,出力の不確かさから雑音の分を差し引いた残りだ,と読める.送る側に選べるのは𝑝だけで,𝑝を動かせば𝐻(𝑌)も𝐻(𝑌 ∣𝑋)も動く.その中で差を最大にする𝑝を探す,というのが次の定義である.
通信路容量
定義 6.1.4(通信路容量). 通信路𝑊の 通信路容量 を,入力分布全体にわたる上限
𝐶(𝑊):=sup𝑝𝐼(𝑝;𝑊)で定める.𝑝はX上の分布全体,すなわち𝑝(𝑥) ≥0かつ∑𝑥𝑝(𝑥) =1を満たす関数の全体を動く.
𝐶(𝑊)は通信路だけで決まる数で,入力分布にも符号にも依存しない.定義の時点では「1 回の使用あたりどれだけ運べるか」という操作的な意味はまだ持っておらず,単に相互情報量の最大化問題の値でしかない.それが誤りなく送れるレートの限界と一致する,というのが 6.3 節以降の内容である.
上限と書いたが,実際には最大値である.それを示すには,𝐼( ⋅ ;𝑊)が入力分布のつくる集合の上で連続であることと,その集合が有界閉集合であることの二つが要る.前者は,一次結合・合成・有限和が連続性を保つという微積分の計算規則と,1.1 節で認めた𝜑(𝑡) = −𝑡log𝑡の[0,∞)上の連続性から下で組み立てる(計算規則のほうは既知とする).後者から最大値の存在を引き出すために,次の一つを証明せずに借りる.
Weierstrass の最大値定理を借りる. 借りるのは「有限次元の実ベクトル空間の空でない有界閉集合の上の実数値連続関数は最大値をとる」という形である.当てる相手は,実ベクトル空間ℝXの部分集合である確率単体
Δ:={𝑝: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章は実ベクトル空間ℝX׈Xの部分集合である制約集合に,第11章はℝの閉区間[0,1]と,確率単体の閉部分集合とに,第15章は銘柄の番号の上の確率単体,すなわちℝ𝑚の部分集合であるポートフォリオ全体に当てる.本書はこの最大値定理を証明しないが,形式化されていないわけではない.定理 6.1.5 の形式化は,Mathlib にある無条件の機械検証済みのこの定理をそのまま呼び出しているからである.「本書で証明しない」ことと「形式化されていない」ことは別である.
定理 6.1.5. 𝑊を通信路(定義 6.1.1)とし,X上の分布全体のなす集合をΔと書く.このとき𝑝 ↦𝐼(𝑝;𝑊)はΔの上で連続であり,𝐼( ⋅ ;𝑊)を最大にする入力分布𝑝∗ ∈Δが存在して𝐶(𝑊) =𝐼(𝑝∗;𝑊)である.
証明. 連続性から示す.(𝑋,𝑌)を結合分布𝑝(𝑥)𝑊(𝑦 ∣𝑥)に従う対とすると,命題 6.1.3 より𝐼(𝑝;𝑊) =𝐻(𝑋) +𝐻(𝑌) −𝐻(𝑋,𝑌)である.三つの項を𝑝の関数として書き下すと
𝐻(𝑋)=∑𝑥𝜑(𝑝(𝑥)),𝐻(𝑌)=∑𝑦𝜑(∑𝑥𝑝(𝑥)𝑊(𝑦∣𝑥)),𝐻(𝑋,𝑌)=∑𝑥,𝑦𝜑(𝑝(𝑥)𝑊(𝑦∣𝑥))となる.どの項も,𝑝の成分の一次結合に𝜑を合成したものの有限和である.一次結合は連続で,𝜑は[0,∞)上で連続だから,𝑝 ↦𝐼(𝑝;𝑊)はΔの上で連続である.
最大値の存在に移る.ΔはℝXの有界閉集合である(有界性は0 ≤𝑝(𝑥) ≤1から,閉であることは定義が等式と不等式で書けていることから従う).よって借用した Weierstrass の最大値定理により,𝐼( ⋅ ;𝑊)はΔの上で最大値をとる.それを与える点を𝑝∗とする.
最後に𝐶(𝑊) =𝐼(𝑝∗;𝑊)を見る.𝐼(𝑝∗;𝑊)は値の集合{𝐼(𝑝;𝑊) :𝑝 ∈Δ}の上界であり,かつその集合の元でもある.上限は上界のうち最小のものだから,この二つから𝐶(𝑊) =𝐼(𝑝∗;𝑊)である.◼
命題 6.1.6. 任意の通信路𝑊に対し𝐶(𝑊) ≥0.
命題 6.1.7. 任意の通信路𝑊に対し𝐶(𝑊) ≤logmin(|X|,|Y|).
証明. 任意の入力分布𝑝をとり,(𝑋,𝑌)を結合分布に従う対とする.命題 6.1.3 より𝐼(𝑝;𝑊) =𝐻(𝑌) −𝐻(𝑌 ∣𝑋)であり,条件付きエントロピーは非負量の平均なので非負(定義 1.2.2 と命題 1.1.4).よって𝐼(𝑝;𝑊) ≤𝐻(𝑌)で,定理 1.1.5 より𝐻(𝑌) ≤log|Y|である.同じことを命題 6.1.3 のもう一つの形𝐼(𝑝;𝑊) =𝐻(𝑋) −𝐻(𝑋 ∣𝑌)で行えば𝐼(𝑝;𝑊) ≤log|X|を得る.𝑝は任意だったので,logmin(|X|,|Y|)は値の集合の上界であり,上限はそれ以下である.◼
読み方は素直である.1 回の使用で運べる情報は,入れ物の大きさを超えられない.送れる記号が|X|通りしかなければ,雑音がまったくなくても 1 回あたりlog|X|より多くは運べないし,受け取る側の区別が|Y|通りしかないなら同じことが出口の側で起きる.命題 6.1.7 は,容量が雑音の少なさだけでなく,入口と出口の狭いほうにも押さえられていることを言っている.
二つの通信路で計算する
例 6.1.8(二元対称通信路). 𝜌 ∈[0,1]を固定する.X =Y ={0,1}とし,
𝑊(1∣0)=𝑊(0∣1)=𝜌,𝑊(0∣0)=𝑊(1∣1)=1−𝜌で定まる通信路を,反転確率 𝜌の 二元対称通信路 という.送ったビットが確率𝜌で裏返り,確率1 −𝜌でそのまま届く通信路である.その容量は
𝐶(𝑊)=1−𝐻𝑏(𝜌)であり,一様入力𝑝(0) =𝑝(1) =1/2で達成される(𝐻𝑏は例 1.1.2 の二値エントロピー).さらに,どの入力分布𝑝についても,𝜋 :=𝑝(1)とおくと
𝐼(𝑝;𝑊)=𝐻𝑏((1−𝜌)𝜋+𝜌(1−𝜋))−𝐻𝑏(𝜌)である.
証明. まず𝐻(𝑌 ∣𝑋)が入力分布によらないことを見る.𝑋 =𝑥を固定すると,𝑌は確率1 −𝜌で𝑥,確率𝜌で𝑥の反対をとるので,その分布のエントロピーは𝑥によらず𝐻𝑏(𝜌)である.定義 1.2.2 はこれを𝑝(𝑥)で平均したものだから,どんな入力分布𝑝に対しても𝐻(𝑌 ∣𝑋) =𝐻𝑏(𝜌)である.
次に𝐼(𝑝;𝑊)を入力分布の式で書く.𝜋 :=𝑝(1)とおくと
𝑞(1)=𝑝(0)𝑊(1∣0)+𝑝(1)𝑊(1∣1)=(1−𝜋)𝜌+𝜋(1−𝜌)であり,𝑞(0) =1 −𝑞(1)である.𝑌は二つの値しかとらないから例 1.1.2 より𝐻(𝑌) =𝐻𝑏(𝑞(1))であり,命題 6.1.3 より𝐼(𝑝;𝑊) =𝐻(𝑌) −𝐻(𝑌 ∣𝑋) =𝐻𝑏(𝑞(1)) −𝐻𝑏(𝜌)となって,主張の一般の式を得る.
上界はここから出る.|Y| =2に定理 1.1.5 を当てると𝐻(𝑌) ≤log2 =1だから𝐼(𝑝;𝑊) ≤1 −𝐻𝑏(𝜌)である.最後にこれが達成されることを見る.𝑝(0) =𝑝(1) =1/2とすると𝜋 =1/2で𝑞(1) =12𝜌 +12(1 −𝜌) =12であり,出力は一様だから例 1.1.2 より𝐻(𝑌) =𝐻𝑏(1/2) =log2 =1で,𝐼(𝑝;𝑊) =1 −𝐻𝑏(𝜌)に等しい.上界が達成されたので,定義 6.1.4 の上限はこの値である.◼
例 6.1.9(二元消失通信路). 𝜂 ∈[0,1]を固定する.X ={0,1},Y ={0,1, ∗}とし,
𝑊(𝑥∣𝑥)=1−𝜂,𝑊(∗∣𝑥)=𝜂,𝑊(1−𝑥∣𝑥)=0(𝑥∈{0,1})で定まる通信路を,消失確率 𝜂の 二元消失通信路 という.送ったビットは確率1 −𝜂でそのまま届き,確率𝜂で「読めなかった」ことを示す記号∗に化ける.裏返って届くことはない.その容量は
𝐶(𝑊)=1−𝜂であり,一様入力𝑝(0) =𝑝(1) =1/2で達成される.さらに,どの入力分布𝑝についても,𝜋 :=𝑝(1)とおくと𝐼(𝑝;𝑊) =(1 −𝜂) 𝐻𝑏(𝜋)である(𝐻𝑏は例 1.1.2 の二値エントロピー).
証明. 入力分布𝑝を任意にとり,𝜋 :=𝑝(1)とおく.二元対称通信路のときと同じく,𝑋 =𝑥を固定したときの𝑌の分布は「確率1 −𝜂で𝑥,確率𝜂で∗」で,そのエントロピーは𝑥によらず𝐻𝑏(𝜂)である.よって𝐻(𝑌 ∣𝑋) =𝐻𝑏(𝜂).
出力分布は𝑞(0) =(1 −𝜋)(1 −𝜂),𝑞(1) =𝜋(1 −𝜂),𝑞( ∗) =𝜂である.定義 1.1.1 に代入し,logの積を和に分けると
𝐻(𝑌)=−(1−𝜂)(𝜋log𝜋+(1−𝜋)log(1−𝜋))−(1−𝜂)log(1−𝜂)−𝜂log𝜂=(1−𝜂)𝐻𝑏(𝜋)+𝐻𝑏(𝜂)となる(確率 0 の項は0log0 =0の約束により両辺で消えるので,分け方に注意は要らない).命題 6.1.3 より
𝐼(𝑝;𝑊)=𝐻(𝑌)−𝐻(𝑌∣𝑋)=(1−𝜂)𝐻𝑏(𝜋)である.例 1.1.2 より𝐻𝑏(𝜋) ≤log2 =1だから𝐼(𝑝;𝑊) ≤1 −𝜂であり,一様入力𝜋 =1/2ではこの値が達成される.上界が達成されたので,定義 6.1.4 の上限はこの値である.◼
二つを数で比べる. 反転確率𝜌 =0.1の二元対称通信路の容量は1 −𝐻𝑏(0.1) ≈0.531ビットである.いっぽう消失確率𝜂 =0.1の二元消失通信路は1 −0.1 =0.9ビットで,こちらのほうがずっと大きい.どちらも10 回に 1 回は届いた記号が信用できないのに,差がこれだけ開くのは,どこが壊れたかを受け取る側が知っているかどうかの違いによる.消失通信路では壊れた位置が∗として見えるので,残りの9割はそのまま信用してよい.対称通信路では壊れた位置が見えないので,届いたビットのどれもが疑わしくなる.反転確率を0.11まで上げると容量は約0.500ビットになり,雑音がないときの 1 ビットのちょうど半分に落ちる.
両端で確かめる. 二元対称通信路で𝜌 =0なら𝐻𝑏(0) =0で容量は 1 ビットである.雑音がないので 1 回に 1 ビット運べる.𝜌 =1でも𝐻𝑏(1) =0なので容量はやはり 1 ビットで,これは「必ず裏返る」通信路が受け取る側で読み替えれば無雑音と同じだからである.いちばん悪いのは𝜌 =1/2で,𝐻𝑏(1/2) =1より容量は 0 になる.このとき出力は入力に依らず一様で,何も運べない.二元消失通信路のほうは𝜂 =0で1 ビット,𝜂 =1で 0 ビットと,消失確率にそのまま比例して減る.
容量を達成する入力分布𝑝∗は定理 6.1.5 で存在だけが分かっている.それが満たさなければならない条件は 6.6 節で与える(逆向き,つまりその条件を満たせば最大化子であるという主張は扱わない).二つの例で一様入力が最大化したのは通信路の対称性のおかげで,一般にはそうならない.
InformationTheory — 形式化検証つき情報理論教科書(レビュー版).数式は MathJax + AMS Euler で事前レンダリング.