東京工業大学 工学院 情報通信系 2019年8月実施 午後 S1
Author
思齐塾, 祭音Myyura (co-authored with GPT 5.6 SOL)
Description
問題 S1
k を 1 以上の整数とし、ガウス分布に従う実数の確率変数 X を考える。すなわち、 x を X の実現値とすると、 X の確率密度関数 pX(x) は
pX(x)=2πσ21exp[−2σ2(x−m)2](S1.1)
となる。ただし、 m は実数の定数で平均値を表し、 σ は正の定数で標準偏差を表すものとする。
以下の問では次の積分を用いてもよい。
∫−∞∞exp(−αx2)dx=απ(S1.2)
なお、 α は任意の正の実数である。
-
∫−∞∞pX(x)dx の値を求めよ。
-
⟨⋅⟩ は集合平均を表し、確率密度関数を用いた平均操作とする。次式で表される ⟨(X−m)2k−1⟩ の値を求めよ。
⟨(X−m)2k−1⟩=∫−∞∞(x−m)2k−1pX(x)dx(S1.3)
-
式 (S1.3) を用いて、 ⟨X⟩ の値を求めよ。
-
式 (S1.2) の両辺を α について k 回微分して
∫−∞∞x2kexp(−αx2)dx=2k∏l=1k(2l−1)π21α−22k+1(S1.4)
を証明せよ。
-
式 (S1.4) を用いて、 ⟨(X−m)2k⟩ の値を求めよ。
-
- の結果を用いて、 ⟨X2⟩ の値を求めよ。
-
⟨X3⟩ の値を求めよ。
题目描述
设 k≥1 为整数,实随机变量 X 服从 Gaussian 分布,密度为
pX(x)=2πσ21exp[−2σ2(x−m)2],(S1.1)
其中 m∈R 是均值参数,σ>0 是标准差。可以使用对任意 α>0 成立的积分
∫−∞∞e−αx2dx=απ.(S1.2)
符号 ⟨⋅⟩ 表示按该概率密度进行的总体平均。
-
求 ∫−∞∞pX(x)dx。
-
计算奇数阶中心矩
⟨(X−m)2k−1⟩=∫−∞∞(x−m)2k−1pX(x)dx.(S1.3)
-
利用 (S1.3) 求 ⟨X⟩。
-
将 (S1.2) 两边关于 α 求导 k 次,证明
∫−∞∞x2ke−αx2dx=2k∏l=1k(2l−1)π1/2α−(2k+1)/2.(S1.4)
-
利用 (S1.4) 求偶数阶中心矩 ⟨(X−m)2k⟩。
-
利用第 5 问求 ⟨X2⟩。
-
求 ⟨X3⟩。
Kai
u=x−m と置換し Gaussian 積分を使うと
∫−∞∞pX(x)dx=2πσ211/(2σ2)π=1.
2),3)
u2k−1e−u2/(2σ2) は奇関数なので E[(X−m)2k−1]=0。特に k=1 より E[X]=m。
Gaussian の減衰により積分内での微分が許され,
(−1)k∫−∞∞x2ke−αx2dx=dαkdk(πα−1/2)=(−1)k2k(2k−1)!!πα−k−1/2.
したがって
∫−∞∞x2ke−αx2dx=2kαk+1/2(2k−1)!!π.
- に α=1/(2σ2) を代入して
E[(X−m)2k]=(2k−1)!!σ2k.
6),7)
Y=X−m とおけば E[Y]=E[Y3]=0, E[Y2]=σ2。よって
E[X2]=m2+σ2,E[X3]=m3+3mσ2.