QuantLib 金融计算——随机过程之通常 Black Scholes 过程

目录

• QuantLib 金融计算——随机过程之通常 Black Scholes 过程
  • 通常 Black Scholes 过程

若是未作特别说明，文中的程序都是 Python3 代码。

QuantLib 金融计算——随机过程之通常 Black Scholes 过程

载入模块

import QuantLib as ql
import pandas as pd
import numpy as np
import seaborn as sn

print(ql.__version__)

1.12

通常 Black Scholes 过程

quantlib-python 中 Black Scholes 框架下常见的几种随机过程均派生自基类 GeneralizedBlackScholesProcess，而 GeneralizedBlackScholesProcess 模拟下列 SDE 描述的一维随机过程：

\[ d \ln S_t = \left( r ( t ) - q ( t ) - \frac { \sigma \left( t , S_t \right)^2 } 2 \right) d t + \sigma d W_t \]

等式使用风险中性漂移而不是通常漂移 \(\mu\)。风险中性利率由股息率 \(q(t)\) 调整，而且相应的扩散项是 \(\sigma\)。

做为基类，GeneralizedBlackScholesProcess 的构造函数为

GeneralizedBlackScholesProcess(x0,
                               dividendTS,
                               riskFreeTS,
                               blackVolTS)

其中：

• x0：QuoteHandle 对象，表示 SDE 的起始值；
• dividendTS：YieldTermStructureHandle 对象，表示股息率的期限结构
• riskFreeTS：YieldTermStructureHandle 对象，表示无风险利率的期限结构
• blackVolTS：BlackVolTermStructureHandle 对象，表示波动率的期限结构

GeneralizedBlackScholesProcess 提供了相应的检查器，返回构造函数接受的关键参数：

• stateVariable；
• dividendYield；
• riskFreeRate；
• blackVolatility

从 StochasticProcess1D 继承来的离散化函数 evolve，描述 SDE 从 \(t\) 到 \(t + \Delta t\) 的变化。

QuantLib 提供了一些具体的派生类，这些类表明众所周知的具体过程，如

• BlackScholesProcess：没有股息率的通常 BS 过程；
• BlackScholesMertonProcess：通常 BS 过程；
• BlackProcess：通常 Black 过程；
• GarmanKohlagenProcess：包含外汇利率的通常 BS 过程

这些派生类在构造和调用方式上大同小异，在下面的例子中，咱们将创建一个具备平坦无风险利率、股息率和波动率期限结构的 Black-Scholes-Merton 过程，并画出模拟结果。

def testingStochasticProcesses1():
    refDate = ql.Date(27, ql.January, 2019)
    riskFreeRate = 0.0321
    dividendRate = 0.0128
    spot = 52.0
    vol = 0.2144
    cal = ql.China()
    dc = ql.ActualActual()

    rdHandle = ql.YieldTermStructureHandle(
        ql.FlatForward(refDate, riskFreeRate, dc))
    rqHandle = ql.YieldTermStructureHandle(
        ql.FlatForward(refDate, dividendRate, dc))

    spotQuote = ql.SimpleQuote(spot)
    spotHandle = ql.QuoteHandle(
        ql.SimpleQuote(spot))

    volHandle = ql.BlackVolTermStructureHandle(
        ql.BlackConstantVol(refDate, cal, vol, dc))

    bsmProcess = ql.BlackScholesMertonProcess(
        spotHandle, rqHandle, rdHandle, volHandle)

    seed = 1234
    unifMt = ql.MersenneTwisterUniformRng(seed)
    bmGauss = ql.BoxMullerMersenneTwisterGaussianRng(unifMt)

    dt = 0.004
    numVals = 250

    bsm = pd.DataFrame()

    for i in range(10):
        bsmt = pd.DataFrame(
            dict(
                t=np.linspace(0, dt * numVals, numVals + 1),
                path=np.nan,
                n='p' + str(i)))

        bsmt.loc[0, 'path'] = spotQuote.value()

        x = spotQuote.value()

        for j in range(1, numVals + 1):
            dw = bmGauss.next().value()
            x = bsmProcess.evolve(bsmt.loc[j, 't'], x, dt, dw)
            bsmt.loc[j, 'path'] = x

        bsm = pd.concat([bsm, bsmt])

    sn.lineplot(
        x='t', y='path',
        data=bsm,
        hue='n', legend=None)


testingStochasticProcesses1()