Когда я пишу функцию плотности вероятности для нормального распределения в Sympy (v. 1.5.1) следующим образом:
pdf_normal = (1 / (sigma*sqrt(2*pi))) * exp(-(1/2) * ((x - mu)/sigma)**2)
на основе такого выражения уравнения (Википедия):
sympy «упрощает» выражение, но значительно усложняет его чтение:
Есть ли способ предотвратить это?
(Это обсуждение не помогает, потому что предлагаемый ответ вообще не касается выражения )
Заранее спасибо!
Вы можете использовать evaluate=False для управления этим, но, скорее всего, это будет оцениваться в последующих операциях:
In [11]: Pow(Mul(sigma, sqrt(2*pi, evaluate=False), evaluate=False), -1, evaluate=False)
Out[11]:
1
─────────
_____
σ⋅╲╱ 2⋅π
Работает ли обман символов для вас?
>>> r2pi = Symbol('\\sqrt{2 \\pi}')
>>> (1 / (sigma*r2pi)) * exp(-(1/2) * ((x - mu)/sigma)**2)
Класс UnevaluatedExpr предотвратит повторную оценку выражения во время последующих операций. С изменениями, предложенными здесь, я думаю, это решит вашу проблему:
>>> from sympy import S
>>> from sympy.abc import *
>>> from sympy import S
>>> usqrt = lambda x: UnevaluatedExpr(sqrt(x, evaluate=False))
>>> (1 / (sigma*usqrt(2*pi))) * exp(((x - mu)/sigma)**2/-2)
exp(-(-mu + x)**2/(2*sigma**2))/(sigma*sqrt(2*pi))
