定義一組隨機數

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 

不過，我想能夠定義說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統計數據目前不支持多個單變量分佈的組合，我們不能一般處理多變量分佈。

如果你願意去數字，那麼你可能會發現感興趣的PyMC3。