python - Black Scholes 看跌期权，分析解决方案实现 - Python

Question

我正在尝试创建一个简单的函数来解决给定股票值(value)数组 x0、某个执行价格 K、无风险利率 r、波动性和到期时间 T 的看跌期权的值(value)。公式在这里找到:

https://en.wikipedia.org/wiki/Black%E2%80%93Scholes_model

我认为这是一个相对直接的功能来实现，但我似乎得到了更大股票价格的错误值。这是我的代码:

def FinalTimeAnalyticalSolution(x0,K,r,vol,T): 

   x0 = np.array(x0)
   N = len(x0) - 1
   t = 0.0

   U = np.empty(N+1)

   U[0] = K * math.exp(-r * (T - t))

   for i in range(1,N+1):

       d1 = (1.0/(vol*(math.sqrt(T - t)))) * (math.log(x0[i]/K) + (r - (vol**2)/2)*(T-t))
       d2 = d1 - vol * math.sqrt(T-t)

       U[i] = norm.cdf(-d2) * K * math.exp(-r * (T - t)) - norm.cdf(-d1) * x0[i]

return U

例如，

x = np.linspace(0,3,7)
U = FinalTimeAnalyticalSolution(x,2.0,0.04,0.2,2.0) 
for i in range(len(x)):
    print x[i]
    print U[i]

给我结果

0.0, 1.84623269277

0.5, 1.34623275338

1.0, 0.847891622826

1.5, 0.404448071531

2.0, 0.139395618544

2.5, 0.0372942572641

3.0, 0.00832440115615

我一直在用期权值(value)计算器在线测试它，比如
http://www.danielsoper.com/fincalc/calc.aspx?id=38我得到以下信息。

0.0, 1.834133,
0.5, 1.334133,
1.0, 0.836973,
1.5, 0.405401,
2.0, 0.15085,
2.5, 0.047112,
3.0, 0.013392,

正如你所看到的，我的解决方案对于小 x 是相对正确的，但随着它的增加，它们越来越远。我真的不知道为什么。我觉得这是我的代码中的一个小错误，但我已经被困了一天。

提前致谢
詹姆士

Answer 1

错误是 d1 公式中的单号:而不是

... (r - (vol**2)/2) ...

它应该是

... (r + (vol**2)/2) ...