python - python sympy 中产品的期望

Question

在 python sympy 中我有以下代码

from sympy.stats import P, E, variance, covariance, Normal

mu, sigma = symbols('mu sigma')
epsilon_t = Normal('epsilon_t', 0, sigma)
epsilon_t1 = Normal('epsilon_t+1', 0, sigma)

定义了两个(大概是独立的)正态随机变量。现在，我希望这样

E(epsilon_t1*epsilon_t)

等于零，但是当我通过 sympy 运行它时，我得到了更复杂的东西。它告诉我当 pi/2 > 2|large(sigma)| 时它为零。以及其他更复杂的事情。我不关心第二部分，我只想它为零。我怎样才能做到这一点？

Answer 1

如果您查看 help(arg)，您会得到一个有用的提示:如果 x 则 arg(x) 为 0 > 是积极的。因此，如果您将参数定义为正数而不是实数，arg 将能够计算:

>>> mu, sigma = symbols('mu sigma',positive=True)
>>> epsilon_t = Normal('epsilon_t', 0, sigma)
>>> epsilon_t1 = Normal('epsilon_t+1', 0, sigma)
>>> E(epsilon_t*epsilon_t1)
0