tecktom.com

Lévy の反転公式

2026年10月15日26分で読めますAI生成教科書『基礎理論』第9章 特性関数の一意性と反転公式
目次
  1. Dirichlet 積分
  2. 反転公式
  3. 密度の反転公式と原子の公式
  4. 確認テスト
  5. 木下の層

AI 生成・木下が未勉強

確率論入門 第9章 §9.1 ★★ 区間の確率を特性関数の積分で書く公式。端点に原子があれば、その確率の半分が数えられる

第8章では、確率変数 XX の分布 PXP_X に次の特性関数を対応させた。

φX(t)=E[eitX]\varphi_X(t)=E[\mathrm{e}^{\mathrm{i} tX}]

そして独立和の特性関数が積になること、0 のまわりの Taylor 展開がモーメントを与えることを示した。

ところが第8章では、「特性関数が等しければ分布も等しい」という主張を証明せずに使った箇所がいくつかある。

  • 積率母関数の一意性(第8章 §8.2)。
  • 特性関数が実数値なら分布が対称であること(第8章 §8.3)。
  • 独立和の分布をその特性関数から読み取る議論。

この章ではその主張を証明する。道具は、特性関数から分布を具体的に復元する公式である。

  • §9.1 では、区間の確率を特性関数の積分で表す Lévy の反転公式を示す。
  • §9.2 では、反転公式から一意性定理を導く。さらに Rd\mathbb{R}^d 値の確率ベクトルについても一意性を証明して、独立性を同時特性関数の分解で判定する基準を得る。
  • §9.3 では、一意性定理と「独立和の特性関数は積」を組み合わせて、第6章・第7章で畳み込みの計算によって示した再生性を短く証明し直す。
  • §9.4 では、第8章のキュムラントを多変数へ広げる。無相関であっても独立とは限らない2つの確率変数を4次のキュムラントで見分ける方法と、その応用である独立成分分析の原理を述べる。

この章で使う前の章の結果を、式の形で書いておく。第8章 §8.3 の特性関数の定義と基本性質から、

φX(t)=E[eitX]=∫Reitx PX(dx),φX(0)=1,∣φX(t)∣≤1,φaX+b(t)=eibtφX(at)(9.1)\varphi_X(t)=E[\mathrm{e}^{\mathrm{i} tX}]=\int_{\mathbb{R}}\mathrm{e}^{\mathrm{i} tx}\,P_X(\mathrm{d} x),\qquad\varphi_X(0)=1,\qquad\left\vert \varphi_X(t)\right\vert \le1,\qquad\varphi_{aX+b}(t)=\mathrm{e}^{\mathrm{i} bt}\varphi_X(at)\tag{9.1}

を使う。X,YX,Y が独立なら

φX+Y(t)=φX(t) φY(t)(9.2)\varphi_{X+Y}(t)=\varphi_X(t)\,\varphi_Y(t)\tag{9.2}

である(第8章 §8.3)。標準正規分布の特性関数が

∫−∞∞eitx e−x2/22π dx=e−t2/2(t∈R)(9.3)\int_{-\infty}^{\infty}\mathrm{e}^{\mathrm{i} tx}\,\frac{\mathrm{e}^{-x^2/2}}{\sqrt{2\pi}}\,\mathrm{d} x=\mathrm{e}^{-t^2/2}\qquad(t\in\mathbb{R})\tag{9.3}

であることも第8章 §8.3 の例題で示した。第5章からは次の2つを借りる。

  • 優収束定理:次の3つが成り立てば E[Yn]→E[Y]E[Y_n]\to E[Y] である。

Yn→Y (各点),∣Yn∣≤W,E[W]<∞Y_n\to Y\ (\text{各点}),\qquad \left\vert Y_n\right\vert \le W,\qquad E[W]<\infty

  • Fubini の定理(第5章 §5.3):有限測度の積の上で可積分な関数は、積分の順序を交換してよい。

第3章 §3.3 からは「分布関数が等しければ分布も等しい」ことを借りる。原子(P(X=x)>0P(X=x)>0 となる点 xx)が高々可算個であることも借りる。

この章の結果は次のように使われる。一意性定理は第13章 §13.1 の Lévy の連続性定理の証明の中で、極限の分布を特定する段階に使う。多変量の一意性は、第6章 §6.2 の多変量正規分布 Nd(μ,Σ)\mathrm N_d(\bm\mu,\Sigma) が平均と Σ\Sigma だけで決まることを保証する。これは第13章 §13.3 の多変量の中心極限定理(Cramér–Wold の方法)の前提になる。

Dirichlet 積分

Fourier 解析では、関数 ff の Fourier 変換と、そこから ff を復元する逆変換を

f^(t)=∫eitxf(x) dx,f(x)=12π∫e−itxf^(t) dt\hat f(t)=\int\mathrm{e}^{\mathrm{i} tx}f(x)\,\mathrm{d} x,\qquad f(x)=\frac{1}{2\pi}\int\mathrm{e}^{-\mathrm{i} tx}\hat f(t)\,\mathrm{d} t

の組で扱う。特性関数は分布の Fourier 変換なので、同じ式で密度を復元できそうに見える。

しかし分布は密度をもつとは限らず、特性関数は可積分とも限らない(たとえば原子をもつ分布の特性関数は ∣t∣→∞\left\vert t\right\vert \to\infty で 0 に近づかない)。そこで、密度の値の代わりに区間の確率を復元する。その計算の中心にあるのが次の積分である。

Si(T)≔∫0Tsin⁡xx dx(T≥0)(9.4)\mathrm{Si}(T)\coloneqq\int_0^T\frac{\sin x}{x}\,\mathrm{d} x\qquad(T\ge0)\tag{9.4}

被積分関数は x=0x=0 で値 1 に連続に延長できる。よって Si\mathrm{Si} は [0,∞)[0,\infty) 上の連続関数である。Si(−T)=−Si(T)\mathrm{Si}(-T)=-\mathrm{Si}(T) と約束して実数全体に広げる。

補題 9.1(Dirichlet 積分)

lim⁡T→∞Si(T)=π2\displaystyle\lim_{T\to\infty}\mathrm{Si}(T)=\frac{\pi}{2} である。また CSi≔sup⁡T∈R∣Si(T)∣<∞C_{\mathrm{Si}}\coloneqq\sup_{T\in\mathbb{R}}\left\vert \mathrm{Si}(T)\right\vert <\infty である。さらに実数 θ\theta と T>0T>0 について

∫0Tsin⁡θtt dt=Si(θT)(9.5)\int_0^T\frac{\sin\theta t}{t}\,\mathrm{d} t=\mathrm{Si}(\theta T)\tag{9.5}

が成り立つ。

証明

x>0x>0 について 1x=∫0∞e−ux du\frac1x=\int_0^\infty\mathrm{e}^{-ux}\,\mathrm{d} u である。これを代入すると

Si(T)=∫0T ⁣ ⁣∫0∞sin⁡x e−ux du dx\mathrm{Si}(T)=\int_0^T\!\!\int_0^\infty\sin x\,\mathrm{e}^{-ux}\,\mathrm{d} u\,\mathrm{d} x

となる。被積分関数の絶対値の積分は

∫0T ⁣ ⁣∫0∞∣sin⁡x∣e−ux du dx=∫0T∣sin⁡x∣x dx≤T\int_0^T\!\!\int_0^\infty\left\vert \sin x\right\vert \mathrm{e}^{-ux}\,\mathrm{d} u\,\mathrm{d} x=\int_0^T\frac{\left\vert \sin x\right\vert }{x}\,\mathrm{d} x\le T

と有限なので、Fubini の定理により積分の順序を交換できる。部分積分を2回行うと

∫0Te−uxsin⁡x dx=1−e−uT(usin⁡T+cos⁡T)1+u2\int_0^T\mathrm{e}^{-ux}\sin x\,\mathrm{d} x=\frac{1-\mathrm{e}^{-uT}(u\sin T+\cos T)}{1+u^2}

となる。右辺を TT で微分すると e−uTsin⁡T\mathrm{e}^{-uT}\sin T になる。また右辺は T=0T=0 で 0 になる。この2点から確かめられる。したがって

Si(T)=∫0∞du1+u2−∫0∞e−uT(usin⁡T+cos⁡T)1+u2 du(9.6)\mathrm{Si}(T)=\int_0^\infty\frac{\mathrm{d} u}{1+u^2}-\int_0^\infty\frac{\mathrm{e}^{-uT}(u\sin T+\cos T)}{1+u^2}\,\mathrm{d} u\tag{9.6}

である。第1項は π2\frac\pi2 である。T≥1T\ge1 では、第2項の被積分関数の絶対値は次のように抑えられる。

e−uT∣usin⁡T+cos⁡T∣1+u2≤e−u(1+u)1+u2≤11+u2\frac{\mathrm{e}^{-uT}\left\vert u\sin T+\cos T\right\vert }{1+u^2}\le\frac{\mathrm{e}^{-u}(1+u)}{1+u^2}\le\frac{1}{1+u^2}

2つ目の不等号は eu≥1+u\mathrm{e}^u\ge1+u による。また被積分関数は各 u>0u>0 で、T→∞T\to\infty のとき 0 に収束する。優収束定理から第2項は 0 に収束し、Si(T)→π2\mathrm{Si}(T)\to\frac\pi2 を得る。

次に有界性を示す。Si\mathrm{Si} は連続で、両側の無限遠で有限の極限 ±π2\pm\frac\pi2 をもつ。

  • 十分大きな RR をとれば、∣T∣≥R\left\vert T\right\vert \ge R の範囲で次が成り立つ。

∣Si(T)∣≤π2+1\left\vert \mathrm{Si}(T)\right\vert \le\frac\pi2+1

  • [−R,R][-R,R] の上では、連続関数として有界である。

よって CSi<∞C_{\mathrm{Si}}<\infty である。

最後に (9.5) を θ\theta の符号で場合分けして示す。

  • θ>0\theta>0 なら、x=θtx=\theta t と置換すれば得られる。
  • θ=0\theta=0 なら、両辺とも 0 である。
  • θ<0\theta<0 なら、正弦関数と Si\mathrm{Si} がともに奇関数であることから従う。

Si\mathrm{Si} は単調ではなく、π2\frac\pi2 のまわりで振動しながら近づく。最大値は T=πT=\pi でとり、CSi=Si(π)≈1.8519C_{\mathrm{Si}}=\mathrm{Si}(\pi)\approx1.8519 である。(9.5) と極限をまとめると、T→∞T\to\infty のとき

1π∫0Tsin⁡θtt dt=1πSi(θT)⟶12sgn⁡θ={12(θ>0)0(θ=0)−12(θ<0)(9.7)\frac1\pi\int_0^T\frac{\sin\theta t}{t}\,\mathrm{d} t=\frac1\pi\mathrm{Si}(\theta T)\longrightarrow\frac12\operatorname{sgn}\theta=\begin{cases}\frac12&(\theta>0)\\0&(\theta=0)\\-\frac12&(\theta<0)\end{cases}\tag{9.7}

となる。θ=0\theta=0 での値が 0 であることが、反転公式の端点に 12\frac12 が現れる理由になる。

反転公式

定理 9.2(Lévy の反転公式)

確率変数 XX の特性関数を φ\varphi とする。実数 a<ba<b について

lim⁡T→∞12π∫−TTe−ita−e−itbit φ(t) dt=P(a<X<b)+12P(X=a)+12P(X=b)(9.8)\lim_{T\to\infty}\frac{1}{2\pi}\int_{-T}^{T}\frac{\mathrm{e}^{-\mathrm{i} ta}-\mathrm{e}^{-\mathrm{i} tb}}{\mathrm{i} t}\,\varphi(t)\,\mathrm{d} t =P(a<X<b)+\frac12P(X=a)+\frac12P(X=b)\tag{9.8}

が成り立つ。とくに a,ba,b が分布関数 FF の連続点なら、右辺は次に等しい。

F(b)−F(a)F(b)-F(a)

被積分関数は t=0t=0 で定義されていないが、

e−ita−e−itbit=∫abe−itu du(9.9)\frac{\mathrm{e}^{-\mathrm{i} ta}-\mathrm{e}^{-\mathrm{i} tb}}{\mathrm{i} t}=\int_a^b\mathrm{e}^{-\mathrm{i} tu}\,\mathrm{d} u\tag{9.9}

と書ける。よって t=0t=0 での値を b−ab-a と定めれば連続で、絶対値は区間の長さ以下である。したがって有限区間 [−T,T][-T,T] の上の積分は問題なく定まる。

左辺の極限を T→∞T\to\infty でとることが必要で、∫−∞∞\int_{-\infty}^\infty に置き換えることは一般にはできない。原子があると被積分関数は可積分でないからである。

証明

左辺の積分を ITI_T と書く。φ(t)=E[eitX]\varphi(t)=E[\mathrm{e}^{\mathrm{i} tX}] を代入すると

IT=12π∫−TTE[eit(X−a)−eit(X−b)it]dtI_T=\frac{1}{2\pi}\int_{-T}^{T}E\Bigl[\frac{\mathrm{e}^{\mathrm{i} t(X-a)}-\mathrm{e}^{\mathrm{i} t(X-b)}}{\mathrm{i} t}\Bigr]\mathrm{d} t

である。中身の絶対値は (9.9) により b−ab-a 以下で、[−T,T]×Ω[-T,T]\times\Omega は有限測度の空間である。よって Fubini の定理により

IT=E[JT(X)],JT(x)≔12π∫−TTeit(x−a)−eit(x−b)it dt(9.10)I_T=E\bigl[J_T(X)\bigr],\qquad J_T(x)\coloneqq\frac{1}{2\pi}\int_{-T}^{T}\frac{\mathrm{e}^{\mathrm{i} t(x-a)}-\mathrm{e}^{\mathrm{i} t(x-b)}}{\mathrm{i} t}\,\mathrm{d} t\tag{9.10}

と期待値と積分を交換できる。

JT(x)J_T(x) を計算する。α=x−a\alpha=x-a、β=x−b\beta=x-b と置く。指数関数を eitα=cos⁡tα+isin⁡tα\mathrm{e}^{\mathrm{i} t\alpha}=\cos t\alpha+\mathrm{i}\sin t\alpha と実部と虚部に分けると

eitα−eitβit=cos⁡tα−cos⁡tβit+sin⁡tα−sin⁡tβt\frac{\mathrm{e}^{\mathrm{i} t\alpha}-\mathrm{e}^{\mathrm{i} t\beta}}{\mathrm{i} t}=\frac{\cos t\alpha-\cos t\beta}{\mathrm{i} t}+\frac{\sin t\alpha-\sin t\beta}{t}

となる。2つの項を別々に扱う。

  • 第1項は tt の奇関数である。∣cos⁡tα−cos⁡tβ∣≤∣t∣ ∣α−β∣\left\vert \cos t\alpha-\cos t\beta\right\vert \le\left\vert t\right\vert \,\left\vert \alpha-\beta\right\vert から有界でもあるので、対称な区間の上の積分は 0 である。
  • 第2項は偶関数なので、[0,T][0,T] の上の積分の2倍になる。

(9.5) を使うと

JT(x)=1π∫0Tsin⁡t(x−a)−sin⁡t(x−b)t dt=1π(Si(T(x−a))−Si(T(x−b)))(9.11)J_T(x)=\frac{1}{\pi}\int_0^T\frac{\sin t(x-a)-\sin t(x-b)}{t}\,\mathrm{d} t=\frac1\pi\Bigl(\mathrm{Si}\bigl(T(x-a)\bigr)-\mathrm{Si}\bigl(T(x-b)\bigr)\Bigr)\tag{9.11}

を得る。

この式から2つのことが分かる。

  1. すべての xx と TT で次の評価が成り立つ。

∣JT(x)∣≤2CSiπ\left\vert J_T(x)\right\vert \le\frac{2C_{\mathrm{Si}}}{\pi} 2. (9.7) から T→∞T\to\infty のとき

JT(x)⟶J(x)≔12sgn⁡(x−a)−12sgn⁡(x−b)={0(x<a または x>b)12(x=a または x=b)1(a<x<b)J_T(x)\longrightarrow J(x)\coloneqq\frac12\operatorname{sgn}(x-a)-\frac12\operatorname{sgn}(x-b) =\begin{cases}0&(x<a\ \text{または}\ x>b)\\[2pt]\frac12&(x=a\ \text{または}\ x=b)\\[2pt]1&(a<x<b)\end{cases}

となる。たとえば x=ax=a では、符号関数の値が次のようになる。

sgn⁡(x−a)=0,sgn⁡(x−b)=−1,J(a)=12\operatorname{sgn}(x-a)=0,\qquad\operatorname{sgn}(x-b)=-1,\qquad J(a)=\tfrac12

JT(X)J_T(X) は定数で抑えられ、各標本点で J(X)J(X) に収束する。よって優収束定理により

lim⁡T→∞IT=E[J(X)]=P(a<X<b)+12P(X=a)+12P(X=b)\lim_{T\to\infty}I_T=E[J(X)]=P(a<X<b)+\frac12P(X=a)+\frac12P(X=b)

を得る。a,ba,b が FF の連続点なら、端点の確率は 0 である(第3章 §3.3)。よって右辺は次に等しい。

P(a<X≤b)=F(b)−F(a)P(a<X\le b)=F(b)-F(a)

図 9.1 に、区間を (a,b)=(0,1)(a,b)=(0,1) としたときの JT(x)J_T(x) を示す。TT が大きくなるにつれて、JTJ_T はこの区間の指示関数に近づく。端点での値は

JT(0)=JT(1)=1πSi(T)J_T(0)=J_T(1)=\frac1\pi\mathrm{Si}(T)

で、T→∞T\to\infty のとき 12\frac12 に近づく。

端点の近くには一定の高さのはみ出しが残る(Fourier 級数の Gibbs 現象と同じ現象である)。このはみ出しは JTJ_T の各点収束を妨げないので、証明には影響しない。

図

図 9.1 反転公式の核 JT(x)=1π(Si(Tx)−Si(T(x−1)))J_T(x)=\frac1\pi\bigl(\mathrm{Si}(Tx)-\mathrm{Si}(T(x-1))\bigr)(a=0a=0、b=1b=1)。破線は T=10T=10、太線は T=60T=60。T→∞T\to\infty で区間 (0,1)(0,1) の指示関数に近づき、端点での値は極限値 12\frac12(点)に近づく。

密度の反転公式と原子の公式

特性関数が可積分なら、極限をとる必要がなくなり、密度そのものを復元できる。

定理 9.3(密度の反転公式)

∫−∞∞∣φ(t)∣ dt<∞\displaystyle\int_{-\infty}^{\infty}\left\vert \varphi(t)\right\vert \,\mathrm{d} t<\infty なら、XX は有界で連続な密度

f(x)=12π∫−∞∞e−itxφ(t) dt(9.12)f(x)=\frac{1}{2\pi}\int_{-\infty}^{\infty}\mathrm{e}^{-\mathrm{i} tx}\varphi(t)\,\mathrm{d} t\tag{9.12}

をもつ。∣f(x)∣≤12π∫∣φ∣ dt\left\vert f(x)\right\vert \le\frac{1}{2\pi}\int\left\vert \varphi\right\vert \,\mathrm{d} t である。

証明

まず原子がないことを示す。(9.8) の被積分関数は、(9.9) から (b−a)∣φ(t)∣(b-a)\left\vert \varphi(t)\right\vert 以下である。これは実数全体の上で可積分なので、優収束定理により極限は全体の積分に置き換えられ、

P(a<X<b)+12P(X=a)+12P(X=b)=12π∫−∞∞e−ita−e−itbit φ(t) dt(9.13)P(a<X<b)+\frac12P(X=a)+\frac12P(X=b)=\frac{1}{2\pi}\int_{-\infty}^{\infty}\frac{\mathrm{e}^{-\mathrm{i} ta}-\mathrm{e}^{-\mathrm{i} tb}}{\mathrm{i} t}\,\varphi(t)\,\mathrm{d} t\tag{9.13}

となる。右辺の絶対値は b−a2π∫∣φ∣ dt\frac{b-a}{2\pi}\int\left\vert \varphi\right\vert \,\mathrm{d} t 以下である。

点 xx と h>0h>0 を取る。区間を (a,b)=(x−h,x+h)(a,b)=(x-h,x+h) とする。左辺は P(X=x)P(X=x) 以上なので

P(X=x)≤2h2π∫−∞∞∣φ(t)∣ dtP(X=x)\le\frac{2h}{2\pi}\int_{-\infty}^{\infty}\left\vert \varphi(t)\right\vert \,\mathrm{d} t

である。h→0h\to0 として P(X=x)=0P(X=x)=0 を得る。原子がないので、(9.13) の左辺は F(b)−F(a)F(b)-F(a) に等しい。

次に (9.12) の ff を調べる。

  • 被積分関数の絶対値は ∣φ∣\left\vert \varphi\right\vert 以下で可積分なので、ff は定まる。同じ評価から ∣f∣≤12π∫∣φ∣\left\vert f\right\vert \le\frac1{2\pi}\int\left\vert \varphi\right\vert である。
  • xn→xx_n\to x なら被積分関数は各 tt で収束する。よって優収束定理から ff は連続である。

(9.9) を (9.13) に代入する。[a,b]×R[a,b]\times\mathbb{R} の上で ∣e−ituφ(t)∣\left\vert \mathrm{e}^{-\mathrm{i} tu}\varphi(t)\right\vert は可積分なので、Fubini の定理を使えて

F(b)−F(a)=12π∫−∞∞(∫abe−itu du)φ(t) dt=∫abf(u) duF(b)-F(a)=\frac{1}{2\pi}\int_{-\infty}^{\infty}\Bigl(\int_a^b\mathrm{e}^{-\mathrm{i} tu}\,\mathrm{d} u\Bigr)\varphi(t)\,\mathrm{d} t=\int_a^b f(u)\,\mathrm{d} u

となる。

左辺は実数で非負である。もし Im⁡f(u0)≠0\operatorname{Im} f(u_0)\neq0 や Re⁡f(u0)<0\operatorname{Re} f(u_0)<0 となる点があったとする。ff の連続性から、その近くの小さな区間 [a,b][a,b] で右辺の虚部が 0 でないか実部が負になり、矛盾する。

よって ff は実数値で非負であり、すべての区間で F(b)−F(a)=∫abfF(b)-F(a)=\int_a^bf が成り立つ。a→−∞a\to-\infty とすれば F(b)=∫−∞bfF(b)=\int_{-\infty}^bf となる。したがって ff は XX の密度である。

逆向きは成り立たない。一様分布 U(−1,1)\mathrm U(-1,1) は密度をもつが、その特性関数 sin⁡tt\frac{\sin t}{t} は可積分でない。実際

∫−∞∞∣sin⁡t∣∣t∣ dt=∞\int_{-\infty}^{\infty}\frac{\left\vert \sin t\right\vert }{\left\vert t\right\vert }\,\mathrm{d} t=\infty

である。定理 9.3 の結論の密度は連続なので、密度に跳びがある分布の特性関数は可積分になりえない。

原子の確率も特性関数から読み取れる。

命題 9.4(原子の公式)

すべての実数 aa について

P(X=a)=lim⁡T→∞12T∫−TTe−itaφ(t) dt.(9.14)P(X=a)=\lim_{T\to\infty}\frac{1}{2T}\int_{-T}^{T}\mathrm{e}^{-\mathrm{i} ta}\varphi(t)\,\mathrm{d} t .\tag{9.14}

証明

被積分関数は絶対値 1 以下なので Fubini の定理が使え、右辺の積分は次に等しい。

E[12T∫−TTeit(X−a) dt]E\Bigl[\frac{1}{2T}\int_{-T}^{T}\mathrm{e}^{\mathrm{i} t(X-a)}\,\mathrm{d} t\Bigr]

中の積分を実数 cc について計算すると

12T∫−TTeitc dt={12T∫−TTcos⁡tc dt=sin⁡TcTc(c≠0)1(c=0)\frac1{2T}\int_{-T}^T\mathrm{e}^{\mathrm{i} tc}\,\mathrm{d} t=\begin{cases}\dfrac{1}{2T}\displaystyle\int_{-T}^T\cos tc\,\mathrm{d} t=\dfrac{\sin Tc}{Tc}&(c\neq0)\\[8pt]1&(c=0)\end{cases}

である。この関数は絶対値 1 以下である。

T→∞T\to\infty のとき、この関数は 1{c=0}\mathbf 1_{\{c=0\}} に収束する。つまり極限は、0 でない cc では 0、c=0c=0 では 1 である。c=X−ac=X-a と置いて優収束定理を使うと、右辺の極限は E[1{X=a}]=P(X=a)E[\mathbf 1_{\{X=a\}}]=P(X=a) になる。

整数値の確率変数では φ\varphi が周期 2π2\pi をもつ。よって (9.14) の平均は1周期の平均に一致し、

P(X=k)=12π∫−ππe−iktφ(t) dtP(X=k)=\frac{1}{2\pi}\int_{-\pi}^{\pi}\mathrm{e}^{-\mathrm{i} kt}\varphi(t)\,\mathrm{d} t

となる。これは第8章の確認テスト 8.3-5 で直接示した式である。

例題 9.1 Laplace 分布から Cauchy 分布の特性関数を得る ★★

オリジナルLaplace 分布(密度 12e−∣x∣\frac12\mathrm{e}^{-\left\vert x\right\vert })の特性関数は 11+t2\frac{1}{1+t^2} である(第8章 §8.3)。

  1. 定理 9.3 を使い、∫−∞∞e−itx1+t2 dt=π e−∣x∣\displaystyle\int_{-\infty}^{\infty}\frac{\mathrm{e}^{-\mathrm{i} tx}}{1+t^2}\,\mathrm{d} t=\pi\,\mathrm{e}^{-\left\vert x\right\vert } を示せ。
  2. (1) から、標準 Cauchy 分布 Cauchy(0,1)\mathrm{Cauchy}(0,1) の特性関数が e−∣t∣\mathrm{e}^{-\left\vert t\right\vert } であることを導け。
  3. Cauchy(m,γ)\mathrm{Cauchy}(m,\gamma) の特性関数を求めよ。
解答
  1. φ(t)=11+t2\varphi(t)=\frac{1}{1+t^2} は ∫∣φ∣ dt=π<∞\int\left\vert \varphi\right\vert \,\mathrm{d} t=\pi<\infty を満たす。定理 9.3 から XX は連続な密度 12π∫e−itxφ(t) dt\frac1{2\pi}\int\mathrm{e}^{-\mathrm{i} tx}\varphi(t)\,\mathrm{d} t をもつ。密度は Lebesgue 測度 Leb\mathrm{Leb} について a.e. の違いを除いて一意で、ここでは両方とも連続なので、すべての xx で

12e−∣x∣=12π∫−∞∞e−itx1+t2 dt\frac12\mathrm{e}^{-\left\vert x\right\vert }=\frac{1}{2\pi}\int_{-\infty}^{\infty}\frac{\mathrm{e}^{-\mathrm{i} tx}}{1+t^2}\,\mathrm{d} t

である(連続な2つの関数が a.e. で等しければ、等しくない点の集合は開集合で測度 0 なので空集合である)。両辺に 2π2\pi を掛ければよい。2. (1) の式で、文字を2つ書き換える。点 xx を −s-s に替える。さらに積分変数を xx と書き直すと

∫−∞∞eisx 1π(1+x2) dx=1π⋅π e−∣−s∣=e−∣s∣\int_{-\infty}^{\infty}\mathrm{e}^{\mathrm{i} sx}\,\frac{1}{\pi(1+x^2)}\,\mathrm{d} x=\frac1\pi\cdot\pi\,\mathrm{e}^{-\left\vert -s\right\vert }=\mathrm{e}^{-\left\vert s\right\vert }

となる。左辺は標準 Cauchy 分布の特性関数 φ(s)\varphi(s) の定義そのものである。留数定理を使わずに e−∣s∣\mathrm{e}^{-\left\vert s\right\vert } が得られた。3. X∼Cauchy(0,1)X\sim\mathrm{Cauchy}(0,1) なら m+γX∼Cauchy(m,γ)m+\gamma X\sim\mathrm{Cauchy}(m,\gamma) である(密度の変換、第7章 §7.1)。(9.1) の最後の式から

φm+γX(t)=eimtφX(γt)=exp⁡(imt−γ∣t∣)\varphi_{m+\gamma X}(t)=\mathrm{e}^{\mathrm{i} mt}\varphi_X(\gamma t)=\exp\bigl(\mathrm{i} mt-\gamma\left\vert t\right\vert \bigr)

である。

Laplace 分布の密度と Cauchy 分布の特性関数はどちらも e−∣⋅∣\mathrm{e}^{-\left\vert \cdot\right\vert } の形で、Laplace 分布の特性関数と Cauchy 分布の密度はどちらも 11+(⋅)2\frac{1}{1+(\cdot)^2} の形である。

例題 9.2 台が有界な特性関数 ★★

オリジナル三角形の密度 g(x)=(1−∣x∣)+g(x)=(1-\left\vert x\right\vert )^+((⋅)+(\cdot)^+ は正の部分)をもつ分布を考える。

  1. gg の特性関数が次であることを示せ。

ψ(t)=2(1−cos⁡t)t2(t≠0),ψ(0)=1\psi(t)=\frac{2(1-\cos t)}{t^2}\quad(t\neq0),\qquad\psi(0)=1 2. ∫−∞∞1−cos⁡tt2 dt=π\displaystyle\int_{-\infty}^{\infty}\frac{1-\cos t}{t^2}\,\mathrm{d} t=\pi を示せ。3. h(x)=1−cos⁡xπx2h(x)=\dfrac{1-\cos x}{\pi x^2} は確率密度で、その特性関数は (1−∣t∣)+(1-\left\vert t\right\vert )^+ であることを示せ。

解答
  1. gg は偶関数なので虚部は消え、

ψ(t)=2∫01(1−x)cos⁡tx dx=2[(1−x)sin⁡txt]01+2t∫01sin⁡tx dx=0+2t⋅1−cos⁡tt\psi(t)=2\int_0^1(1-x)\cos tx\,\mathrm{d} x=2\Bigl[\frac{(1-x)\sin tx}{t}\Bigr]_0^1+\frac2t\int_0^1\sin tx\,\mathrm{d} x=0+\frac{2}{t}\cdot\frac{1-\cos t}{t}

となる(部分積分)。t→0t\to0 では 1−cos⁡t∼t2/21-\cos t\sim t^2/2 である。よって ψ(t)→1=ψ(0)\psi(t)\to1=\psi(0) である。2. 0≤ψ(t)≤min⁡(1,4/t2)0\le\psi(t)\le\min(1,4/t^2) なので ψ\psi は可積分である。定理 9.3 の密度の公式 (9.12) を x=0x=0 で使うと、gg は連続なので(例題 9.1 (1) と同じ理由で)すべての点で一致し、

1=g(0)=12π∫−∞∞ψ(t) dt=1π∫−∞∞1−cos⁡tt2 dt1=g(0)=\frac{1}{2\pi}\int_{-\infty}^{\infty}\psi(t)\,\mathrm{d} t=\frac1\pi\int_{-\infty}^{\infty}\frac{1-\cos t}{t^2}\,\mathrm{d} t

を得る。3. h≥0h\ge0 で、(2) から ∫h=1\int h=1 なので確率密度である。(9.12) を一般の xx で書くと

(1−∣x∣)+=12π∫−∞∞e−itx 2(1−cos⁡t)t2 dt=∫−∞∞e−itx h(t) dt(1-\left\vert x\right\vert )^+=\frac{1}{2\pi}\int_{-\infty}^{\infty}\mathrm{e}^{-\mathrm{i} tx}\,\frac{2(1-\cos t)}{t^2}\,\mathrm{d} t=\int_{-\infty}^{\infty}\mathrm{e}^{-\mathrm{i} tx}\,h(t)\,\mathrm{d} t

である。点 xx を −s-s と書く。さらに積分変数の名前を yy に替えると、右辺は

∫−∞∞eisyh(y) dy\int_{-\infty}^{\infty}\mathrm{e}^{\mathrm{i} sy}h(y)\,\mathrm{d} y

となり、密度 hh の特性関数そのものになる。よって hh の特性関数は (1−∣s∣)+(1-\left\vert s\right\vert )^+ である。

特性関数 (1−∣t∣)+(1-\left\vert t\right\vert )^+ は ∣t∣≥1\left\vert t\right\vert \ge1 で 0 になる。一方で、hh の裾は x−2x^{-2} の程度なので、平均は存在しない。

∫−∞∞∣x∣h(x) dx=∞\int_{-\infty}^{\infty}\left\vert x\right\vert h(x)\,\mathrm{d} x=\infty

特性関数が t=0t=0 で尖っていて微分できないことと、平均が存在しないことが対応している。

確認テスト

オリジナル

  1. 9.1-1 定数 X=cX=c の特性関数は eitc\mathrm{e}^{\mathrm{i} tc} である。a=c<ba=c<b のとき、(9.8) の左辺の積分を Si 関数で書け。また、その極限が 12\frac12 になることを確かめよ。
  2. 9.1-2 φ(t)=e−∣t∣\varphi(t)=\mathrm{e}^{-\left\vert t\right\vert } に密度の反転公式 (9.12) を使う。次の積分を直接計算して、右辺を得よ。

12π∫−∞∞e−itxe−∣t∣ dt=1π(1+x2)\frac1{2\pi}\int_{-\infty}^{\infty}\mathrm{e}^{-\mathrm{i} tx}\mathrm{e}^{-\left\vert t\right\vert }\,\mathrm{d} t=\frac{1}{\pi(1+x^2)} 3. 9.1-3 次の混合分布を考える。確率 12\frac12 で X=0X=0 とする。残りの確率で X∼N(0,1)X\sim\mathrm N(0,1) とする。その特性関数は

φ(t)=12+12e−t2/2\varphi(t)=\frac12+\frac12\mathrm{e}^{-t^2/2}

である。原子の公式 (9.14) を使い、P(X=0)P(X=0) と P(X=1)P(X=1) を求めよ。4. 9.1-4 XX の原子を x1,x2,…x_1,x_2,\dots とする。次を示せ。

lim⁡T→∞12T∫−TT∣φX(t)∣2 dt=∑kP(X=xk)2\lim_{T\to\infty}\frac{1}{2T}\int_{-T}^{T}\left\vert \varphi_X(t)\right\vert ^2\,\mathrm{d} t=\sum_kP(X=x_k)^2

(ヒント:XX と独立で同分布な YY をとる。その差の特性関数を考える。)

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『基礎理論』の記事一覧へ