3.5 Shannon–McMillan–Breiman 定理

第2章の漸近等分配性(定理 2.1.4)は,i.i.d. 情報源で 1 文字あたりの驚きがに近づく,という主張だった.本節はそれをエルゴード的な定常情報源に広げる.近づく先はではなくエントロピーレートで,収束は概収束(確率 1 の事象の上で各点の収束)である.この一般化を Shannon–McMillan–Breiman 定理という.

まず量の名前を用意する.

定義 3.5.1(経験エントロピー). 情報源に対し,ブロックの確率をと書く.に対し

ˆ𝐻𝑛:=1𝑛log𝑝(𝑋𝑛)

を長さ経験エントロピー という.

式は定義 2.1.1 とまったく同じで,i.i.d. の仮定を外して同じ量を使う,という宣言である.違いは使える道具のほうにある.あちらでは独立性からが周辺分布の積にほどけ(補題 2.1.3),大数の法則にかけられる形になった.記憶のある情報源ではこのほどきが使えない.ほどく代わりに何をするかが,本節の中身である.

確率のブロックが実現する事象は,各で高々有限個の確率の事象の和だから確率である.したがって実現したブロックの確率は確率 1 で正で,は確率 1 で定まる有限の値をとる.

形式化: 経験エントロピーにあたる定義 blockLogAvg (ソース)

定理 3.5.2(Shannon–McMillan–Breiman). をエルゴード的な定常情報源とする.確率 1 で

ˆ𝐻𝑛=1𝑛log𝑝(𝑋𝑛)𝐻(X)(𝑛).

証明は上下からの挟み撃ちで行う.上からは記憶を有限の深さで打ち切った近似を,下からは無限の過去まで使った近似を当てる.どちらの側でも働くのは「確率分布どうしの比の期待値は1 を超えない」という一手だけなので,それを先に補題として切り出しておく.証明本体は補題 3.5.7補題 3.5.8 で,二つがそろったところで定理 3.5.2 の証明に戻る.

借りる道具

3.4 節の Birkhoff の個別エルゴード定理に加えて,下界(補題 3.5.8)の証明でだけ次の二つを借りる.上界(補題 3.5.7)はどちらも使わない.

両側への拡張. 定常情報源に対し,負の時刻まで延ばした列で,定常であり,かつ非負の時刻の部分の同時分布が元の情報源と一致するものが存在する.元がエルゴード的なら,延ばした列もシフトについてエルゴード的である.これは第15章 15.7 節 が借り直す.そちらが当てるのもアルファベットが株価比の有限集合になるだけの定常列なので,射程は広がらない.

無限の過去による条件付け. 両側に延ばした列について,有限の過去で条件付けた確率で確率 1 で収束する.その極限をと書く.さらに期待値も収束し,

𝔼[log𝑝(𝑋0𝑋1,𝑋2,)]=lim𝑘𝐻(𝑋0𝑋1,,𝑋𝑘)

が成り立つ.これは,増えていく情報で条件付けた確率が,その全体で条件付けた確率に確率 1 で収束するという Lévy のマルチンゲール収束定理と,期待値の側でその収束を許す優収束定理から得られる.

この極限について,以下では有限個の条件付けと同じ二つの操作を認めて使う.第一に,各時刻の条件付き確率が同じように定まり,その積が「無限の過去を与えたときのの条件付き確率」に等しいこと.これが条件付き確率のチェイン則である.第二に,の実現ごとにの分布が定まり,期待値が「を固定してとる期待値」と「についてとる期待値」の二段に分かれること.有限個の条件付けなら第1章の道具から出るが,無限の過去についてはここで認める.これは第15章 15.7 節 も引く.そちらが節の終わりで借りる二つは,無限の過去のもとでのの条件付き分布をこの極限として読み,その筋書きは,各期の条件付き評価を期ぶんの積へ渡す段が上の二つの操作の外に出ることを,足りないものとして名指す.当てる対象は同じく株価比の有限集合に値をとる定常列なので,射程は広がらない.

以上はどれも本書では証明しないが,3.4 節の Birkhoff の定理と同じく,無条件の機械検証済みの定理として形式化されている.

形式化: 両側への拡張 μZ,非負の時刻での一致 μZ_nat_proj_eq,シフトの測度保存性measurePreserving_shiftZ,エルゴード性 ergodic_shiftZ (ソース),Lévy の収束condProbPast_tendsto_condProbInfty とその極限 condProbInfty (ソース),無限の過去で条件付けた負対数の期待値 integral_pmfLogCondInfty_eq_entropyRate (ソース)

形式化上の注記. 期待値の等式は,形式化では右辺を条件付きエントロピーの極限ではなくエントロピーレートそのものとして述べている.本文が補題 3.5.8 の冒頭で定常性と定理 3.2.6 を経由して書き換える 2 段が,あちらでは 1 本の宣言に畳み込まれている.

比は指数の尺度では見えない

補題 3.5.3. 非負の確率変数の列が任意のを満たすなら,確率 1 で

limsup𝑛1𝑛logΛ𝑛0

である(のときはと読む).

証明. とおく.だから Markov の不等式が使えてである.したがって任意のに対し

Pr[𝑛𝑁𝐴𝑛]𝑛𝑁1𝑛2

であり,が収束するので右辺はに向かう.左辺はについて非増加だから,「が無限回起こる」事象の確率はである.

すなわち確率 1 で,十分大きなすべてのについて,したがって

1𝑛logΛ𝑛<2log𝑛𝑛

である.右辺はに収束するので,上極限は以下である.

補題 3.5.3 が言っているのは,期待値が 1 で抑えられた比はより速くは大きくならず,という尺度ではそれがに潰れる,ということである.第2章の議論が指数の肩をの精度でしか見ていなかったのと同じ粗さで,この粗さのおかげで「多項式ぶんのずれ」は最初から見えない.

記憶を文字で打ち切る

記憶の深さを文字に制限した近似分布を作る.

定義 3.5.4(次打ち切り近似). 定常情報源𝑘 0に対し

𝑞𝑘(𝑥𝑛):=𝑝(𝑥𝑘)𝑛1𝑖=𝑘𝑝(𝑥𝑖𝑥𝑖𝑘,,𝑥𝑖1)

と定める.ここでは,定常性により時刻に依らない量を表す(𝑝(𝑥0) :=1.条件にあたる長さの並びの確率がのときはと約束する).

これは「直前の文字しか覚えていない情報源」が同じブロックに与える確率である.なら 1 文字ずつ独立に出す情報源,ならマルコフ情報源(3.3 節)にあたる.真の分布と違って総和はちょうど 1 にならないが,超えることはない.

補題 3.5.5. 定義 3.5.4は,任意のに対しを満たす.

証明. 後ろの文字から順に和をとる.最後の文字について和をとると,条件にあたるの確率が正なら条件付き確率の和はちょうど 1 で,確率がなら約束により和はである.いずれにせよ 1 以下だから,全体はまでの部分で抑えられる.同じことをまで繰り返すとが残る.

補題 3.5.6. をエルゴード的な定常情報源とする.各に対し,確率 1 で

1𝑛log𝑞𝑘(𝑋𝑛)𝐻(𝑋𝑘𝑋𝑘)(𝑛).

証明. 実現した並びの確率は確率 1 で正だから,であり,またの各について長さの並びの確率も正である.後者から各段の条件付き確率が正なので,確率 1 でである.対数をとって

1𝑛log𝑞𝑘(𝑋𝑛)=1𝑛log𝑝(𝑋𝑘)+1𝑛𝑛1𝑖=𝑘(log𝑝(𝑋𝑖𝑋𝑖𝑘,,𝑋𝑖1))

と分かれる.第 1 項はに依らない確率変数をで割ったものだからに収束する.

第 2 項に 3.4 節の Birkhoff の定理を当てる.先頭文字だけで決まる関数

𝑓𝑘(𝑥):=log𝑝(𝑥𝑘𝑥0,,𝑥𝑘1)

をとると,である(条件付き確率を時刻に依らない量にとったのは定義 3.5.4 の側なので,ここは定義の書き換えにすぎない).と置き換えると第 2 項は

1𝑛𝑛𝑘1𝑗=0𝑓𝑘(𝜎𝑗𝑋)=𝑛𝑘𝑛1𝑛𝑘𝑛𝑘1𝑗=0𝑓𝑘(𝜎𝑗𝑋)

になる.であり,その期待値は定義 1.2.2 よりで,これは定理 1.2.4定理 1.1.5 より以下である.期待値が有限だからは確率 1 で有限な値をとり(確率の並びの上でになりうるが,そこは期待値に効かない),3.4 節で借りた定理が要求する形を満たす.Birkhoff の定理より第 2 因子は確率 1 でに収束し,と合わせて主張を得る.

形式化: negLogQk_div_tendsto_condEntropyTail (ソース)

上界

補題 3.5.7. をエルゴード的な定常情報源とする.確率 1 で

証明. を固定する.確率 1 でだから

Λ𝑛:=𝑞𝑘(𝑋𝑛)𝑝(𝑋𝑛)(𝑛>𝑘)

は確率 1 で定まる非負の確率変数である.補題 3.5.5 より,その期待値は

𝔼[Λ𝑛]=𝑥𝑛:𝑝(𝑥𝑛)>0𝑝(𝑥𝑛)𝑞𝑘(𝑥𝑛)𝑝(𝑥𝑛)=𝑥𝑛:𝑝(𝑥𝑛)>0𝑞𝑘(𝑥𝑛)1

である.補題 3.5.3 より確率 1 で,すなわち

1𝑛logΛ𝑛=ˆ𝐻𝑛(1𝑛log𝑞𝑘(𝑋𝑛))

と書き直してである.補題 3.5.6 より引かれている側はに収束するから

limsup𝑛ˆ𝐻𝑛𝐻(𝑋𝑘𝑋𝑘).

これが各について確率 1 で成り立つ.可算個の確率 1 の事象の共通部分もまた確率 1 だから,確率 1 ですべてのについて同時に成り立つ.そこでとすれば,定理 3.2.6より右辺はに収束する.

形式化: algoet_cover_limsup_bound (ソース)

下界

補題 3.5.8. をエルゴード的な定常情報源とする.確率 1 で

証明. 借りた両側への拡張により,以下ではに延ばしたエルゴード的な定常列とみなす.

𝑔(𝑥):=log𝑝(𝑥0𝑥1,𝑥2,)

とおくと

𝔼[𝑔(𝑋)]=lim𝑘𝐻(𝑋0𝑋1,,𝑋𝑘)=lim𝑘𝐻(𝑋𝑘𝑋𝑘)=𝐻(X)

である.最初の等号は借りた期待値の収束,真ん中の等号は両側列の定常性で時計をずらしただけ,最後の等号は定理 3.2.6 である.の期待値が有限なのだから,とくには確率 1 で有限な値をとり,3.4 節で借りた定理を当てられる.

比の期待値を抑える. 無限の過去をと書き,

𝑞(𝑋𝑛):=𝑛1𝑖=0𝑝(𝑋𝑖𝑋𝑖1,𝑋𝑖2,)

とおく.借りた条件付き確率のチェイン則より,これはを与えたときのの条件付き確率にほかならない.とおくと,を固定した条件付き期待値は

𝔼[Λ𝑛P]=𝑥𝑛𝑝(𝑥𝑛P)𝑝(𝑥𝑛)𝑝(𝑥𝑛P)=𝑥𝑛:𝑝(𝑥𝑛P)>0𝑝(𝑥𝑛)1

である.両辺の期待値をとってを得る.

補題 3.5.3 を当てる. だから,確率 1 で

limsup𝑛((1𝑛log𝑞(𝑋𝑛))ˆ𝐻𝑛)0.

一方であり,両側に延ばした列もエルゴード的な定常列だから,Birkhoff の定理よりこれは確率 1 でに収束する.𝑎𝑛 :=ˆ𝐻𝑛と書くとだから

liminf𝑛𝑎𝑛lim𝑛𝑏𝑛limsup𝑛(𝑏𝑛𝑎𝑛)𝐻(X)

である.

最後に,主張はだけで決まる量についてのものであり,両側に延ばした列の非負時刻の部分の分布は元の情報源と一致するから,元の情報源についても成り立つ.

形式化: algoet_cover_liminf_bound (ソース)

証明(定理 3.5.2). 補題 3.5.7補題 3.5.8 より,確率 1 で

𝐻(X)liminf𝑛ˆ𝐻𝑛limsup𝑛ˆ𝐻𝑛𝐻(X)

である.したがって極限が存在し,その値はである.

形式化: shannon_mcmillan_breiman (ソース)

二つの近似で挟む. 上からは記憶を文字で打ち切ったを,下からは無限の過去まで使ったを当てた.は真の情報源より予測が下手で(記憶が浅い),は上手である(の外にある情報まで使う).上からの評価はで止まり,まで下りてくる.下からの評価のほうは,Lévy の収束のおかげで最初からちょうどに着く.

同じ一手が両側で効く. 上界ではを,下界では補題 3.5.3 に入れた.使ったのはどちらも「確率分布どうしの比の期待値は 1 を超えない」だけである.比が多項式ぶんしか増えないことが分かれば,の尺度ではそれが消える.証明の仕事は,比の期待値が 1 以下になるように分母と分子を選ぶところに集中している.

なぜ無限の過去に出ていくのか. だけの関数なので,分母に置いて下界にも使えそうに見える.しかしの期待値は 1 で抑えられない.が極端に小さい並びで比が爆発しうるからである.を与えたときのの条件付き確率そのものなので,を固定してについて和をとると,分母がそこで割り戻されて残るのは真の確率の総和になり,1 を超えない.効いているのは大小関係ではなく,分母がそれ自身の確率で重みづけされているという正規化である.負の時刻に出ていく理由はここにあり,両側への拡張はそのためだけに借りている.

期待値のレベルでは何が起きているか

命題 3.5.9. 任意の情報源とに対しである.

証明. 定義 1.1.1 よりである.両辺をで割れば,右辺は定義 3.5.1の期待値である.

形式化: expected_blockLogAvg_eq (ソース)

命題 3.5.9 は定義を書き直しただけである.それでも役に立つのは,定理 3.2.6 と合わせるとが直ちに出るからである.つまり「平均すればに近づく」ところまでは,本節の道具を一つも使わずに言える.定理 3.5.2 が言い足しているのは「1 本 1 本の実現でもそうなる」という一点であり,Birkhoff の定理も両側への拡張も,その一点のためだけに要る.

第2章の符号化定理はどこまで生き延びるか

系 3.5.10. とし,定常情報源に対して定義 2.3.1 のブロック情報源符号を考える.

  1. 情報源がさらにエルゴード的でならば,かつを満たす符号の族が存在する.
  2. 逆に,符号の族がを満たし,レートの列が上に有界ならばである.

証明. 主張 1 の達成可能性を示す.定義 2.2.1 の典型集合を,に,をブロックの確率に読み替えて同じ式で定める.定理 3.5.2 は概収束を述べているので,2.1 節で借りた「概収束は確率収束を含む」より確率収束も従い,定理 2.2.3 の証明は定理 2.1.4 の確率収束しか使っていないから,そのまま通る.定理 2.2.4 は典型集合の定義の書き直しにすぎず,定理 2.2.5定理 2.2.3定理 2.2.4 と「確率の総和が 1」しか使わない.したがって定理 2.3.2 の符号の構成がそのまま働く.

主張 2 の逆に移る.定理 2.3.4 の証明が i.i.d. を使っているのは,補題 2.3.3と書き換える 1 箇所だけである.そこをと読み替えると残りはそのまま通って,各

𝐻𝑛𝑛𝑅𝑛+𝛿𝑛,𝛿𝑛0

を得る.あとはを言えばよい.補題 3.2.4 より数列は非増加であり,定理 3.2.6 よりその極限がだから,各項は以上である.補題 3.2.3 よりはその相加平均だから,やはり以上である.よってとなり,下極限をとれば主張を得る.

系 3.5.10 は,情報源符号化定理(定理 2.3.6)が i.i.d. の仮定なしに成り立つ,と言っている.置き換わったのは 1 文字のエントロピーがエントロピーレートになったところだけで,圧縮の限界を決めているのが「1 文字あたりの不確かさ」であることは変わらない.

エルゴード性がどこで効くか. 逆定理 (2) の証明はエルゴード性を使っていない.定常でありさえすればが言えるからである.エルゴード性が要るのは達成可能性 (1) のほう,つまり定理 3.5.2 を経由して典型集合を作るところである.例 3.4.2 がそこを壊す.最初のコインを固定すればその先は i.i.d. で,混合の確率のうち自分の側の項が支配するから,経験エントロピーは定理 2.1.4 のとおり分岐ごとの値に落ち着く.すなわち表ならビット,裏ならビットである.ただ一つの値に収束するという定理 3.5.2 の結論が,ここでは成り立っていない.いっぽうエントロピーレートは,命題 3.5.9から二つの値の平均ビットになる.したがってにとった符号は,表の側の実現を覆いきれない.長さのブロックが通りあるのに,用意した符号語は個しかないためである.誤り確率はに落ちず,近くに残る.

形式化上の注記. 系 3.5.10 に対応する単独の宣言は形式化されていない.エルゴード的な定常情報源に対する典型集合とブロック情報源符号は,第2章の i.i.d. 版とは別に立てる必要があり,そこまでは形式化されていないためである.本文の証明が第2章の主張を「そのまま通る」と書いているのは紙の上での議論で,機械検証されているのは定理 3.5.2 までである.

第2章はエントロピーを,分布から計算する量・長い系列が見せる振る舞い・圧縮率の下限,という三つの顔で捉えた.本章はその三つがそのままエントロピーレートに引き継がれることを見たことになる.第6章では情報源ではなく通信路(信号を受け取って別の信号を返すもの)に目を移し,そこを通せる情報の量にも同じように限界があることを見る.

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