8.1 ガウス通信路

第6章の通信路は,入力も出力も有限個の記号だった.遷移確率を並べた表がその全体で,容量は有限個の数の上での最大化問題だった.しかし実際に信号を運ぶ物理的な経路は,電圧や電磁波の振幅といった実数を運ぶ.そこに乗る雑音も実数値で,しかも多くの小さな擾乱が重なった結果として生じる.本章が扱うのは,出力が入力にそういう雑音を足したものになる通信路であり,雑音の分布はガウス分布にとる.理由は二つある.一つはいま述べた成り立ちで,多くの小さな擾乱の和がガウス分布に近づくことは中心極限定理として知られている.もう一つは,散らばりの大きさを固定すると微分エントロピーが上から抑えられるという第7章 定理 7.2.4 の性質が,8.2 節で容量の上界を与えるところまでそのまま効くことである.

第7章で作った量は,そのまま当てられる.出力の散らばりの広さを微分エントロピーで測り,入力を知ってもなお残る散らばりを差し引けば,通信路が運んだ情報になる.第6章 命題 6.1.3という読み方が,に取り替えてそのまま生きる.ただし差し引くべき量,すなわち条件付き微分エントロピーは,第7章が用意しなかった.本節はまずそれを定義し,ガウス通信路ではその値が入力によらない定数になることを見る.

本章ではの底をにとる.8.3 節以降での形の量を扱うので,第6章と同じく指数と対数の底をそろえておくと式が読みやすい.第7章の等式・不等式は底をより大きくとるかぎりそのまま成り立つ(7.1 節)ので,底を決めても引き継ぎに支障はない.

記号を二つ断っておく.は本章を通じて信号の電力を表す実数である.第3章は同じ字で遷移確率を,第10章は分布を表しており,どちらも裸のを書くので,引数の形では見分けられない.電力をと書くのはこの分野の標準なので改名しない.本章に現れる裸のはつねに電力である.添字の付いた系(定義 8.3.1)は誤り確率を,系(8.6 節)は各通信路に配る電力を表す.は雑音の分散を表す実数で,こちらは第2章の個数第4章の個数と違って引数をとらない.

通信路の定義

定義 8.1.1(ガウス通信路). とする.ガウス通信路 とは,入力に対して出力の分布を定義 7.2.1)と定める対応である.をこの通信路の 雑音の分散 と呼ぶ.長さの入力に対する出力の結合密度は

𝑛1𝑖=0𝑔𝑥𝑖,𝑁(𝑦𝑖)

とする(定義 7.2.1 の記号).入力に密度をもつ実数値確率変数を与えたとき,入力と出力の対結合密度

𝑓(𝑥)𝑔𝑥,𝑁(𝑦)((𝑥,𝑦)2)

と定め,この出力 と呼ぶ.出力から入力を引いた雑音 と呼び,入力と出力の関係を

𝑌=𝑋+𝑍

と書く.

定義 7.2.1 の密度を書き下すとである.すなわち,入力を送ったときの出力の分布は,平均分散のガウス分布をだけ平行移動したものにほかならない.ここでのは,定義 8.1.1 が通信路への入力として固定した一つの実数であって,確率変数を観測して得た値ではない.したがってこの読み下しは,確率の事象で条件付けることを含まない(その条件付けを本書がどう避けるかは,次に条件付き微分エントロピーを定義するところで述べる).雑音の側で読むと,実数を入れたときのの密度はで,を含まない.という式そのものはを移項しただけのもので,それだけでは何も言っていない.内容があるのは,どの入力を選んでも雑音の分布がのままで,出力はそこへ入力を足しただけのものになる,といういま見たことのほうである.8.3 節8.4 節で符号を扱うときも,長さの符号語に対する出力はその符号語に雑音の列を足したものになる.長さの入力に対する結合密度を積の形にとったのは,各時刻の雑音がほかの時刻の入出力に依存しないということで,定義 6.1.1 が積の形で書いた記憶のなさと同じ条件である.入力と出力の対の結合密度をと定めたのも,定義 6.1.1 が結合分布をと定めたのと同じ形である.

形式化: awgnChannel (ソース)

形式化上の注記. 形式化の awgnChannel は,入力に平均分散のガウス分布を返す核であり,定義 8.1.1回の使用の部分にあたる.長さへの延長は,一般の核を成分ごとの積測度に延ばす Channel.toBlock (InformationTheory/Shannon/BlockwiseChannel/Definition.lean) と,成分ごとに雑音の分散を変えられる parallelGaussianChannel (InformationTheory/Shannon/ParallelGaussian/Basic.lean) の二つが与える.ただしガウス通信路の符号を扱う段では,そのどちらも経由せず,誤り確率を定める errorProbAt (InformationTheory/Shannon/ChannelCoding/Basic.lean) の中で積測度が直に組まれる.入力と出力の対の結合密度にあたる部分も,形式化では密度ではなく,入力の分布と核から合成される結合分布として扱われる.

入力に何の制限も置かなければ何が起きるかを,見当として見ておく.出力の散らばりを抑えているのは雑音だけなので,入力の振れ幅を大きくとれば,離れた二つの入力は出力を見るだけで区別がつくようになる.二つに限らない.雑音の散らばりにくらべて十分に離れた点を個並べれば,どの点を送ったかがほぼ確実に読み取れてビットが運べ,はいくらでも大きくとれる.ここまでは見当であって,メッセージがその個数だけ区別できることを示したわけではない.入力に制限を置く理由としてはこれで足りるので,そこで送る信号の電力を抑える.を電力と呼ぶのは,が電圧や電流の振幅を表すとき,単位時間あたりに消費されるエネルギーがに比例するからである.回の使用については入力の乗の平均以下にすること,長さのブロックについては各符号語の乗和を以下にすることが,本章で置く制約である(前者を 8.2 節,後者を 8.3 節で使う).どちらも「回の使用あたりに使える電力がまで」という同じ制限を,一つの確率変数と一本の符号語のそれぞれについて書いたものである.

条件付き微分エントロピー

第7章 7.4 節は,条件付き微分エントロピーを用意しなかった.その定義にを知ったときのの密度が要り,条件付ける相手が確率の事象になるからである.差の形で書けば,この困難を通らずに定義できる.離散のときの定理 1.2.3 のチェイン則)を,そのまま定義に据えればよい.

定義 8.1.2(条件付き微分エントロピー). を実数値確率変数とし,が結合密度をもつとする.がどちらも定まるとき,を与えたときの条件付き微分エントロピー

(𝑌𝑋):=(𝑋,𝑌)(𝑋)

で定め,このとき定まる という.

定義 8.1.2 が測っているのは,二つをまとめて見たときの散らばりから,だけを見たときの散らばりを差し引いた残りである.を知ったあとにに残る散らばり,と読める.が負になりうる(例 7.1.3)ので,も負になりうる.定義 8.1.2 が要求しているのは対が結合密度をもつことである.8.3 節で符号を扱うときはこの前提が満たされないことを見るので,そこでは別の形の道具を使う.

形式化上の注記. 形式化には condDifferentialEntropy (InformationTheory/Shannon/EPI/G2/ConvEntropyMonotone.lean) があるが,定義 8.1.2 とは形が違う.あちらは条件付き分布を測度論の道具(condDistrib)で取り出し,その微分エントロピーを条件の分布で平均したもので,本文が避けた「を知ったときのの分布」を正面から使っている.しかも平均をとる積分は,被積分関数が可積分でないときにを返す約束のものである.二つが等しいことを述べる宣言は在庫にないから,同じ値になると無条件には言えない.

命題 8.1.3. を実数値確率変数とし,が結合密度をもつとする.(𝑋,𝑌)(𝑋)がいずれも定まるならば,は定まり

𝐼(𝑋;𝑌)=(𝑌)(𝑌𝑋)0

である.

証明. 仮定より定義 8.1.2 の条件が満たされ,が定まる.同じ仮定は定理 7.5.6 の仮定でもあるから,は定まってである.右辺のにほかならない.

形式化上の注記. 命題 8.1.3 にあたる単独の宣言はない.differentialEntropy_sub_condDifferentialEntropy_eq_toReal_klDiv (InformationTheory/Shannon/EPI/G2/ConvEntropyMonotone.lean) が同じ形の等式を与えるが,そこで差し引かれているのは condDifferentialEntropy,すなわち条件付き分布を測度論の道具で取り出して平均した量であり,本文の定義 8.1.2 とは対象が違う.

命題 8.1.3 は,第6章 命題 6.1.3を微分エントロピーで書き直したものである.読み方も同じで,運べた情報とは出力の散らばりから雑音の分を差し引いた残りである.離散のときと違うのは,もそれ単独では散らばりの大きさとして読めない(7.1 節)ことで,読めるのは差のほうだけである.

雑音の分は入力によらない

ガウス通信路では,差し引くべきが入力によらない一つの定数になる.雑音の分布が入力によって平行移動するだけで,形が変わらないからである.

証明のために道具を二つ借りる.一つめは 累次積分 で,としての座標を二つの組に分けるとき,上の実数値関数が非負であるかであるならば,全体での積分は,片方の組について先に積分してからもう片方の組について積分したものに等しく,が非負のときは両辺が同時にになることを許す,という積分の事実である.で書けば

2𝑢(𝑥,𝑦)𝑑𝑥𝑑𝑦=(𝑢(𝑥,𝑦)𝑑𝑦)𝑑𝑥

である.当てる対象は,定義 8.1.1 の結合密度に対数や次以下の多項式を掛けた関数と,その絶対値と,8.3 節に現れるガウス密度の積と,8.4 節に現れる立方体の上の定数関数である.依存するのは命題 8.1.4 と,8.2 節定理 8.2.2 と,8.3 節補題 8.3.3 と,8.4 節命題 8.4.1 の証明である.二つめは学部の確率で扱う ガウス分布の畳み込み で,と実数について,どの実数に対しても

𝑔𝜇,𝜎21(𝑥)𝑔𝑥,𝜎22(𝑦)𝑑𝑥=𝑔𝜇,𝜎21+𝜎22(𝑦)

が成り立つ,という事実である(独立なガウス分布に従う二つの確率変数の和がふたたびガウス分布に従い,分散が足し算になることの,密度による言い換えである).当てる対象は定義 7.2.1 のガウス密度だけで,依存するのは例 8.1.6 の証明と,同じ事実を 9.4 節で借り直している第9章命題 9.4.4命題 9.4.6 の証明と,9.4 節の宣言を確率変数の形のまま引き直している第16章例 16.7.9命題 16.8.1 の証明である.ただしあちらが借りるのは確率変数の形で,和だけでなく差についても同じことをあわせて借りるので,借りる形は 9.4 節の宣言のほうが広い.第16章が引くのは和についてだけなので,そのぶんは 9.4 節より狭い.

命題 8.1.4. とし,密度をもつ実数値確率変数をガウス通信路(定義 8.1.1)に与えて出力を得るとする.が定まるならばは定まり

(𝑋,𝑌)=(𝑋)+12log(2𝜋𝑒𝑁),(𝑌𝑋)=12log(2𝜋𝑒𝑁)

である.

証明. 筋は行で書ける.結合密度の対数がの和に分かれ,第項の積分がによらない,というだけである.以下は,その計算が許されることを確かめながら進める.

定義 8.1.1 よりの結合密度はである.これが上の密度であることを見る.非負であり,借りた累次積分と,7.2 節で借りたガウス分布の基本性質から得られるにより

2𝐹(𝑥,𝑦)𝑑𝑥𝑑𝑦=𝑓(𝑥)(𝑔𝑥,𝑁(𝑦)𝑑𝑦)𝑑𝑥=𝑓(𝑥)𝑑𝑥=1

だからである.

対数を二つに分ける.の点ではであり,は実軸のすべての点で正だから,そこでである.の点では定義 7.1.2 の約束によりである.

の項の積分が入力によらないことを見る.定理 7.2.2に当てると,は有限であり

𝑔𝑥,𝑁(𝑦)log𝑔𝑥,𝑁(𝑦)𝑑𝑦=12log(2𝜋𝑒𝑁)

である.右辺はを含まない.また 7.3 節で借りた置換積分を𝑐 =1𝑏 =𝑥)という置き換えで使うと,だから

𝑔𝑥,𝑁(𝑦)log𝑔𝑥,𝑁(𝑦)𝑑𝑦=𝑔0,𝑁(𝑢)log𝑔0,𝑁(𝑢)𝑑𝑢=:𝐴

であり,このによらない有限の値である.

が定まることを見る.だから,借りた累次積分を非負関数に当てて

2𝐹log𝐹𝑓(𝑥)log𝑓(𝑥)(𝑔𝑥,𝑁(𝑦)𝑑𝑦)𝑑𝑥+𝑓(𝑥)(𝑔𝑥,𝑁(𝑦)log𝑔𝑥,𝑁(𝑦)𝑑𝑦)𝑑𝑥

を得る.右辺の第項はで,が定まることから有限,第項はで有限である.よっては定まる.

値を求める.いま見たとおり積分は絶対収束するから,借りた累次積分をに当てて二つの項に分けてよく

(𝑋,𝑌)=𝑓(𝑥)log𝑓(𝑥)(𝑔𝑥,𝑁(𝑦)𝑑𝑦)𝑑𝑥+𝑓(𝑥)(𝑔𝑥,𝑁(𝑦)log𝑔𝑥,𝑁(𝑦)𝑑𝑦)𝑑𝑥

である.第項はに等しく,第項は上で見た値によりに等しい.符号を戻すとである.最後に定義 8.1.2 を当てるとを得る.

形式化上の注記. 命題 8.1.4 の第の等式にあたる内容は,形式化では入力を一つ固定した形で述べられている.differentialEntropy_awgnChannel_apply_eq_noise (InformationTheory/Shannon/AWGN/MIBridge.lean) が,入力に対する出力の分布の微分エントロピーが雑音の分布のそれに等しいことを,によらず与える.それを入力の分布で平均した形が awgn_cond_entropy_eq_noise_entropy_of_const (InformationTheory/Shannon/AWGN/MIBridge.lean) で,こちらは入力をガウス分布にとった場合の宣言である.入力の分布に条件を置かない形の宣言も在庫にあり,condDifferentialEntropy_indep_add_eq (InformationTheory/Shannon/EPI/G2/ConvEntropyMonotone.lean) が,独立な二つの実数値確率変数の和について,一方で条件付けた微分エントロピーが他方の微分エントロピーに等しいことを与える.雑音を和の第項,入力を第項,係数をにとれば,命題 8.1.4 の第の等式にあたる内容になる.いずれも条件付き微分エントロピーを入力ごとの微分エントロピーの平均として書いており,結合と周辺の差で定めた定義 8.1.2 とは対象が違う.第の等式にあたる宣言は在庫にない.

系 8.1.5. とし,密度をもつ実数値確率変数をガウス通信路(定義 8.1.1)に与えて出力を得るとする.がどちらも定まるならば,は定まり

𝐼(𝑋;𝑌)=(𝑌)12log(2𝜋𝑒𝑁)

である.

証明. 命題 8.1.4 よりが定まるから,命題 8.1.3 の仮定がすべて満たされ,である.命題 8.1.4 の第の等式でその第項を置き換えればよい.

系 8.1.5 が本章のこれ以降を動かす.運べる情報を大きくしたければ,を大きくするほかに手はない.差し引かれる量は通信路の雑音だけで決まっていて,入力の選び方では動かないからである.そしてをどこまで大きくできるかは,第7章がすでに答えている.出力の分散が抑えられているかぎり微分エントロピーには上限があり(定理 7.2.4),その上限をガウス分布が達成する(系 7.2.5).電力制約が出力の分散を抑え,最大エントロピー性が上限を与える,という二段構えが 8.2 節の内容である.

例 8.1.6(ガウス入力). 𝑁 >0とし,に従うとしてガウス通信路(定義 8.1.1)に与え,出力を得るとする.このときに従い,は定まって

𝐼(𝑋;𝑌)=12log(1+𝑃𝑁)

である.

証明. の密度はだから,定義 8.1.1 よりの結合密度はである.7.4 節で借りた周辺密度の存在により,その第周辺密度(定義 7.4.1)はの密度であり,借りたガウス分布の畳み込みによりそれは

𝑔0,𝑃(𝑥)𝑔𝑥,𝑁(𝑦)𝑑𝑥=𝑔0,𝑃+𝑁(𝑦)

である.よってに従う.

定理 7.2.2 よりは定まってである.系 8.1.5 より

𝐼(𝑋;𝑌)=12log(2𝜋𝑒(𝑃+𝑁))12log(2𝜋𝑒𝑁)=12log𝑃+𝑁𝑁

であり,右辺はである.

形式化上の注記. 例 8.1.6 にあたる内容は形式化されているが,相互情報量が本文と別の量である.mutualInfoOfChannel_gaussianInput_closed_form' (InformationTheory/Shannon/AWGN/MIClosedForm.lean) が,ガウス入力に対する相互情報量がに等しいことを,のもとで与える(通信路の核が可測であるという正則性の前提も付くが,これは別の宣言で解消される).出力がふたたびガウス分布になることは,核と入力の合成が畳み込みに等しいという isAwgnBindEqConv (InformationTheory/Shannon/AWGN/BindConvolution.lean) と,そこから出力の分布を読む awgn_output_gaussian_of_bind_eq_conv (InformationTheory/Shannon/AWGN/MutualInfoBridge.lean) の合成で得られる.ただしそこでの相互情報量は,通信路の核と入力の分布から測度論の言葉で組み立てたもので,定義 7.5.4 の密度で書いた相互情報量とは別の量である.

幅で読む. 例 8.1.6 の値は,をその分布の実効的な幅と読む物差し(7.1 節)を当てると読み下せる.底がなので,出力の実効的な幅はであり,入力を知ったあとに残る散らばりの幅はである.系 8.1.5 の右辺はこの二つの幅の比の対数にほかならず,比はである.平方根がを,二つの分散の比が「1 +」を出している.測っているのは雑音の幅の中に信号の幅がいくつ入るかではなく,出力全体の幅が雑音の幅の何倍かなので,電力をに近づければ比はに近づき,運べる量もに近づく.同じ比の読み方が 8.4 節でもういちど,こんどは体積の比として現れる.

比で決まる. 例 8.1.6 の値はの比だけで決まっていて,電力と雑音を同じ倍率で大きくしても変わらない.この比信号対雑音比 と呼ぶ.目盛を取り替えても相互情報量が動かないこと(系 7.5.7)が,ここに具体的な形で現れている.入力も雑音も同じ倍率で引き伸ばせば,だけ増える(系 7.3.4)が,雑音の分散がになるので差し引くも同じだけ増え,系 8.1.5 の差は動かない.

規模感. 底をにとってならビット,ならビット,なら約ビット,なら約ビットである.信号対雑音比を倍にしても運べる量が倍になるわけではなく,増えるのは高々ビットである(からビット,からビットの増加であり,比が大きいほどこの値に近づく).比の対数で効くので,効くのは掛け算ではなく足し算のほうである.第6章の二元対称通信路が反転確率で約ビットだった(例 6.1.8)のと並べると,信号対雑音比のガウス通信路がおよそそれと同じくらいの通信路にあたる.

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