跳到主要内容

東京工業大学 工学院 情報通信系 2020年8月実施 概率统计

Author

思齐塾, 祭音Myyura

Description

問題 S1

kk を 1 以上の整数とし、ガウス分布に従う実数の確率変数 XX を考える。すなわち、 xxXX の実現値とすると、 XX の確率密度関数 pX(x)p_X(x)

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

となる。ただし、 mm は実数の定数で平均値を表し、 σ\sigma は正の定数で標準偏差を表すものとする。

以下の問では次の積分を用いてもよい。

exp(αx2)dx=πα(S1.2)\int_{-\infty}^{\infty} \exp(-\alpha x^2) dx = \sqrt{\frac{\pi}{\alpha}} \tag{S1.2}

なお、 α\alpha は任意の正の実数である。


  1. pX(x)dx\int_{-\infty}^{\infty} p_X(x) dx の値を求めよ。

  2. \langle \cdot \rangle は集合平均を表し、確率密度関数を用いた平均操作とする。次式で表される (Xm)2k1\langle (X - m)^{2k-1} \rangle の値を求めよ。

(Xm)2k1=(xm)2k1pX(x)dx(S1.3)\langle (X - m)^{2k-1} \rangle = \int_{-\infty}^{\infty} (x - m)^{2k-1} p_X(x) dx \tag{S1.3}
  1. 式 (S1.3) を用いて、 X\langle X \rangle の値を求めよ。

  2. 式 (S1.2) の両辺を α\alpha について kk 回微分して

x2kexp(αx2)dx=l=1k(2l1)2kπ12α2k+12(S1.4)\int_{-\infty}^{\infty} x^{2k} \exp(-\alpha x^2) dx = \frac{\prod_{l=1}^{k} (2l - 1)}{2^k} \pi^{\frac{1}{2}} \alpha^{-\frac{2k+1}{2}} \tag{S1.4}

を証明せよ。

  1. 式 (S1.4) を用いて、 (Xm)2k\langle (X - m)^{2k} \rangle の値を求めよ。

    1. の結果を用いて、 X2\langle X^2 \rangle の値を求めよ。
  2. X3\langle X^3 \rangle の値を求めよ。

题目描述

k1k\geq1 为整数,实随机变量 XX 服从 Gaussian 分布,密度为

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

其中 mRm\in\mathbb R 是均值参数,σ>0\sigma>0 是标准差。可以使用对任意 α>0\alpha>0 成立的积分

eαx2dx=πα.(S1.2)\int_{-\infty}^{\infty}e^{-\alpha x^2}\,dx =\sqrt{\frac\pi\alpha}. \tag{S1.2}

符号 \langle\cdot\rangle 表示按该概率密度进行的总体平均。

  1. pX(x)dx\int_{-\infty}^{\infty}p_X(x)\,dx
  2. 计算奇数阶中心矩
(Xm)2k1=(xm)2k1pX(x)dx.(S1.3)\left\langle(X-m)^{2k-1}\right\rangle =\int_{-\infty}^{\infty} (x-m)^{2k-1}p_X(x)\,dx. \tag{S1.3}
  1. 利用 (S1.3) 求 X\langle X\rangle
  2. 将 (S1.2) 两边关于 α\alpha 求导 kk 次,证明
x2keαx2dx=l=1k(2l1)2kπ1/2α(2k+1)/2.(S1.4)\int_{-\infty}^{\infty}x^{2k}e^{-\alpha x^2}\,dx = \frac{\prod_{l=1}^k(2l-1)}{2^k} \pi^{1/2}\alpha^{-(2k+1)/2}. \tag{S1.4}
  1. 利用 (S1.4) 求偶数阶中心矩 (Xm)2k\langle(X-m)^{2k}\rangle
  2. 利用第 5 问求 X2\langle X^2\rangle
  3. X3\langle X^3\rangle

Kai

問題 S1

1)

u=(xm)/(2σ)u=(x-m)/(\sqrt2\sigma) と変数変換し,式 (S1.2) を用いると

pX(x)dx=1πeu2du=1.\int_{-\infty}^{\infty}p_X(x)\,dx =\frac1{\sqrt\pi}\int_{-\infty}^{\infty}e^{-u^2}\,du =\boxed{1}.

2)

u=xmu=x-m とおけば, u2k1eu2/(2σ2)u^{2k-1}e^{-u^2/(2\sigma^2)} は奇関数であり,積分区間は原点について対称である。よって

(Xm)2k1=0.\boxed{\langle(X-m)^{2k-1}\rangle=0}.

3)

  1. k=1k=1 とすると
0=Xm=Xm,0=\langle X-m\rangle=\langle X\rangle-m,

したがって X=m\boxed{\langle X\rangle=m} である。

4)

式 (S1.2) の左辺を積分記号下で kk 回微分すると

dkdαkeαx2dx=(1)kx2keαx2dx.\frac{d^k}{d\alpha^k}\int_{-\infty}^{\infty}e^{-\alpha x^2}\,dx =(-1)^k\int_{-\infty}^{\infty}x^{2k}e^{-\alpha x^2}\,dx.

一方,右辺について

dkdαk(πα1/2)=(1)kl=1k(2l1)2kπα(2k+1)/2.\frac{d^k}{d\alpha^k}\left(\sqrt\pi\,\alpha^{-1/2}\right) =(-1)^k\frac{\prod_{l=1}^{k}(2l-1)}{2^k}\sqrt\pi\, \alpha^{-(2k+1)/2}.

両辺の (1)k(-1)^k を消去すれば

x2keαx2dx=l=1k(2l1)2kπ1/2α(2k+1)/2\boxed{ \int_{-\infty}^{\infty}x^{2k}e^{-\alpha x^2}\,dx =\frac{\prod_{l=1}^{k}(2l-1)}{2^k}\pi^{1/2}\alpha^{-(2k+1)/2}}

となり,式 (S1.4) が示された。

5)

u=xmu=x-m , α=1/(2σ2)\alpha=1/(2\sigma^2) として式 (S1.4) を使うと

(Xm)2k=12πσ2u2keu2/(2σ2)du=(l=1k(2l1))σ2k.\begin{aligned} \langle(X-m)^{2k}\rangle &=\frac1{\sqrt{2\pi\sigma^2}} \int_{-\infty}^{\infty}u^{2k}e^{-u^2/(2\sigma^2)}\,du\\ &=\left(\prod_{l=1}^{k}(2l-1)\right)\sigma^{2k}. \end{aligned}

したがって

(Xm)2k=(2k1)!!σ2k.\boxed{\langle(X-m)^{2k}\rangle=(2k-1)!!\,\sigma^{2k}}.

6)

  1. k=1k=1 とすれば (Xm)2=σ2\langle(X-m)^2\rangle=\sigma^2 である。3) の X=m\langle X\rangle=m を用いて展開すると
σ2=X22mX+m2=X2m2.\sigma^2=\langle X^2\rangle-2m\langle X\rangle+m^2 =\langle X^2\rangle-m^2.

ゆえに

X2=m2+σ2.\boxed{\langle X^2\rangle=m^2+\sigma^2}.

7)

  1. k=2k=2 とすると (Xm)3=0\langle(X-m)^3\rangle=0 である。展開し,3), 6) の結果を代入すると
0=X33m(m2+σ2)+3m2mm3.0=\langle X^3\rangle-3m(m^2+\sigma^2)+3m^2m-m^3.

したがって

X3=m3+3mσ2.\boxed{\langle X^3\rangle=m^3+3m\sigma^2}.