UTF8gbsn

假如x,y∼U(0,1)x,y\sim U(0,1)x,yU(0,1), 也就是说x,yx,yx,y都是(0,1)(0,1)(0,1)上的均匀分布.
问题是:z=∣x−y∣z=|x-y|z=xy是什么分布?

采样并列图表

import numpy as np
import matplotlib.pyplot as plt

N = 500000

data = []
for k in range(N):
    x1 = np.random.rand()
    x2 = np.random.rand()
    d = np.abs(x1-x2)
    data.append(d)

n, bins, patches = plt.hist(data, 50)

plt.title('Histogram some')
plt.grid(True)
plt.show()

在这里插入图片描述

可见这z=∣x−y∣z=|x-y|z=xy的分布是一个三角形分布.

分析

通过上面的分析, 我们可以确定z=∣x−y∣z=|x-y|z=xy的概率密度函数为 p(z)={2(1−z),z∈(0,1)0,(otherwise)p(z)=\left\{ \begin{aligned} 2(1-z), z\in (0,1)\\ 0, \quad (otherwise) \end{aligned} \right.p(z)={2(1z),z(0,1)0,(otherwise)

但是如何分析得到这个概率密度函数呢? 我的方法是画出z=∣x−y∣z=|x-y|z=xy的函数图像.

在这里插入图片描述

注意上面的图像是个3D图像. 那么我们如何计算概率密度? 首先计算,
图像表面积是多少. ∫0122(1−z)dz=2\int_{0}^{1}{2\sqrt{2}(1-z)dz}=\sqrt{2}0122 (1z)dz=2

高度为z的点构成两条直线, 这两条直线的长度为22(1−z)2\sqrt{2}(1-z)22 (1z),
有此可见,点zzz处的概率密度为 22(1−z)/2=2(1−z)2\sqrt{2}(1-z)/\sqrt{2}=2(1-z)22 (1z)/2 =2(1z).
于是原式得证. 这里其实就是要除以一个表面积的总和而已.
加粗样式

Logo

有“AI”的1024 = 2048,欢迎大家加入2048 AI社区

更多推荐