泊松分布、泊松过程、泊松点过程

西方无朔 2018-03-15

展开全文

1.泊松分布

泊松分布是二项分布的极限分布，假设有一列二项分布B(n,p_n)，均值为 $λ$ ,即 $lim_{n \to \infty} n p_{n} = λ > 0$ ，对任何非负整数k(即发生k次的概率)有 $lim_{n \to \infty} b (k; n, p_{n}) = lim_{n \to \infty} C_{n}^{k} p_{n}^{k} (1 - p_{n})^{n - k} = e^{- λ} \frac{λ^{k}}{k!}$ 。
证明：
$\begin{aligned} (1) & C_{n}^{k} p_{n}^{k} (1 - p_{n})^{n - k} & = \frac{n!}{k! (n - k)!} p_{n}^{k} (1 - p_{n})^{n - k} \\ (2) & = \frac{1 \times 2 \times 3 \times . . . \times n}{k! \times 1 \times 2 \times 3... \times (n - k) \times n^{k}} (1 - p_{n})^{- k} (n p_{n})^{k} (1 - p_{n})^{n} \\ (3) & = \frac{n \times (n - 1) \times (n - 2) \times . . . \times (n - k + 1)}{k! \times n^{k}} (1 - p_{n})^{- k} (n p_{n})^{k} (1 - p_{n})^{n} \\ (4) & = \frac{1}{k!} (1 - \frac{1}{n}) (1 - \frac{2}{n}) . . . (1 - \frac{k - 1}{n}) (1 - p_{n})^{- k} (n p_{n})^{k} (1 - p_{n})^{n} \end{aligned}$
注意到 $lim_{n \to \infty} (n p_{n})^{k} = λ^{k}$ ，和 $lim_{n \to \infty} (1 - p_{n})^{n} = e^{- λ}$ 。

定理证毕。

泊松分布是二项分布的极限分布，当n很大，p很小时，二项分布就可以近似地看成时参数 $λ = n p$ 的泊松分布。

2.泊松过程

实验结果满足泊松分布的实验即为泊松过程。

3.泊松点过程

泊松点过程其实和泊松过程并无区别。只是在我初接触的时候不自觉的把它当成一个二维的撒点过程。所以我想更多人会把这个术语当做是如何在二维平面撒满足泊松分布点的方法。放心，这里也是介绍方法的。

3.1一维的撒点方法

3.1.1算法1

我们注意到，在齐次泊松过程中，两次事件的距离是满足均值为 $\frac{1}{λ}$ 的指数分布。

(0) 初始化 t = 0;

(1) 取一个满足均匀分布u~U(0,1)的随机数u;

(2) $t = t - \frac{1}{λ} l o g (u)$ ;

(3) 生成一个点t；

(4) 返回(1)。

3.1.2算法2

假设在固定的时长 $[0, t_{0}]$ ，事件发生次数为 $N (t_{0}) = n ，事件发生的时间$ T_1,T_2,...,T_n $(排序过的) 满足均匀分布。 (1) 生成随机变量 n 满足泊松分布$ n\sim Poisson(\lambda t);

(2) 如果n=0,结束退出。否则，独立地生成 $u_{1} \sim U (0, 1), u_{2} \sim U (0, 1), . . ., u_{n} \sim U (0, 1)$ ，然后将得到的 $u_{i}$ 进行排序（按升序）得到 $u_{(1)}, u_{(2)}, u_{(3)}, . . ., u_{(n)}$

(3) $t_{i} = t_{0} u_{(i)}, i = 1, 2, . . . ., n$

(4) 生成 $t_{i}$ 。

3.2二维的撒点方法

二维撒点满足的泊松分布最关键的唯一和一维的不同点在于泊松分布参数由 $λ t$ 变为 $λ \times A (R)$ ,A(R)是区域R的容量(对更高维也可以用的定义)。在半径为r>0的的圆平面更容易完成满足泊松分布的撒点过程。

3.2.1算法1

3.2.1.1定理1

假设 $(R_{1}, θ_{1}), (R_{2}, θ_{2}), . . ., (R_{N}, θ_{N})$ 是极坐标系中代表N>0的事件的位置，分布满足在圆 $C \equiv (x, y) : x^{2} + y^{2} \leq r^{2}$ 内的齐次泊松过程。设N=n,排好序的事件的半径 $R_{(1)}, R_{(2)}, . . ., R_{(n)}$ 符合密度函数为 $f (z) = \frac{2 z}{r^{2}}, z \in [0, r]$ 的分布， $θ_{1}, θ_{2}, . . ., θ_{n}$ 独立同分布于均匀分布 $U [0, 2 π]$ 并且和 $R_{1}, R_{2}, . . ., R_{n}$ 独立。
下面就提出利用齐次泊松过程在半径为r的圆上生成N=n个点，满足参数为 $π r^{2} λ$ 的泊松分布。生成的点的极坐标半径满足上述定理。

3.2.1.2算法实现

(0) $n \sim P o i s s o n (π r^{2} λ)$ 。如果n=0，结束退出。否则，独立地生成 $u_{1} \sim U (0, 1), u_{2} \sim U (0, 1), . . ., u_{n} \sim U (0, 1)$

(1) $R_{1} \leftarrow r \sqrt{u_{1}}, R_{2} \leftarrow r {\sqrt{u}}_{2}, . . . R_{n} \leftarrow r {\sqrt{u}}_{n}$ 。

(2) 对 $R_{1}, R_{2}, . . ., R_{n}$ 排序（升序）得到 $R_{(1)}, R_{(2)}, . . ., R_{(n)}$ 。

(3) 独立地生成 $u_{n + 1} \sim U (0, 1), u_{n + 2} \sim U (0, 1), . . ., u_{2 n} \sim U (0, 1)$ 。

(4) $θ_{1} \leftarrow 2 π u_{n + 1}, θ_{2} \leftarrow 2 π u_{n + 2}, . . ., θ_{n} \leftarrow 2 π u_{2 n}$ 。

(5) 生成点的坐标 $(R_{(1)}, θ_{1}), (R_{(2)}, θ_{2}), . . ., (R_{(n)}, θ_{n})$ 。

3.2.2算法2

在更加不规则的区域内完成。如:感兴趣的区域为 $C \equiv (x, y) : x \geq 0, y \geq 0, y \leq f (x), x \leq T$ .其中排过序的 $X_{(1)}, X_{(2)}, . . ., X_{(N)}$ 是事件的横坐标，那么 $\int_{X_{(} i - 1)}^{X_{(} i)} f (x) d x, X_{(0)} = 0, i = 1, 2, . . ., 独立同分布于均值为 \frac{1}{λ} 的指数分布$ 。其相应的纵坐标 $Y_{i}$ 服从均匀分布 $U (0, f (X_{(i)}))$