1、# 基于Black - Scholes 公式的期权定价公式def bsm娄多骋能_call_value(s,,k,T,r,sigma):from math import log,sqrt,expfro罪焐芡拂m scipy import statss=float(s)d1=(log(s/k)+(r+0.5*sigma**2)*T)/(sqrt(T))d2=(log(s/k)+(r-0.5*sigma**2)*T)/(sqrt(T))value=(s*stats/norm.cdf(d1,0.0,1.0))-K*exp(-r*T)*stats.norm.cdf(d2,0.0,1.0)return value当前价s行权价k上面的函数只是计算隐含波动率所需的基本函数。当然，还需要单独的期权报价，也就是基于VSTOXX指数的欧式看涨期权，以及生成隐含波动率的代码，下面就用IPython来完成。比如把参考日期t=0，指数的收盘价为v=7，6639对于无风险短期利率，假定年率r=0,01这样就可以利用print '期权价格 : %.4f' % bsm_call_value(s,,k,T,r,sigma):
2、另一种写法def bsm_call_value(spot, strike, maturity, r, vol）:from math import log,sqrt,expfrom scipy import statsd1 = (log(spot/strike) + (r + 0.5 * vol *vol) * maturity) / vol / sqrt(maturity)d2 = d1 - vol * sqrt(maturity)price = spot * norm.cdf(d1) - strike * exp(-r*maturity) * norm.cdf(d2)return price当前价spot行权价strike到期期限maturity无风险利率r波动率vol后面会说到迭代所有的到期日并绘制图形
