python - 定义一组随机数

Question

我正在尝试使用 sympy.stats 模块来帮助计算表达式的期望值。在最简单的形式中，我的问题如下:

我知道如何定义两个正态分布的随机变量并计算期望值

from sympy.stats import Normal, E
import sympy as sp
sigma = sp.symbols('sigma', positive=True)

x1 = Normal("x1", 0, sigma)
x2 = Normal("x2", 0, sigma)
print E(x1 * x2)

Out: 0

但是，我希望能够定义 say X，它是 N 的向量(或列表)正态分布的随机变量。问题是我需要 N 来保持任意性我无法创建列表或任何东西。

有人对如何进行有任何建议吗？

我需要这样做的原因是我将对 X 的组件进行求和，其中一些将相关，而另一些则不会。

Answer 1

我认为您正在寻找类似以下的内容？

mu = MatrixSymbol('mu', n, 1)
Sigma = MatrixSymbol('Sigma', n, n)

X = Normal('X', mu, sigma)  # A multivariate normal distribution

如果是这样，那么您目前不走运。 SymPy stats 目前只支持几个单变量分布的组合，我们一般不处理多变量分布。

如果您愿意使用数字，那么您可能会对 PyMC3 感兴趣。