7.1 微分エントロピー
第1章から第6章まで,エントロピーの引数に置いた確率変数はつねに有限個の値をとるものだった.エントロピーの定義(定義 1.1.1)は記号ごとの確率を足し合わせる形をしていて,記号が有限個であることをその形に織り込んでいる.しかし実際に測る量の多くは実数である.電圧,到着時刻,位置.これらを情報理論の言葉で扱うには,和を積分に置き換えた量が要る.本章はその量を定義し,離散のエントロピーと何が同じで何が違うのかを見る.
置き換えそのものは素直である.確率𝑝(𝑥)を密度𝑓(𝑥)に,和∑𝑥を積分∫𝑑𝑥に取り替えればよい.ただし,そうしてできた量は離散のエントロピーの単なる連続版ではない.本節の終わりで見るように,実数を有限の精度で読み取った値のエントロピーは精度を上げると発散し,新しい量はその発散する部分を差し引いた残りにあたる.差し引いた残りなので負にもなる.
定義
定義 7.1.1(確率密度関数). 𝑛 ≥1とする.𝑓 :ℝ𝑛 →[0,∞)が積分できて∫ℝ𝑛𝑓(𝑥) 𝑑𝑥 =1を満たすとき,𝑓をℝ𝑛上の 密度 という.ℝ𝑛に値をとる確率変数𝑋とℝ𝑛上の密度𝑓について,任意の直方体𝑅 =[𝑎1,𝑏1] ×⋯ ×[𝑎𝑛,𝑏𝑛]で
Pr[𝑋∈𝑅]=∫𝑅𝑓(𝑥)𝑑𝑥が成り立つとき,𝑓を𝑋の 確率密度関数 といい,本書では以下,単に𝑋の密度と呼ぶ.𝑛 ≥2のときこれを 結合密度 とも呼ぶ.
定義 7.1.1 が積分を置いたのはℝ𝑛の全体と直方体の上だけだが,以下では{𝑥 :𝑓(𝑥) >0}のように条件で定めた集合𝐴の上での積分も使う.これは,𝐴の外で0とおいた関数のℝ𝑛全体での積分と読む.本章に現れる集合はすべてこの読み方ができるものとし,積分の線形性,非負な関数についての単調性(𝐴 ⊆𝐵ならば∫𝐴 ≤∫𝐵),および絶対収束する積分を項ごとに分けてよいことを,微積分の計算規則として既知とする.同じ資格で,有界閉区間の上の連続関数はその絶対値も積分できること,部分積分と,そこに現れる境界項の極限(多項式と減衰する指数関数の積が𝑥 →∞で0に近づくこと,𝑡2log𝑡が𝑡 ↓0で0に近づくこと),log𝑡の導関数が(log𝑒)/𝑡であること,および有限な値に収束する広義積分の扱いも既知とする.当てる対象は,本章に現れる密度と,それに多項式や対数を掛けて作る関数,および例 7.1.5 の証明でそれらを上下から抑えるのに使うlog𝑥とloglog𝑥である.第8章もこれらを同じ資格で使う.そちらでの当てる対象は,第8章に現れる密度とそれに多項式や対数を掛けて作る関数,および 8.4 節の命題 8.4.1 に現れる球と立方体の上での定数関数である.第16章 16.7 節も同じ資格でこれを借り直すが,そちらの一覧は本節のものと同じではない.本節が挙げた規則のうち引き継ぐのは積分の線形性・log𝑡の導関数・部分積分・𝑡2log𝑡の極限・広義積分の扱いの五つで,そこに合成関数の微分・指数関数の導関数・微積分の基本定理を加えている.そちらでの当てる対象は,16.7 節と 16.8 節に現れる密度と,それに多項式や対数を掛けて作る関数のほかに,16.7 節が雑音の分散を変数として書く1次式の対数とその逆数を含む.
確率と密度をつなぐ事実も二つ借りる.一つめは 期待値の積分表示 で,ℝ𝑛に値をとり密度𝑓をもつ確率変数𝑋と実数値関数𝑢について𝔼[𝑢(𝑋)] =∫ℝ𝑛𝑢(𝑥)𝑓(𝑥) 𝑑𝑥であり,一方が絶対収束すればもう一方も絶対収束する,という事実である.当てる対象は,𝑢が密度の対数(密度が正の点でその対数,そのほかの点で0ととったもの)である場合と,𝑢が2次以下の多項式である場合だけである.二つめは 体積0の集合の確率 で,𝑛次元の体積が0である集合の上に,密度をもつ確率変数は確率を置かない,という事実である.当てる対象は,7.4 節の終わりに出てくる平面の対角線と,第8章 8.2 節が電力0の制約を扱うところに出てくる実軸の1点,および 8.3 節の初めに出てくる「第1座標の組が有限個の点のどれかに等しい」という集合の三つである.どちらも定義 7.1.1 からその場で読み取れる一歩ではない.定義 7.1.1 が確率と積分を結びつけたのは,直方体の上でだけだからである.
定義 7.1.2(微分エントロピー). 𝑛 ≥1とし,𝑋をℝ𝑛に値をとり密度𝑓をもつ確率変数とする.∫ℝ𝑛|𝑓(𝑥)log𝑓(𝑥)| 𝑑𝑥 <∞であるとき,𝑋の 微分エントロピー を
ℎ(𝑋):=−∫ℝ𝑛𝑓(𝑥)log𝑓(𝑥)𝑑𝑥で定め,このときℎ(𝑋)が 定まる という.𝑓(𝑥) =0の点では定義 1.1.1 と同じく𝑓(𝑥)log𝑓(𝑥) =0と約束する.𝑋 =(𝑋1,…,𝑋𝑛)と成分で書くときはℎ(𝑋1,…,𝑋𝑛)とも書き,𝑛 ≥2のときこれを 結合微分エントロピー と呼ぶ.
式の形は,離散のときとまったく同じである.借りた期待値の積分表示によりℎ(𝑋) =𝔼[ −log𝑓(𝑋)]と読めて,−log𝑓(𝑥)を「その点で密度が薄いこと」の目盛とし,それを𝑋自身の分布で平均している.離散のエントロピー𝐻(𝑋) =𝔼[ −log𝑝(𝑋)]と並べると,中身の𝑝が𝑓に替わっただけである.
絶対値をつけた積分が有限であることを要求したのは,ℎ(𝑋)を一つの実数として確定させるためである.離散のときは項が有限個だったので値はつねに実数だったが,積分では正の部分と負の部分の一方または両方が発散しうる.裾の重い密度ではそれが実際に起きて,ℎ(𝑋)を定めない場合がある.例 7.1.5 がその一つである.
替わったことの帰結は小さくない.確率𝑝(𝑥)はつねに[0,1]に入るので−log𝑝(𝑥) ≥0であり,平均をとった𝐻(𝑋)も非負だった(命題 1.1.4).密度𝑓(𝑥)にはその上限がなく,1を超えることができる.そこでは−log𝑓(𝑥)が負になり,ℎ(𝑋)も負になりうる.文字を𝐻ではなくℎにするのはこのためで,二つは同じ量の離散版と連続版ではない.
ℎ(𝑋)が𝑋の値でなく密度だけで決まるのは,離散のときと同じである.なお同じ𝑋の密度は一通りには決まらない.ある密度の値を有限個の点で書き換えても,どの直方体の上での積分も変わらないので,得られる関数はやはり同じ𝑋の密度である.二つの密度が食い違えるのはこの種の,積分を変えない範囲に限られ,したがってℎ(𝑋)の値は密度の選び方によらない.この事実は測度論を使って示されるもので,本書では証明しない.依存するのは,密度を「その一つ」として選んでから値を計算する箇所,すなわち命題 7.3.2 と命題 7.3.3 の証明,補題 7.4.2 を経由する定理 7.4.3 と系 7.4.4,および定義 7.5.4・命題 7.5.5・系 7.5.7 である.第8章と第9章も,微分エントロピーの値を密度から計算するたびに同じ事実に依存する(第8章 命題 8.1.4 と第9章 補題 9.4.2 がその形である).
𝑛 =1が基本で,7.4 節に入るまで𝑛 ≥2の形で述べる主張は補題 7.2.3 だけである.𝑛 ≥2は,複数の観測をまとめて一つの確率変数と見るときに要る.
例
例 7.1.3(一様分布). 𝑤 >0とし,𝑋が長さ𝑤の区間の上の一様分布に従うとする.このときℎ(𝑋)は定まり,ℎ(𝑋) =log𝑤である.とくに𝑤 <1ならばℎ(𝑋) <0であり,𝑤 =1ならばℎ(𝑋) =0である.
証明. 区間を[𝑎,𝑎 +𝑤]とすると,密度は[𝑎,𝑎 +𝑤]の上で𝑓(𝑥) =1/𝑤,その外で0である.𝑓(𝑥)log𝑓(𝑥)は[𝑎,𝑎 +𝑤]の上で定数1𝑤log1𝑤,その外で0だから,その絶対値の積分は有限でℎ(𝑋)は定まる.値は
ℎ(𝑋)=−∫𝑎+𝑤𝑎1𝑤log1𝑤𝑑𝑥=−log1𝑤=log𝑤である.logは引数が1で0,1より小さいところで負だから,後半が従う.◼
離散のエントロピーが負にならなかった(命題 1.1.4)のに対し,微分エントロピーは負になる.区間の幅を半分にするとℎはlog2だけ減り,幅をいくらでも小さくとればℎはいくらでも小さくなる.「不確かさの大きさ」という読み方をそのまま持ち込めないことが,最小の例で見えている.
次の例にはlog𝑒という因子が現れる.logの底を𝑒にとればこれは1だが,本書は底を固定していないので落とせない.本章に現れる等式・不等式は,底を1より大きくとるかぎり,そのままの形で成り立つ.底を𝑏から𝑏′に取り替えると,どの対数も一斉にlog𝑏′𝑏倍になるので,ℎの値もlog𝑒も同じ倍率で変わるからである.
例 7.1.4(指数分布). 𝑚 >0とし,𝑋が平均𝑚の指数分布,すなわち密度
𝑓(𝑥)=⎧{
{⎨{
{⎩1𝑚𝑒−𝑥/𝑚(𝑥≥0),0(𝑥<0)をもつとする.このときℎ(𝑋)は定まり,ℎ(𝑋) =log(𝑒 𝑚)である.
証明. 𝑥 ≥0でlog𝑓(𝑥) = −log𝑚 −𝑥𝑚log𝑒だから
𝑓(𝑥)log𝑓(𝑥)=−1𝑚𝑒−𝑥/𝑚(log𝑚+𝑥𝑚log𝑒)である.右辺の絶対値は𝑥 →∞で指数的に0に近づくので積分は絶対収束し,ℎ(𝑋)は定まる.値は
ℎ(𝑋)=log𝑚∫∞01𝑚𝑒−𝑥/𝑚𝑑𝑥+log𝑒𝑚∫∞0𝑥1𝑚𝑒−𝑥/𝑚𝑑𝑥となる.第 1 の積分は密度の全積分だから1である.第 2 の積分は部分積分で
∫∞0𝑥1𝑚𝑒−𝑥/𝑚𝑑𝑥=[−𝑥𝑒−𝑥/𝑚]∞0+∫∞0𝑒−𝑥/𝑚𝑑𝑥=𝑚となる.よってℎ(𝑋) =log𝑚 +log𝑒 =log(𝑒 𝑚)である.◼
幅で読む. 例 7.1.3 と例 7.1.4 では,どちらも微分エントロピーが長さの対数の形になった.一様分布では区間の長さそのものが,指数分布では平均の𝑒倍がその長さである.底を2にとると,この形は2ℎ(𝑋)をその分布の実効的な幅と読む物差しになる.幅𝑤の一様分布では2ℎ(𝑋) =𝑤で区間の幅そのもの,平均𝑚の指数分布では2ℎ(𝑋) =𝑒 𝑚である.長さの対数だから,長さの単位を変えれば値も変わる.引き伸ばして幅を何倍かにすればℎはその倍率の対数だけ増えるはずで,実際そうなることは 7.3 節で示す.
最後に,ℎが定まらない例を見ておく.定義 7.1.2 が絶対値をつけた積分の有限性を要求したのは,こういう密度を外すためである.
例 7.1.5(裾の重い密度). 𝑐 :=(log𝑒)2とおき,𝑥 ≥𝑒で𝑓(𝑥) :=𝑐𝑥(log𝑥)2,𝑥 <𝑒で𝑓(𝑥) :=0と定める.このとき𝑓はℝ上の密度であり,𝑓を密度にもつ実数値確率変数𝑋について∫ℝ|𝑓(𝑥)log𝑓(𝑥)| 𝑑𝑥 =∞である.すなわちℎ(𝑋)は定まらない.
証明. 𝑓が密度であることを見る.𝑓は非負である.log𝑡の導関数が(log𝑒)/𝑡であることから,−1/log𝑥の導関数は(log𝑒)/(𝑥(log𝑥)2)である.𝑥 →∞で−1/log𝑥 →0であり,𝑥 =𝑒での値は−1/log𝑒だから
∫∞𝑒𝑑𝑥𝑥(log𝑥)2=1log𝑒[−1log𝑥]∞𝑒=1(log𝑒)2である.よって∫ℝ𝑓 =𝑐/(log𝑒)2 =1であり,𝑓はℝ上の密度である.
絶対値をつけた積分が発散することを見る.𝑥 ≥𝑒では−log𝑓(𝑥) = −log𝑐 +log𝑥 +2loglog𝑥であり,𝑥 →∞でlog𝑥もloglog𝑥も+∞に発散するから,𝑥が十分大きいところでは−log𝑓(𝑥) ≥12log𝑥である.そこでは|𝑓(𝑥)log𝑓(𝑥)| =𝑓(𝑥)( −log𝑓(𝑥)) ≥𝑐/(2𝑥log𝑥)が成り立つ.ところがloglog𝑥の導関数が(log𝑒)2/(𝑥log𝑥)であることから∫∞𝑑𝑥𝑥log𝑥は発散する.よって∫ℝ|𝑓log𝑓| =∞である.◼
離散のエントロピーとの関係
微分エントロピーが離散のエントロピーの何にあたるのかは,𝑋を有限の精度で読み取ったときのエントロピーを調べればわかる.実数を幅𝛿の区間に切って,どの区間に入ったかだけを記録すれば,得られるのは有限個の値をとる確率変数であり,そのエントロピーは定義 1.1.1 で測れる.精度を上げる(𝛿を小さくする)と何が起きるかが,次の定理である.
証明が使う一般の道具を二つ借りる.一つめは 積分の平均値の定理,すなわち「有界閉区間の上で連続な実数値関数について,その区間の上での積分は,区間の中のある1点での関数の値に区間の長さを掛けたものに等しい」という微積分の定理である.二つめは リーマン和の収束,すなわち「有界閉区間の上で連続な実数値関数について,区間を𝐽等分して各小区間から1点ずつ選んで作ったリーマン和は,𝐽 →∞でその関数の積分に収束する」という事実である.どちらも当てる相手は有界閉区間の上の連続な実数値関数だけで,依存するのは定理 7.1.6 の証明だけである.ほかに,𝜑(𝑡) = −𝑡log𝑡が[0,∞)の全体で連続であること(補題 1.1.8)を使う.
次の定理で区間に付ける名前𝐼𝑖は,定義 1.3.1 の相互情報量𝐼とは関係がない.
定理 7.1.6(量子化した確率変数のエントロピー). 𝑎 <𝑏を実数とし,𝑋を密度𝑓をもつ実数値確率変数で,𝑓は[𝑎,𝑏]の外で0であり[𝑎,𝑏]の上で連続であるとする.各整数𝐽 ≥1について𝛿𝐽 :=(𝑏 −𝑎)/𝐽とおき,[𝑎,𝑏]を
𝐼𝑖=[𝑎+(𝑖−1)𝛿𝐽,𝑎+𝑖𝛿𝐽)(1≤𝑖≤𝐽−1),𝐼𝐽=[𝑎+(𝐽−1)𝛿𝐽,𝑏]と分ける.̃𝑋𝐽 ∈{1,…,𝐽}を,𝑋 ∈𝐼𝑖のときその番号𝑖,𝑋 ∉[𝑎,𝑏]のとき1と定める.このときℎ(𝑋)は定まり,
𝐻(̃𝑋𝐽)+log𝛿𝐽⟶ℎ(𝑋)(𝐽→∞)である.
仮定は,密度が有界閉区間[𝑎,𝑏]の外で0になることを求めている.例 7.1.4 の指数分布も,7.2 節で扱うガウス分布もこれを満たさないので,定理 7.1.6 はそのままでは当たらない.この仮定を外したときに同じ極限が成り立つかどうかは,本書では扱わない.
証明. はじめにℎ(𝑋)が定まることを見る.第1章の記号にならって𝜑(𝑡) := −𝑡log𝑡と書く.第1章が認めたとおり𝜑は[0,∞)の上で連続だから𝜑 ∘𝑓は[𝑎,𝑏]の上で連続で,有界閉区間の上の連続関数はその絶対値も積分できる.𝑓は[𝑎,𝑏]の外で0で,そこでは𝜑(𝑓(𝑥)) =𝜑(0) =0だから,∫ℝ|𝑓log𝑓| =∫𝑏𝑎|𝜑(𝑓(𝑥))| 𝑑𝑥 <∞であり,ℎ(𝑋) =∫𝑏𝑎𝜑(𝑓(𝑥)) 𝑑𝑥である.
̃𝑋𝐽の分布を求める.定義 7.1.1 を[𝑎,𝑏]に当てるとPr[𝑋 ∈[𝑎,𝑏]] =∫𝑏𝑎𝑓 =1だから,𝑋 ∉[𝑎,𝑏]は確率0でしか起きず,̃𝑋𝐽 =𝑖の確率は𝑝𝑖 :=Pr[𝑋 ∈𝐼𝑖]に等しい(𝑖 =1については確率0の場合が足されるだけである).次に𝑝𝑖を,密度の1点での値で書き直す.1点からなる閉区間に定義 7.1.1 を当てると,𝑋がその1点をとる確率は0である.したがって𝑝𝑖は,𝐼𝑖に右端を加えた閉区間――𝐼𝑖に𝑋が入る確率に等しく,定義 7.1.1 よりそれは――𝐼𝑖の上での𝑓の積分である.𝑓は――𝐼𝑖の上で連続だから,積分の平均値の定理により,ある𝑡𝑖 ∈――𝐼𝑖について𝑝𝑖 =𝑓(𝑡𝑖) 𝛿𝐽が成り立つ.
エントロピーを二つの部分に分ける.𝑠 >0と𝑟 ≥0について𝜑(𝑠𝑟) =𝑠 𝜑(𝑟) −𝑠 𝑟log𝑠である(𝑟 >0では両辺とも−𝑠𝑟log(𝑠𝑟)を展開したもの,𝑟 =0では両辺0).𝑠 =𝛿𝐽,𝑟 =𝑓(𝑡𝑖)ととると𝜑(𝑝𝑖) =𝛿𝐽 𝜑(𝑓(𝑡𝑖)) −𝑓(𝑡𝑖) 𝛿𝐽log𝛿𝐽である.̃𝑋𝐽は{1,…,𝐽}に値をとりPr[̃𝑋𝐽 =𝑖] =𝑝𝑖だから,定義 1.1.1 より
𝐻(̃𝑋𝐽)=𝐽∑𝑖=1𝜑(𝑝𝑖)=𝐽∑𝑖=1𝛿𝐽𝜑(𝑓(𝑡𝑖))−log𝛿𝐽𝐽∑𝑖=1𝑓(𝑡𝑖)𝛿𝐽である.∑𝑖𝑓(𝑡𝑖)𝛿𝐽 =∑𝑖𝑝𝑖 =1だから,移項して
𝐻(̃𝑋𝐽)+log𝛿𝐽=𝐽∑𝑖=1𝛿𝐽𝜑(𝑓(𝑡𝑖))を得る.
右辺の極限をとる.右辺は,[𝑎,𝑏]を𝐽等分して各小区間から点𝑡𝑖を選んで作った𝜑 ∘𝑓のリーマン和にほかならない.𝜑 ∘𝑓は[𝑎,𝑏]の上で連続だから,リーマン和の収束により,𝐽 →∞でこの和は∫𝑏𝑎𝜑(𝑓(𝑥)) 𝑑𝑥に収束する.はじめに見たとおりこれはℎ(𝑋)である.◼
定理 7.1.6 が,微分エントロピーと離散のエントロピーの関係を与えている.𝐽を大きくすると𝛿𝐽 →0なのでlog𝛿𝐽 → −∞であり,したがって𝐻(̃𝑋𝐽)自身は+∞に発散する.実数を無限の精度で読み取れば無限の情報が要る,ということである.微分エントロピーはその発散する部分−log𝛿𝐽を差し引いた残りで,精度を上げるたびに機械的に増えるぶんを除いた「散らばりの広さ」だけを測っている.
桁数で読む. 区間の両端𝑎と𝑏がどちらも2進の有限小数である場合を,底を2にとって読んでみる.整数𝑘を十分大きくとれば𝛿𝐽 =2−𝑘となる𝐽があり,そのとき𝐼𝑖による分割は,𝑋を小数点以下𝑘桁(2進)まで切り捨てて読む操作と一致する.定理 7.1.6 をこの𝐽に沿って読むと𝐻(̃𝑋𝐽) −𝑘 →ℎ(𝑋)(𝑘 →∞)である.𝑘桁まで読み取った値のエントロピーが,𝑘ビットよりℎ(𝑋)だけ多い値に近づいていく,ということである.各𝑘で等しくなると言っているのではない.ℎ(𝑋)が負でありうる(例 7.1.3)のは,𝑘桁の読み取りが𝑘ビットに満たない情報しか運ばない場合があるからである.𝑋が[0,1/2]の上の一様分布に従うときがそれで,例 7.1.3 よりℎ(𝑋) = −1ビットであり,𝑘桁まで読み取った値は2𝑘−1個の値を等しい確率でとるから,そのエントロピーは例 1.1.3 より𝑘 −1ビットである.
この読み方からは,微分エントロピーを単独で「情報量」と読めない理由も見える.𝐻(̃𝑋𝐽)自身は発散し,ℎ(𝑋)は例 7.1.3 のように負にもなる.ℎ(𝑋)が働くのは,−log𝛿𝐽という共通の項を差し引いたあとの比較としてである.同じ区間[𝑎,𝑏]の上で連続でその外で密度が0である二つの確率変数に定理 7.1.6 を当てて引き算すると,−log𝛿𝐽が打ち消し合い,読み取った値のエントロピーの差が微分エントロピーの差に収束する.7.5 節で相対エントロピーと相互情報量を定義すると,そこには差だけが現れる.離散のときと同じ形の式が,そのまま成り立つことになる.
InformationTheory — 形式化検証つき情報理論教科書(レビュー版).数式は MathJax + AMS Euler で事前レンダリング.