Python中具有指数分布的马尔可夫链

Question

我是Python中随机数生成的新手。

我理解马尔可夫链的概念，但我无法将其转换为有用的代码：

Markov Chain - P(X_m+1 = j|X_m = i, X_m-1 = x_m-1,.....,X_0 = i_0) = P(X_m+1 = j |X_m = i)

例如，我如何使用Python生成根据指数分布：

P(x) = Ne^{−|x|}

其中N是归一化因子？

Answer 1

我想你要找的是

import numpy.random as npr    
p_x = npr.exponential(N,t)

其中N是缩放因子的逆，t是要生成的随机数。您必须根据链的顺序提前计算参数(在您的情况下，为1(。