随机变量函数的联合概率律

第9节中,我们较为详细地讨论了求随机变量函数的个体概率律的问题。接下来自然要考虑求作为函数产生的多个随机变量的联合概率律的问题。原则上,这个问题与之前考虑的问题并无不同。然而,细节更为复杂。因此,在本节中,我们仅给出一个常用公式,用于求 n 个随机变量 Y 1 , Y 2 , , Y n 的联合概率密度函数,这些随机变量是由 n 个联合连续的随机变量 X 1 , X 2 , , X n 的函数产生的:我们仅考虑函数 g 1 ( x 1 , x 2 , , x n ) g 2 ( x 1 , x 2 , , x n ) , g n ( x 1 , x 2 , , x n ) 在所有点 ( x 1 , x 2 , , x n ) 处具有连续的一阶偏导数,且使得雅可比行列式

J ( x 1 , x 2 , , x n ) = | g 1 x 1 g 1 x 2 g 1 x n g 2 x 1 g 2 x 2 g 2 x n g n x 1 g n x 2 g n x n | 0 在所有点 ( x 1 , x 2 , , x n ) 处成立。设 C 为使得 n 个方程\begin{align} ended with \end{aligned}\begin{align} y_{1}&=g_{1}\left(x_{1}, x_{2}, \ldots, x_{n}\right),\\ y_{2}&=g_{2}\left(x_{1}, x_{2}, \ldots, x_{n}\right),\\ &\vdots\\ y_{n}&=g_{n}\left(x_{1}, x_{2}, \ldots, x_{n}\right) \end{aligned}\tag{10.3}\]至少有一个解 ( x 1 , x 2 , , x n ) 的点 ( y 1 , y 2 , , y n ) 的集合。那么,(10.3)中的方程组恰好有一个解,我们将其记为\begin{align} ended with \end{aligned}\begin{align} x_{1}&=g^{-1}\left(y_{1}, y_{2}, \ldots, y_{n}\right),\\ x_{2}&=g^{-1}\left(y_{1}, y_{2}, \ldots, y_{n}\right),\\ &\vdots\\ x_{n}&=g^{-1}\left(y_{1}, y_{2}, \ldots, y_{n}\right). \end{aligned} \tag{10.4} \]

如果 X 1 , X 2 , , X n 是联合连续的随机变量,其联合概率密度函数在 ( x 1 , x 2 , , x n ) 空间中除有限个点外均连续,那么由(10.1)定义的随机变量 Y 1 , Y 2 , , Y n 是联合连续的,其联合概率密度函数由下式给出 如果 ( y 1 , y 2 , , y n ) 属于 C ,且 x 1 , x 2 , , x n 由(10.4)给出;对于不属于 C ( y 1 , y 2 , , y n )  

应当注意,(10.5)是(8.18)的推广。我们留给读者去阐述(8.22)的类似推广。

我们省略随机变量 Y 1 , Y 2 , , Y n 是联合连续且具有联合概率密度的证明。我们概述一下由(10.5)给出的联合概率密度函数公式的证明。可以证明,对于任意实数 u 1 , u 2 , , u n (10.6)式右边的概率等于

 

其中 D n = { ( x 1 , x 2 , , x n ) : u 1 g 1 ( x 1 , x 2 , , x n ) u 1 + h 1 , , u n g n ( x 1 , x 2 , , x n ) u n + h n } .  

现在,如果 ( u 1 , u 2 , , u n ) 不属于 C ,那么对于足够小的 h 1 , h 2 , , h n 值, D n 中没有点 ( x 1 , x 2 , , x n ) ,且(10.7)中的概率为0。由于(10.6)中取极限的量对于足够小的 h 1 , h 2 , , h n 值为0,因此对于不在 C 中的 ( u 1 , u 2 , , u n ) ,有 f Y 1 , Y 2 , , Y n ( u 1 , u 2 , , u n ) = 0 。这样就证明了(10.5)。为了证明(10.5),我们使用著名的多重积分变量替换公式(参见 R. Courant, 微积分学, Interscience, New York, 1937, Vol II, p. 253, 或 T. Apostol, 数学分析, Addison-Wesley, Reading, Massachusetts, 1957, p. 271)将(10.7)式右边的积分变换为积分

用(10.8)中的积分替换(10.6)式右边的概率,然后取(10.6)中所示的极限,我们最终得到(10.5)。

例10A。设 X 1 X 2 是联合连续的随机变量。令 U 1 = X 1 + X 2 , U 2 = X 1 X 2 。对于任意实数 u 1 u 2 ,证明

 

 

g 1 ( x 1 , x 2 ) = x 1 + x 2 g 2 ( x 1 , x 2 ) = x 1 x 2 。方程 u 1 = x 1 + x 2 u 2 = x 1 x 2 显然有解 x 1 = ( u 1 + u 2 ) / 2 x 2 = ( u 1 u 2 ) / 2 。雅可比行列式 J 由下式给出

 

J = | g 1 x 1 g 1 x 2 g 2 x 1 g 2 x 2 | = | 1 1 1 1 | = 2.  

鉴于这些事实,(10.9)是(10.5)的直接结果。

用完全相同的方法可以建立以下结果:

例10B。设 X 1 X 2 是联合连续的随机变量。令 那么对于任意实数 r α ,满足 r 0 0 α 2 π ,有 应当注意,我们立即从(10.11)得到由(9.13 ) 给出的 f R ( r ) 的公式,因为

 

例10C坐标轴旋转。设 X 1 X 2 是联合分布的随机变量。对于区间 0 α 2 π 中的某个角度 α ,令 那么 为了说明(10.14)的用法,考虑两个联合正态分布的随机变量,其联合概率密度函数由(9.31)给出,且 m 1 = m n 2 = 0 。那么 其中

从(10.15)可以看出,通过坐标轴旋转从联合正态分布的随机变量 X 1 X 2 得到的两个随机变量 Y 1 Y 2 是联合正态分布的。此外,如果旋转角度 α 选择得使

 

那么 B = 0 ,且 Y 1 Y 2 是独立正态分布的。因此通过适当的坐标轴旋转,两个联合正态分布的随机变量可以变换为两个独立正态分布的随机变量

理论习题

10.1。设 X 1 X 2 是独立的随机变量,每个都服从参数为 λ 的指数分布。证明随机变量 X 1 + X 2 X 1 / X 2 是独立的。

10.2。设 X 1 X 2 是独立的随机变量,每个都服从参数为 m = 0 σ > 0 的正态分布。证明 X 1 2 + X 2 2 X 1 / X 2 是独立的。

10.3。设 X 1 X 2 是独立的随机变量,分别服从自由度为 n 1 n 2 χ 2 分布。证明 X 1 + X 2 X 1 / X 2 是独立的。

10.4。设 X 1 , X 2 X 3 是独立同分布的正态随机变量。令 X ¯ = ( X 1 + X 2 + X 3 ) / 3 S = ( X 1 X ¯ ) 2 + ( X 2 X ¯ ) 2 + ( X 3 X ¯ ) 2 。证明 X ¯ S 是独立的。

10.5生成正态分布随机变量的随机样本。设 U 1 , U 2 是独立的随机变量,每个都在区间0到1上均匀分布。证明随机变量是独立的随机变量,每个都服从均值为0、方差为1的正态分布。(关于此结果的讨论,参见 G. E. P. Box 和 Mervin E. Muller, “A note on the generation of random normal deviates,” 数理统计年鉴, Vol. 29 (1958), pp. 610-611。)

习题

10.1。设 X 1 X 2 是独立的随机变量,每个都服从参数为 λ = 1 2 的指数分布。求 Y 1 Y 2 的联合概率密度函数,其中 (i) Y 1 = X 1 + X 2 , Y 2 = X 1 X 2 ,(ii) Y 1 = 最大值 ( X 1 , X 2 ) , Y 2 = ( X 1 , X 2 )

 

答案

(i) 1 8 e 1 / 2 y 1 如果 y 1 > 0 , | y 2 | y 1 ; 0 否则;(ii) 1 2 e 1 / 2 ( y 1 + y 2 ) 如果 0 y 2 < y 1 y 1 0 ; 0 否则。

 

10.2。设 X 1 X 2 具有由下式给出的联合概率密度函数如果否则 ( R , θ ) 的联合概率密度函数,其中 R = X 1 2 + X 2 2 θ = tan 1 X 2 / X 1 。证明并解释为什么 R 2 是均匀分布的而 R 不是。

10.3。设 X Y 是独立的随机变量,每个都在区间0到1上均匀分布。求随机变量 R θ 的个体和联合概率密度函数,其中 R = X 2 + Y 2 θ = tan 1 Y X

 

Answer

f R , 0 ( r , α ) = r if 0 < r cos α , r sin α < 1 ; 0 otherwise; f R ( r ) = π 2 r for 0 < r 1 , ( 2 csc 1 r π 2 ) for 1 r 2 ; 0 otherwise; f 0 ( θ ) = 1 2 sec 2 θ for 0 θ π 4 ; 1 2 csc 2 θ for π 4 θ π 2 ; 0 otherwise.

 

10.4 . Two voltages X ( t ) and Y ( t ) are independently and normally distributed with parameters m = 0 and σ = 1 . These are combined to give two new voltages, U ( t ) = X ( t ) + Y ( t ) and V ( t ) = X ( t ) Y ( t ) . Find the joint probability density function of U ( t ) and V ( t ) . Are U ( t ) and V ( t ) independent? Find P [ U ( t ) > 0 , V ( t ) < 0 ] .