8.3 逆定理

8.2 節は,相互情報量の最大化問題の値でしかない.第6章 6.1 節と同じで,この段階では「回の使用あたりビット運べる」と言う資格がない.本節はまず,ガウス通信路を使って送るとはどういう操作かを符号として定義し,そのうえで,どんな符号の族を設計しても達成レートがを超えられないことを示す.本の符号については,誤り確率を大きく許せばレートはいくらでも上げられるので,超えられないのは誤り確率をに近づける族が残せるレートのほうである.

論証の形は第6章 6.4 節と同じである.一様なメッセージのもつ情報量を,ファノの不等式で「通信路が運んだぶん」と「誤りで説明されるぶん」に分ける.運んだぶんは,系 8.1.5 と同じ形で出力の微分エントロピーから雑音の分を引いたものになり,を電力制約で抑えればが出る.抑えるのに使う道具は二つで,微分エントロピーの劣加法性(定理 7.4.3)と,分散を固定したときの最大エントロピー性(定理 7.2.4)である.道具の対応も見ておく.第6章では,記憶のなさがの右辺を各時刻の相互情報量の和に落とした(定理 6.4.5).本節でその役をするのは,出力の側を成分ごとにほどく劣加法性と,雑音の側の成分が独立なのでが和になること(系 7.4.4)の二つである.記憶のなさは,長さの出力の結合密度を積の形にとった定義 8.1.1 に,すでに書き込まれている.

条件付けについて,先に断っておくことがある.8.1 節が条件付き微分エントロピーの定義で手間をかけたのは,条件付ける相手が実数値の入力で,という確率の事象で条件付けることになるからだった.本節が条件付ける相手はメッセージで,これは有限個の値を確率ずつでとるから,その困難には当たらない.避けたのは連続な側で条件付けることであって,離散な側で条件付けることではない.

そのうえで,一つだけ第6章と事情が違うところがある.符号語は有限個しかないので,入力がとる値はの有限個の点である.したがって対は,の中で「第座標の組がその有限個の点のどれかに等しい」という集合の上に確率を置く.この集合の次元の体積がであることを微積分の計算規則として既知とすると,7.1 節で借りた体積の集合の確率により,には定義 7.1.1 の意味の結合密度がない.したがって定義 8.1.2 の条件付き微分エントロピーは当てられず,ファノの不等式を通す段も第6章のままでは書けない.そこはひとまとめに借用する.

借りる道具

観測が実数値のときのファノの不等式を,単発逆定理(実数値の観測) という名前で証明せずに借りる.主張の形は次のとおりである.をとり,上の一様分布に従う確率変数,に値をとる確率変数とする.上の密度で,を与えたときのの条件付き分布がを密度にもつものとし,を密度にもつ確率変数の微分エントロピーをと書く.が密度をもってが定まり,各も定まるとする.からへの復号器が,どのについてもを集合の上で積分できるものであるならば,について

log𝑀(𝑉)1𝑀𝑀𝑚=1𝑚+𝐻𝑏(𝑃𝑒)+𝑃𝑒log(𝑀1)

が成り立つ(例 1.1.2 の二値エントロピー).当てる対象は,定義 8.3.1 のブロック通信路符号をガウス通信路に通して得られる組だけであり,依存するのは定理 8.3.5 の証明だけである.

これは定理 6.4.1 の観測が実数値の場合にあたる.右辺の読み方を書いておく.を与えるごとにの分布はに定まるから,その微分エントロピーで平均した第項は,メッセージを知ったあとにに残る散らばりを測っている.したがって右辺の第項と第項の差は,命題 8.1.3と同じ形をしている(メッセージは有限個の値をとるので対に結合密度はなく,定義 8.1.2 の意味の条件付き微分エントロピーそのものは書けない).借りるのは,この差が定理 6.4.1の役を果たすということである.

本書がこれを証明しないのは,観測が実数値のときにを定めるのに,第7章が避けた「確率の事象で条件付けた分布」の道具が要るからである.形式化は測度論の言葉で相互情報量を直接扱えるので第6章と同じ道を通れるが,本書はその道具を用意していないので,ここでは結論の形だけを借りる.借りたからといって,そこから出る不等式が保証の外に出るわけではない.定理 8.3.5 の不等式そのものは無条件に機械検証されている.ただし形式化はこの一歩を通らず,別のルートで同じ不等式に達している.

電力制約つきブロック通信路符号

定義 8.3.1(電力制約つきブロック通信路符号). 𝑁 >0𝑃 >0𝑀 1とする.雑音の分散のガウス通信路(定義 8.1.1)に対する長さ,電力電力制約つきブロック通信路符号 とは,符号化写像と復号器の組であって,各について 電力制約

𝑛1𝑖=0(𝑐(𝑚)𝑖)2𝑛𝑃

を満たし,どのとどのについても,密度を集合の上で積分できるものをいう.その レートで定める.メッセージを送ったときの 誤り確率,およびそれをについて最悪値と平均で束ねた 最大誤り確率平均誤り確率

𝑃𝑒(𝑚):={𝑦𝑛:𝑑(𝑦𝑛)𝑚}𝑛1𝑖=0𝑔𝑐(𝑚)𝑖,𝑁(𝑦𝑖)𝑑𝑦𝑛,𝑃𝑒,max:=max1𝑚𝑀𝑃𝑒(𝑚),¯𝑃𝑒:=1𝑀𝑀𝑚=1𝑃𝑒(𝑚)

で定める.

第6章定義 6.2.1 との違いは三つである.符号語がではなくに住むこと,電力制約が付いたこと,そして誤りの事象の上で密度を積分できることを求めていることである.三つめは,出力が実数値になったせいで誤り確率が和ではなく積分で定まることから来る条件で,第6章では誤りの事象が有限集合だったので要らなかった.7.1 節が積分の線形性や単調性を微積分の計算規則として既知としたのと同じ資格で,積分できる集合どうしの有限個の和と差もまた積分できる集合であることを,本章では既知とする(使うのは系 8.4.2 の証明だけである).電力制約をで書いたのは,回の使用あたりの電力をまでに抑える,という意味である.レートの読み方も第6章と同じで,通りのメッセージを区別するのに要るビットを回の使用で割った,回あたりの取り分である.

形式化: AwgnCode (ソース)

形式化上の注記. 形式化の AwgnCode は符号化写像・復号器・電力制約を一つにまとめた構造で,定義 8.3.1 の「密度を積分できる」にあたる条件は,復号器が可測であることとして書かれている.誤り確率は errorProbAt を核 awgnChannel に当てたもので,本文の積分にあたる.最大誤り確率を束ねた宣言はなく,第6章 6.2 節と同じく「すべてのメッセージについて誤り確率が小さい」という形で書かれる.

定義 8.3.2(達成レート). 𝑁 >0とする.長さ,電力のブロック通信路符号を第項とする族達成可能 であるとは,最大誤り確率がを満たし,かつレートの列が上に有界なことをいう.このときをその族の 達成レート と呼び,達成レート全体の集合をと書く.

定義 6.2.2 をそのままガウス通信路に移した定義である.電力制約は符号の側に入っているので,達成レートの定義そのものは第6章と変わらない.

形式化上の注記. 達成レートに対応する宣言はガウス通信路の側にはない.第6章 6.2 節で触れたとおり,通信路の側に達成レートの集合そのものを表す宣言がまだ用意されていないためである.

出力の微分エントロピーを抑える

符号語を有限個並べて一様に選ぶと,出力の分布はガウス分布の混合になる.その微分エントロピーが定まることと,各成分の乗の平均がどんな値になるかを,先に確かめておく.その値に電力制約を当てて上から抑えるのは,補題 8.3.4 のほうである.

補題 8.3.3(ガウス分布の混合). 𝑁 >0𝑀 1とし,とする.

¯𝑞(𝑦):=1𝑀𝑀𝑚=1𝑘1𝑗=0𝑔𝑎𝑚,𝑗,𝑁(𝑦𝑗)(𝑦𝑘)

上の密度である.を密度にもつ値の確率変数をとすると,は定まり,各については密度をもち,は定まって

𝔼[𝑉2𝑗]=𝑁+1𝑀𝑀𝑚=1𝑎2𝑚,𝑗

である.

証明. この補題があとで使う内容は一つで,混合の乗平均が,混ぜる前の各分布の乗平均の平均になる,ということである.残りはすべて,が密度であること,その周辺密度が何になるか,微分エントロピーが定まること,の確認にあたる.

が密度であることを見る.非負であり,7.4 節で借りた独立と積の密度により各上の密度だから,その平均であるの全積分もである.

周辺密度を求める.8.1 節で借りた累次積分により,を第変数以外について積分すると,各因子の全積分がだからが残る.よって定義 7.4.1 の第周辺密度はであり,7.4 節で借りた周辺密度の存在によりこれはの密度である.

乗の平均を求める.定理 8.2.2 の証明と同じ計算で,7.2 節で借りたガウス分布の基本性質からである.よって 7.1 節で借りた期待値の積分表示より

𝔼[𝑉2𝑗]=1𝑀𝑀𝑚=1𝑔𝑎𝑚,𝑗,𝑁(𝑡)𝑡2𝑑𝑡=1𝑀𝑀𝑚=1(𝑁+𝑎2𝑚,𝑗)

であり,これが主張の値である.同じ計算をのかわりにについて行うと,もやはり有限である.

が定まることを見る.のすべての点で正である.上からは,各因子がを満たすからであり,したがってである.下からは,和の一項だけを残してであり,定義 7.2.1 の密度の対数をとると

log¯𝑞(𝑦)log𝑀+𝑘2log(2𝜋𝑁)+log𝑒2𝑁𝑘1𝑗=0(𝑦𝑗𝑎1,𝑗)2

である.二つを合わせると,右辺の第項を絶対値に替えた次以下の多項式

log𝑀+𝑘2log(2𝜋𝑁)+log𝑒2𝑁𝑘1𝑗=0(𝑦𝑗𝑎1,𝑗)2

を上から抑える.いま見たとおり各は有限だから,7.1 節で借りた期待値の積分表示より,これにを掛けた積分は有限である.よっても有限で,は定まる.の場合を各の密度に当てれば,が定まることも同じ議論から従う.

補題 8.3.4(出力の微分エントロピーの上界). 𝑁 >0𝑃 >0𝑀 1とし,長さ,電力のブロック通信路符号(定義 8.3.1)をとる.上の一様分布に従うメッセージの符号語をガウス通信路に通して得られる出力をとすると,は定まり

(𝑌𝑛)𝑛2log(2𝜋𝑒(𝑃+𝑁))

である.

証明. 出力の密度を求める.を与えたときのの分布は,定義 8.1.1 より密度をもつ.個の値を等確率でとるから,どの直方体についても

Pr[𝑌𝑛𝐽]=1𝑀𝑀𝑚=1𝐽𝑛1𝑖=0𝑔𝑐(𝑚)𝑖,𝑁(𝑦𝑖)𝑑𝑦𝑛

であり,右辺は補題 8.3.3ととったもののの上での積分である.よって定義 7.1.1 の意味での密度はこのであり,補題 8.3.3 よりと各は定まって

𝔼[𝑌2𝑖]=𝑁+1𝑀𝑀𝑚=1(𝑐(𝑚)𝑖)2

である.以下この値をと書く.よりである.

成分ごとに抑える.補題 8.3.3 よりは密度をもち,は定まる.が有限だからよりは平均をもち,その分散はである.よって定理 7.2.4 よりである.

和を抑える.定理 7.4.3 の劣加法性よりである.補題 1.1.9 を,凹関数補題 1.1.6 による),重み,点ととって当てると

1𝑛𝑛1𝑖=0log(2𝜋𝑒𝜎2𝑖)log(1𝑛𝑛1𝑖=02𝜋𝑒𝜎2𝑖)

である.電力制約より

1𝑛𝑛1𝑖=0𝜎2𝑖=𝑁+1𝑛𝑀𝑀𝑚=1𝑛1𝑖=0(𝑐(𝑚)𝑖)2𝑁+1𝑛𝑀𝑀𝑛𝑃=𝑁+𝑃

である.補題 8.2.5 よりは単調だから,上の二つを合わせて

(𝑌𝑛)𝑛21𝑛𝑛1𝑖=0log(2𝜋𝑒𝜎2𝑖)𝑛2log(2𝜋𝑒(𝑃+𝑁))

を得る.

補題 8.3.4 の二つの不等号は,別のことを言っている.劣加法性のほうは,回の使用をまとめて見たときの散らばりが,回ずつ見たときの和を超えないことである.捨てているのは,出力の成分のあいだの依存で,符号語が時刻をまたいで作る相関はそこに現れる.補題 1.1.9 のほうは,各時刻の分散を平均で置き換えて上から抑えている段で,成分ごとの偏りをここで捨てている.どちらもの向きが逆定理に必要な向きにそろっている.の凹性が上から抑える向きに効くのは,8.2 節例 8.2.6 で「薄く広げるほうがを大きくする」ことを見たのと同じ形である.

逆定理

定理 8.3.5(ガウス通信路の逆定理). 𝑁 >0𝑃 >0𝑀 2とし,長さ,電力のブロック通信路符号(定義 8.3.1)をとる.メッセージが上の一様分布に従うとしたときの平均誤り確率をとすると

log𝑀𝑛2log(1+𝑃𝑁)+𝐻𝑏(¯𝑃𝑒)+¯𝑃𝑒log(𝑀1)

である(例 1.1.2 の二値エントロピー).

証明. 借りた単発逆定理を,𝑘 :=𝑛𝑉 :=𝑌𝑛ととって当てる.仮定を確かめる.を与えたときのの分布がを密度にもつことは定義 8.1.1 そのものであり,が密度をもってが定まることは補題 8.3.4 で見た.復号器が各についてを集合の上で積分できることは,定義 8.3.1の場合として求めている条件である.

を求める.7.4 節で借りた独立と積の密度により,を密度にもつ確率変数の成分は独立で,第成分はを密度にもつ.定理 7.2.2 より各成分の微分エントロピーは定まってだから,系 7.4.4 より

𝑚=𝑛1𝑖=012log(2𝜋𝑒𝑁)=𝑛2log(2𝜋𝑒𝑁)

である.これはによらないので,も同じ値である.

誤り確率を読み替える.メッセージが一様分布に従うときは,メッセージごとの誤り確率を一様な重みで平均したものだから,定義 8.3.1 の平均誤り確率に一致する.

以上を借りた単発逆定理に入れると

log𝑀(𝑌𝑛)𝑛2log(2𝜋𝑒𝑁)+𝐻𝑏(¯𝑃𝑒)+¯𝑃𝑒log(𝑀1)

である.補題 8.3.4 よりだから,第項と第項の差は

𝑛2log(2𝜋𝑒(𝑃+𝑁))𝑛2log(2𝜋𝑒𝑁)=𝑛2log𝑃+𝑁𝑁

以下であり,右辺はである.

形式化: awgn_converse (ソース)

形式化上の注記. 形式化の awgn_converse は,本文と同じ不等式を無条件に与える.仮定は本文と一つずつ対応していて,符号が AwgnCode であることが定義 8.3.1 の電力制約と可測性に,メッセージ数が以上であることが本文のに,長さが正であることがに,雑音の分散がでないことが本文のに,平均誤り確率の等式が本文のにあたる.本文に対して真に増えている条件は,核が可測であるという正則性のもの一つだけで,これは同じ名前空間の isAwgnChannelMeasurable (InformationTheory/Shannon/AWGN/ChannelMeasurability.lean) が与えるので,実際には条件になっていない.形式化はを自然対数にとるので単位はナットであり,本文が底をにとったのと単位の名前だけが違う.

証明のルートは本文と食い違う.本文は,実数値の観測についてのファノの不等式を単発逆定理としてまとめて借り,そこから出るを微分エントロピーの劣加法性(定理 7.4.3)で成分ごとにほどいた.形式化は劣加法性をまったく通らない.離散のファノの不等式をメッセージと復号結果に当て,データ処理不等式でメッセージと出力の相互情報量に移し,相互情報量のチェイン則と記憶のなさで各時刻の相互情報量の和に落とし,各時刻の相互情報量を出力の微分エントロピーと雑音の微分エントロピーの差に書き換え(awgn_per_letter_mi_bridge_genuine (InformationTheory/Shannon/AWGN/Converse.lean)),最大エントロピー性と Jensen の不等式で抑える.保証されているのは定理 8.3.5 の不等式そのもので,そこに至る手順ではない.

左辺が誰のものかを見る. 第6章 定理 6.4.1 のときと同じで,左辺のは符号を設計する側が決める量である.右辺の第項が通信路と電力制約で決まる量,残る項が誤りで説明のつく量である.両辺をで割ってから誤り確率をに近づけると,レートが上に有界でありさえすれば後ろの項は消える(そこは系 8.3.6 で確かめる).するとをすべてが支えなければならず,レートがを超えられないという主張になる.

弱逆定理

系 8.3.6(弱逆定理). 𝑁 >0とする.長さ,電力のブロック通信路符号を第項とする族定義 8.3.2 の意味で達成可能なら,その達成レートは以下である.すなわちのどの元もを超えない.

証明. 族を一つとり,とおく.定義 8.3.2 よりであり,レートの列は上に有界で,よりである.

メッセージに一様分布を入れる. 定義 8.3.2 の達成可能性は最大誤り確率で書かれているので,メッセージにどんな分布を与えても平均誤り確率は以下である.そこで各について,メッセージが上の一様分布に従うとして議論する.これは符号に置く追加の仮定ではなく,評価のためにこちらが選ぶ分布である.このとき定義 8.3.1 の平均誤り確率以下だからである.

とおく.これは非負であり,の有界性から,補題 6.4.7 よりである.以下,どのでもであることを示せばよい.

の項を片付ける. そのようなではである.よりであり,補題 8.2.5 よりである.よってが成り立つ.残るではである.

定理 8.3.5 を当てる. について定理 8.3.5 より

log𝑀𝑛𝑛2log(1+𝑃𝑁)+𝐻𝑏(¯𝑃(𝑛)𝑒)+¯𝑃(𝑛)𝑒log(𝑀𝑛1)

である.よりだから,で割って

𝑅𝑛12log(1+𝑃𝑁)+𝐻𝑏(¯𝑃(𝑛)𝑒)𝑛+¯𝑃(𝑛)𝑒𝑅𝑛=12log(1+𝑃𝑁)+𝛿𝑛

を得る.

どのでもでありだから,両辺で下極限をとってを得る.

形式化上の注記. 系 8.3.6 に対応する単独の宣言はない.達成レートの集合そのものに対応する宣言が通信路の側に用意されていないので,この主張は形式化の外にある.機械検証が及んでいるのは,各についての不等式である定理 8.3.5 に紐付けた宣言までで,補正項が消える段(第6章 補題 6.4.7 にあたる部分)は形式化されていない.

系 8.3.6 が言えたので,あとは逆向き,すなわち未満のレートが実際に達成できることを示せば,8.2 節で最大化問題の値として定めたに操作的な意味が付く.それが 8.4 節の内容である.

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