10.4 Lagrange 乗数の存在

10.2 節10.3 節の主張は,どれも「そのの Gibbs 分布がモーメント制約を満たすならば」という条件のもとで述べられている.順序を確かめておきたい.先にあるのは制約の値のほうで,はそれに合うように後から探すものである.探せば見つかるのか.見つかるとして,1 つに決まるのか.この問いに答えないかぎり,最大エントロピー定理は「もし解があれば」という条件つきの主張にとどまる.実行可能な分布はいくらでもあるのに Gibbs 分布がその中に 1 つも入っていないという制約があるかどうかが,分からないままだからである.

これまでにが実際に見つかった例は 2 つある.例 10.2.8 の二値の場合は式で解け,例 10.1.4 のサイコロは数値で解けた.だがそれは 2 つの例でそうだったというだけで,一般に何が起きるのかは何も言っていない.

本節の道具は微分である.10.3 節が主役に据えた対数分配関数を 1 回微分すると特徴関数の平均が出て(命題 10.4.1),もう 1 回微分すると分散が出る(命題 10.4.2).特徴関数が 1 つのときは,この二つだけで存在と一意性がともに出る(定理 10.4.3).以下,は空でない有限アルファベット,はその上の特徴関数,10.2 節で置いた記法である.

分配関数を微分する

微分の計算規則を借りる. 借りるのは微積分の計算規則で,次の 5 つである.有限個の微分可能な関数の和・積・商(分母が 0 でないところ)・合成はふたたび微分可能で,導関数は和の法則・積の法則・商の法則・合成関数の微分の法則で与えられること.指数関数と対数関数の導関数がであること.1 次式の導関数がであること.2 階の導関数についても同じ規則を繰り返し使えること.そして,区間の上で 2 階導関数がつねに 0 以上なら,関数はその区間で凸であること.最後の 1 つは,1.1 節が狭義凹の判定に借りたのと同じ微積分の判定法であり,符号を逆にした形で使う.当てる相手は,有限アルファベット上の和として書かれたとその商,および直線に沿ってを1 変数関数に制限したものである.この借用に依存するのは命題 10.4.1命題 10.4.2 だけで,以降はこの 2 つの結論だけを使う.本書はこの 5 つを証明しないが,形式化されていないわけではない.実際,どれも Mathlib にある無条件の機械検証済みの定理として形式化されている.

形式化上の注記(本節共通). 本節の主張(命題 10.4.1命題 10.4.2定理 10.4.3系 10.4.4例 10.4.5例 10.4.6)には,対応する形式化された宣言がない.本書はこれらを証明するが,そのままの形では機械検証されていない.本節に形式化ポインタが 1 つも現れないのはそのためである.なお命題 10.4.1の内容は,有限アルファベットではなく一般の確率測度と傾き変換の言葉で,特徴関数が1 つの場合に限って形式化されている.本節の記号との橋渡しがないので,ここでは紐付けない.

命題 10.4.1. を空でない有限アルファベット,をその上の特徴関数とする.定義 10.3.1 の対数分配関数は各変数について偏微分可能で,任意のに対して

𝜕𝜓𝜕𝜆𝑖(𝜆)=𝔼𝑝𝜆[𝑓𝑖]=𝑥X𝑝𝜆(𝑥)𝑓𝑖(𝑥)

が成り立つ.

証明. を固定し,第成分だけを動かす.各についての 1 次式であり,そのによる導関数はである.合成関数の微分とより

𝜕𝜕𝜆𝑖exp(𝜆,𝑓(𝑥))=𝑓𝑖(𝑥)exp(𝜆,𝑓(𝑥))

となる.は有限だからはこれらの有限和であり,項別に微分できて

𝜕𝑍𝜕𝜆𝑖(𝜆)=𝑥𝑓𝑖(𝑥)exp(𝜆,𝑓(𝑥))

である.指数関数の値は正では空でないからであり,に合成関数の微分とを当てられて

𝜕𝜓𝜕𝜆𝑖(𝜆)=1𝑍(𝜆)𝑥𝑓𝑖(𝑥)exp(𝜆,𝑓(𝑥))=𝑥𝑓𝑖(𝑥)𝑝𝜆(𝑥)

を得る.最後の等号は定義 10.2.1 による.

微分するとモーメントが出る. 分配関数を微分すると特徴関数の平均が出てくる,というこの形が本節を動かす.という 1 つの関数が,ごとの Gibbs 分布のもとでの平均を全部覚えていることになる.が母関数(微分するたびに分布の量を吐き出す関数)として使われるのは,この性質による.

を探す問題も,これで書き換わる.モーメント制約は,命題 10.4.1 により

𝜕𝜓𝜕𝜆𝑖(𝜆)=𝑐𝑖(𝑖=1,,𝑘)

と同じことである.すなわちを探すとは,の勾配が指定された値になる点を探すことにほかならない.勾配が指定の値をとる点を探す問題では,その関数が凸かどうかが効いてくる.次の命題がそれを与える.

命題 10.4.2. を空でない有限アルファベット,をその上の特徴関数とし,に対して上の関数で定める.任意のに対し𝑡 )とおくと,は 2 回微分可能で,すべてのに対して

𝑔(𝑡)=𝔼𝑝𝜆+𝑡𝑣[(𝑣𝑓)2](𝔼𝑝𝜆+𝑡𝑣[𝑣𝑓])2

が成り立つ.右辺はのもとでのの分散であり,つねに0 以上である.したがって上の凸関数である.さらに,となるのは上で定数のとき,かつそのときに限る.とくにととるとの平行移動になり,は 2 回微分可能でである.

証明. 𝑎(𝑥) :=𝜆,𝑓(𝑥)とおくとである.そこで

𝑆(𝑡):=𝑍(𝜆+𝑡𝑣)=𝑥exp(𝑎(𝑥)+𝑡𝑏(𝑥))

と書く.各項はの 1 次式と指数関数の合成だから,合成関数の微分によりについて何度でも微分でき,が有限だから項別に微分できて

𝑆(𝑡)=𝑥𝑏(𝑥)exp(𝑎(𝑥)+𝑡𝑏(𝑥)),𝑆(𝑡)=𝑥𝑏(𝑥)2exp(𝑎(𝑥)+𝑡𝑏(𝑥))

である.指数関数の値は正では空でないからであり,に合成関数の微分と商の法則を当てて

𝑔(𝑡)=𝑆(𝑡)𝑆(𝑡),𝑔(𝑡)=𝑆(𝑡)𝑆(𝑡)𝑆(𝑡)2𝑆(𝑡)2=𝑆(𝑡)𝑆(𝑡)(𝑆(𝑡)𝑆(𝑡))2

を得る.定義 10.2.1 よりだから,𝑆(𝑡)/𝑆(𝑡) =𝔼𝑝𝜆+𝑡𝑣[𝑏]であり,主張の等式が出る.

右辺が 0 以上であることを見る.𝑤 :=𝑝𝜆+𝑡𝑣とおくと,を使って展開すれば

𝔼𝑤[𝑏2]¯𝑏2=𝑥𝑤(𝑥)(𝑏(𝑥)¯𝑏)2

である.各項は非負だから和も非負である.しかも和が 0 になるのは全項が 0 のとき,かつそのときに限り,命題 10.2.2 よりだから,それは各であること,すなわち上で定数であることと同値である.

凸性を見る.上の関数が凸であるとは,任意の 2 点を結ぶ線分の上で凸だということであり,線分はを選んでと書けるから,すべてのについてが凸であることと同じである.いまが示せたので,借用した「2 階導関数が 0 以上なら凸」によりは凸であり,したがっても凸である.

最後にの場合を見る.このときだから,逆にと書ける.は 1 次式でその導関数は1 だから,借用した合成関数の微分と 1 次式の導関数により,は 2 回微分可能でである.とおけばを得る.

分散が非負だから凸である. 凸性の理由がこれ以上ないほど短い形で出た.2 階微分は分散で,分散は 2 乗の平均だから負になりようがない.逆に,凸性が崩れかける(2 階微分が 0 になる)のは分散が消えるときで,それは上で定数になる向きに限られる.特徴関数どうしが一次従属だったり,特徴関数の一次結合が定数関数になったりすると,そういう向きが現れる.定数関数を加えたが一次独立なら,はどの向きにも狭義に曲がっている.

のときにこれをどう使うかを言っておく.特徴関数が 1 つで,それが定数関数でなければ,上の狭義性によりはどの点でも正である.すると,命題 10.4.1 によれば特徴関数の平均であるは狭義単調増加になる.勾配が指定の値をちょうど 1 回とる,という次の定理の結論は,凸性をこの単調性の形で使って出す.

特徴関数が 1 つのときの存在と一意性

中間値の定理と,導関数から単調性を読む規則を借りる. 借りるのは 4 つである.一つめは中間値の定理で,形は「有界閉区間の上の実数値連続関数は,両端の値のあいだの値をすべてとる」.二つめは「微分可能な関数は連続である」.三つめは「区間の上で導関数がつねに正なら,関数はその区間で狭義単調増加である」.四つめは指数関数の基本性質で,が狭義単調増加であること,であること,および有限個の項からなる和と商の極限が項ごとの極限から定まることである.当てる相手はいずれもの全体で定義された 1 変数の実数値関数で,具体的にはとその有界閉区間への制限,記号を固定したときの,および指数関数の値の有限個からなる和と商である.一つめから三つめに依存するのは定理 10.4.3 だけ,四つめに依存するのは定理 10.4.3例 10.4.5例 10.4.6 だけで,ほかの主張はこれらの結論しか使わない.本書はこの 4 つを証明しないが,形式化されていないわけではない.実際,どれもMathlib にある無条件の機械検証済みの定理として形式化されている.

定理 10.4.3(制約に整合する Lagrange 乗数の存在と一意性). を空でない有限アルファベット,とし,をその上の特徴関数とする.

𝑓min:=min𝑥X𝑓1(𝑥),𝑓max:=max𝑥X𝑓1(𝑥)

とおく.制約の値を満たすならば,Gibbs 分布がモーメント制約を満たすような実数が,ただ 1 つ存在する.

証明. とおく.示すことは,が値をちょうど 1 回とることである.

は定数関数でない. 仮定よりである.が有限だから最小値と最大値はどちらもの値で,それが異なるということはが定数関数でないということである.

は狭義単調増加で連続である. なので偏微分は普通の微分であり,命題 10.4.1 よりである.命題 10.4.2にとって使うと,であり,は 2 回微分可能でのもとでのの分散に等しい.すなわち

𝑚(𝜆)=𝜓(𝜆)=𝔼𝑝𝜆[𝑓21](𝔼𝑝𝜆[𝑓1])2

である.は定数関数でないから,命題 10.4.2 の狭義性の主張によりがすべてので成り立つ.借用した「導関数がつねに正なら狭義単調増加」によりの上で狭義単調増加であり,微分可能だから連続でもある.

両端の極限を計算する. とおく.は空でない有限集合だからも空でなく,である.定義 10.2.1 によればで割ったものである.この分子と分母をで割ると

𝑚(𝜆)=𝑥𝑓1(𝑥)exp(𝜆(𝑓1(𝑥)𝑓max))𝑥exp(𝜆(𝑓1(𝑥)𝑓max))

である.の項では指数の中身が 0 なのでの値は 1,の項ではだからで指数の中身はに向かい,の値は 0 に向かう.項は有限個だから,分子はに,分母はに向かう.よって

lim𝜆+𝑚(𝜆)=𝑓max

である.を与える点の集合について同じ計算をすればを得る.

中間値の定理を当てる. だから,での極限より,となる実数がある.同じくより,となる実数がある.は狭義単調増加だからである(ならとなって矛盾する).は閉区間の上で連続だから,借用した中間値の定理によりを満たすがこの区間にある.これで存在が言えた.

一意性は狭義単調性から出る.ならばだから,両方がに等しいことはない.

系 10.4.4. を空でない有限アルファベット,とし,をその上の特徴関数とする.

𝑓min:=min𝑥X𝑓1(𝑥),𝑓max:=max𝑥X𝑓1(𝑥)

とおく.制約の値を満たすとし,定理 10.4.3 が与えるただ 1 つの実数をとする.このときは実行可能集合の上でエントロピーを最大にする唯一の分布であり,その最大値は

𝐻(𝑝𝜆)=𝜓(𝜆)𝜆𝑐1

である.

証明. 定理 10.4.3 により,はモーメント制約を満たす.とくには実行可能である.この定理 10.2.5 を当てると,実行可能な任意の分布についてであり,は最大にする分布である.同じ定理 10.2.6 を当てると,実行可能なを満たすのはのとき,かつそのときに限るから,最大にする分布はこれ以外にない.値については,同じ定理 10.3.3 を当ててを得る.ではである.

形式化上の注記. 系 10.4.4 に対応する単独の宣言はない.定理 10.2.5定理 10.2.6定理 10.3.3 の形式化を合わせれば,を与えたときの最大性・一意性・値までは出る.そのの存在と一意性を与える定理 10.4.3 の側が形式化されていないので,系全体としては機械検証に届かない.

条件つきだった主張が,条件から解放される. 10.2 節10.3 節の主張はどれも「そのの Gibbs 分布が制約を満たすならば」という条件つきだった.系 10.4.4 はその条件を,制約の値がの最小値と最大値のあいだに真にあるという,制約の側だけで確かめられる条件に置き換えている.を先に見つけておかなくても,最大化分布が存在すること,それが Gibbs 分布であること,1 つしかないこと,そしてその値までが言える.

例 10.1.4 のサイコロがまさにこの場合である.についてだから仮定はぴったり満たされ,10.2 節で数値的に求めたは確かに存在して,しかもそれ以外にない.あの数値解は当て推量ではなく,系 10.4.4 の保証する唯一解だった.

には踏み込まない. 特徴関数が 2 つ以上あるときに同じことが言えるかは,本書では扱わない.理由は道具立てにある.で効いたのは,が実数上の 1 変数関数であって,狭義単調増加なら中間値の定理がそのまま当たるという点だった.では,個の偏微分を並べたの勾配はからへの写像になり,「両端の値のあいだ」という言い方がそもそもできない.が存在するの範囲は,点𝑥 X)の凸包の言葉で述べることになり,その特徴づけと証明には多変数の凸解析の道具が要る.本書はその道具を借りていないので,この場合は証明しない.形式化もされていない.

例で確かめる

例 10.4.5(二値の場合の乗数). X ={0,1}𝑘 =1とする.このとき任意の実数に対して

𝔼𝑝𝜆[𝑓1]=𝑒𝜆1+𝑒𝜆=11+𝑒𝜆

であり,この値はについて狭義単調増加で,で 0 に,で 1 に向かう.制約の値をとすると,定理 10.4.3 が与える唯一のである.

証明. 𝑓1(0) =0だから,定義 10.2.1 よりであり,

𝔼𝑝𝜆[𝑓1]=0𝑝𝜆(0)+1𝑝𝜆(1)=𝑒𝜆1+𝑒𝜆

である.分子と分母をで割れば 2 つめの表示になる.借用した指数関数の性質によりについて狭義単調減少で正だから,は狭義単調増加である.またではなので値は 0 に向かい,ではなので値は 1 に向かう.

𝑓min =0だから,に対して定理 10.4.3 の仮定が満たされ,制約を満たすがただ 1 つ存在する.いっぽう例 10.2.8がその制約を満たすことを示している.一意性により,これが定理 10.4.3 の与えるである.

定理 10.4.3 の道具立てが,そのまま目に見える. 例 10.4.5 の式は,定理 10.4.3 の証明に出てきたを最小の場合に書き下したものである.狭義単調増加であること,両端の極限がであること,そして値がその 2 つには決して届かないこと,すなわちが 0 にも 1 にもならないことが式のまま読める.制約の値が開区間の中にあるという仮定が何を排除しているのかも,これで見当がつく.端点を指定したらどうなるのかが,次の例である.

例 10.4.6(制約の値が端にあるとき). を空でない有限アルファベット,とし,特徴関数は定数関数でないとする.とおき,とおいて,制約の値をにとる.このとき次の 4 つが成り立つ.

  1. 上の分布が実行可能であることと,の外のすべての点でが 0 をとることは同値である.
  2. どの実数についても,は実行可能でない.
  3. 実行可能集合の上でエントロピーを最大にする分布は上の一様分布ただ 1 つで,その値はである.
  4. のとき,に,で 0 に収束する.

証明.

  1. すべてのであり,だから

    𝑓max𝔼𝑃[𝑓1]=𝑥𝑃(𝑥)(𝑓max𝑓1(𝑥))

    であって,右辺の各項は非負である.が実行可能であること,すなわち左辺が 0 であることは,全項が 0 であることと同値である.全項が 0 であるとは各またはが成り立つことであり,これは「ならば」と同じことである.

  2. 命題 10.2.2 よりはすべての点で正である.は定数関数でないからに属さない点があり,そこでである.主張 1 よりは実行可能でない.

  3. 主張 1 より,実行可能な分布の外で 0 をとる.よってであり,に制限したものは上の分布である.の項はエントロピーの和に寄与しない(定義 1.1.1 の約束)から,はその制限のエントロピーに等しい.を空でない有限アルファベットとみて定理 1.1.5 を当てるとであり,等号は上で一様なとき,かつそのときに限る.上の一様分布はの外で 0 をとるので,主張 1 より実行可能である.

  4. 定理 10.4.3 の証明と同じく,定義 10.2.1 の分子と分母をで割ると

    𝑝𝜆(𝑥)=exp(𝜆(𝑓1(𝑥)𝑓max))𝑦exp(𝜆(𝑓1(𝑦)𝑓max))

    である.のとき,分母はに向かい,分子はならつねに 1,なら 0 に向かう.

Gibbs の形は最大化分布を全部は覆わない. 例 10.4.6 が言っているのは,最大化分布は存在するのに,それが Gibbs 分布ではない,ということである.10.2 節の定理が言ったのは「制約を満たすがあれば,それが最大化分布」ということであって,その逆向きの「最大化分布はかならず Gibbs 分布である」ということではなかった.その差が現れるのがここである.最大化分布は Gibbs 分布の極限としては届く.届かないのは,が有限の値であるかぎり,ということである.

定理 10.4.3 が開区間を仮定した理由も,これで見える.制約の値が端に来ると実行可能集合が痩せて,全点で正の分布を 1 つも含まなくなる.はつねに全点で正だから(命題 10.2.2),その中に入りようがない.開区間という仮定は,実行可能集合が全点で正の分布を含むだけの余裕をもつことを保証している.

第10章が示したことを並べておく.制約が分布について線形であるかぎり(定義 10.1.1),最大エントロピー分布の形はに比例するものに絞られる(定義 10.2.1定理 10.2.5).その根拠は相対エントロピーの非負性ただ 1 つで(定理 1.6.1),凸性も微分も要らなかった.最大値そのものは対数分配関数を使ってという閉じた式で書ける(定理 10.3.3).最後に残った問いは,制約に整合するがあるのかどうかだった.特徴関数が 1 つのときは,を 2 回微分するだけで答えが出て(定理 10.4.3),条件つきだった主張が制約の値だけで確かめられる形になった(系 10.4.4).「知っていることは反映し,知らないことは決めつけない」という 10.1 節の原則を 1 つ置いただけで,指数型分布族という具体的な分布の形と,その形を決めるパラメータの求め方までが出てくる.

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