跳到主要内容

東京工業大学 環境・社会理工学院 土木・環境工学系 2020年8月実施 午前 基礎科目 問題4

Author

思齐塾, 祭音Myyura

Description

  1. 互いに独立なふたつの正規確率変数 Xi(i=1,2)X_i (i = 1, 2) を考え,これらの確率変数の平均および分散をそれぞれ, μXi,σXi2(i=1,2)\mu_{X_i}, \sigma_{X_i}^2 (i = 1, 2) とする。このとき, Λ=α1X1+α2X2\Lambda = \alpha_1 X_1 + \alpha_2 X_2 によって定義される確率変数 Λ\Lambda は,以下のような平均 μΛ\mu_\Lambda および分散 σΛ2\sigma_\Lambda^2 をもつ正規確率変数となる。ただし, αi(i=1,2)\alpha_i (i = 1, 2) はゼロではない任意の実数である。
μΛ=α1μX1+α2μX2(1)\mu_\Lambda = \alpha_1 \mu_{X_1} + \alpha_2 \mu_{X_2} \tag{1}
σΛ2=α12σX12+α22σX22(2)\sigma_\Lambda^2 = \alpha_1^2 \sigma_{X_1}^2 + \alpha_2^2 \sigma_{X_2}^2 \tag{2}

このことを用いて以下の問いに答えよ。

なお,任意の確率変数 XX が平均 μ\mu ,分散 σ2\sigma^2 の正規分布に従うとき, XX の確率密度関数は次式で表わされることを用いてよい。

fX(x)=12πσ2exp[12((xμ)2σ2)](3)f_X(x) = \frac{1}{\sqrt{2\pi\sigma^2}} \exp\left[ -\frac{1}{2}\left( \frac{(x - \mu)^2}{\sigma^2} \right) \right] \tag{3}

また,解答にあたっては次ページの表 2,表 3 に示す近似値を正しい値とみなして計算に用いること。

(1) NN 個 ( NN は任意の自然数) の互いに独立な正規確率変数 Xi(i=1,,N)X_i (i = 1, \ldots, N) を考え,これらの確率変数の平均および分散はそれぞれ, μXi,σXi2(i=1,,N)\mu_{X_i}, \sigma_{X_i}^2 (i = 1, \ldots, N) とする。このとき,次式で表わされる確率変数 ZZ の確率密度関数を求めよ。

Z=i=1NαiXi(4)Z = \sum_{i=1}^{N} \alpha_i X_i \tag{4}

ただし, αi\alpha_i は任意の実数で,少なくともひとつ以上の ii について αi0\alpha_i \neq 0 である。

(2) A 市から B 市へ移動する場合,図 1 に示すように P 町を経由するルート a と Q 村を経由するルート b がある。A → P, P → B, A → Q, Q → B の各区間の所要時間は正規分布に従うと仮定することができて,それらの平均および分散は表 1 に示す値であることがわかっている。

表 1 各区間の所要時間の平均および分散

ルート区間平均 [時間]分散 [時間^2]
aA → P32
P → B74
bA → Q51
Q → B62

このとき,A 市から B 市へ移動する場合について,以下の問いに答えよ。

問 1. ルート a の所要時間が 8 時間 12 分以内となる確率を求めよ。

問 2. ルート b の所要時間が 10 時間 9 分以上,12 時間 42 分以下となる確率を求めよ。

問 3. ルート a のほうがルート b よりも早く B 市に到着する確率を求めよ。

题目描述

考虑两个相互独立的正态随机变量 XiX_ii=1,2i=1,2),其均值和方差分别为 μXi\mu_{X_i}σXi2\sigma_{X_i}^2。对任意非零实数 α1,α2\alpha_1,\alpha_2,随机变量

Λ=α1X1+α2X2\Lambda=\alpha_1X_1+\alpha_2X_2

仍服从正态分布,其均值与方差为

μΛ=α1μX1+α2μX2,(1)\mu_\Lambda =\alpha_1\mu_{X_1}+\alpha_2\mu_{X_2}, \tag{1}
σΛ2=α12σX12+α22σX22.(2)\sigma_\Lambda^2 =\alpha_1^2\sigma_{X_1}^2+\alpha_2^2\sigma_{X_2}^2. \tag{2}

并可使用如下事实:若随机变量 XX 服从均值 μ\mu、方差 σ2\sigma^2 的正态分布,则其概率密度为

fX(x)=12πσ2exp ⁣[12(xμ)2σ2].(3)f_X(x) =\frac{1}{\sqrt{2\pi\sigma^2}} \exp\!\left[ -\frac12\frac{(x-\mu)^2}{\sigma^2} \right]. \tag{3}

作答时,题面要求把下一页表 2、表 3 中给出的近似值视为准确值使用。回答下列问题。

  1. XiX_ii=1,,Ni=1,\ldots,N)是 NN 个相互独立的正态随机变量,其中 NN 为任意自然数,XiX_i 的均值和方差分别为 μXi\mu_{X_i}σXi2\sigma_{X_i}^2。对任意实数 αi\alpha_i,且至少有一个 αi0\alpha_i\neq0,定义

    Z=i=1NαiXi.(4)Z=\sum_{i=1}^{N}\alpha_iX_i. \tag{4}

    求随机变量 ZZ 的概率密度函数。

  2. 从 A 市到 B 市有两条路线。图 1 的含义是:路线 a 依次经过 A \to P、P \to B,路线 b 依次经过 A \to Q、Q \to B。各区段所需时间均可假定服从正态分布,其均值与方差如下。

    路线区段均值(小时)方差(小时2^2
    aA \to P32
    aP \to B74
    bA \to Q51
    bQ \to B62

    回答以下三问。

    1. 求路线 a 的总行程时间不超过 88 小时 1212 分的概率。

    2. 求路线 b 的总行程时间在 1010 小时 99 分以上且在 1212 小时 4242 分以下的概率。

    3. 求走路线 a 比走路线 b 更早到达 B 市的概率。

図表の補足

公式問題PDF、4–5ページ

AからPまたはQを経由してBへ向かう二経路

表2:指定された平方根の近似値

kk123456789101112
k\sqrt{k}11.41.722.22.42.62.833.23.33.5

表3:Φ(s)=P(N(0,1)s)\Phi(s)=P(N(0,1)\le s) の近似値

ssΦ(s)\Phi(s)ssΦ(s)\Phi(s)ssΦ(s)\Phi(s)ssΦ(s)\Phi(s)ssΦ(s)\Phi(s)ssΦ(s)\Phi(s)
0.000.5000.200.5790.400.6550.600.7260.800.7881.000.841
0.010.5040.210.5830.410.6590.610.7290.810.7911.010.844
0.020.5080.220.5870.420.6630.620.7320.820.7941.020.846
0.030.5120.230.5910.430.6660.630.7360.830.7971.030.848
0.040.5160.240.5950.440.6700.640.7390.840.8001.040.851
0.050.5200.250.5990.450.6740.650.7420.850.8021.050.853
0.060.5240.260.6030.460.6770.660.7450.860.8051.060.855
0.070.5280.270.6060.470.6810.670.7490.870.8081.070.858
0.080.5320.280.6100.480.6840.680.7520.880.8111.080.860
0.090.5360.290.6140.490.6880.690.7550.890.8131.090.862
0.100.5400.300.6180.500.6910.700.7580.900.8161.100.864
0.110.5440.310.6220.510.6950.710.7610.910.8191.110.867
0.120.5480.320.6260.520.6980.720.7640.920.8211.120.869
0.130.5520.330.6290.530.7020.730.7670.930.8241.130.871
0.140.5560.340.6330.540.7050.740.7700.940.8261.140.873
0.150.5600.350.6370.550.7090.750.7730.950.8291.150.875
0.160.5640.360.6410.560.7120.760.7760.960.8311.160.877
0.170.5670.370.6440.570.7160.770.7790.970.8341.170.879
0.180.5710.380.6480.580.7190.780.7820.980.8361.180.881
0.190.5750.390.6520.590.7220.790.7850.990.8391.190.883

原題は区間時間の独立性を個別には再記していない。以下では冒頭の独立性を各区間にも適用する。そうでなければ共分散情報が不足する。

Kai

4. 独立な正規確率変数の線形結合

以下では,問題冒頭の独立性を各区間の所要時間にも適用し,四つの区間時間は互いに独立であるという通常の出題意図で計算する。

(1)

二変数について与えられた事実を帰納的に用いると,

Z=i=1NαiXiZ=\sum_{i=1}^{N}\alpha_iX_i

も正規分布に従い,その平均と分散は

μZ=i=1NαiμXi,σZ2=i=1Nαi2σXi2\mu_Z=\sum_{i=1}^{N}\alpha_i\mu_{X_i},\qquad \sigma_Z^2=\sum_{i=1}^{N}\alpha_i^2\sigma_{X_i}^2

である。少なくとも一つの αi\alpha_i が非零なので σZ2>0\sigma_Z^2>0 であり,求める密度は

fZ(z)=12πi=1Nαi2σXi2exp ⁣[(zi=1NαiμXi)22i=1Nαi2σXi2].\boxed{ f_Z(z)=\frac{1}{\sqrt{2\pi\sum_{i=1}^{N}\alpha_i^2\sigma_{X_i}^2}} \exp\!\left[-\frac{\left(z-\sum_{i=1}^{N}\alpha_i\mu_{X_i}\right)^2}{2\sum_{i=1}^{N}\alpha_i^2\sigma_{X_i}^2}\right] }.

(2) 問1

ルート a の総所要時間を TaT_a とすると

TaN(3+7,2+4)=N(10,6).T_a\sim N(3+7,\,2+4)=N(10,6).

8時間12分は 8.28.2 時間であり、表2に従って 6=2.4\sqrt6=2.4 を用いる。標準正規分布の分布関数を Φ\Phi と書けば

P(Ta8.2)=Φ ⁣(8.2106)=Φ(1.8/2.4)=Φ(0.75)=10.773=0.227.P(T_a\le8.2)=\Phi\!\left(\frac{8.2-10}{\sqrt6}\right) =\Phi(-1.8/2.4)=\Phi(-0.75) =1-0.773=\boxed{0.227}.

(2) 問2

ルート b の総所要時間を TbT_b とすると

TbN(5+6,1+2)=N(11,3).T_b\sim N(5+6,\,1+2)=N(11,3).

表2の 3=1.7\sqrt3=1.7 を使用する。10時間9分は 10.1510.15 時間,12時間42分は 12.712.7 時間であるから

P(10.15Tb12.7)=Φ ⁣(12.7113)Φ ⁣(10.15113)=Φ(1)Φ(0.5)=0.841(10.691)=0.532.\begin{aligned} P(10.15\le T_b\le12.7) &=\Phi\!\left(\frac{12.7-11}{\sqrt3}\right)-\Phi\!\left(\frac{10.15-11}{\sqrt3}\right)\\ &=\Phi(1)-\Phi(-0.5)\\ &=0.841-(1-0.691)\\ &=\boxed{0.532}. \end{aligned}

(2) 問3

独立性から差 D=TaTbD=T_a-T_b も正規分布に従い,

DN(1011,6+3)=N(1,9).D\sim N(10-11,\,6+3)=N(-1,9).

したがって

P(Ta<Tb)=P(D<0)=Φ ⁣(0(1)3)=Φ(1/3).P(T_a<T_b)=P(D<0)=\Phi\!\left(\frac{0-(-1)}{3}\right) =\Phi(1/3).

表3では 1/31/3 を小数第2位の 0.330.33 に丸めて参照すると 0.629\boxed{0.629} となる。原題は表の刻みの間の扱いを指定していないため、Φ(0.33)=0.629\Phi(0.33)=0.629Φ(0.34)=0.633\Phi(0.34)=0.633 を線形補間する場合は約 0.63030.6303 である。分布関数の式 Φ(1/3)\Phi(1/3) は共通である。

参考として、指定された表を使わずに高精度で計算すると問1–3はそれぞれ約 0.2312,0.5250,0.63060.2312,0.5250,0.6306 になる。上では「表の近似値を正しい値とみなす」という設問の指定に従った。