2.2 典型集合
定理 2.1.4 は「実際に出る系列の確率は2−𝑛𝐻近辺にある」と言っていた.この「近辺にある」を満たす系列だけを集めた集合を作り,それがどれくらい大きいかを数えるのが本節である.数え上げの結果(約2𝑛𝐻個)が,次節のデータ圧縮の限界をそのまま与える.
以下,周辺分布は全アルファベット上で正,すなわち𝑝(𝑎) >0(𝑎 ∈X)と仮定する.これは情報源への制約ではなく記法の整理である.𝑝(𝑎) =0の文字は𝑋𝑖が決してとらないのだから,Xからあらかじめ除いておけばよい.
定義
定義 2.2.1(典型集合). 𝜀 >0と𝑛 ≥1に対し,長さ𝑛の 典型集合 を
𝑇(𝑛)𝜀:={𝑥∈X𝑛:∣−1𝑛log𝑝(𝑥)−𝐻(𝑋)∣<𝜀}で定める.𝑇(𝑛)𝜀の元を 典型系列 という.
定義は経験エントロピーの条件をそのまま集合に書き直したものである.𝑥が典型であるとは,その確率が
2−𝑛(𝐻(𝑋)+𝜀)<𝑝(𝑥)<2−𝑛(𝐻(𝑋)−𝜀)の範囲にあることにほかならない(両辺のlogをとって−1𝑛倍すれば定義に戻る).つまり典型集合とは「確率がこの幅に収まっている系列の集まり」であって,系列の見た目(どの文字がどこに並んでいるか)には一切触れていない.
例 2.2.2(偏ったコインの典型集合). 例 2.1.2 の情報源(𝑝(1) =0.9,𝑝(0) =0.1,𝐻(𝑋) ≈0.469)で,長さ𝑛のブロック𝑥に1が𝑞𝑛個含まれるとする.すると
−1𝑛log𝑝(𝑥)=−𝑞log0.9−(1−𝑞)log0.1≈0.152𝑞+3.322(1−𝑞)であり,これが𝐻(𝑋) ≈0.469に近いことと,割合𝑞が0.9に近いことが同値になる(左辺は𝑞の一次式で,傾きの絶対値は3.322 −0.152 =3.17だから,|𝑞 −0.9| <𝜀/3.17が典型性の条件である).
見落としやすい点が一つある.単独でいちばん確率の高いブロックは1ばかりが並んだ(1,1,…,1)で,その確率は0.9𝑛である.しかしこの系列の割合は𝑞 =1なので,𝜀が小さければ典型集合に入らない.
例 2.2.2 の最後の注意は,典型集合の性質を一つに凝縮している.典型集合は「起こりやすい系列の集まり」ではない.個々の系列としては,1ばかりの系列のほうが典型系列より確率が高い.それでも次に見るように,ブロックが典型集合に落ちる確率は 1 に近づく.確率のやや低い系列が,数の多さで圧倒するからである.典型集合が捉えているのは「1 個あたりの確率」ではなく「確率×個数」の釣り合いのほうである.
性質1:典型集合に入る確率は 1 に近づく
定理 2.2.3. 任意の𝜀 >0に対しPr[𝑋𝑛 ∈𝑇(𝑛)𝜀] →1(𝑛 →∞).
証明. 定義 2.2.1 より,事象{𝑋𝑛 ∈𝑇(𝑛)𝜀}は{ |ˆ𝐻𝑛 −𝐻(𝑋)| <𝜀 }そのものである.定理 2.1.4 の確率収束は,その余事象の確率Pr[ |ˆ𝐻𝑛 −𝐻(𝑋)| ≥𝜀 ]が 0 に収束すると述べているから,もとの事象の確率は 1 に収束する.◼
性質2:典型系列の確率はどれもほぼ等しい
定理 2.2.4. 𝜀 >0とする.𝑥 ∈𝑇(𝑛)𝜀ならば
2−𝑛(𝐻(𝑋)+𝜀)≤𝑝(𝑥)≤2−𝑛(𝐻(𝑋)−𝜀).
証明. 𝑥 ∈𝑇(𝑛)𝜀は定義より∣ −1𝑛log𝑝(𝑥) −𝐻(𝑋)∣ <𝜀,すなわち
𝑛(𝐻(𝑋)−𝜀)<−log𝑝(𝑥)<𝑛(𝐻(𝑋)+𝜀)と同値である.各辺に−1を掛けて不等号の向きを反転させ,2を底とする指数をとれば主張を得る(実際には狭義の不等号で成り立つが,以降は等号込みの形で足りる).◼
定理 2.2.4 が「漸近等分配性」という名前の由来である.典型系列はどれも確率が2−𝑛𝐻の同じ幅の中に収まっている.すなわち全確率のほぼ全部が,たがいにほぼ等確率な系列のあいだで分け合われ,系列ごとの差は指数の肩の𝜀のぶんしか残らない.𝑛が大きいとき2±𝑛𝜀という因子は決して小さくないが,𝜀はいくらでも小さくとれるので,1𝑛log𝑝(𝑥)という尺度で見るかぎり典型系列は互いに区別がつかない.
性質3:典型集合の大きさは約2𝑛𝐻
次の評価は上下で形が揃っていない.上界は確率の総和が 1 だという勘定だけで出るのでどの𝑛でも成り立つが,下界は「ブロックがほぼ確実に典型集合に入る」という漸近的な事実(定理 2.2.3)を使うので,𝑛が大きいときにしか言えない.しかもその確率は 1 に近づくだけで 1 には届かないから,取りこぼしのぶんの因子1 −𝜂が残る.
定理 2.2.5. 𝜀 >0とする.任意の𝑛に対し
∣𝑇(𝑛)𝜀∣≤2𝑛(𝐻(𝑋)+𝜀).また任意の𝜂 ∈(0,1)に対し,𝑛が十分大きければ
∣𝑇(𝑛)𝜀∣≥(1−𝜂)2𝑛(𝐻(𝑋)−𝜀).
証明. 上界を示す.全系列にわたる確率の総和は 1 だから,典型系列だけの和はそれ以下である.定理 2.2.4 の下界を各項に当てると
1=∑𝑥∈X𝑛𝑝(𝑥)≥∑𝑥∈𝑇(𝑛)𝜀𝑝(𝑥)≥∣𝑇(𝑛)𝜀∣⋅2−𝑛(𝐻(𝑋)+𝜀).両辺に2𝑛(𝐻(𝑋)+𝜀)を掛ければ上界を得る.
下界に移る.定理 2.2.3 より,𝑛を十分大きくとればPr[𝑋𝑛 ∈𝑇(𝑛)𝜀] ≥1 −𝜂にできる.この確率は典型系列の確率の和であり,定理 2.2.4 の上界を各項に当てると
1−𝜂≤∑𝑥∈𝑇(𝑛)𝜀𝑝(𝑥)≤∣𝑇(𝑛)𝜀∣⋅2−𝑛(𝐻(𝑋)−𝜀).両辺に2𝑛(𝐻(𝑋)−𝜀)を掛ければ下界を得る.◼
証明の要点は一つだけである.確率の和を「個数× 1 個あたりの確率」で挟むと,確率の側は定理 2.2.4 が押さえてくれるので,残った個数の側が決まる.上界と下界で使う向きが逆になっているのに注意したい.個数を上から抑えるには 1 個あたりの確率を下から,下から抑えるには上から評価する必要があり,定理 2.2.4 が両側評価であることが両方に効いている.
規模感. 全系列は|X|𝑛 =2𝑛log|X|個ある.定理 1.1.5 により𝐻(𝑋) ≤log|X|であり,等号は一様分布のときに限る.分布が一様でなければ肩の差は正なので,𝜀をそれより小さくとれば,典型集合の割合はおよそ2−𝑛(log|X|−𝐻(𝑋)−𝜀)で指数的に 0 へ向かう(𝜀 →0の極限が2−𝑛(log|X|−𝐻(𝑋))である).例 2.2.2 の偏ったコインなら,全系列2𝑛個に対し典型集合は約20.469𝑛個で,割合は2−0.531𝑛となる.𝑛 =100なら10−16程度である.それでもブロックはほぼ確実にその中に落ちる.定理 2.2.3 と定理 2.2.5 は,この一見矛盾する二つの事実を同時に述べている.
要素数を符号長に読み替える. 典型集合の要素数が約2𝑛𝐻なら,その各元に高々2𝑛(𝐻+𝜀)通りの通し番号を振れる.番号を書き下すのに必要なのは𝑛(𝐻 +𝜀)ビット,1 文字あたり𝐻 +𝜀ビットである.典型集合に入らない系列は捨ててよい.定理 2.2.3 がその確率を 0 に押し込むからである.これが次節の圧縮方式そのものであり,定理 2.2.5 の上界がそのままレートの上界になる.
InformationTheory — 形式化検証つき情報理論教科書(レビュー版).数式は MathJax + AMS Euler で事前レンダリング.