python sympy中对产品的期望

Posted 2023-02-24

【中文标题】python sympy中对产品的期望

【英文标题】：Expectation of a product in python sympy

【发布时间】：2020-02-04 11:12:57

【问题描述】：

在 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|larg(sigma)| 为零还有一些更复杂的东西。我不在乎第二部分，我只希望它为零。我怎样才能做到这一点？

【参考方案1】：

如果您查看help(arg)，您会得到一个有用的提示：如果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