随机变量的变换

                     

贡献者: addis

预备知识 概率分布函数

1. 求新变量的分布函数

   我们先来讨论这样一个问题:令两个随机变量 x1,x2 间有函数关系 x1=g(x2),若已知 x1 的分布函数为 f1(x1),求 x2 的分布函数 f2(x2)

   将两个概率分布写成微分形式,有

(1)dP=f1(x1)dx1(2)dP=f2(x2)dx2 .
若将式 1 中的 x1 替换成 g(x2)dx1 替换成 g(x2)dx2,有
(3)dP=f1[g(x2)]g(x2)dx2 .
对比式 2 ,得
(4)f2(x)=f1[g(x)]g(x) ,
这样,就求出了 x2 的分布函数 f2(x)

例 1 线性变换

   若随机变量 x 的分布为 f(x),那么 y=ax 的分布为 g(y)=f(y/a)/a。平均值变为 aμ,方差变为 a2σ2

例 2 

   已知 f1(x1)=3x12x1[0,1]x2=x12,求 x2 的分布函数。

   用 x2 表示 x1x1=x2,代入 dP=f1(x1)dx1,得

(5)dP=3x22d(x2)=32x2dx2 .
所以 x2 的分布函数为 f2(x2)=3x2/2

2. 求两变量的关系

预备知识 可分离变量的微分方程

   另一个常见的问题是已知 x1x2 的分布函数 f1(x1),f2(x2),求两个随机变量需要满足的函数关系。

   对比式 1 式 2 可得一个已分离变量的微分方程

(6)f1(x1)dx1=f2(x2)dx2 .
将方程两边积分即可得到两变量所满足的函数关系
(7)F1(x1)=F2(x2)+C .
其中函数 F1,F2 分别是函数 f1,f2 的一个原函数,待定常数 C 通常可以由 x1x2 的取值范围确定。

   这个问题的一个常见应用是在程序中生成指定分布函数的随机变量。在许多编程语言中,随机数生成器只能生成一个从 0 到 1 均匀分布的随机变量(即 f(x)=1),若我们需要一个其他分布的随机变量,就可以使用以上方法。

例 3 

   已知随机变量 x1  (x1[0,1]) 的分布函数为 f1(x1)=1,求函数关系 x2=g(x1) 使得 x2 的分布函数为 f2(x2)=2x2  (x2[0,1])

   将 f1,f2 代入式 6 并两边积分得

(8)x1=x22+C .
由于 x1x2 的区间关系得 C=0,所以有 x2=x1

例 4 

   给出两个随机变量 ξ1,ξ2  (ξ1,ξ2[0,1]),分布函数均为 f(ξi)=1,用 ξ1,ξ2 表示某随机点的极坐标 (r,θ) 使得该点在单位圆内均匀随机分布。

   要使随机点在单位圆内随机分布,θ 显然应该在 [0,2π] 间均匀随机分布,所以令 θ=2πξ2 即可。要决定 r 的分布函数,我们把单位圆划分为许多小圆环,随机点出现在某圆环内的概率等于该圆环的面积比单位圆的面积,即

(9)dP=2πrdrπ=2rdr .
所以 r 的分布函数为 2r。令 rξ1 间存在函数关系,由式 6
(10)1dξ1=2rdr .
两边积分得 ξ1=r2ξ1=0r=0,所以积分常数为零),即 r=ξ1。所以,随机点的极坐标取 (ξ1,2πξ2) 即可。


致读者: 小时百科一直以来坚持所有内容免费无广告,这导致我们处于严重的亏损状态。 长此以往很可能会最终导致我们不得不选择大量广告以及内容付费等。 因此,我们请求广大读者热心打赏 ,使网站得以健康发展。 如果看到这条信息的每位读者能慷慨打赏 20 元,我们一周就能脱离亏损, 并在接下来的一年里向所有读者继续免费提供优质内容。 但遗憾的是只有不到 1% 的读者愿意捐款, 他们的付出帮助了 99% 的读者免费获取知识, 我们在此表示感谢。

                     

友情链接: 超理论坛 | ©小时科技 保留一切权利