是否可以用sympy创建符号多元正态分布？

Question

可以使用给定的均值和方差值创建多元法线：

from sympy.stats import Normal
from sympy import *

mu = Matrix([1,0])
Sigma = Matrix([[1,0],[0,1]])
X = Normal("X", mu, Sigma)

我想对符号矩阵做同样的事情：

n=5
mu = MatrixSymbol("mu", n, 1)
Sigma = MatrixSymbol("Sigma", n, n)
X = Normal("X", mu, Sigma)

导致TypeError: Mix of Matrix and Scalar symbols

这有可能吗？

Answer 1

如果您使用的是sympy 1.4，我认为您需要使用：

from sympy.stats.joint_rv_types import MultivariateNormalDistribution

但是您上面的代码仅适用于开发版本，即即将发布的版本1.5。