1.1 エントロピー
確率変数𝑋を観測する前,その結果にはどれだけの「不確かさ」があるだろうか.あるいは観測したとき,私たちはどれだけの「情報」を得るのだろうか.エントロピーは,この素朴な問いに一つの定量的な答えを与える量である.
良い不確かさの尺度に何を期待するか,条件を並べてみよう.確実な事象(ある値を確率 1 でとる)は不確かさ 0 であってほしい.とりうる値が多く,かつそれらが均等に起こりやすいほど不確かさは大きいはずだ.そして独立な二つの実験を同時に行ったときの不確かさは,それぞれの不確かさの和であってほしい(加法性).
この最後の要請が,式の形をほとんど決めてしまう.まず,確率𝑝で起きる一つの事象を観測したときに得られる情報の量を𝑠(𝑝)と書こう.めったに起きないことが起きたときほど驚きは大きいので𝑠は減少関数で,確実に起きること(𝑝 =1)には驚きがないから𝑠(1) =0である.ここに加法性を課す.独立な二つの実験で確率𝑝の事象と確率𝑞の事象が同時に起きたとき,その同時事象の確率は𝑝𝑞であり,得た情報は二つの情報の和であってほしい:
𝑠(𝑝𝑞)=𝑠(𝑝)+𝑠(𝑞).「掛け算を足し算に変える」というこの性質を持つ単調な関数は,定数倍を除いて対数しかない.これは Cauchy の関数方程式として知られる古典的な事実で,本書では証明しない(単調性は上で𝑠に課した減少性がそのまま与える).よって𝑠(𝑝) = −log𝑝と書ける(0 <𝑝 ≤1で−log𝑝 ≥0となるよう符号を選んだ.定数倍の自由度は対数の底の選び方に吸収される).
𝑠(𝑝) = −log𝑝は「その値が出たときの驚き」であって,まだ観測前の不確かさではない.観測前にわかっているのは分布だけなので,実際に出る値ごとの驚きを,その値が出る確率で平均する.これが次のエントロピーである.なお Shannon は,いま並べた三つの要請(確実な事象で 0・値が多く均等なほど大・加法性)に連続性を加えると,尺度が定数倍を除いてこの形に限ることを示した(Shannon の公理的特徴づけ).
定義
定義 1.1.1(エントロピー). 有限アルファベットX上に分布する離散確率変数𝑋(分布𝑝(𝑥))の エントロピー を
𝐻(𝑋)=−∑𝑥∈X𝑝(𝑥)log𝑝(𝑥)で定める.𝑝(𝑥) =0の項は𝑡log𝑡 →0(𝑡 ↓0)にならい0log0 =0と約束する.
はじめに見たとおり,−log𝑝(𝑥)は「𝑥が出たときの驚き」であり,𝐻(𝑋)はその驚きを𝑋自身の分布で平均した 平均驚き量 である.∑𝑥𝑝(𝑥) ⋅(驚き)という形をしているだけで,新しい要素は何もない.分布が一点に集中していれば出る値は決まっていて驚きが起きず𝐻 =0,分布が平らに広がっているほどどの値が出てもそこそこ驚くので𝐻は大きくなる.後で「符号長」として読み直すこともできる.確率𝑝(𝑥)の記号におよそ−log𝑝(𝑥)文字ぶんの符号語を割り当てるのが最適で,そのとき平均符号長がちょうど𝐻(𝑋)になる.なぜそれが最適なのかは第2章の情報源符号化定理で示す.𝑞を信じて符号化する人なら,同じ読み方で払う平均符号長は−∑𝑥𝑝(𝑥)log𝑞(𝑥)になる.いまは先の見通しとして聞き流してよい.
エントロピーは𝑋の値そのものではなく分布𝑝のみに依存する量であることに注意したい.𝑋の値を別の名前に読み替えても(たとえば目の数1,…,6を色の名前に置き換えても)𝐻(𝑋)は変わらない.記号𝐻(𝑋)は慣用だが,より正確には𝐻(𝑝)と書くべき汎関数である.
例
例 1.1.2(ベルヌーイ分布). 𝑋が確率𝑝で1,確率1 −𝑝で0をとるとき,
𝐻(𝑋)=−𝑝log𝑝−(1−𝑝)log(1−𝑝)=:𝐻𝑏(𝑝).この𝐻𝑏を 二値エントロピー関数 と呼ぶ.𝐻𝑏(0) =𝐻𝑏(1) =0(結果が確定していて不確かさがない)であり,𝑝 =1/2で最大値log2をとる(最も予測しづらい公平なコイン).𝐻𝑏は[0,1]上で凹かつ𝑝 =1/2を軸に左右対称である.二値エントロピーは 1.10 節のファノの不等式でふたたび主役として現れる.
例 1.1.3(一様分布). 𝑋がX(要素数|X| =𝑀)上で一様,すなわち各𝑥で𝑝(𝑥) =1/𝑀のとき,
𝐻(𝑋)=−∑𝑥1𝑀log1𝑀=log𝑀.どの値も等しく確率1/𝑀なので驚きはどれもlog𝑀で一定であり,平均をとってもlog𝑀のままである.底を 2 にとればlog2𝑀ビットになる.これは𝑀 =2𝑘個の等確率な結果を区別するのに𝑘ビット必要,という素朴な勘定と一致する.あとで見るように(定理 1.1.5),これは要素数𝑀のアルファベット上で達成可能なエントロピーの最大値 である.
以降,単位が要る箇所ではlogの底を 2 にとり,log2𝑀の値を ビット と呼ぶ.底を1より大きい範囲で変えても以下のすべての等式・不等式はそのまま成り立つ(単位の名前だけが変わる).
規模感. 公平なコイン 1 枚はlog2(1 ビット),公平なサイコロ 1 個はlog6(約 2.58 ビット),英字 26 文字が一様ならlog26(約 4.70 ビット).実際の英文は文字の出現が偏り,さらに前後の文字に依存するので,1 文字あたりの不確かさはこれよりずっと小さい.その「偏りと依存でどれだけ下がるか」を測るのが本章の残りと第3章の主題である.
性質1:非負性
命題 1.1.4. 定義 1.1.1 の設定(有限アルファベットX上の𝑋)で𝐻(𝑋) ≥0.等号は𝑋がある一点に確率 1 で集中するときに限り成り立つ.
証明. 各𝑥について0 ≤𝑝(𝑥) ≤1だからlog𝑝(𝑥) ≤0,したがって−𝑝(𝑥)log𝑝(𝑥) ≥0.和をとっても非負である(𝑝(𝑥) =0の項は約束により 0).
等号𝐻(𝑋) =0は,非負項の和が 0 であることだから,すべての𝑥で−𝑝(𝑥)log𝑝(𝑥) =0を要する.0 <𝑝(𝑥) <1の𝑥があれば−𝑝(𝑥)log𝑝(𝑥) >0となり矛盾.ゆえに各𝑝(𝑥)は0または1であり,総和が 1 という条件と合わせると,ちょうど一つの𝑥で𝑝(𝑥) =1.◼
命題 1.1.4 は「不確かさが負になることはなく,不確かさ 0 とは結果が確定していること」という当然の事実を述べている.当然に見えるが,これは「驚きは負にならない」(0 ≤𝑝(𝑥) ≤1なので−log𝑝(𝑥) ≥0)ことの言い換えであり,平均をとる前の各項がすでに非負である,という構造から来ている.逆に𝐻(𝑋) =0なら「どの値が出ても驚かない」=「出る値が決まっている」ということになる.
性質2:一様分布が最大化する(𝐻(𝑋) ≤log|X|)
定理 1.1.5. |X| =𝑀のとき𝐻(𝑋) ≤log𝑀.等号は𝑋がX上で一様分布のときに限り成り立つ.
この定理は「アルファベットの大きさが不確かさの上限を決め,その上限は均等に散らばったときに達成される」という主張である.偏りは予測の手がかりであり,手がかりがあるぶんだけ不確かさは減る.だから偏りのない一様分布がいちばん予測しづらい,と読める.例 1.1.2 の二値エントロピー𝐻𝑏が𝑝 =1/2で最大になったのは,この定理の𝑀 =2の場合にほかならない.
これは情報理論で繰り返し使われる基本的な上界で,たとえば「𝑀個の記号で送れる情報量は 1 記号あたり高々log𝑀」という符号化の素朴な限界を裏づける.裏返せば,log𝑀と実際の𝐻(𝑋)の差が「偏りのぶんだけ圧縮できる余地」であり,この差をきちんと測る量が 1.6 節で出てくる.
本書では,この定理を予備知識なしの数学的帰納法で証明する.必要な解析的事実は凹性(下に張り出さず,弦が常にグラフの下に来ること)だけで,Jensen の不等式のような道具は既知として持ち込まず,その有限版を帰納法でその場で組み立てる.
使うのは次の二つの関数についての事実で,どれも本節で証明する.
𝜑(𝑡)=−𝑡log𝑡(𝑡≥0, 𝜑(0):=0),log𝑡(𝑡>0).証明のために微積分から二つの判定法を借りる.一つめは「開区間の上で 2 階微分が負なら,その区間で狭義凹である」,二つめは「微分可能な狭義凹関数のグラフは,接点を除いて接線より真に下にある」である.当てる相手はどちらも上の二つの関数だけで,依存するのは 補題 1.1.6・補題 1.1.7・補題 1.1.8 の証明だけである.初等関数の微分の計算規則と,連続関数の積が連続であることも既知とする.対数そのものについては,底を1より大きくとることと,log1 =0,log𝑠𝑎 =𝑎log𝑠(𝑠 >0)という定義から直ちに出る規則だけを使う.本節が借りるのはこれだけで,logと𝜑の性質はすべてそこから出す.
補題 1.1.6(対数の狭義凹性). logは(0,∞)の上で狭義凹である.
証明. log𝑡の 2 階微分は−(log𝑒)/𝑡2である.底は1より大きいのでlog𝑒 >0であり,𝑡 >0では 2 階微分が負である.借りた一つめの判定法から狭義凹性が従う.◻
補題 1.1.7(対数の接線不等式). どの𝑡 >0についても
log𝑡≤(𝑡−1)log𝑒が成り立ち,等号は𝑡 =1に限る.
証明. log𝑡の導関数は(log𝑒)/𝑡だから,𝑡 =1での接線は𝑦 =(𝑡 −1)log𝑒である.補題 1.1.6 と,借りた二つめの判定法(微分可能な狭義凹関数のグラフは接点を除いて接線より真に下にある)を合わせると,接点𝑡 =1を除いて真の不等号が成り立つ.◻
底を𝑒にとればlog𝑒 =1なので右辺は𝑡 −1になるが,単位が要る箇所でとる底2ではlog𝑒 =1.4426…で,因子を落とすと不等式は偽になる(𝑡 =1.5でlog21.5 =0.585… >0.5).底を1より大きくとるかぎりlog𝑒 >0なので,この不等式に非負の重みを掛けて足すという使い方をするかぎり,因子は結論の向きを変えない.この一本は 1.6 節以降で繰り返し使う.命題 1.1.4 の証明で使ったlogの符号もここから出る.0 <𝑡 <1なら𝑡 −1 <0とlog𝑒 >0からlog𝑡 <0である.
補題 1.1.8(𝜑の連続性と狭義凹性). 𝜑(𝑡) = −𝑡log𝑡(𝑡 ≥0,𝜑(0) :=0)は[0,∞)の全体で連続であり,[0,∞)の上で狭義凹である(端点0を含む弦についても不等号は狭義である).
証明. 連続性から見る.𝑡 >0では連続関数の積だから連続である.𝑡 =0での連続性を示す.0 <𝑡 <1では 補題 1.1.7 よりlog𝑡 ≤(𝑡 −1)log𝑒 <0だから𝜑(𝑡) >0である.補題 1.1.7 を𝑡−1/2 >1に当てると−12log𝑡 ≤(𝑡−1/2 −1)log𝑒 ≤𝑡−1/2log𝑒であり,両辺に2𝑡 >0を掛けて
0<𝜑(𝑡)=−𝑡log𝑡≤2√𝑡log𝑒(0<𝑡<1)を得る.右辺は𝑡 ↓0で0に収束するから,挟み撃ちにより𝜑(𝑡) →0 =𝜑(0)である.
(0,∞)での狭義凹性に移る.𝜑の 2 階微分は−(log𝑒)/𝑡であり,log𝑒 >0と𝑡 >0から負である.借りた一つめの判定法による.
端点0を含む弦に移る.𝜑は𝑡 =0で微分可能でないので,この弦は判定法の外にある.𝑏 >0を固定し,弦との差𝑔(𝑡) :=𝜑(𝑡) −𝜑(0) −𝑡(𝜑(𝑏) −𝜑(0))/𝑏を[0,𝑏]の上で考える.𝑔はいま示した連続性から連続で,𝜑から一次式を引いただけだから(0,𝑏)で狭義凹であり,𝑔(0) =𝑔(𝑏) =0である.𝑡 ∈(0,𝑏)と0 <𝜀 <𝑡に対し,[𝜀,𝑏]上の凹性から𝑔(𝑡)は𝑔(𝜀)と𝑔(𝑏)を結ぶ線分の𝑡での値以上であり,その値は𝜀 ↓0で0に収束するから,𝑔 ≥0が[0,𝑏]の全体で成り立つ.極限が保つのは広義の不等号だけなので,狭義性はここから改めて出す.𝑐 ∈(0,𝑏)に対して0 <𝑐1 <𝑐 <𝑐2 <𝑏と𝑐 =𝜆𝑐1 +(1 −𝜆)𝑐2(0 <𝜆 <1)をとれば,𝑔の狭義凹性から𝑔(𝑐) >𝜆 𝑔(𝑐1) +(1 −𝜆) 𝑔(𝑐2) ≥0となり,端点0を含む弦についても狭義の不等式が成り立つ.◻
端点込みの狭義性は,定理 1.1.5 の等号条件で確率 0 の記号を排除するのに実際に効く.補題 1.1.8 の凹性を 2 点の形で書くと,𝑎,𝑏 ∈[0,∞)と𝜆 ∈[0,1]に対して
𝜑(𝜆𝑎+(1−𝜆)𝑏)≥𝜆𝜑(𝑎)+(1−𝜆)𝜑(𝑏),(∗)であり,𝑎 ≠𝑏かつ0 <𝜆 <1のとき不等号は狭義である.
補題 1.1.9(有限 Jensen の不等式・帰納法版). 𝜑を区間𝐼上の凹関数とする.有限個の点𝑡1,…,𝑡𝑛 ∈𝐼と重み𝑤1,…,𝑤𝑛 ≥0(∑𝑖𝑤𝑖 =1)に対して
𝑛∑𝑖=1𝑤𝑖𝜑(𝑡𝑖)≤𝜑(𝑛∑𝑖=1𝑤𝑖𝑡𝑖).さらに𝜑が狭義凹なら,等号は正の重みをもつ点𝑡𝑖がすべて互いに等しいとき,かつそのときに限る.
証明(𝑛についての数学的帰納法). 𝑛 =1のときは𝑤1 =1なので両辺とも𝜑(𝑡1)であり,等号が成り立って条件も自明である.
𝑛 =2のときはまさに凹性の 2 点形( ∗)そのものである(𝜆 =𝑤1,1 −𝜆 =𝑤2).狭義凹なら𝑡1 ≠𝑡2かつ両重み正で狭義不等号になる.
以下,𝑛点で補題が成り立つと仮定して𝑛 +1点を示す.重み𝑤1,…,𝑤𝑛+1 ≥0(和 1)が与えられたとき,もし𝑤𝑛+1 =1なら他は 0 で両辺𝜑(𝑡𝑛+1)となり成立.そうでなければ𝑊 :=1 −𝑤𝑛+1 =∑𝑛𝑖=1𝑤𝑖 >0とおき,𝑣𝑖 :=𝑤𝑖/𝑊(𝑖 =1,…,𝑛,∑𝑖𝑣𝑖 =1)と正規化する.内側の凸結合𝑠 :=∑𝑛𝑖=1𝑣𝑖𝑡𝑖は𝐼に属する(区間は凸結合で閉じている).全体の重心は
𝑛+1∑𝑖=1𝑤𝑖𝑡𝑖=𝑊𝑠+𝑤𝑛+1𝑡𝑛+1と 2 点𝑠,𝑡𝑛+1の凸結合(重み𝑊,𝑤𝑛+1)に書ける.ここで 2 点凹性( ∗)と𝑛点への帰納法の仮定 𝜑(𝑠) ≥∑𝑛𝑖=1𝑣𝑖𝜑(𝑡𝑖)を順に使うと
𝜑(𝑛+1∑𝑖=1𝑤𝑖𝑡𝑖)≥𝑊𝜑(𝑠)+𝑤𝑛+1𝜑(𝑡𝑛+1)≥𝑊𝑛∑𝑖=1𝑣𝑖𝜑(𝑡𝑖)+𝑤𝑛+1𝜑(𝑡𝑛+1)=𝑛+1∑𝑖=1𝑤𝑖𝜑(𝑡𝑖),最後の等号は𝑊𝑣𝑖 =𝑤𝑖による.これで𝑛 +1点でも不等式が成り立つ.
等号条件(狭義凹). 上の連鎖で 2 箇所の≥がいずれも等号になるときのみ全体が等号になる.
𝑤𝑛+1 =0のときは,( ∗)の重みの一方が 0 なので第 1 の≥はつねに等号であり,条件は第 2 の不等式だけから決まる.帰納法の仮定の等号条件により,それは「正の重みをもつ𝑡𝑖 (𝑖 ≤𝑛)が互いに等しい」(𝑡𝑛+1は重み 0 なので条件に入らない)であって,主張のとおりである.
𝑤𝑛+1 >0のときを見る.第 1 の( ∗)が等号になるのは,𝜑が狭義凹で両重み𝑊,𝑤𝑛+1がともに正だから,𝑠 =𝑡𝑛+1のとき,かつそのときに限る.第 2 の不等式が等号になるのは,帰納法の仮定により正の重みをもつ𝑡𝑖 (𝑖 ≤𝑛)が互いに等しいとき,かつそのときに限る.後者が成り立つとき,その共通の値を𝑐と書けば𝑠 =∑𝑖𝑣𝑖𝑡𝑖 =𝑐である(𝑣𝑖は𝑖 ≤𝑛上の重みで総和 1,正の重みをもつ点はすべて𝑐).したがって前者の条件𝑠 =𝑡𝑛+1は𝑐 =𝑡𝑛+1と同じことになり,二つを合わせると「正の重みをもつ点が𝑡𝑛+1を含めてすべて一致する」となる.◻
証明(定理 1.1.5). エントロピーは定義から𝜑の和である:
𝐻(𝑋)=−∑𝑥𝑝(𝑥)log𝑝(𝑥)=∑𝑥𝜑(𝑝(𝑥)).補題 1.1.9 を,重み𝑤𝑥 =1/𝑀(一様),点𝑡𝑥 =𝑝(𝑥)ととって適用する.各𝑝(𝑥) ≥0は𝜑の定義域[0,∞)に入るので,質量 0 の記号を場合分けして除く必要もない.∑𝑥(1/𝑀)𝑝(𝑥) =1/𝑀だから
1𝑀∑𝑥𝜑(𝑝(𝑥))≤𝜑(∑𝑥1𝑀𝑝(𝑥))=𝜑(1𝑀)=−1𝑀log1𝑀=1𝑀log𝑀.両辺を𝑀倍すれば𝐻(𝑋) ≤log𝑀.等号は𝜑が狭義凹なので,補題の等号条件により,正の重み(ここではすべての𝑥で𝑤𝑥 =1/𝑀 >0)をもつ点𝑝(𝑥)が𝑥によらず一定のとき,すなわち∑𝑥𝑝(𝑥) =1と合わせて𝑝(𝑥) =1/𝑀,つまり𝑋が一様分布のときに限る.◼
InformationTheory — 形式化検証つき情報理論教科書(レビュー版).数式は MathJax + AMS Euler で事前レンダリング.