Можно ли создать символическое многомерное нормальное распределение с 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.