我给出了这两个方程,a(α)和v都有给定值. 我知道这些问题不能用classic 的方法来解决,而是迭代地解决,所以我想用scipy和scipy来计算Phi(分布).
我的主要目标是计算k1和k2(在代码中用x和y替换).
(我用x替换K1,用y替换K2,用a替换阿尔法.)
from sympy import symbols, Eq, solve, exp
from scipy.stats import norm
x, y = symbols('x y')
v = 2.2526 # given
a = 0.1 # given
eq1 = (x + (1/v)) * exp((-(1/2)) * (x**2)) - (-y + (1/v)) * exp((-(1/2)) * (y**2))
eq2 = (((1 - norm.cdf(-1/v))**-1) * (norm.cdf(y) - norm.cdf(-x)) + a - 1)
sol = solve((eq1, eq2), (x, y))
print(sol)
如果我运行这个,我会得到
TypeError("无法确定关系的真值")
这个错误是什么意思?
发帖主题:Re:Kolibrios
x, y = symbols('x y', real=True)
引向
TypeError:输入类型不支持ufunc‘isnan’,并且 无法安全地将输入强制为任何受支持的类型 选角规则是"安全"