我正在尝试根据两个均匀随机变量创建一个混合随机变量:
import sympy as sp
import sympy.stats as ss
x = sp.symbols('x', real=True, positive=True)
mixed = ss.Uniform('X', 10, 20) * sp.Rational(1, 3) + ss.Uniform('Y', 15, 30) * sp.Rational(2, 3)
d = ss.density(mixed)(x).simplify()
我已经能够使用两个Normal随机变量拼凑出一个类似的示例,因此我认为使用Uniform随机变量w这样的事情将是微不足道的。
import sympy as sp
import sympy.stats as ss
x = sp.symbols('x', real=True, positive=True)
mixed = ss.Normal('X', 1, 1) * sp.Rational(1, 3) + ss.Normal('Y', 2, 1) * sp.Rational(2, 3)
d = ss.density(mixed)(x).simplify()
display(d.subs(x, 1.5))