東京工業大学 工学院 情報通信系 2020年8月実施 概率统计
Author
思齐塾, 祭音Myyura
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
問題 S1
u=(x−m)/(2σ) と変数変換し,式 (S1.2) を用いると
∫−∞∞pX(x)dx=π1∫−∞∞e−u2du=1.
u=x−m とおけば, u2k−1e−u2/(2σ2) は奇関数であり,積分区間は原点について対称である。よって
⟨(X−m)2k−1⟩=0.
- で k=1 とすると
0=⟨X−m⟩=⟨X⟩−m,
したがって ⟨X⟩=m である。
式 (S1.2) の左辺を積分記号下で k 回微分すると
dαkdk∫−∞∞e−αx2dx=(−1)k∫−∞∞x2ke−αx2dx.
一方,右辺について
dαkdk(πα−1/2)=(−1)k2k∏l=1k(2l−1)πα−(2k+1)/2.
両辺の (−1)k を消去すれば
∫−∞∞x2ke−αx2dx=2k∏l=1k(2l−1)π1/2α−(2k+1)/2
となり,式 (S1.4) が示された。
u=x−m , α=1/(2σ2) として式 (S1.4) を使うと
⟨(X−m)2k⟩=2πσ21∫−∞∞u2ke−u2/(2σ2)du=(l=1∏k(2l−1))σ2k.
したがって
⟨(X−m)2k⟩=(2k−1)!!σ2k.
- で k=1 とすれば ⟨(X−m)2⟩=σ2 である。3) の ⟨X⟩=m を用いて展開すると
σ2=⟨X2⟩−2m⟨X⟩+m2=⟨X2⟩−m2.
ゆえに
⟨X2⟩=m2+σ2.
- で k=2 とすると ⟨(X−m)3⟩=0 である。展開し,3), 6) の結果を代入すると
0=⟨X3⟩−3m(m2+σ2)+3m2m−m3.
したがって
⟨X3⟩=m3+3mσ2.