7.2 ガウス分布の微分エントロピー

例 7.1.3例 7.1.4 は,どちらも分布を先に決めて微分エントロピーを計算するものだった.本節が扱うのはその逆向きの問いである.散らばりの大きさをある値に固定したとき,微分エントロピーをいちばん大きくする分布は何か.離散のときの答えは一様分布だった(定理 1.1.5).そこで散らばりを抑えていたのはアルファベットの大きさだが,実軸の上ではその役目を果たすものがない.実軸の全体で一定の値をとる関数は,積分がにならないので定義 7.1.1 の密度になりえないからである.代わりに分散を固定すると,ガウス分布がその最大を与える.何を固定するかは一通りではない.本節が分散を選ぶのは,ガウス通信路の章で扱う,送る信号の電力に制約のかかる通信路のためである.そこで置く制約は送る信号の乗の平均を抑えるもので,平均がの信号では,それが分散を抑えることと同じになる.

本節は,ガウス分布の微分エントロピーがであることを示し,続いて分散を固定した分布のうちでこれが最大であることを示す.二つめの証明の要は,比べる相手の密度を一つ持ち出して上から抑える,という形の補題である.この補題は 7.4 節7.5 節でも同じ形で使う.

借りる道具は二つある.一つめは学部の確率で扱う ガウス分布の基本性質,すなわち「定義 7.2.1 の関数は実軸の上での積分がであり,それを密度とする確率変数の平均は,分散はである」ことである.当てる対象は定義 7.2.1 のガウス分布だけで,依存するのは定理 7.2.2 の証明,定理 7.2.4 の証明,系 7.2.5 の証明,7.5 節例 7.5.8 の証明,および第8章命題 8.1.4定理 8.2.2系 8.2.3補題 8.3.3 の証明と,8.2 節定義 8.2.1 の直後で上限をとる範囲が空でないことを確かめる箇所,8.4 節が球の勘定のために雑音の乗和の平均を求める箇所,同じ性質を 9.4 節が確率変数の言葉で借り直している第9章命題 9.4.4命題 9.4.6 の証明,および同じ性質をそのまま引いている第16章系 16.7.6命題 16.7.8 の証明である(第16章が引くのは,実軸の上での積分がであることと,という形で書いた平均と分散で,射程はここから広がらない).二つめは 零集合の上での積分,すなわち「非負な関数と実数値関数について,でありであるならばである」という積分の事実である.当てる対象は,が本章に現れる密度で,が密度から「その密度が正の点でその対数,そのほかの点で」として作った関数(の第成分に当てたものを含む)である場合だけである.依存するのは補題 7.2.3 の証明,7.4 節補題 7.4.2 の証明,および 7.5 節定理 7.5.2 の証明である.

ガウス分布の微分エントロピー

定義 7.2.1(ガウス分布). を実数,とする.密度

𝑔𝜇,𝜎2(𝑥)=12𝜋𝜎2exp((𝑥𝜇)22𝜎2)(𝑥)

をもつ実数値確率変数の分布を,平均,分散ガウス分布(正規分布)といい,と書く.

定理 7.2.2(ガウス分布の微分エントロピー). を実数,とし,に従うとする.このときは定まり

(𝑋)=12log(2𝜋𝑒𝜎2)

である.

証明. 定義 7.2.1 の密度の対数をとると

log𝑔𝜇,𝜎2(𝑥)=12log(2𝜋𝜎2)+(𝑥𝜇)22𝜎2log𝑒

である.右辺は次式であり,借りたガウス分布の基本性質よりだから,の積分は有限では定まる.値は,同じ性質から得られるを使って

(𝑋)=𝑔𝜇,𝜎2(𝑥)(log𝑔𝜇,𝜎2(𝑥))𝑑𝑥=12log(2𝜋𝜎2)+log𝑒2𝜎2𝜎2

となる.右辺はである.

形式化: differentialEntropy_gaussianReal,標準ガウス分布の場合 differentialEntropy_gaussianReal_std,「定まる」にあたる積分可能性 integrable_density_log_density_of_gaussian (ソース)

値が平均によらないのが目につく.ガウス分布の形はで左右に動くだけで幅が変わらないから,散らばりの広さを測る量は動かない.平行移動で微分エントロピーが変わらないことは,ガウス分布に限らず成り立つ(命題 7.3.2).一方には依存していて,倍するとだけ増える.例 7.1.3 の一様分布で幅を倍したときと同じ増え方である.𝜎2 =1/(2𝜋𝑒),すなわちのときになり,それより分散が小さければである.底をにとるとの微分エントロピーはビットである.7.1 節の幅の読みに乗せればで,標準偏差のおよそ倍がこの分布の実効的な幅にあたる.一様分布の幅,指数分布の平均の倍と,同じ物差しの上に並ぶ.

分散を固定したときの最大性

比べる相手の密度による上界を,先に用意する.相手の密度を一つ選ぶごとにという数が決まり,それがを上から抑える,という形の補題である.どのを選ぶかはこちらの自由なので,計算しやすいを選べば,その値がそのままの上界になる.

補題 7.2.3(比べる相手の密度による上界). とし,に値をとり密度をもつ確率変数,上の密度とする.とおく.であり,が定まり,であるならば

(𝑋)𝑆𝑓(𝑥)log𝑞(𝑥)𝑑𝑥

である.

証明. の各点でかつだから,補題 1.1.7 の対数不等式に当てられて

𝑓(𝑥)log𝑞(𝑥)𝑓(𝑥)𝑓(𝑥)(𝑞(𝑥)𝑓(𝑥)1)log𝑒=(𝑞(𝑥)𝑓(𝑥))log𝑒

が成り立つ.

両辺をの上で積分する.左辺の被積分関数はと分かれ,前者の絶対値の積分は仮定より有限,後者の絶対値の積分はが定まることから有限だから,左辺の積分は絶対収束してに等しい.右辺の積分は,が密度だからであり,仮定よりで,なので,以下である.よって

𝑆𝑓(𝑥)log𝑞(𝑥)𝑑𝑥𝑆𝑓(𝑥)log𝑓(𝑥)𝑑𝑥0

である.

左辺の第 2 項をに直す.の外ではのどちらかであり,の点ではである.残る「かつ」の点の全体をとすると,である.借りた零集合の上での積分を,と「の点で,そのほかの点で」という関数に当てるとである.よってであり,上の不等式はとなる.移項すれば主張を得る.

形式化上の注記. の場合の補題 7.2.3 は,differentialEntropy_le_cross_entropy (InformationTheory/Shannon/EPI/InfiniteVariance/Truncation/Construction.lean) が単独で検証している.仮定は本文と一つずつ対応していて,の分布が Lebesgue 測度について絶対連続であることが「が密度をもつ」に,比べる相手の分布についての同じ絶対連続性が「が密度である」に,の分布が比べる相手の分布について絶対連続であることがに,密度にを合成したものが積分できることが「が定まる」に,比べる相手の対数密度がの分布について積分できることがにあたる.一般のに対応する単独の宣言はない.

補題 7.2.3 が言っているのは,仮定を満たすをどう選んでもを下回れない,ということである.に課したのはの二つである.として自身をとればどちらも成り立ち,両辺は等しくなるので,はこの量が届きうる下端にあたる.離散のときに 1.6 節で見た「思い込んで符号化したときの平均符号長」という読み方は,そのままの形では通らない.も負になりうるので(例 7.1.3),どちらも符号長として読めないからである.読み方が戻るのは二つの差のほうで,そこは 7.5 節で扱う.を一般にとってあるのは,7.4 節7.5 節でも同じ補題を使うからである.

仮定は,になる点にが確率を置いていない,ということである.これが外れると不等式は破れる.を長さの区間の上の一様分布の密度,をその区間の左半分の上の一様分布の密度にとると,は左半分でである.このとき例 7.1.3 よりだが,の上で値をとるのでであり,は成り立たない.

定理 7.2.4(ガウス分布の最大エントロピー性). とし,を密度をもつ実数値確率変数とする.が定まり,が平均をもち,その分散が定まって以下であるならば

(𝑋)12log(2𝜋𝑒𝜎2)

である.

証明. の密度を,平均をとし,補題 7.2.3定義 7.2.1)ととって当てる.比べる相手の平均をの平均に合わせるのは,こちらで選べる設定である.

補題 7.2.3 の仮定を確かめる.借りたガウス分布の基本性質より上の密度である.これは実軸のすべての点で正だからであり,である.また定義 7.2.1 の密度の対数をとるとであり,これは次式で,仮定よりは有限だから,7.1 節で借りた期待値の積分表示よりである.

補題 7.2.3 を当てる.の外ではなのでであり

(𝑋)𝑓(𝑥)log𝑔𝜇,𝜎2(𝑥)𝑑𝑥=12log(2𝜋𝜎2)+log𝑒2𝜎2𝔼[(𝑋𝜇)2]

を得る.の分散で,仮定より以下だから,右辺は以下である.

形式化: differentialEntropy_le_gaussian_of_variance_le (ソース)

形式化上の注記. 形式化の仮定は四つある.分布が Lebesgue 測度について絶対連続であることが本文の「密度をもつ」に,密度にを合成したものが積分できることが「が定まる」に,ある中心からの乗が積分できてその積分が以下であることが「分散が定まって以下である」にあたり,四つめはその中心がの平均に等しいことで,これが「が平均をもち」にあたる.本文が借りたガウス分布の基本性質のうち,分散がであることは gaussianReal_variance_eq (InformationTheory/Shannon/DifferentialEntropy.lean) として形式化されている.

系 7.2.5. を実数,とし,に従うとする.このとき定理 7.2.4 の仮定を満たし,その不等式は等号で成り立つ.

証明. 借りたガウス分布の基本性質よりの平均は,分散はであり,定理 7.2.2 よりは定まる.よって定理 7.2.4 の仮定を満たす.その値は定理 7.2.2 よりで,定理 7.2.4 の右辺に等しい.

系 7.2.5 は,定理 7.2.4 の上界が達成されることを言っている.その逆向き,すなわち定理 7.2.4 で等号が成り立つのは分散がちょうどのガウス分布に限る,ということを本書は述べない.仮定が分散以下なので,等号を言うには分散がちょうどであることまで込みで言うことになる.密度どうしの一致は,確率の集合の上での食い違いを許す形でしか言えず,そこを正確に述べるには測度論が要るからである.定理 7.2.4 のこの等号条件は形式化もされていない.

定理 7.2.4 は,離散のときの定理 1.1.5 に対応する主張である.どちらも散らばりの大きさを固定してエントロピーの最大を求めているが,何を固定するかが違う.離散ではとりうる値の個数を固定し,最大を与えたのは一様分布だった.実軸の上では分散を固定し,ガウス分布がその最大を与える.ガウス分布が出てくる理由は,補題 7.2.3 の使い方にある.上界は,が固定した量の次式であれば,その固定した量の平均だけで決まる.定義 7.2.1 の密度はまさにそういうで,次式だから,上界はの分散だけで決まり,分散を以下に固定したどの分布も同じ一つの値で抑えられる.しかも自身がその固定を満たすので,抑えた値がそのまま達成される(系 7.2.5).離散での「いちばん平らな分布が最大にする」という読み方は,ここには持ち込めない.ガウス分布の密度は平らではないし,実軸の上に平らな密度が置けないことが本節の出発点だった.固定するものを変えれば最大を与える分布も変わり,制約を与えたときにどの分布がエントロピーを最大にするかを,有限アルファベットの上で一般に扱うのが第10章である.

例 7.2.6(同じ分散をもつ一様分布との比較). とする.長さの区間の上の一様分布は分散がであり,その微分エントロピーはである.これは,どの実数についてもの微分エントロピーよりだけ小さく,底をにとればこの差は約ビットである.

証明. とし,長さの区間の上の一様分布の分散を求める.区間を平行移動すると平均も同じだけ動くので,分散は変わらない.そこで区間をにとってよく,そのとき平均はで,密度は区間の上でだから,7.1 節で借りた期待値の積分表示より分散は

𝑤/2𝑤/2𝑥21𝑤𝑑𝑥=1𝑤23(𝑤2)3=𝑤212

である.これがに等しいのはのときで,そのとき例 7.1.3 より微分エントロピーはである.

差をとる.定理 7.2.2 より

12log(2𝜋𝑒𝜎2)12log(12𝜎2)=12log2𝜋𝑒12=12log𝜋𝑒6

であり,だからこの値は正である.底をにとるとである.

例 7.2.6定理 7.2.4 の一つの実例である.分散という一つの数だけをそろえて比べると,実軸の全体に裾を引くガウス分布のほうが,有限の区間に押し込めた一様分布より広がっている.差がによらない定数になる理由は次節でわかる.どちらの分布もを掛けて作れるので,掛け算のぶんだけ微分エントロピーが同じだけ増えるからである(命題 7.3.3).

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