跳到主要内容

法政大学 理工学研究科 システム理工学専攻 経営システム系 2025年8月実施 概率统计

Author

思齐塾, 祭音Myyura

Description

ある工場で生産している製品はまず2個のモジュールA, Bを別々に製造した後、これら2個を組み立てて最終的に製品が完成する。モジュールAとBを製造するためにかかる時間はそれぞれ XA,XBX_A, X_B (単位: 時間)であり、これらの確率密度関数は

fXA(xA)=αeαxA(xA>0),fXB(xB)=βeβxB(xB>0)f_{X_A}(x_A) = \alpha e^{-\alpha x_A} \quad (x_A > 0), \quad f_{X_B}(x_B) = \beta e^{-\beta x_B} \quad (x_B > 0)

で与えられている。また、製造されたモジュールAとBを組み立てて製品を完成させるために要する時間 YY の確率密度関数は

gY(y)=γ2yeγy(y>0)g_Y(y) = \gamma^2 y e^{-\gamma y} \quad (y > 0)

である。 XA,XB,YX_A, X_B, Y は全て独立で、2個のモジュールの製造を同時に開始し、両方の製造が終了するまで組み立てを開始することはできない。モジュールの製造が終了したらすぐに組み立てを開始し、モジュールの製造完了から組み立て開始までの間の時間はゼロと見なすことができる。このとき、以下の問いに答えよ。

(1) モジュールAとBの2個の製造を同時に開始してから2個の製造が完了するまでに要する時間 WW の確率分布を必要な文字を用いて表せ。

(2) E[XA]=E[XB]=E[Y]=2E[X_A] = E[X_B] = E[Y] = 2 のとき、 α,β,γ\alpha, \beta, \gamma の値を求めよ。

(3) (2)で求めた α,β,γ\alpha, \beta, \gamma の下、2個のモジュールの製造を同時に開始してからこれらを組み立てて製品が完成するまでの時間 TT の確率密度関数を示せ。

题目描述

某工厂生产一种产品:先分别制造模块 A、B,再将两个模块组装成最终产品。制造模块 A、B 所需时间分别为 XA,XBX_A,X_B(单位:小时),其概率密度函数为

fXA(xA)=αeαxA(xA>0),f_{X_A}(x_A)=\alpha e^{-\alpha x_A}\qquad(x_A>0),
fXB(xB)=βeβxB(xB>0).f_{X_B}(x_B)=\beta e^{-\beta x_B}\qquad(x_B>0).

将制成的模块 A、B 组装成产品所需时间 YY 的概率密度函数为

gY(y)=γ2yeγy(y>0).g_Y(y)=\gamma^2y e^{-\gamma y}\qquad(y>0).

XA,XB,YX_A,X_B,Y 相互独立。两个模块同时开始制造,只有二者都制造完毕后才能开始组装;两个模块均完成后立即组装,制造完成与开始组装之间的等待时间视为零。回答下列问题。

(1)设 WW 为从两个模块同时开始制造到两个模块均制造完成所需的时间,用必要的参数表示 WW 的概率分布。

(2)当

E[XA]=E[XB]=E[Y]=2E[X_A]=E[X_B]=E[Y]=2

时,求 α,β,γ\alpha,\beta,\gamma

(3)使用(2)求得的 α,β,γ\alpha,\beta,\gamma,设 TT 为从两个模块同时开始制造到组装完成、最终产品制成所需的总时间,求 TT 的概率密度函数。

Kai

以下の (1) の分布関数・密度は w0w\le0 では 00 とする。また、(3) の密度の計算は t>0t>0 に対するものであり、 t0t\le0 では fT(t)=0f_T(t)=0 である。

(1) W=max(XA,XB)W = \max(X_A, X_B) である。 WwW \le w となるのは XAwX_A \le w かつ XBwX_B \le w のときであるから、 WW の累積分布関数は

FW(w)=P(Ww)=P(XAw,XBw)=P(XAw)P(XBw)(XA,XB は独立)=(1eαw)(1eβw)(w>0)\begin{aligned} F_W(w) &= P(W \le w) = P(X_A \le w, X_B \le w) \\ &= P(X_A \le w) P(X_B \le w) \quad (X_A, X_B \text{ は独立}) \\ &= (1 - e^{-\alpha w})(1 - e^{-\beta w}) \quad (w > 0) \end{aligned}

よって、 WW の確率密度関数は

fW(w)=ddwFW(w)=αeαw(1eβw)+βeβw(1eαw)=αeαw+βeβw(α+β)e(α+β)w(w>0)f_W(w) = \frac{d}{dw} F_W(w) = \alpha e^{-\alpha w}(1 - e^{-\beta w}) + \beta e^{-\beta w}(1 - e^{-\alpha w}) = \alpha e^{-\alpha w} + \beta e^{-\beta w} - (\alpha + \beta)e^{-(\alpha + \beta)w} \quad (w > 0)

(2) E[XA]=0xAαeαxAdxA=1α=2E[X_A] = \int_0^{\infty} x_A \alpha e^{-\alpha x_A} dx_A = \frac{1}{\alpha} = 2 より、 α=12\alpha = \frac{1}{2} である。

同様に、 E[XB]=1β=2E[X_B] = \frac{1}{\beta} = 2 より、 β=12\beta = \frac{1}{2} である。

E[Y]=0yγ2yeγydy=2γ=2E[Y] = \int_0^{\infty} y \gamma^2 y e^{-\gamma y} dy = \frac{2}{\gamma} = 2 より、 γ=1\gamma = 1 である。

(3) T=W+YT = W + Y である。 WWYY は独立であるから、 TT の確率密度関数は

fT(t)=0tfW(w)gY(tw)dwf_T(t) = \int_0^t f_W(w) g_Y(t - w) dw

ここで、 α=β=12,γ=1\alpha = \beta = \frac{1}{2}, \gamma = 1 より、

fW(w)=12e12w+12e12wew=e12wew(w>0)f_W(w) = \frac{1}{2} e^{-\frac{1}{2}w} + \frac{1}{2} e^{-\frac{1}{2}w} - e^{-w} = e^{-\frac{1}{2}w} - e^{-w} \quad (w > 0)
gY(y)=yey(y>0)g_Y(y) = y e^{-y} \quad (y > 0)

よって、

fT(t)=0t(e12wew)(tw)e(tw)dw=0te12w(tw)e(tw)dw0tew(tw)e(tw)dw=et0t(tw)e12wdwet0t(tw)dw=et[2(tw)e12w0t+20te12wdw]et[tw12w20t]=et[2te0+4e12w0t]et[t212t2]=et[2t+4(e12t1)]et[12t2]=(12t22t4)et+4e12t\begin{aligned} f_T(t) &= \int_0^t (e^{-\frac{1}{2}w} - e^{-w}) (t - w) e^{-(t - w)} dw \\ &= \int_0^t e^{-\frac{1}{2}w} (t - w) e^{-(t - w)} dw - \int_0^t e^{-w} (t - w) e^{-(t - w)} dw \\ &= e^{-t} \int_0^t (t - w) e^{\frac{1}{2}w} dw - e^{-t} \int_0^t (t - w) dw \\ &= e^{-t} \left[ 2(t-w)e^{\frac{1}{2}w} \Big|_0^t + 2 \int_0^t e^{\frac{1}{2}w} dw \right] - e^{-t} \left[ tw - \frac{1}{2}w^2 \Big|_0^t \right] \\ &= e^{-t} \left[ -2te^0 + 4e^{\frac{1}{2}w} \Big|_0^t \right] - e^{-t} \left[ t^2 - \frac{1}{2}t^2 \right] \\ &= e^{-t} \left[ -2t + 4(e^{\frac{1}{2}t} - 1) \right] - e^{-t} \left[ \frac{1}{2}t^2 \right] \\ &= (- \frac{1}{2}t^2 - 2t - 4)e^{-t} + 4e^{-\frac{1}{2}t} \end{aligned}

(tw)ew/2dw=(tw)2ew/2(1)2ew/2dw=2(tw)ew/2+4ew/2\int (t-w)e^{w/2}dw = (t-w)2e^{w/2} - \int (-1)2e^{w/2}dw = 2(t-w)e^{w/2} + 4e^{w/2} fT(t)=et(2t+4et/24)et(t2/2)f_T(t) = e^{-t}(-2t+4e^{t/2}-4) - e^{-t}(t^2/2) fT(t)=et(4et/2(t2/2)2t4)f_T(t) = e^{-t}(4e^{t/2}-(t^2/2)-2t-4)