跳到主要内容

東京大学 新領域創成科学研究科 複雑理工学専攻 2019年8月実施 専門基礎科目 第4問

Author

Miyake, 祭音Myyura

Description

Independent summary

Source: University of Tokyo, AY2020 examination, Problem 4, pp. 10–11.

All journey components below are independent normal variables. Times are in minutes.

RouteComponentMeanStandard deviation
1Walking152
1Train A254
1Bus A255
2Bicycle254
2Train B158
2Bus B158

Use fX(x)=e(xμ)2/(2σ2)/2πσ2f_X(x)=e^{-(x-\mu)^2/(2\sigma^2)}/\sqrt{2\pi\sigma^2} and MX(t)=E[etX]M_X(t)=E[e^{tX}], with σ>0\sigma>0. The supplied standard-normal upper-tail table is:

x0x_01.001.281.441.651.962.332.58
P(X0x0)P(X_0\geq x_0)16%10%7.5%5.0%2.5%1.0%0.5%
  1. Find α,β\alpha,\beta such that αX+β\alpha X+\beta is standard normal.
  2. For independent XiN(μi,σi2)X_i\sim N(\mu_i,\sigma_i^2), find A,B,CA,B,C in MaX1+bX2(t)=eAt2+Bt+CM_{aX_1+bX_2}(t)=e^{At^2+Bt+C}.
  3. Select the interval containing P(T175)P(T_1\geq75).
  4. Find the departure time on route 1 for arrival by 10 a.m. with probability 99%.
  5. Select the interval containing P(T1T225)P(T_1-T_2\geq25).

Interval choices: 10–16%, 7.5–10%, 5.0–7.5%, 2.5–5.0%, 1.0–2.5%, 0.5–1.0%, 0.0–0.5%. You may use 21.4\sqrt2\simeq1.4, 31.7\sqrt3\simeq1.7, 52.2\sqrt5\simeq2.2, 72.6\sqrt7\simeq2.6.

题目描述

所有行程分量均为相互独立的正态随机变量,时间单位为分钟,标准差为正。

路线分量均值标准差
1步行152
1电车A254
1公交A255
2自行车254
2电车B158
2公交B158
  1. 考虑均值为 μ\mu、标准差为 σ\sigma 的随机变量 XX 的线性标准化:求实数 α,β\alpha,\beta,使 αX+β\alpha X+\beta 的均值为 0、方差为 1,并分别讨论 α>0\alpha>0α<0\alpha<0

  2. 先由积分推导正态变量

    XN(μ,σ2)X\sim N(\mu,\sigma^2)

    的矩母函数,再对相互独立的

    X1N(μ1,σ12),X2N(μ2,σ22)X_1\sim N(\mu_1,\sigma_1^2),\qquad X_2\sim N(\mu_2,\sigma_2^2)

    aX1+bX2aX_1+bX_2 的矩母函数,并将其写成 exp(At2+Bt+C)\exp(At^2+Bt+C) 的形式以确定 A,B,CA,B,C

  3. 给定“通学路线 1”的时间分布为 N(65,45)N(65,45),把超过 75 的概率标准化为标准正态右尾 P(X01.49)P(X_0\ge1.49),并用给定表格界定在 5.0%5.0\%7.5%7.5\% 之间。根据下方所列正态右尾表选择概率区间。

  4. 求通过路线1到达学校的概率为99%、到达时限为上午10时所对应的最晚出发时间。

  5. 使用相互独立的两条通学路线时间

    P1N(65,45),P2N(55,144),P_1\sim N(65,45),\qquad P_2\sim N(55,144),

    从而令 ΔT=P1P2N(10,189)\Delta T=P_1-P_2\sim N(10,189),并估计 P(ΔT>25)P(\Delta T>25)10%10\%16%16\% 之间。从下列概率区间中选择。

标准正态右尾表:

x0x_01.001.281.441.651.962.332.58
P(X0x0)P(X_0\ge x_0)16%10%7.5%5.0%2.5%1.0%0.5%

概率区间选项为10–16%、7.5–10%、5.0–7.5%、2.5–5.0%、1.0–2.5%、0.5–1.0%、0.0–0.5%。可以使用 21.4,31.7,52.2,72.6\sqrt2\simeq1.4,\sqrt3\simeq1.7,\sqrt5\simeq2.2,\sqrt7\simeq2.6

Kai

(問1)

与えられた条件より、

0=αμ+β,    1=α2σ2\begin{aligned} 0 = \alpha \mu + \beta , \ \ \ \ 1 = \alpha^2 \sigma^2 \end{aligned}

であるから、

α>0\alpha \gt 0 とすると、

α=1σ,    β=αμ=μσ\begin{aligned} \alpha = \frac{1}{\sigma} , \ \ \ \ \beta = - \alpha \mu = - \frac{\mu}{\sigma} \end{aligned}

であり、

α<0\alpha \lt 0 とすると、

α=1σ,    β=αμ=μσ\begin{aligned} \alpha = - \frac{1}{\sigma} , \ \ \ \ \beta = - \alpha \mu = \frac{\mu}{\sigma} \end{aligned}

である。

(問2)

まず、 平均 μ\mu 、標準偏差 σ>0\sigma>0 の正規確率変数 XX の モーメント母関数 MX(t)M_X(t) は、次のように求められる:

MX(t)=EX[exp(tX)]=12πσ2exp((xμ)22σ2)exp(tx)dx=exp(μt+σ2t22)2πσ2exp({x(μ+σ2t)}22σ2)dx=exp(μt+σ2t22).\begin{aligned} M_X(t) &= E_X \left[ \exp (tX) \right] \\ &= \frac{1}{\sqrt{2 \pi \sigma^2}} \int_{-\infty}^\infty \exp \left(- \frac{(x-\mu)^2}{2 \sigma^2} \right) \exp \left( tx \right) dx \\ &= \frac{\exp \left( \mu t + \frac{\sigma^2 t^2}{2} \right) } {\sqrt{2 \pi \sigma^2}} \int_{-\infty}^\infty \exp \left(- \frac {\left\{ x - (\mu + \sigma^2 t) \right\}^2 }{2 \sigma^2} \right) dx \\ &= \exp \left( \mu t + \frac{\sigma^2 t^2}{2} \right) . \end{aligned}

そこで、 X1X_1X2X_2 が独立であることを考慮して、

MaX1+bX2(t)=EX1[exp(taX1)]EX2[exp(tbX2)]=exp(μ1ta+σ12t2a22)exp(μ2tb+σ22t2b22)=exp(a2σ12+b2σ222t2+(aμ1+bμ2)t)\begin{aligned} M_{aX_1+bX_2}(t) &= E_{X_1} \left[ \exp (taX_1) \right] E_{X_2} \left[ \exp (tbX_2) \right] \\ &= \exp \left( \mu_1 ta + \frac{\sigma_1^2 t^2 a^2}{2} \right) \exp \left( \mu_2 tb + \frac{\sigma_2^2 t^2 b^2}{2} \right) \\ &= \exp \left( \frac{a^2 \sigma_1^2 + b^2 \sigma_2^2 }{2} t^2 + (a \mu_1 + b \mu_2 ) t \right) \end{aligned}

であるから、

A=a2σ12+b2σ222,    B=aμ1+bμ2,    C=0\begin{aligned} A = \frac{a^2 \sigma_1^2 + b^2 \sigma_2^2 }{2} , \ \ \ \ B = a \mu_1 + b \mu_2 , \ \ \ \ C = 0 \end{aligned}

を得る。

(問3)

通学路1 N(65,45)\sim N(65, 45)

Pr(t6545756545)Pr(X01.49)\Pr \big(\frac{t-65}{\sqrt{45}} \ge \frac{75-65}{\sqrt{45}} \big) \simeq \Pr (X_0 \ge 1.49)
 1.441.491.655.0%Pr(X01.49)7.5%\begin{aligned} &\because \ 1.44 \le 1.49 \le 1.65 \\ &\therefore 5.0\% \le \Pr (X_0 \ge 1.49) \le 7.5\% \end{aligned}

(問4)

出発から午前10時までの時間を LL 分とすると、

Pr(T1L)=0.99,L65452.33.\Pr(T_1\leq L)=0.99,\qquad \frac{L-65}{\sqrt{45}}\simeq2.33.

従って L65+2.334580.63L\simeq65+2.33\sqrt{45}\simeq80.63 分であり、出発時刻は午前8時39分ごろである。分単位で少なくとも99%の到着確率を確保するなら、午前8時39分に出発する。

(問5)

通学路1 P1P_1 N(65,45)\sim N(65, 45) と 通学路2 P2P_2 N(55,144)\sim N(55, 144) より

ΔT=P1P2N(10,189)\Delta T = P_1 - P_2 \sim N(10, 189)

であるから、

Pr(ΔT>25)=Pr(ΔT10189>15189)      1<15189<1.28      10%<Pr(ΔT>25)<16%\begin{aligned} &\Pr (\Delta T > 25) = \Pr \big(\frac{\Delta T-10}{\sqrt{189}} > \frac{15}{\sqrt{189}} \big) \\ &\because \ \ \ \ \ \ 1 < \frac{15}{\sqrt{189}} < 1.28 \\ &\therefore \ \ \ \ \ \ 10\% < \Pr(\Delta T > 25) < 16\% \end{aligned}

を得る。