绘制Vasicek模型。只能将大小为-1的数组转换为python标量

Posted 2021-04-10

[尝试绘制Vasiceks投资组合损失分布图，当我想为每个时间步长在0和1之间建立x的数字数组时，我遇到了问题。

x.astype(int) 
x = np.arange(0.000001, 1, 0.000001)
a1 = math.sqrt((1-rho)/rho)
a2 = -1/(2*rho)*((math.sqrt(1-rho)*norm.ppf(x)-norm.ppf(p)))**2
a3 = ((1/2)*norm.ppf(x))**2

lossdist = a1*math.exp(a2+a3))  

以上是我的模型，但是当我运行它时，它进入了lossdist函数并暂停，生成

TypeError: only size-1 arrays can be converted to Python scalars

我看了一个关于上一个问题的答案，但没有任何效果，尝试使用x.astype(int)和np.vectorise解决方案。我该如何解决？

答案

math.exp仅适用于标量输入，但是您要传递数组作为输入。您应该使用numpy.exp，因为它接受一个数组作为输入。

参考：https://docs.scipy.org/doc/numpy/reference/generated/numpy.exp.html