Sympy计算集合概率测度时调用P传入Interval触发TypeError如何解决

问题原因 #

Python原生的in运算符会强制要求运算结果为布尔值，而Sympy中随机变量与Interval做in判断时返回的是符号逻辑表达式，无法直接转为布尔值，因此触发类型错误。

解决方案 #

有两种可行的写法：

• 直接传入区间对应的不等式组合给P函数：

from sympy.stats import Normal, P
n = Normal("n", 0, 1)
result = P((n >= -1) & (n <= 1))

• 如果需要使用Interval对象，调用Interval的contains方法生成符号逻辑表达式传入：

from sympy.stats import Normal, P
from sympy import Interval
n = Normal("n", 0, 1)
result = P(Interval(-1, 1).contains(n))

两种写法都可以得到标准正态分布在±1σ区间的概率结果：erf(√2/2)，转换为数值约为0.6827。

内容的提问来源于stack exchange，提问作者Wojciech Kulma