跳到主要内容

東京工業大学 工学院 情報通信系 2019年8月実施 S1 正規分布のモーメント

Author

祭音Myyura (co-authored with GPT 5.6 SOL)

Description

k1k\ge1 を整数,mm を実定数,σ>0\sigma>0 とし,XX の確率密度を

pX(x)=12πσ2exp[(xm)22σ2]p_X(x)=\frac1{\sqrt{2\pi\sigma^2}}\exp\left[-\frac{(x-m)^2}{2\sigma^2}\right]

とする。α>0\alpha>0 に対する公式 eαx2dx=π/α\int_{-\infty}^\infty e^{-\alpha x^2}dx=\sqrt{\pi/\alpha} を用いてよい。

  1. pX(x)dx\int_{-\infty}^\infty p_X(x)dx を求めよ。
  2. E[(Xm)2k1]E[(X-m)^{2k-1}] を求めよ。
  3. E[X]E[X] を求めよ。
  4. Gaussian 積分を α\alphakk 回微分して x2keαx2dx\int_{-\infty}^\infty x^{2k}e^{-\alpha x^2}dx の公式を示せ。
  5. E[(Xm)2k]E[(X-m)^{2k}] を求めよ。
  6. E[X2]E[X^2] を求めよ。
  7. E[X3]E[X^3] を求めよ。

题目描述

由高斯积分及参数微分计算正态分布的归一化、各阶中心矩,以及一至三阶原点矩。

Kai

1)

u=xmu=x-m と置換し Gaussian 積分を使うと

pX(x)dx=12πσ2π1/(2σ2)=1.\int_{-\infty}^\infty p_X(x)dx =\frac1{\sqrt{2\pi\sigma^2}}\sqrt{\frac\pi{1/(2\sigma^2)}}=\boxed1.

2),3)

u2k1eu2/(2σ2)u^{2k-1}e^{-u^2/(2\sigma^2)} は奇関数なので E[(Xm)2k1]=0\boxed{E[(X-m)^{2k-1}]=0}。特に k=1k=1 より E[X]=m\boxed{E[X]=m}

4)

Gaussian の減衰により積分内での微分が許され,

(1)kx2keαx2dx=dkdαk(πα1/2)=(1)k(2k1)!!2kπαk1/2.(-1)^k\int_{-\infty}^\infty x^{2k}e^{-\alpha x^2}dx =\frac{d^k}{d\alpha^k}(\sqrt\pi\alpha^{-1/2}) =(-1)^k\frac{(2k-1)!!}{2^k}\sqrt\pi\alpha^{-k-1/2}.

したがって

x2keαx2dx=(2k1)!!π2kαk+1/2.\boxed{\int_{-\infty}^\infty x^{2k}e^{-\alpha x^2}dx =\frac{(2k-1)!!\sqrt\pi}{2^k\alpha^{k+1/2}}}.

5)

  1. α=1/(2σ2)\alpha=1/(2\sigma^2) を代入して
E[(Xm)2k]=(2k1)!!σ2k.\boxed{E[(X-m)^{2k}]=(2k-1)!!\,\sigma^{2k}}.

6),7)

Y=XmY=X-m とおけば E[Y]=E[Y3]=0E[Y]=E[Y^3]=0, E[Y2]=σ2E[Y^2]=\sigma^2。よって

E[X2]=m2+σ2,E[X3]=m3+3mσ2.\boxed{E[X^2]=m^2+\sigma^2},\qquad \boxed{E[X^3]=m^3+3m\sigma^2}.