東京大学 新領域創成科学研究科 人間環境学専攻 2023年8月実施 概率统计
Author
思齐塾 , 祭音Myyura
Description
確率変数の変数変換に関する以下の問いに答えよ。
X X X を標準正規分布に従う確率変数とする。ここで、 Y = − 3 X + 2 Y = -3X + 2 Y = − 3 X + 2 で得られる確率変数 Y Y Y を導入し、 Y Y Y が従う確率密度関数を f Y ( y ) f_Y(y) f Y ( y ) とする。このとき、 f Y ( y ) f_Y(y) f Y ( y ) を求めよ。また、 Y Y Y の期待値と分散を求めよ。ただし、確率変数 X X X の従う確率密度関数を f X ( x ) f_X(x) f X ( x ) とすると、
f X ( x ) = 1 2 π e − x 2 2 f_X(x) = \frac{1}{\sqrt{2\pi}} e^{-\frac{x^2}{2}} f X ( x ) = 2 π 1 e − 2 x 2
である。
Y 1 , Y 2 Y_1, Y_2 Y 1 , Y 2 を独立な標準正規分布に従う確率変数とする。 U = 2 Y 1 + Y 2 , V = Y 1 − 2 Y 2 U = 2Y_1 + Y_2, V = Y_1 - 2Y_2 U = 2 Y 1 + Y 2 , V = Y 1 − 2 Y 2 を導入し、2次元確率変数 ( U , V ) (U, V) ( U , V ) の従う2次元同時確率密度関数を g ( u , v ) g(u, v) g ( u , v ) とする。このとき、 g ( u , v ) g(u, v) g ( u , v ) を求めよ。また、確率変数 V V V の期待値と分散を求めよ。
题目描述
回答以下关于随机变量变换的问题:
设随机变量 X X X 服从标准正态分布,其概率密度为
f X ( x ) = 1 2 π e − x 2 / 2 . f_X(x)=\frac{1}{\sqrt{2\pi}}e^{-x^2/2}. f X ( x ) = 2 π 1 e − x 2 /2 .
定义 Y = − 3 X + 2 Y=-3X+2 Y = − 3 X + 2 ,并以 f Y ( y ) f_Y(y) f Y ( y ) 表示 Y Y Y 的概率密度。求 f Y ( y ) f_Y(y) f Y ( y ) ,以及 Y Y Y 的期望和方差。
设 Y 1 , Y 2 Y_1,Y_2 Y 1 , Y 2 是相互独立且均服从标准正态分布的随机变量,定义
U = 2 Y 1 + Y 2 , V = Y 1 − 2 Y 2 . U=2Y_1+Y_2,\qquad V=Y_1-2Y_2. U = 2 Y 1 + Y 2 , V = Y 1 − 2 Y 2 .
以 g ( u , v ) g(u,v) g ( u , v ) 表示二维随机变量 ( U , V ) (U,V) ( U , V ) 的联合概率密度。求 g ( u , v ) g(u,v) g ( u , v ) ,以及随机变量 V V V 的期望和方差。
Kai
Solution:
Since Y = − 3 X + 2 Y = -3X + 2 Y = − 3 X + 2 , we have X = 2 − Y 3 X = \frac{2 - Y}{3} X = 3 2 − Y . The Jacobian of this transformation is given by d X d Y = − 1 3 \frac{dX}{dY} = -\frac{1}{3} d Y d X = − 3 1 . Thus, the probability density function of Y Y Y is given by
f Y ( y ) = f X ( x ) ∣ d X d Y ∣ = 1 2 π e − ( 2 − y 3 ) 2 2 ∣ − 1 3 ∣ = 1 3 2 π e − ( 2 − y ) 2 18 = 1 3 2 π e − ( y − 2 ) 2 18 f_Y(y) = f_X(x) \left| \frac{dX}{dY} \right| = \frac{1}{\sqrt{2\pi}} e^{-\frac{(\frac{2-y}{3})^2}{2}} \left| -\frac{1}{3} \right| = \frac{1}{3\sqrt{2\pi}} e^{-\frac{(2-y)^2}{18}} = \frac{1}{3\sqrt{2\pi}} e^{-\frac{(y-2)^2}{18}} f Y ( y ) = f X ( x ) d Y d X = 2 π 1 e − 2 ( 3 2 − y ) 2 − 3 1 = 3 2 π 1 e − 18 ( 2 − y ) 2 = 3 2 π 1 e − 18 ( y − 2 ) 2
Therefore, Y ∼ N ( 2 , 9 ) Y \sim N(2, 9) Y ∼ N ( 2 , 9 ) . The expected value of Y Y Y is E [ Y ] = 2 E[Y] = 2 E [ Y ] = 2 and the variance of Y Y Y is V a r ( Y ) = 9 Var(Y) = 9 Va r ( Y ) = 9 .
Solution:
Since U = 2 Y 1 + Y 2 U = 2Y_1 + Y_2 U = 2 Y 1 + Y 2 and V = Y 1 − 2 Y 2 V = Y_1 - 2Y_2 V = Y 1 − 2 Y 2 , we can express Y 1 Y_1 Y 1 and Y 2 Y_2 Y 2 in terms of U U U and V V V as follows:
2 V = 2 Y 1 − 4 Y 2 2V = 2Y_1 - 4Y_2 2 V = 2 Y 1 − 4 Y 2 , so U − 2 V = 5 Y 2 U - 2V = 5Y_2 U − 2 V = 5 Y 2 and Y 2 = U − 2 V 5 Y_2 = \frac{U - 2V}{5} Y 2 = 5 U − 2 V .
V = Y 1 − 2 Y 2 V = Y_1 - 2Y_2 V = Y 1 − 2 Y 2 , so Y 1 = V + 2 Y 2 = V + 2 ( U − 2 V 5 ) = 5 V + 2 U − 4 V 5 = 2 U + V 5 Y_1 = V + 2Y_2 = V + 2(\frac{U - 2V}{5}) = \frac{5V + 2U - 4V}{5} = \frac{2U + V}{5} Y 1 = V + 2 Y 2 = V + 2 ( 5 U − 2 V ) = 5 5 V + 2 U − 4 V = 5 2 U + V .
The Jacobian of this transformation is
J = ∣ ∂ Y 1 ∂ U ∂ Y 1 ∂ V ∂ Y 2 ∂ U ∂ Y 2 ∂ V ∣ = ∣ 2 5 1 5 1 5 − 2 5 ∣ = − 4 25 − 1 25 = − 5 25 = − 1 5 J = \begin{vmatrix} \frac{\partial Y_1}{\partial U} & \frac{\partial Y_1}{\partial V} \\ \frac{\partial Y_2}{\partial U} & \frac{\partial Y_2}{\partial V} \end{vmatrix} = \begin{vmatrix} \frac{2}{5} & \frac{1}{5} \\ \frac{1}{5} & -\frac{2}{5} \end{vmatrix} = -\frac{4}{25} - \frac{1}{25} = -\frac{5}{25} = -\frac{1}{5} J = ∂ U ∂ Y 1 ∂ U ∂ Y 2 ∂ V ∂ Y 1 ∂ V ∂ Y 2 = 5 2 5 1 5 1 − 5 2 = − 25 4 − 25 1 = − 25 5 = − 5 1
The joint density of Y 1 Y_1 Y 1 and Y 2 Y_2 Y 2 is given by f Y 1 , Y 2 ( y 1 , y 2 ) = f Y 1 ( y 1 ) f Y 2 ( y 2 ) f_{Y_1, Y_2}(y_1, y_2) = f_{Y_1}(y_1) f_{Y_2}(y_2) f Y 1 , Y 2 ( y 1 , y 2 ) = f Y 1 ( y 1 ) f Y 2 ( y 2 ) since Y 1 Y_1 Y 1 and Y 2 Y_2 Y 2 are independent.
Since Y 1 Y_1 Y 1 and Y 2 Y_2 Y 2 are standard normal, f Y 1 ( y 1 ) = 1 2 π e − y 1 2 2 f_{Y_1}(y_1) = \frac{1}{\sqrt{2\pi}} e^{-\frac{y_1^2}{2}} f Y 1 ( y 1 ) = 2 π 1 e − 2 y 1 2 and f Y 2 ( y 2 ) = 1 2 π e − y 2 2 2 f_{Y_2}(y_2) = \frac{1}{\sqrt{2\pi}} e^{-\frac{y_2^2}{2}} f Y 2 ( y 2 ) = 2 π 1 e − 2 y 2 2 .
Thus, f Y 1 , Y 2 ( y 1 , y 2 ) = 1 2 π e − y 1 2 + y 2 2 2 f_{Y_1, Y_2}(y_1, y_2) = \frac{1}{2\pi} e^{-\frac{y_1^2 + y_2^2}{2}} f Y 1 , Y 2 ( y 1 , y 2 ) = 2 π 1 e − 2 y 1 2 + y 2 2 .
g ( u , v ) = f Y 1 , Y 2 ( y 1 ( u , v ) , y 2 ( u , v ) ) ∣ J ∣ = 1 2 π e − ( 2 u + v 5 ) 2 + ( u − 2 v 5 ) 2 2 ∣ − 1 5 ∣ = 1 10 π e − 4 u 2 + 4 u v + v 2 + u 2 − 4 u v + 4 v 2 50 = 1 10 π e − 5 u 2 + 5 v 2 50 = 1 10 π e − u 2 + v 2 10 g(u, v) = f_{Y_1, Y_2}(y_1(u, v), y_2(u, v)) |J| = \frac{1}{2\pi} e^{-\frac{(\frac{2u+v}{5})^2 + (\frac{u-2v}{5})^2}{2}} \left| -\frac{1}{5} \right| = \frac{1}{10\pi} e^{-\frac{4u^2 + 4uv + v^2 + u^2 - 4uv + 4v^2}{50}} = \frac{1}{10\pi} e^{-\frac{5u^2 + 5v^2}{50}} = \frac{1}{10\pi} e^{-\frac{u^2 + v^2}{10}} g ( u , v ) = f Y 1 , Y 2 ( y 1 ( u , v ) , y 2 ( u , v )) ∣ J ∣ = 2 π 1 e − 2 ( 5 2 u + v ) 2 + ( 5 u − 2 v ) 2 − 5 1 = 10 π 1 e − 50 4 u 2 + 4 uv + v 2 + u 2 − 4 uv + 4 v 2 = 10 π 1 e − 50 5 u 2 + 5 v 2 = 10 π 1 e − 10 u 2 + v 2
To find the expected value and variance of V V V , we need to find the marginal distribution of V V V .
f V ( v ) = ∫ − ∞ ∞ g ( u , v ) d u = ∫ − ∞ ∞ 1 10 π e − u 2 + v 2 10 d u = 1 10 π e − v 2 10 ∫ − ∞ ∞ e − u 2 10 d u = 1 10 π e − v 2 10 10 π = 1 10 π e − v 2 10 f_V(v) = \int_{-\infty}^{\infty} g(u, v) du = \int_{-\infty}^{\infty} \frac{1}{10\pi} e^{-\frac{u^2 + v^2}{10}} du = \frac{1}{10\pi} e^{-\frac{v^2}{10}} \int_{-\infty}^{\infty} e^{-\frac{u^2}{10}} du = \frac{1}{10\pi} e^{-\frac{v^2}{10}} \sqrt{10\pi} = \frac{1}{\sqrt{10\pi}} e^{-\frac{v^2}{10}} f V ( v ) = ∫ − ∞ ∞ g ( u , v ) d u = ∫ − ∞ ∞ 10 π 1 e − 10 u 2 + v 2 d u = 10 π 1 e − 10 v 2 ∫ − ∞ ∞ e − 10 u 2 d u = 10 π 1 e − 10 v 2 10 π = 10 π 1 e − 10 v 2
Thus, V ∼ N ( 0 , 5 ) V \sim N(0, 5) V ∼ N ( 0 , 5 ) , so E [ V ] = 0 E[V] = 0 E [ V ] = 0 and V a r ( V ) = 5 Var(V) = 5 Va r ( V ) = 5 .