跳到主要内容

電気通信大学 情報理工学研究科 情報学専攻 2022年8月実施 選択問題 確率・オペレーションズリサーチ

Author​

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

Description​

この科目では問1を必ず解き、問2と問3のいずれか一方を選ぶ。以下では両方を扱う。

問1​

XX の密度を

fX(x)=(1−x2)I(0≤x<2)f_X(x)=\left(1-\frac{x}{2}\right)I(0\le x<2)

とする。平均と分散を求めよ。また、XX と独立で P(Y=1)=P(Y=−1)=1/2P(Y=1)=P(Y=-1)=1/2 である YY に対し、Z=XYZ=XY の密度と分散を求めよ。

問2​

  1. 独立な Xi∼Poisson⁡(λ)X_i\sim\operatorname{Poisson}(\lambda) の和 Sn=∑i=1nXiS_n=\sum_{i=1}^nX_i の確率関数を求め、λ=5,n=20\lambda=5,n=20 のとき P(S20>120)P(S_{20}>120) を求めよ。また、SnS_n が初めて正の整数 KK を超える時刻 NN の確率関数を求めよ。
  2. 密度をもつ正の確率変数 XX の平均を μ\mu とする。X>μ+δX>\mu+\delta(δ>0\delta>0)なら a>1a>1 枚、そうでなければ 00 枚のメダルが返るゲームを考え、各回は独立とする。
    1. P(X≥μ+δ)P(X\ge\mu+\delta) の上限を求めよ。
    2. μ=δ=2\mu=\delta=2 のとき、10 回で返るメダル枚数の期待値の上限を求めよ。
    3. さらに V[X]=σ2V[X]=\sigma^2 が既知の場合の確率の上限を求めよ。
    4. μ=δ=2,σ=1\mu=\delta=2,\sigma=1 のとき、10 回で返る枚数の期待値の上限を求めよ。

問3​

次の線形計画問題を Simplex 法で解け。

maximize5x1+9x2+7x3,subject tox1+3x2+2x3≤10,3x1+4x2+2x3≤12,2x1+x2+2x3≤8,x1,x2,x3≥0.\begin{array}{ll} \text{maximize}&5x_1+9x_2+7x_3,\\ \text{subject to} &x_1+3x_2+2x_3\le10,\\ &3x_1+4x_2+2x_3\le12,\\ &2x_1+x_2+2x_3\le8,\\ &x_1,x_2,x_3\ge0. \end{array}

スラック変数を x4,x5,x6x_4,x_5,x_6 とし、初期実行基底解、0 回目・1 回目の表、2 回目の実行基底解を求め、最適になるまで計算を続けよ。

题目描述​

必答部分考查连续分布及随机变量乘积;选答部分分别考查 Poisson 分布、停止时刻与概率不等式,或线性规划的单纯形法。这里给出全部答案。

Kai​

問1​

1-1, 1-2​

E[X]=∫02x(1−x2)dx=23,E[X2]=∫02x2(1−x2)dx=23.\begin{aligned} E[X] &=\int_0^2x\left(1-\frac{x}{2}\right)dx =\frac23,\\ E[X^2] &=\int_0^2x^2\left(1-\frac{x}{2}\right)dx =\frac23. \end{aligned}

したがって、

E[X]=23,V[X]=23−49=29.\boxed{E[X]=\frac23},\qquad \boxed{V[X]=\frac23-\frac49=\frac29}.

1-3​

z>0z>0 では Y=1Y=1、z<0z<0 では Y=−1Y=-1 の場合を用いれば、

fZ(z)={12−∣z∣4,∣z∣<2,0,∣z∣≥2.\boxed{ f_Z(z)= \begin{cases} \dfrac12-\dfrac{|z|}{4},&|z|<2,\\[2mm] 0,&|z|\ge2. \end{cases}}

1-4​

E[Y]=0E[Y]=0, Y2=1Y^2=1 より

E[Z]=0,E[Z2]=E[X2]E[Y2]=23.E[Z]=0,\qquad E[Z^2]=E[X^2]E[Y^2]=\frac23.

ゆえに、

V[Z]=23.\boxed{V[Z]=\frac23}.

問2 A​

2-1​

Poisson 分布の再生性より Sn∼Poisson⁡(nλ)S_n\sim\operatorname{Poisson}(n\lambda) である。したがって、

pn(s)=P(Sn=s)=e−nλ(nλ)ss!(s=0,1,…).\boxed{ p_n(s)=P(S_n=s) =e^{-n\lambda}\frac{(n\lambda)^s}{s!} \quad(s=0,1,\ldots)}.

2-2​

S20∼Poisson⁡(100)S_{20}\sim\operatorname{Poisson}(100) なので、

P(S20>120)=1−∑s=0120e−100100ss!≃0.02267.\boxed{ P(S_{20}>120) =1-\sum_{s=0}^{120}e^{-100}\frac{100^s}{s!} \simeq0.02267}.

標準正規近似では

P(S20>120)≃1−Φ(2)=0.0228P(S_{20}>120)\simeq1-\Phi(2)=0.0228

となる。

2-3​

原卷 PDF 5 ページ 。KK は正の整数とする。S0=0S_0=0、p0(s)=1{s=0}p_0(s)=\mathbf1_{\{s=0\}} と定めれば、N=nN=n は Sn−1≤K<SnS_{n-1}\le K<S_n と同値である。したがって n≥1n\ge1 に対し、

P(N=n)=P(Sn−1≤K)−P(Sn≤K)=∑s=0K{pn−1(s)−pn(s)}.\begin{aligned} P(N=n) &=P(S_{n-1}\le K)-P(S_n\le K)\\ &=\boxed{\sum_{s=0}^{K}\{p_{n-1}(s)-p_n(s)\}}. \end{aligned}

問2 B​

2-4​

Markov の不等式より、

P(X≥μ+δ)≤μμ+δ.\boxed{ P(X\ge\mu+\delta)\le\frac{\mu}{\mu+\delta}}.

2-5​

μ=δ=2\mu=\delta=2 では一回の当選確率は 1/21/2 以下である。よって、10 回で返るメダル数の期待値は

10a⋅12=5a 以下.\boxed{10a\cdot\frac12=5a\text{ 以下}}.

2-6​

Chebyshev の不等式より、

P(X≥μ+δ)≤P(∣X−μ∣≥δ)≤σ2δ2.\begin{aligned} P(X\ge\mu+\delta) &\le P(|X-\mu|\ge\delta)\\ &\le\boxed{\frac{\sigma^2}{\delta^2}}. \end{aligned}

2-7​

μ=δ=2,σ=1\mu=\delta=2,\sigma=1 では当選確率は 1/41/4 以下である。したがって、

10a⋅14=5a2 以下.\boxed{10a\cdot\frac14=\frac{5a}{2}\text{ 以下}}.

問3​

3-1​

スラック変数 x4,x5,x6x_4,x_5,x_6 を導入する。最初の実行基底解は

(x1,x2,x3,x4,x5,x6)=(0,0,0,10,12,8).\boxed{(x_1,x_2,x_3,x_4,x_5,x_6)=(0,0,0,10,12,8)}.

3-2​

0 回目の Simplex 表は次のとおりである。

基底x1x_1x2x_2x3x_3x4x_4x5x_5x6x_6定数項
x4x_413210010
x5x_534201012
x6x_62120018
cj−zjc_j-z_j5970000

x2x_2 を入れ、最小比 12/4=312/4=3 を与える x5x_5 を出す。

3-3​

1 回目の表は

基底x1x_1x2x_2x3x_3x4x_4x5x_5x6x_6定数項
x4x_4−5/4-5/401/21/21−3/4-3/401
x2x_23/43/411/21/201/41/403
x6x_65/45/403/23/20−1/4-1/415
cj−zjc_j-z_j−7/4-7/405/25/20−9/4-9/4027

である。2 回目では x3x_3 が入り x4x_4 が出る。更新された実行基底解は

(x1,x2,x3,x4,x5,x6)=(0,2,2,0,0,2).\boxed{(x_1,x_2,x_3,x_4,x_5,x_6)=(0,2,2,0,0,2)}.

3-4​

2 回目の表は

基底x1x_1x2x_2x3x_3x4x_4x5x_5x6x_6定数項
x3x_3−5/2-5/2012−3/2-3/202
x2x_2210−1-1102
x6x_6500−3-3212
cj−zjc_j-z_j9/29/200−5-53/23/2032

である。x1x_1 を入れ、x6x_6 を出すと、3 回目の表は

基底x1x_1x2x_2x3x_3x4x_4x5x_5x6x_6定数項
x3x_30011/21/2−1/2-1/21/21/23
x2x_20101/51/51/51/5−2/5-2/56/56/5
x1x_1100−3/5-3/52/52/51/51/52/52/5
cj−zjc_j-z_j000−23/10-23/10−3/10-3/10−9/10-9/10169/5169/5

となる。すべての被約費用が非正なので最適であり、

(x1,x2,x3)=(25,65,3),zmax⁡=1695.\boxed{ (x_1,x_2,x_3)=\left(\frac25,\frac65,3\right), \qquad z_{\max}=\frac{169}{5}}.