連続型確率変数の特性関数
確率空間\(\left( \Omega ,\mathcal{F},P\right) \)に加えて絶対連続型の確率変数\begin{equation*}X:\Omega \rightarrow \mathbb{R} \end{equation*}が与えられているものとします。さらに、確率変数\(X\)の確率分布が確率密度関数\begin{equation*}f_{X}:\mathbb{R} \rightarrow \mathbb{R} \end{equation*}によって記述されているものとします。つまり、\begin{equation*}
\forall B\in \mathcal{B}\left( \mathbb{R} \right) :P\left( X\in B\right) =\int_{B}f_{X}\left( x\right) dx
\end{equation*}が成り立つということです。
実数\(t\in \mathbb{R} \)を任意に選んだ上で、変数\(x\in \mathbb{R} \)に関する複素数値関数\begin{equation*}e^{itx}:\mathbb{R} \rightarrow \mathbb{C} \end{equation*}を定義します。ただし、\(i\in \mathbb{C} \)は虚数単位です。関数\(e^{itx}\)の定義域は\(\mathbb{R} \)であるため、これと先の確率変数\(X\)を用いることにより、それぞれの標本点\(\omega \in \Omega \)に対して、以下の複素数\begin{equation*}e^{itX\left( \omega \right) }\in \mathbb{C} \end{equation*}を値として定める新たな複素数値関数\begin{equation*}
e^{itX}:\Omega \rightarrow \mathbb{C} \end{equation*}が定義可能です。つまり、\(e^{itX}\)は\(X\)と\(e^{itx}\)の合成関数です。この関数\(e^{itX}\)は複素数を値としてとり得ますが、これもまた確率変数とみなした上で、その期待値を、\begin{equation}E\left( e^{itX}\right) =\int_{-\infty }^{+\infty }e^{itx}\cdot f_{X}\left(
x\right) dx \quad \cdots (1)
\end{equation}と定義します。
オイラーの等式より、任意の\(\omega \in \Omega \)に対して、\begin{equation*}e^{itX\left( \omega \right) }=\cos \left( tX\left( \omega \right) \right)
+i\sin \left( tX\left( \omega \right) \right)
\end{equation*}となるため、\begin{equation}
e^{itX}=\cos \left( tX\right) +i\sin \left( tX\right) \quad \cdots (2)
\end{equation}を得ます。したがって、\begin{eqnarray*}
E\left( e^{itX}\right) &=&\int_{-\infty }^{+\infty }e^{itx}\cdot
f_{X}\left( x\right) dx\quad \because \left( 1\right) \\
&=&\int_{-\infty }^{+\infty }\left[ \cos \left( tx\right) +i\sin \left(
tx\right) \right] \cdot f_{X}\left( x\right) dx\quad \because \left(
2\right) \\
&=&\int_{-\infty }^{+\infty }\cos \left( tx\right) \cdot f_{X}\left(
x\right) dx+i\int_{-\infty }^{+\infty }\sin \left( tx\right) \cdot
f_{X}\left( x\right) dx \\
&=&E\left( \cos \left( tX\right) \right) +iE\left( \sin \left( tX\right)
\right) \quad \because \text{LOTUS}
\end{eqnarray*}を得ます。正弦関数および余弦関数は\(\left[ -1,1\right] \)上の値をとる有界関数であるため\(E\left( \cos \left(tX\right) \right) \)および\(E\left( \sin \left( tX\right)\right) \)は有限な実数として定まり、したがって\(E\left( \cos \left( tX\right) \right) +iE\left( \sin \left(tX\right) \right) \)すなわち\(E\left( e^{itX}\right) \)は1つの複素数として定まることが保証されます。任意の\(t\in \mathbb{R} \)について同様の議論が成立します。
このような事情を踏まえると、絶対連続型の確率変数\(X:\Omega \rightarrow \mathbb{R} \)と確率密度関数\(f_{X}:\mathbb{R} \rightarrow \mathbb{R} \)が与えられた場合には、それぞれの\(t\in \mathbb{R} \)に対して、\begin{equation*}\varphi _{X}\left( t\right) =E\left( e^{itX}\right)
\end{equation*}を値として定める複素数値関数\begin{equation*}
\varphi _{X}:\mathbb{R} \rightarrow \mathbb{C} \end{equation*}が定義可能です。これを確率変数\(X\)の特性関数(characteristic function of \(X\))と呼びます。
\begin{array}{cl}
1 & \left( if\ x\in X\left( \Omega \right) \right) \\
0 & \left( otherwise\right)
\end{array}\right.
\end{equation*}を定めるものとします。特性関数\(\varphi _{X}:\mathbb{R} \rightarrow \mathbb{C} \)がそれぞれの\(t\in \mathbb{R} \backslash \left\{ 0\right\} \)に対して定める値は、\begin{eqnarray*}\varphi _{X}\left( t\right) &=&E\left( e^{itX}\right) \\
&=&\int_{-\infty }^{+\infty }e^{itx}\cdot f_{X}\left( x\right) dx \\
&=&\int_{-\infty }^{+\infty }\cos \left( tx\right) \cdot f_{X}\left(
x\right) dx+i\int_{-\infty }^{+\infty }\sin \left( tx\right) \cdot
f_{X}\left( x\right) dx \\
&=&\int_{0}^{1}\cos \left( tx\right) dx+i\int_{0}^{1}\sin \left( tx\right) dx
\\
&=&\left[ \frac{1}{t}\sin \left( tx\right) \right] _{x=0}^{1}+i\left[ -\frac{1}{t}\cos \left( tx\right) \right] _{x=0}^{1} \\
&=&\left[ \frac{1}{t}\sin \left( t\right) -\frac{1}{t}\sin \left( 0\right) \right] +i\left[ -\frac{1}{t}\cos \left( t\right) +\frac{1}{t}\cos \left(
0\right) \right] \\
&=&\frac{1}{t}\sin \left( t\right) -\frac{1}{t}\sin \left( 0\right) -\frac{i}{t}\cos \left( t\right) +\frac{i}{t}\cos \left( 0\right) \\
&=&\frac{1}{it}\left[ i\sin \left( t\right) -i\sin \left( 0\right)
-i^{2}\cos \left( t\right) +i^{2}\cos \left( 0\right) \right] \\
&=&\frac{1}{it}\left[ i\sin \left( t\right) -i\sin \left( 0\right) +\cos
\left( t\right) -\cos \left( 0\right) \right] \\
&=&\frac{1}{it}\left\{ \left[ \cos \left( t\right) +i\sin \left( t\right) \right] -\left[ \cos \left( 0\right) +i\sin \left( 0\right) \right] \right\}
\\
&=&\frac{e^{it}-e^{i0}}{it} \\
&=&\frac{e^{it}-1}{it}
\end{eqnarray*}となります。\(t=0\)の場合には、\begin{eqnarray*}\varphi _{X}\left( 0\right) &=&E\left( e^{0}\right) \\
&=&E\left( 1\right) \\
&=&1
\end{eqnarray*}となります。
特性関数とモーメントの関係
特性関数はモーメントと深い関係があります。モーメントについて簡単に復習した上で、特性関数との関係を以下で解説します。
絶対連続型の確率変数\(X:\Omega \rightarrow \mathbb{R} \)の確率分布が確率密度関数\(f_{X}:\mathbb{R} \rightarrow \mathbb{R} \)として記述されているものとします。自然数\(m\in \mathbb{N} \)を任意に選んだとき、それぞれの\(\omega \in \Omega \)に対して、\begin{equation*}X^{m}\left( \omega \right) =\left[ X\left( \omega \right) \right] ^{m}
\end{equation*}を定める新たな確率変数\(X^{m}:\Omega \rightarrow \mathbb{R} \)が定義可能です。この確率変数\(X^{m}\)の期待値は、\begin{equation*}E\left( X^{m}\right) =\int_{-\infty }^{+\infty }x^{m}f_{X}\left( x\right)
dx\quad \because \text{LOTUS}
\end{equation*}となりますが、これを確率変数\(X\)の\(m\)次のモーメントと呼びます。
確率変数\(X\)の特性関数\(\varphi _{X}:\mathbb{R} \rightarrow \mathbb{C} \)が与えられた状況を想定します。さらに、自然数\(m\in \mathbb{N} \)について、確率変数\(X\)の\(m\)次のモーメント\(E\left(X^{m}\right) \)が有限な実数として定まる場合には、それは特性関数\(\varphi _{X}\)の点\(0\)における\(m\)次の微分係数と\(\frac{1}{i^{m}}\)の積と一致することが保証されます。つまり、以下の関係\begin{equation*}E\left( X^{m}\right) =\frac{1}{i^{m}}\cdot \left. \frac{d^{m}\varphi
_{X}\left( t\right) }{dt^{m}}\right\vert _{t=0}
\end{equation*}が成り立つということです。
_{X}\left( t\right) }{dt^{m}}\right\vert _{t=0}
\end{equation*}が成り立つ。
\begin{array}{cl}
1 & \left( if\ x\in X\left( \Omega \right) \right) \\
0 & \left( otherwise\right)
\end{array}\right.
\end{equation*}を定めるものとします。先に示したように、\(X\)の特性関数\(\varphi _{X}:\mathbb{R} \rightarrow \mathbb{C} \)はそれぞれの\(t\in \mathbb{R} \)に対して、\begin{equation*}\varphi _{X}\left( t\right) =\left\{
\begin{array}{cl}
\dfrac{e^{it}-1}{it} & \left( if\ t\not=0\right) \\
1 & \left( if\ t=0\right)
\end{array}\right.
\end{equation*}を定めます。先の命題より、\(m\in \mathbb{N} \)を任意に選んだとき、以下の関係\begin{equation}E\left( X^{m}\right) =\frac{1}{i^{m}}\cdot \left. \frac{d^{m}\varphi
_{X}\left( t\right) }{dt^{m}}\right\vert _{t=0} \quad \cdots (1)
\end{equation}が成立します。複素指数関数の定義より、任意の\(z\in \mathbb{C} \)について、\begin{equation*}e^{z}=1+z+\frac{z^{2}}{2}+\frac{z^{3}}{3!}+\cdots
\end{equation*}が成り立ちます。\(t\in \mathbb{R} \backslash \left\{ 0\right\} \)を任意に選んだとき\(it\in \mathbb{C} \)であるため、\begin{equation*}e^{it}=1+it+\frac{\left( it\right) ^{2}}{2}+\frac{\left( it\right) ^{3}}{3!}+\cdots
\end{equation*}を得ます。したがって、\begin{eqnarray*}
\dfrac{e^{it}-1}{it} &=&\frac{1+it+\frac{\left( it\right) ^{2}}{2}+\frac{\left( it\right) ^{3}}{3!}+\cdots -1}{it} \\
&=&\frac{it+\frac{\left( it\right) ^{2}}{2}+\frac{\left( it\right) ^{3}}{3!}+\cdots }{it} \\
&=&1+\frac{it}{2}+\frac{\left( it\right) ^{2}}{3!}+\frac{\left( it\right)
^{3}}{4!}+\cdots
\end{eqnarray*}を得ます。すると、\begin{eqnarray*}
\frac{d}{dt}\left( \dfrac{e^{it}-1}{it}\right) &=&\frac{i}{2}+\frac{2i^{2}}{3!}t+\frac{3i^{3}}{4!}t^{2}+\cdots \\
\frac{d^{2}}{dt^{2}}\left( \dfrac{e^{it}-1}{it}\right) &=&\frac{2\cdot
1\cdot i^{2}}{3!}+\frac{3\cdot 2\cdot i^{3}}{4!}t+\frac{4\cdot 3\cdot i^{4}}{5!}t^{2}+\cdots \\
\frac{d^{3}}{dt^{3}}\left( \dfrac{e^{it}-1}{it}\right) &=&\frac{3\cdot
2\cdot 1\cdot i^{3}}{4!}t+\frac{4\cdot 3\cdot 2\cdot i^{4}}{5!}t^{2}+\frac{5\cdot 4\cdot 3\cdot i^{5}}{6!}t^{3}+\cdots
\end{eqnarray*}などとなり、一般的には、\begin{equation*}
\frac{d^{m}}{dt^{m}}\left( \dfrac{e^{it}-1}{it}\right) =\frac{m!\cdot i^{m}}{\left( m+1\right) !}+\frac{\left( m+1\right) m\times \cdots \times 2\times
i^{m+1}}{\left( m+2\right) !}t+\cdots
\end{equation*}となるため、\begin{eqnarray*}
\left. \frac{d^{m}}{dt^{m}}\left( \dfrac{e^{it}-1}{it}\right) \right\vert
_{t=0} &=&\frac{m!\cdot i^{m}}{\left( m+1\right) !} \\
&=&\frac{i^{m}}{m+1}
\end{eqnarray*}を得ます。これと\(\left(1\right) \)より、\begin{eqnarray*}E\left( X^{m}\right) &=&\frac{1}{i^{m}}\cdot \frac{i^{m}}{m+1} \\
&=&\frac{1}{m+1}
\end{eqnarray*}を得ます。他方で、定義にもとづいてモーメントを求めると、\begin{eqnarray*}
E\left( X^{m}\right) &=&\int_{-\infty }^{+\infty }x^{m}f_{X}\left( x\right)
dx\quad \because \text{LOTUS} \\
&=&\int_{0}^{1}x^{m}\cdot 1dx \\
&=&\int_{0}^{1}x^{m}dx \\
&=&\left[ \frac{x^{m+1}}{m+1}\right] _{0}^{1} \\
&=&\frac{1^{m+1}}{m+1}-\frac{0^{m+1}}{m+1} \\
&=&\frac{1}{m+1}
\end{eqnarray*}となりますが、これは先の結果と整合的です。
\end{equation*}が成り立ちます。つまり、\(X\)の原点まわりの\(1\)次のモーメントは\(X\)の期待値と一致します。したがって、\(E\left(X^{1}\right) \)が有限な実数として定まる場合、先の命題より、\begin{equation*}E\left( X^{1}\right) =\frac{1}{i}\cdot \left. \frac{d\varphi _{X}\left(
t\right) }{dt}\right\vert _{t=0}
\end{equation*}が成り立つため、\begin{equation*}
E\left( X\right) =\frac{1}{i}\cdot \left. \frac{d\varphi _{X}\left( t\right)
}{dt}\right\vert _{t=0}
\end{equation*}を得ます。特性関数から期待値を導出できるということです。
\end{equation*}が成り立ちます。また、\(X\)の期待値とモーメントの間には以下の関係\begin{equation*}E\left( X\right) =E\left( X^{1}\right)
\end{equation*}が成り立ちます。\(E\left(X^{1}\right) \)および\(E\left( X^{2}\right) \)が有限な実数として定まる場合、先の命題より、\begin{eqnarray*}E\left( X^{1}\right) &=&\frac{1}{i}\cdot \left. \frac{d\varphi _{X}\left(
t\right) }{dt}\right\vert _{t=0} \\
E\left( X^{2}\right) &=&\frac{1}{i^{2}}\cdot \left. \frac{d^{2}\varphi
_{X}\left( t\right) }{dt^{2}}\right\vert _{t=0}=-\left. \frac{d^{2}\varphi
_{X}\left( t\right) }{dt^{2}}\right\vert _{t=0}
\end{eqnarray*}を得るため、\begin{eqnarray*}
\mathrm{Var}\left( X\right) &=&-\left. \frac{d^{2}\varphi _{X}\left( t\right)
}{dt^{2}}\right\vert _{t=0}-\left[ \frac{1}{i}\cdot \left. \frac{d\varphi
_{X}\left( t\right) }{dt}\right\vert _{t=0}\right] ^{2} \\
&=&-\left. \frac{d^{2}\varphi _{X}\left( t\right) }{dt^{2}}\right\vert
_{t=0}+\left[ \left. \frac{d\varphi _{X}\left( t\right) }{dt}\right\vert
_{t=0}\right] ^{2}
\end{eqnarray*}を得ます。特性関数から分散を導出できるということです。
特性関数による確率分布の特徴づけ
確率変数の確率分布は特性関数によって特徴づけられます。つまり、2つの確率変数が同一の特性関数を持つことと、それらの確率変数が同一の確率分布にしたがうことは必要十分です。順番に解説します。
2つの絶対連続型の確率変数\(X,Y:\Omega \rightarrow \mathbb{R} \)の確率分布が確率密度関数\(f_{X},f_{Y}:\mathbb{R} \rightarrow \mathbb{R} \)としてそれぞれ表現されているものとします。加えて、これらの確率変数\(X,Y\)の特性関数が\(\varphi _{X},\varphi _{Y}:\mathbb{R} \rightarrow \mathbb{C} \)であるものとします。このとき、\(f_{X}\)と\(f_{Y}\)が一致することと\(\varphi _{X}\)と\(\varphi _{Y}\)が一致することは必要十分になります。
&&\left( b\right) \ \forall t\in \mathbb{R} :\varphi _{X}\left( t\right) =\varphi _{Y}\left( t\right)
\end{eqnarray*}は必要十分である。
離散型確率変数の確率分布は分布関数によって表現することもできるため、上の命題より以下を得ます。
&&\left( b\right) \ \forall t\in \mathbb{R} :\varphi _{X}\left( t\right) =\varphi _{Y}\left( t\right)
\end{eqnarray*}は必要十分である。
多くの場合、2つの確率変数の確率分布が一致することを示すことよりも、それらの特性関数が一致することを示すことの方が容易です。上の命題は、2つの確率変数が同一の確率分布にしたがうことを示す代わりに、それらの特性関数が一致することを示してもよいことを保証します。
連続型確率変数との合成関数の特性関数
絶対連続型の確率変数\(X:\Omega \rightarrow \mathbb{R} \)の確率分布が確率密度関数\(f_{X}:\mathbb{R} \rightarrow \mathbb{R} \)によって記述されているものとします。連続関数\(g:\mathbb{R} \rightarrow \mathbb{R} \)を任意に選んだ上で、それぞれの\(\omega \in \Omega \)に対して、\begin{eqnarray*}Y\left( \omega \right) &=&\left( g\circ X\right) \left( \omega \right) \\
&=&g\left( X\left( \omega \right) \right) \quad \because \text{合成関数の定義}
\end{eqnarray*}を定める新たな関数\begin{equation*}
Y:\Omega \rightarrow \mathbb{R} \end{equation*}を定義します。確率変数と連続関数の合成関数は確率変数であるため\(Y\)は確率変数です。\(Y\)の特性関数\(\varphi_{Y}:\mathbb{R} \rightarrow \mathbb{C} \)はそれぞれの\(t\in \mathbb{R} \)に対して、\begin{equation*}\varphi _{Y}\left( t\right) =E\left( e^{itg\left( X\right) }\right)
\end{equation*}を定めます。
\end{equation*}を定める確率変数\(Y:\Omega\rightarrow \mathbb{R} \)を定義する。\(Y\)の特性関数\(\varphi _{Y}:\mathbb{R} \rightarrow \mathbb{C} \)はそれぞれの\(t\in \mathbb{R} \)に対して、\begin{equation*}\varphi _{Y}\left( t\right) =E\left( e^{itg\left( X\right) }\right)
\end{equation*}を定める。
連続型確率変数の定数倍の特性関数
絶対連続型の確率変数\(X:\Omega \rightarrow \mathbb{R} \)と実数\(c\in \mathbb{R} \)が与えられたとき、それぞれの\(\omega \in \Omega \)に対して、\begin{equation*}\left( cX\right) \left( \omega \right) =cX\left( \omega \right)
\end{equation*}を定める新たな関数\begin{equation*}
cX:\Omega \rightarrow \mathbb{R} \end{equation*}が定義可能です。確率変数の定数倍は確率変数であるため\(cX\)は確率変数です。さらに、\(X\)の特性関数\(\varphi _{X}:\mathbb{R} \rightarrow \mathbb{C} \)と\(Y\)の特性関数\(\varphi _{cX}:\mathbb{R} \rightarrow \mathbb{C} \)に関して、任意の\(t\in \mathbb{R} \)について以下の関係\begin{equation*}\varphi _{cX}\left( t\right) =\varphi _{X}\left( ct\right)
\end{equation*}が成り立ちます。
\end{equation*}が成り立つ。
連続型確率変数の1次スケーリングの特性関数
連続型の確率変数\(X:\Omega\rightarrow \mathbb{R} \)と実数\(c,d\in \mathbb{R} \)が与えられたとき、それぞれの\(\omega \in \Omega \)に対して、\begin{equation*}\left( cX+d\right) \left( \omega \right) =cX\left( \omega \right) +d
\end{equation*}を定める新たな関数\begin{equation*}
cX+d:\Omega \rightarrow \mathbb{R} \end{equation*}が定義可能です。確率変数の定数倍は確率変数であり、定数関数は確率変数であり、確率変数どうしの和は確率変数であるため\(cX+d\)は確率変数です。さらに、\(X\)の特性関数\(\varphi _{X}:\mathbb{R} \rightarrow \mathbb{C} \)と\(cX+d\)の特性関数\(\varphi _{cX+d}:\mathbb{R} \rightarrow \mathbb{C} \)に関して、任意の\(t\in \mathbb{R} \)について以下の関係\begin{equation*}\varphi _{cX+d}\left( t\right) =e^{idt}\cdot \varphi _{X}\left( ct\right)
\end{equation*}が成り立ちます。
\end{equation*}が成り立つ。
演習問題
\begin{array}{cl}
\frac{1}{b-a} & \left( if\ x\in X\left( \Omega \right) \right) \\
0 & \left( otherwise\right)
\end{array}\right.
\end{equation*}を定めるものとします。\(X\)の特性関数を求めてください。
プレミアム会員専用コンテンツです
【ログイン】【会員登録】