广东工业大学学报  2024, Vol. 41Issue (1): 127-134.  DOI: 10.12052/gdutxb.220185.
0

引用本文 

孙有发, 彭文彦. 基于前景理论框架和Heston模型的行为期权定价[J]. 广东工业大学学报, 2024, 41(1): 127-134. DOI: 10.12052/gdutxb.220185.
Sun You-fa, Peng Wen-yan. Behavioral Option Pricing under Prospect Theory Framework and Heston Model[J]. JOURNAL OF GUANGDONG UNIVERSITY OF TECHNOLOGY, 2024, 41(1): 127-134. DOI: 10.12052/gdutxb.220185.

基金项目:

国家自然科学基金资助项目(72271064, 71771058);广东省自然科学基金资助项目(2022A1515011125, 2017A030313400)

作者简介:

孙有发(1976–),男,教授,博士,主要研究方向为计算金融、计算实验金融、金融风险管理,E-mail:y.f.sun@gdut.edu.cn

通信作者

彭文彦(1998–),男,硕士研究生,主要研究方向为计算金融,E-mail:pwy13100614752@126.com

文章历史

收稿日期:2022-12-07
基于前景理论框架和Heston模型的行为期权定价
孙有发, 彭文彦    
广东工业大学 经济学院, 广东 广州 510520
摘要: 行为期权定价是当前国际金融领域的热门研究主题之一。虽然随机波动率模型已成为国际衍生品定价领域的标准模型,但该模型对短到期期权(尤其是虚值期权)的定价仍不准确,其原因之一是传统的期权定价方法忽略了现实市场中的非理性心理和行为因素。针对上述问题,本文运用前景理论期权定价框架,引入价值函数来刻画投资者面对收益与损失的前景价值判断,引入决策权重函数来修正Heston随机波动率模型刻画的资产价格路径的概率密度函数,将期权合约签订与交割的现金流视为分散的心理账户情形,在市场均衡条件下推导出Heston模型下欧氏行为期权的定价公式。上证50ETF期权的实证结果表明:考虑了前景理论的Heston随机波动率模型,能显著地提升短到期虚值期权的定价准确度;参数校正结果发现,定价性能的提升要归因于Heston模型中纳入的表征非理性心理与情绪的行为参数;相对而言,投资者对实值期权的风险态度偏中性,因此行为参数对其定价精度的提升有限。
关键词: 行为期权定价    前景理论    心理账户    Heston模型    上证50ETF期权    
Behavioral Option Pricing under Prospect Theory Framework and Heston Model
Sun You-fa, Peng Wen-yan    
School of Economics, Guangdong University of Technology, Guangzhou 510520, China
Abstract: Behavioral option pricing has been one of the hottest frontiers in the area of international finance. Stochastic volatility model has become the de facto standard model of international derivatives pricing, but it is not accurate enough in the pricing of short-term options, especially for OTM. One of the reasons is that the traditional option pricing methods ignore the irrational psychological and behavioral factors in the real market. To solve this problem,prospect theory is brought into the traditional option pricing framework, the different value judgments of investors facing gains and losses are described by introducing the subjective value function, the subjective decision weight function is used to modify the probability density function of the asset price path described by Heston model, and then the cash flow in the two periods of signing and executing European option contracts is regarded as two separate mental accounts. Prices of European options under Heston model are derived under the condition of market equilibrium. The empirical results of SSE 50ETF option show that the Heston stochastic volatility model considering the prospect theory can significantly improve the pricing accuracy of short maturity OTM option. The model parameter correction results show that the improvement of pricing performance is due to the behavioral parameters representing irrational psychology and emotion included in Heston model. Relatively speaking, investors' risk attitude towards ITM options is nearly neutral, so the improvement of behavior parameters on its pricing accuracy is limited.
Key words: behavioral option pricing    prospect theory    mental account    Heston model    SSE 50ETF option    

主流的理性期权定价模型 (Rational Option Pricing Model, ROPM) 的相关文献通常从特定的标的资产价格动态过程出发,应用风险中性定价理论计算欧氏期权价格。为了提升模型的定价性能,学界开发了丰富的标的资产价格和波动率模型,如Black-Scholes (BS) 模型[1]、跳变模型[2-3]和随机波动率模型等[4-5]。然而,这些模型对现实期权产品,尤其是短到期虚值 (Out of The Money, OTM) 期权的定价能力远远不足,其中一个重要原因是BOPM期权定价框架没有考虑投资者的非理性心理与行为因素[6-8],随着期权到期日的临近,这些非理性心理与行为因素会对期权价格产生显著影响。

近年基于行为金融学相关理论与方法而兴起的行为期权定价模型(Behavior Option Pricing Model, BOPM) 则为期权定价研究提供了一个新的视角。已有的BOPM大致可分为2个流派:一是引入非理性因子到标的资产价格动态过程中,再应用ROPM的相关理论推导出期权定价公式[9-14];二是通过行为金融学理论来刻画投资者的价值判断、风险态度和决策行为,再根据市场交易机制和均衡定价原理求解期权价格[15-23]。流派一的研究只是将非理性心理与行为纳入标的资产价格模型中,本质上仍属于传统的期权定价方法;流派二的研究则遵从行为金融学说中的有限理性、有限认知、有限控制力和有限自利等,然后依据期权市场交易机制,推理、推导或计算相应的期权价格,属于正统的行为期权定价。本研究属于第二个流派。

在实践应用层面,当前行为期权定价模型尚未在金融领域有大量应用。不过,Baele等[18]认为:“在未来,行为金融学有很大的潜力成为资产定价的主导方法”。此外,Bar-Gill[24]的研究表明,许多法律规则能用期权理论解释,行为期权定价模型可用来指导制定法律政策。

作为行为金融学的核心理论,前景理论能较好地解释“资产溢价之谜”[25]、OTM期权低回报率、高方差溢价[18]等市场异象。同时,前景理论框架下的期权定价方法较传统理性框架下的期权定价方法呈现出显著优势。在BOPM领域,Versluis[26]率先提出了基于前景理论[27-28]框架的行为期权定价框架。该文假设标的资产价格过程服从几何布朗运动 (Geometric Brownian Motion, GBM) 模型;利用决策权重函数对标的资产价格路径的概率密度函数进行了非线性化修正;并将欧氏期权合约签订与执行两时段的现金流视为分离或聚合的两种心理账户情形[29-30],并计算其前景价值;最后依据市场出清条件,构建这2种情形下的行为期权定价模型。Versluis[26]的实证结果表明,与传统理性框架下BS模型相比,两种期权定价模型的定价精度都有所提高;尤其是现金流聚合情形下的定价模型,拥有与Heston模型相当的定价精度;只不过,该情形下缺乏期权解析公式,计算复杂度相对更高。Pena和Alemanni[31]、Nardon和 Pianca[32-33]进一步拓展了Versluis[26]的研究:前者选取了3个月到期的看涨期权作为研究样本,发现基于经典GBM模型的行为期权定价模型的拟合精度高于Heston模型;后者分别构建了买方和卖方两个角度下的行为期权定价模型,并通过数值实验考察不同的权重函数形式对期权价格的影响。

不难发现,上述文献大多假设标的资产价格过程服从GBM模型。但由于GBM模型假设标的资产价格服从对数正态分布且波动率为常数,与现实资产价格经验统计事实严重不符,导致基于其计算出的期权价格仍然不够准确[34]

随机波动率模型是刻画现实资产价格行为最为准确的模型门类之一,已成为“复杂金融产品定价的事实标准模型”[35]。Heston模型则是其中的代表模型,它不仅能捕捉到资产收益率分布的“尖峰”“厚尾”“有偏”等特性,且兼具模型结构简单、可计算性较强等优势,被广泛应用于金融衍生品定价。

已有的研究表明,Heston模型难以刻画金融市场剧烈波动期间“陡峭的波动率倾斜”[36]。因此,本文利用前景理论期权定价框架结合Heston模型构建了基于Heston随机波动率模型的行为期权定价方法。该方法在一定程度上克服了Heston模型的这个缺点。

假定标的资产价格过程服从Heston模型,并应用费曼卡兹定理推导其特征函数的偏微分方程;再通过反演公式获得相应的概率密度函数和累积分布函数,并将其嵌入Nardon和 Pianca[32-33]的前景理论期权定价框架;基于期权合约签订与执行的现金流视为分离的心理账户情形,分别推导买方和卖方角度的Heston行为期权定价公式。由于Heston模型的概率密度函数和分布函数较为复杂,本文用级数展开式对其进行逼近,以此提高模型的计算性。

上证50ETF期权的定价结果表明,本文提出的行为期权定价方法能够显著提高短到期OTM期权的定价精度,其结果明显优于传统理性框架下的Heston模型的期权定价结果以及基于GBM模型的行为期权定价结果。模型参数校正结果表明,对于OTM期权,投资者非理性情绪突出,在期权定价模型中纳入非理性情绪与行为因素,可以极大提升OTM期权定价精度;而对于实值 (In The Money, ITM) 期权,投资者基本呈现风险中性态度,行为参数的引入只能轻微提升模型定价精度。

本研究构建了基于随机波动率模型的行为期权定价方法,验证了行为金融学说中的前景理论对于提升期权定价性能的价值,并给出了具体的解析定价公式。

相较于传统理性期权定价模型而言,本文的模型参数蕴含了投资者情绪、风险态度等更为丰富的市场交易信息,且定价更为准确,市场监控者、期权做市商、期权交易者等都能从中挖掘丰富的信号以监控市场、调节价格、对冲风险等。因此,本文提出的模型具备较强的学术意义和实践意义。

1 Heston行为期权定价模型 1.1 前景理论和累积前景理论

Kahneman和Tversky提出的前景理论(Prospect Theory) [27]的核心观点为投资者的决策行为不遵循期望效用最大化原则,而是取决于期望效用与参照点的差距:当期望效用高于参照点时,投资者表现出风险厌恶;当期望效用低于参照点时,投资者表现出风险偏好。因此,他们创造性地用价值函数$ v $代替效用函数,来描述投资者面对收益与损失的价值判断。

假设未来状态的有限集合用前景对$ \left(\left(\mathrm{\Delta }{x}_{1},{p}_{1}\right)\right. , \left(\mathrm{\Delta }{x}_{2},{p}_{2}\right) ,\left.\left(\mathrm{\Delta }{x}_{3},{p}_{3}\right), \cdots, \left(\mathrm{\Delta }{x}_{n},{p}_{n}\right) \right) $来表示。其中,$\mathrm{\Delta }{x}_{i}(1\le i\le n)$是第$i$种相对于参照点$ {x}^{*} $定义的前景,若${x_i}$表示第$i$种前景的绝对值,则$ \mathrm{\Delta }{x}_{i}={x}_{i}-{x}^{*} $$ {p}_{i} $是其客观概率,$ v\left(\mathrm{\Delta }{x}_{i}\right) $是第$i$种前景给投资者带来的前景价值,$ {\pi }_{i}=w\left({p}_{i}\right) $是投资者对第$ i $种前景的决策权重函数。决策权重函数是对客观概率的非线性化修正,其现实依据是:投资者常常表现出对小概率事件的过分重视,对大概率事件又重视不足。

根据前景理论,投资者面临收益和损失表现出不同的风险偏好,因此分别用$ {v}^{+} $$ {v}^{-} $表示收益区间($ 0\le i\le n $) 和损失区间($ -m\le i\le 0 $) 的价值函数,相应的决策权重函数为$ {w}^{+} $$ {w}^{-} $,未来状态汇总的前景价值为

$\begin{split} V=&{\displaystyle\sum\limits _{i=-m}^{n}}{\pi }_{i}v\left(\mathrm{\Delta }{x}_{i}\right) ={\displaystyle\sum\limits _{i=-m}^{0}}{w}^-\left({p}_{i}\right) {v}^-\left(\mathrm{\Delta }{x}_{i}\right) +\\ &{\displaystyle\sum\limits_{i=0}^{n}}{w}^+\left({p}_{i}\right) {v}^+\left(\mathrm{\Delta }{x}_{i}\right) \end{split}$ (1)

在前景理论的基础上,Kahneman和Tversky进一步提出了能计算连续前景价值的累积前景理论(Cumulative Prospect Theory) [28],在$ -m\le i\le n $区间,第$ i $种前景的决策权重函数是累积决策权重函数的差分:

$\begin{split} {\pi }_{i}^-=&{w}^-\left({p}_{1}+\cdots +{p}_{i}\right) -{w}^-\left({p}_{2}+\cdots +{p}_{i-1}\right), \\ &2\leqslant i\leqslant k; \\ {\pi }_{i}^+=&{w}^+\left({p}_{i}+\cdots +{p}_{n}\right) -{w}^+\left({p}_{i+1}+\cdots +{p}_{n}\right),\\ &k+1\leqslant i\leqslant n-1 \end{split}$ (2)

Davies和Satchell[37]则给出了连续情形下,累积前景价值的表达式为

$\begin{split} V=&{\displaystyle\int }_{-\mathrm{\infty }}^{0}{\varPsi }^-\left[F\left(x\right) \right]f\left(x\right) {v}^-\left(x\right) \mathrm{d}x+ \\ &{\displaystyle\int }_{0}^{+\mathrm{\infty }}{\varPsi }^+\left[1-F\left(x\right) \right]f\left(x\right) {v}^+\left(x\right) \mathrm{d}x\end{split}$ (3)

式中:$ \varPsi =\dfrac{\mathrm{d}w\left(p\right) }{p} $为决策权重函数的导数,$ F $$ f $分别为累积分布函数和概率密度函数。

1.1.1 关于价值函数

实证研究中,价值函数通常采取式 (4) 的形式:

$ v\left(x\right) =\left\{ \begin{array}{l} -\lambda (-x{) }^{\beta },x < 0\\ {x}^{\alpha },x\geqslant 0\end{array} \right. $ (4)

式中:$ \alpha $$ \beta $分别衡量投资者面对收益和损失的风险偏好,$ \lambda $为损失厌恶系数,衡量投资者对待损失的敏感程度。基于Tversky和Kahneman的参数设置($ \alpha =\beta =0.88,\lambda =2.25 $),图1为价值函数图像。

图 1 价值函数图像 Figure 1 Figure of value function
1.1.2 关于决策权重函数

决策权重函数${w}{\pm }\left(p\left(x\right) \right)$在不同文献的设定形式略有差异,经典的2种函数形式可见Kahneman、Tversky[28]和Prelec[38]的研究。

$ {w}_{\mathrm{K}\mathrm{a}\mathrm{h}\mathrm{n}\mathrm{e}\mathrm{m}\mathrm{a}\mathrm{n}}\left(p\right) =\frac{{p}^{\gamma }}{({p}^{\gamma }+(1-p{) }^{\gamma }{) }^{\textstyle\frac{1}{\gamma }}} $ (5)
$ {w}_{\mathrm{P}\mathrm{r}\mathrm{e}\mathrm{l}\mathrm{e}\mathrm{c}}\left(p\right) ={\mathrm{e}}^{-(-\mathrm{ln}p{) }^{\gamma }} $ (6)

相较于Kahneman、Tversky[28],Prelec[38]决策权重函数的形式更为简洁,可计算性更好,因此本文选取Prelec[38]的决策权重函数进行研究。图2展示了$ \gamma $取值0.7,0.8,0.9,1时Prelec[38]决策权重函数值与客观概率的关系。当$ \gamma < 1 $时,决策权重函数呈现倒“S”形,低概率事件被过分重视,中、高概率事件重视不足,且$ \gamma $值越小,图像弯曲程度越大。

图 2 Prelec决策权重函数图像 Figure 2 Figure of Prelec’s probability weighting function
1.2 Heston随机波动率模型

为了计算式 (3) 中标的资产的概率密度函数和累积分布函数,需要提前设定标的资产价格动态过程。考虑到Heston模型的优良特性,本文使用Heston模型刻画标的资产价格运动。

$ \left\{\begin{array}{l}{\rm{d}}{S}_{t}=\mu {S}_{t}{\rm{d}}t+\sqrt{{v}_{t}}{S}_{t}{\rm{d}}{W}_{1,t}\\ {\rm{d}}{v}_{t}=\kappa (\theta -{v}_{t}) {\rm{d}}t+\sigma \sqrt{{v}_{t}}{\rm{d}}{W}_{2,t}\\ {\rm{d}}{W}_{1,t}\times {\rm{d}}{W}_{2,t}=\rho {\rm{d}}t\end{array}\right. $ (7)

式中:$ {S}_{t} $$ t $时刻的标的资产价格,$ \mu $为资产收益率,$ {v}_{t} $$t$时刻资产收益率的方差,$ \theta $为资产收益率方差的长期均值,$ \kappa $为方差均值的回复速度,$ \sigma $为资产收益率方差的波动系数,$ {W}_{1,t} $$ {W}_{2,t} $为标准的一维布朗运动,$ \,\rho $为相关系数。

$ {t}_{0},T $分别为期权合约签订日和期权合约到期日,$ \tau =T-{t}_{0} $为期权剩余交易天数,$ K $为期权行权价;为下文便于表述,定义$ v={v}_{{t}_{0}},X=\mathrm{ln}({S}_{{t}_{0}}/K) $;在风险中性测度下,$ y=\mathrm{ln}({S}_{T}/K) $的条件特征函数为

$ {g}_{y}\left(\varphi ,\tau ;X,v,{t}_{0}\right) ={\mathrm{E}}^{\mathrm{Q}}\left[{\mathrm{e}}^{\mathrm{i}\varphi y}|{X}_{{t}_{0}}=X,{v}_{{t}_{0}}=v\right]$ (8)

式中:$ {\mathrm{E}}^{\mathrm{Q}}[\cdot ] $为风险中性测度下的条件期望,$ \mathrm{i}=\sqrt{-1} $是虚数单位,$ \varphi $为频率。

通过Feynman-Kac定理推导$ {g}_{y} $的偏微分方程为

$\begin{split} &\frac{1}{2}v\frac{{\partial }^{2}g}{\partial {X}^{2}}+\rho \sigma v\frac{{\partial }^{2}g}{\partial X\partial v}+\frac{1}{2}{\sigma }^{2}v\frac{{\partial }^{2}g}{\partial {v}^{2}}+\\&\left(r-\frac{1}{2}v\right) \frac{\partial g}{\partial X}+\kappa \left(\theta -v\right) \frac{\partial g}{\partial v}-\frac{\partial g}{\partial \tau }=0 \end{split}$ (9)

边界条件为:$ {g}_{y}(\varphi ,0;y,{v}_{T},T) ={\mathrm{e}}^{\mathrm{i}\varphi y} $

假设特征函数的形式为

$ {g}_{y}\left(\varphi ,\tau ;X,v,{t}_{0}\right) ={\mathrm{e}}^{\mathrm{i}\varphi X+A\left(\tau \right) +B\left(\tau \right) v}$ (10)

将式(10) 代入式(9) ,并结合边界条件$ A\left(0\right) =0, B\left(0\right) =0 $,计算可得

$ B\left(\tau \right) =\frac{d-b}{2a}\left(\frac{1-{\mathrm{e}}^{\tau d}}{1+\dfrac{d-b}{d+b}{\mathrm{e}}^{\tau d}}\right) $ (11)
$ A\left(\tau \right) =\mathrm{i}r\tau \varphi +\frac{m}{a}\mathrm{ln}(\frac{2d}{h}) +\frac{m}{2a}(d-b) \tau $ (12)

式中:$a=\dfrac{1}{2}{\sigma }^{2};b\;=\;\rho \sigma \mathrm{i}\varphi \;-\;\kappa \;;\;c\;=\;\dfrac{1}{2}\mathrm{i}\varphi (\mathrm{i}\varphi -1) ;d\;=\; - \sqrt{{b}^{2}-4ac};h=b+d-b{\mathrm{e}}^{d\tau }+d{\mathrm{e}}^{d\tau };m=\theta \kappa$

根据Gil-Pelaez[39]的反演公式,可得$ y=\mathrm{ln}({S}_{T}/K) $的概率密度函数和累积分布函数为

$ {f}_{y}=\frac{1}{\mathrm{\pi }}{\displaystyle\int }_{0}^{+\mathrm{\infty }}{R}({\mathrm{e}}^{-\mathrm{i}\varphi y}{g}_{y}\left(\varphi \right) ) \mathrm{d}\varphi $ (13)
$ F\left(y\right) =\frac{1}{2}-\frac{1}{\mathrm{\pi }}{\displaystyle\int }_{0}^{+\mathrm{\infty }}{I}(\frac{{\mathrm{e}}^{-\mathrm{i}\varphi y}{g}_{y}\left(\varphi \right) }{\varphi }) \mathrm{d}\varphi $ (14)

式中:${R}$${I}$分别为函数的实部和虚部。

显然,式(13) 和式(14) 均涉及到无穷积分,不可避免地增加了计算负担。本文参考Witkovsky[40]的思路,应用级数展开法计算式(13) 和式(14) 的近似逼近值:

$ \begin{split}{f}_{y}\approx &\frac{{\delta }_{t}}{\mathrm{\pi }}(\frac{1}{2}+\displaystyle\sum\nolimits_{j=1}^{N-1}\mathrm{cos}\left({\varphi }_{j}y\right) {g}_{y}\left({\varphi }_{j}\right) +\\ &\frac{1}{2}\mathrm{cos}\left({\varphi }_{N}y\right) {g}_{y}\left({\varphi }_{N}\right) \end{split}$ (15)
$ \begin{split}{F}_{y}\approx& \frac{1}{2}-\frac{{\delta }_{t}}{\mathrm{\pi }}(-\frac{y}{2}+\displaystyle\sum\nolimits _{j=1}^{N-1}\mathrm{sin}\left({\varphi }_{j}y\right) \frac{{g}_{y}\left({\varphi }_{j}\right) }{{\varphi }_{j}}+\\ &\frac{1}{2}\mathrm{sin}\left({\varphi }_{N}y\right) \frac{{g}_{y}\left({\varphi }_{N}\right) }{{\varphi }_{N}}\end{split}$ (16)

式中:$ N $为一个足够大的常数,将$\varphi$的截断区间$ (0,L) $分割为大小为$ {\delta }_{t}=L/N $$ N $个子区间,${\varphi }_{j}$为子区间的第$ j $个子节点,$ L $的大小由6-sigma法则确定。

综上,结合式(3)、式(4)、式(6)、式(15)和式(16)可计算期权投资者的累积前景价值。

1.3 Heston行为期权定价公式

Thaler提出的心理账户理论[29](Mental Account) 认为,投资者在计算损益时,倾向于将收益与损失放置到互相分离的心理账户里分别计算其结果。具体到期权市场而言,期权卖方 (买方) 可能将合约签订时权利金的收入 (支出) 、行权时的报酬 (payoff) 支出 (收入) 划分到不同的心理账户。因此结合心理账户理论与累积前景理论,卖方(writer) 视角的看涨期权前景价值为

$ \begin{split} {V}^{\mathrm{w}}=&{v}^+\left({c}^{\mathrm{w}}{\mathrm{e}}^{rT}\right) +\\[-2pt] &{\displaystyle\int }_{0}^{+\mathrm{\infty }}{\psi }^-\left(1-F\left(y\right) \right) f\left(y\right) {v}^-\left[K\right(1-{\mathrm{e}}^{y}\left) \right]\mathrm{d}y \end{split} $ (17)

式中:$ {c}^{\mathrm{w}} $是看涨期权的卖出价,等式右边第一项为权利金的时间价值,第二项为合约到期时的预期支付。

在市场均衡时,前景价值为0,计算得到看涨期权卖方价格的表达式为

$ {c}^{\mathrm{w}}={\mathrm{e}}^{-rT}(\lambda {\displaystyle\int }_{0}^{+\mathrm{\infty }}{\psi }^-\left(1-F\left(y\right) \right) f\left(y\right) \times (K\left({\mathrm{e}}^{y}-1\right) {) }^{\beta }\mathrm{d}y{) }^{\tfrac{1}{\alpha }} $ (18)

同理,买方(holder) 的看涨期权前景价值为

$ \begin{split}{V}^{\mathrm{h}}=&{v}^-\left(-{c}^{\mathrm{h}}{\mathrm{e}}^{rT}\right) +\\ &{\displaystyle\int }_{0}^{+\mathrm{\infty }}{\psi }^+\left(1-F\left(y\right) \right) f\left(y\right) {v}^-\left(K\left({\mathrm{e}}^{y}-1\right) \right) \mathrm{d}y \end{split}$ (19)

式中:$ {c}^{\mathrm{h}} $是看涨期权的买入价。令式(19) 等于0,可得看涨期权买方价格的表达式为

$ {c}^{\mathrm{h}}={\mathrm{e}}^{-rT}\left(\frac{1}{\lambda }{\displaystyle\int }_{0}^{+\mathrm{\infty }}{\psi }^+\left(1-F\left(y\right) \right) f\left(y\right) (K\left({\mathrm{e}}^{y}-1\right) {) }^{\alpha }\mathrm{d}y\right)^{\tfrac{1}{\beta }} $ (20)

同理,可得看跌期权卖方价格和买方价格的表达式为

$ {p}^{\mathrm{w}} = {\mathrm{e}}^{-rT}(\lambda {\displaystyle\int }_{-\mathrm{\infty }}^{0}{\psi }^-\left(F\left(y\right) \right) f\left(y\right) \times (K\left(1-{\mathrm{e}}^{y}\right) {) }^{\beta }\mathrm{d}y{) }^{\textstyle\frac{1}{\alpha }} $ (21)
$ {p}^{\mathrm{h}} = {\mathrm{e}}^{-rT}(\frac{1}{\lambda }{\displaystyle\int }_{-\mathrm{\infty }}^{0}{\psi }^+\left(1-F\left(y\right) \right) f\left(y\right) \times (K\left(1-{\mathrm{e}}^{y}\right) {) }^{\alpha }\mathrm{d}y{) }^{\textstyle\frac{1}{\beta }} $ (22)

至此,本文在市场均衡条件下推导出Heston模型下欧氏行为期权(看涨和看跌) 的买方和卖方价格。

2 数值定价

本节运用数值定价方法展示Heston行为期权定价公式的性能,并与理性框架下Heston模型定价结果作对比分析。

设定Heston模型的参数:$ S = 3.5 $$ T = 0.1 $$ r = 0.03 $$ {v}_{0}=0.05 $$ \theta =0.1 $$ \rho =-0.5 $$ \kappa =0.03 $$ \sigma =0.3 $

Tversky-Kahneman (TK) 情绪下行为参数为:$ \alpha =\beta =0.88 $$ \lambda =2.25 $$ {\gamma }^{+}=0.61 $$ {\gamma }^{-}=0.69 $

温和情绪下行为参数为:$ \alpha =\beta =0.988 $$ \lambda =0.9 $$ {\gamma }^{+}=0.961 $$ {\gamma }^{-}=0.969 $

对于理性框架下的Heston期权定价模型,由于其不考虑投资者的心理与行为因素,在本文框架中相当于设置了行为参数:$ \alpha =\beta =\lambda =\gamma =1 $

数值定价结果如表1所示,其中,$ {c}_{\mathrm{H}} $$ {c}_{\mathrm{H}\mathrm{B}}^{\mathrm{h}} $$ {c}_{\mathrm{H}\mathrm{B}}^{\mathrm{w}} $分别代表理性框架下Heston模型的期权定价结果、Heston行为期权买方价格、Heston行为期权卖方价格。

表 1 数值定价结果 Table 1 Numerical pricing results

表1的数值结果可以看出,无论在TK情绪或温和情绪下,对于每个行权价,买方视角下的期权价格都低于卖方视角下的期权价格,理性框架下的Heston模型定价结果则位于二者之间。由此可见,考虑了投资者情绪特征的Heston行为期权定价模型能较好地解释现实中的买卖价差情况。

但TK情绪下的行为期权定价结果与理性框架下Heston模型偏差很大,温和情绪下的定价结果则较为相近。这表明Tversky-Kahneman情绪下行为参数可能是一组实验性参数,目的只是为了解释前景理论,并没有真正用于期权定价实践。

3 上证50ETF期权的实证分析 3.1 数据选取

本文选取到期时间为2022年3月23日的上证50ETF认购期权数据进行实证研究,样本区间为2022年1月4日至2022年3月15日,以检验Heston行为期权定价模型的有效性。由于到期期限过短的期权流动性较差、价格波动大,因此在定价实验中剔除到期期限小于7天的期权数据。考虑到Heston模型对中、长到期的期权定价已经较为准确,本文仅关注短到期期权数据。本文的期权数据来源于Wind终端,无风险利率采用一年期Shibor,数据来源于http://www.shibor.org/

3.2 模型参数估计

参照金融业界的做法,对实值期权$ (S/K > 1) $和虚值期权$ (S/K < 1) $分别进行参数估计。临近到期时,利用预测日前一天的上证50ETF期权数据所得的模型参数,来估计预测日当天的期权价格,并选用MATLAB工具箱里的非线性最小二乘优化Isqnonlin求解器估计参数值。篇幅所限,仅展示的Heston行为期权定价模型买方视角的参数估计结果。模型的待估参数为

$ \varTheta =\left\{{v}_{0},\theta ,\rho ,\kappa ,\sigma ,\alpha ,\beta ,\lambda ,\gamma \right\} $ (23)

通过求解如下的最优化问题得到Heston行为期权定价模型的参数估计结果:

$ \mathrm{min}(\displaystyle\sum \left|({C}^{\mathrm{m}\mathrm{o}\mathrm{d}\mathrm{e}\mathrm{l}}-{C}^{\mathrm{m}\mathrm{a}\mathrm{r}\mathrm{k}\mathrm{e}\mathrm{t}}) /{C}^{\mathrm{m}\mathrm{a}\mathrm{r}\mathrm{k}\mathrm{e}\mathrm{t}}) \right|) $ (24)

表2表3分别展示了实值期权数据和虚值期权数据参数校正后的描述性统计结果。可以发现,无论是实值期权,还是虚值期权,由它们校准出来Heston模型基础参数结果大体一致,这是因为Heston基础模型刻画的是同时间段、相同标的资产价格的行为;但扩展的行为参数校准结果是有较大差异的。基于实值期权数据校正得到的行为参数在1附近波动,且波动较小,反映出对于实值期权,投资者基本呈现风险中性态度;而基于虚值期权数据校正得到的行为参数偏离1较远,整体波动较大,反映出投资者情绪态度变化较大,也进一步应证了将实值期权数据和虚值期权数据分开校准的合理性。

表 2 实值期权数据的参数描述性统计结果 Table 2 Descriptive statistic results of ITM's parameters
表 3 虚值期权数据的参数描述性结果 Table 3 Descriptive statistic results of OTM's parameters
3.3 模型定价精度对比

本节展示了Heston行为期权定价模型同传统理性框架下Heston期权定价模型和GBM行为期权定价模型在定价性能的比较结果。以平均绝对百分比误差 (Mean Absolute Percentage Error, MAPE) 来评判模型的整体定价精度。

$ {\rm{MAPE}}=\frac{1}{N}{\displaystyle\sum\nolimits _{i=1}^{N}}\left(\frac{\left|{C}_{i}^{\mathrm{m}\mathrm{o}\mathrm{d}\mathrm{e}\mathrm{l}}-{C}_{i}^{\mathrm{m}\mathrm{a}\mathrm{r}\mathrm{k}\mathrm{e}\mathrm{t}}\right|}{{C}_{i}^{\mathrm{m}\mathrm{a}\mathrm{r}\mathrm{k}\mathrm{e}\mathrm{t}}}\right) $ (25)

对于实值期权数据和虚值期权数据,分别进行参数校准后计算出期权价格,模型的整体定价精度如表4所示:对于实值期权,3个模型均有着较好的定价精度,Heston行为期权定价模型的定价精度相较于传统理性框架下Heston模型提升较小,结合模型参数校正结果,说明行为参数的引入对ITM期权定价精度提升有限;而对于虚值期权,Heston行为期权定价模型的定价精度则要远远强于GBM行为期权定价模型和传统理性框架下的Heston模型,不仅说明了将GBM模型替换为Heston模型的合理性,也表明了在OTM期权定价中考虑投资者非理性情绪与行为的必要性。

表 4 模型定价精度对比 Table 4 Model pricing accuracy comparison
4 结论

本文结合Heston随机波动率模型和前景理论期权定价框架,将期权合约签订时权利金的支付、行权时报酬 (payoff) 的支付划分到不同的心理账户下,依据市场均衡条件,推导了Heston行为期权定价公式。基于上证50ETF期权的实证结果发现,相较于传统理性框架下的Heston模型,Heston行为期权定价模型显著提升短到期OTM期权的定价精度,而对ITM期权的定价精度则提升较小;这是因为对于OTM期权,投资者风险态度、情绪变化较为突出,而行为参数的引入能较好地捕捉这些特征,而对于ITM期权,投资者风险态度偏中性,行为参数的引入对模型性能提升有限。

本文所得到的欧氏期权的良好定价性能,一方面归功于Heston随机波动率模型刻画标的资产价格行为的准确性,另一方面也有赖于行为参数的引入,使模型能捕捉到市场情绪的变化,更好地刻画现实市场。

本文的研究结果表明:基于前景理论的期权定价框架可作为一种风险中性期权定价方法的替代。未来可进一步研究更多样化的心理账户情形,在标的资产价格运动建模中加入非理性因子等。

参考文献
[1]
BLACK F, SCHOLES M. The pricing of options and corporate liabilities[J]. Journal of Political Economy, 1973, 81(3): 637-654. DOI: 10.1086/260062.
[2]
MERTON R C. Option pricing when underlying stock returns are discontinuous[J]. Journal of Financial Economics, 1976, 3(1-2): 125-144. DOI: 10.1016/0304-405X(76)90022-2.
[3]
BATES D S. Jumps and stochastic volatility: exchange rate processes implicit in deutsche mark options[J]. The Review of Financial Studies, 1996, 9(1): 69-107. DOI: 10.1093/rfs/9.1.69.
[4]
HULL J. An analysis of the bias in option pricing caused by a stochastic volatility[J]. Advances in Futures and Options Research, 1988, 3: 29-61.
[5]
HESTON S L. A closed-form solution for options with stochastic volatility with applications to bond and currency options[J]. The Review of Financial Studies, 1993, 6(2): 327-343. DOI: 10.1093/rfs/6.2.327.
[6]
GARLEANU N, PEDERSEN L H, POTESHMAN A M. Demand-based option pricing[J]. The Review of Financial Studies, 2008, 22(10): 4259-4299.
[7]
JONES C S, SHEMESH J. Option mispricing around nontrading periods[J]. The Journal of Finance, 2018, 73(2): 861-900. DOI: 10.1111/jofi.12603.
[8]
HAN B. Investor sentiment and option prices[J]. The Review of Financial Studies, 2008, 21(1): 387-414. DOI: 10.1093/rfs/hhm071.
[9]
SIDDIQI H. Anchoring-adjusted option pricing models[J]. Journal of Behavioral Finance, 2019, 20(2): 139-153. DOI: 10.1080/15427560.2018.1492922.
[10]
孙有发, 姚宇航, 邱梓杰, 等. 基于特征函数局部结构微扰法的行为期权定价研究[J]. 系统工程理论与实践, 2022, 42(12): 3247-3264.
SUN Y F, YAO Y H, QIU Z J, et al. Behavioral option pricing method based on perturbation in the local structure of characteristic function[J]. Systems Engineering-Theory & Practice, 2022, 42(12): 3247-3264. DOI: 10.12011/SETP2021-1685.
[11]
JANG B G, KIM C, KIM K T, et al. Psychological barriers and option pricing[J]. Journal of Futures Markets, 2015, 35(1): 52-74.
[12]
SONG S, WANG Y. Pricing double barrier options under a volatility regime-switching model with psychological barriers[J]. Review of Derivatives Research, 2017, 20: 255-280. DOI: 10.1007/s11147-017-9130-x.
[13]
SONG S, WANG G, WANG Y. Pricing European options under a diffusion model with psychological barriers and leverage effect[J]. The European Journal of Finance, 2020, 26(12): 1184-1206. DOI: 10.1080/1351847X.2020.1725083.
[14]
孙有发, 邱梓杰, 姚宇航, 等. 基于深度学习算法的行为期权定价[J]. 系统管理学报, 2021, 30(4): 697-708.
SUN Y F, QIU Z J, YAO Y H, et al. Behavioral option pricing based on deep-learning algorithm[J]. Journal of Systems & Management, 2021, 30(4): 697-708.
[15]
RACHEV S, STOYANOV S, FABOZZI F J. Behavioral finance option pricing formulas consistent with rational dynamic asset pricing[EB/OL]. arXiv:1710. 03205(2017-10-09) [2022-11-30]. https: //doi. org/10.48550/arXiv. 1710. 03205.
[16]
RACHEV S, FABOZZI F J, RACHEV B, et al. Option pricing with greed and fear factor: The rational finance approach[EB/OL]. arXiv: 1709.08134(2017-09-24) [2022-11-30]. https: //doi. org /10. 48550/arXiv. 1709.08134.
[17]
DAI J, SHIRVANI A, FABOZZI F J. Rational finance approach to behavioral option pricing[EB/OL]. arXiv: 2005.05310(2020-05-10) [2022-11-30]. https: //doi. org/10.48550/arXiv. 2005. 05310.
[18]
BAELE L, DRIESSEN J, EBERT S, et al. Cumulative prospect theory, option returns, and the variance premium[J]. The Review of Financial Studies, 2019, 32(9): 3667-3723. DOI: 10.1093/rfs/hhy127.
[19]
姜继娇, 杨乃定. 基于两心理账户BPT的复合实物期权定价模型[J]. 管理科学学报, 2008, 11(1): 89-94.
JIANG J J, YANG N D. Compound real option pricing model based on BPT with two mental accounts[J]. Journal of Management Sciences in China, 2008, 11(1): 89-94. DOI: 10.3321/j.issn:1007-9807.2008.01.010.
[20]
SHEFRIN H, STATMAN M. Behavioral aspects of the design and marketing of financial products[J]. Financial Management, 1993, 22(2): 123-134. DOI: 10.2307/3665864.
[21]
SHILLER R J. Human behavior and the efficiency of the financial system[J]. Handbook of Macroeconomics, 1999, 1: 1305-1340.
[22]
BREUER W, PERST A. Retail banking and behavioral financial engineering: the case of structured products[J]. Journal of Banking & Finance, 2007, 31(3): 827-844.
[23]
SUZUKI K, SHIMOKAWA T, MISAWA T. Agent-based approach to option pricing anomalies[J]. IEEE Transactions on Evolutionary Computation, 2009, 13(5): 959-972. DOI: 10.1109/TEVC.2008.2011745.
[24]
BAR-GILL O. Pricing legal options: a behavioral perspective[J]. Review of Law & Economics, 2005, 1(2): 204-240.
[25]
BENARTZI S, THALER R H. Myopic loss aversion and the equity premium puzzle[J]. The Quarterly Journal of Economics, 1995, 110(1): 73-92. DOI: 10.2307/2118511.
[26]
VERSLUIS C, LEHNERT T, WOLFF C C P. A cumulative prospect theory approach to option pricing[EB/OL]. (2010-11-30) [2022-11-30]. https: //ssrn. com/ abstract=1717015.
[27]
KAHNEMAN D, TVERSKY A. Prospect theory: an analysis of decision under risk[J]. Econometrica, 1979, 47(2): 363-391.
[28]
TVERSKY A, KAHNEMAN D. Advances in prospect theory: cumulative representation of uncertainty[J]. Journal of Risk and Uncertainty, 1992, 5: 297-323. DOI: 10.1007/BF00122574.
[29]
THALER R H. Mental accounting and consumer choice[J]. Marketing Science, 1985, 4(3): 199-214. DOI: 10.1287/mksc.4.3.199.
[30]
THALER R H. Mental accounting matters[J]. Journal of Behavioral Decision Making, 1999, 12(3): 183-206. DOI: 10.1002/(SICI)1099-0771(199909)12:3<183::AID-BDM318>3.0.CO;2-F.
[31]
PENA A, ALEMANNI B, ZANOTTI G. On the role of behavioral finance in the pricing of financial derivatives: the case of the S&P 500[EB/OL]. (2010-08-26) [2022-11-30]. https://ssrn.com/abstract=1798607.
[32]
NARDON M, PIANCA P. A behavioral approach to the pricing of European options[M] //Mathematical and Statistical Methods for Actuarial Sciences and Finance. Switzerland: Springer International Publishing, 2014: 219-230.
[33]
NARDON M, PIANCA P. European option pricing under cumulative prospect theory with constant relative sensitivity probability weighting functions[J]. Computational Management Science, 2019, 16: 249-274. DOI: 10.1007/s10287-018-0324-y.
[34]
吴鑫育, 李心丹, 马超群. 基于随机波动率模型的上证50ETF期权定价研究[J]. 数理统计与管理, 2019, 38(1): 115-131.
WU X Y, LI X D, MA C Q. The pricing of SSE 50 ETF options based on stochastic volatility model[J]. Journal of Applied Statistics and Management, 2019, 38(1): 115-131. DOI: 10.13860/j.cnki.sltj.20181031-002.
[35]
孙有发, 吴碧云, 郭婷, 等. 非仿射随机波动率模型下的50ETF期权定价: 基于Fourier-Cosine方法[J]. 系统工程理论与实践, 2020, 40(4): 888-904.
SUN Y F, WU B Y, GUO T, et al. Fourier-cosine based option pricing for SSE 50ETF under non-affine stochastic volatility model[J]. Systems Engineering-Theory & Practice, 2020, 40(4): 888-904. DOI: 10.12011/1000-6788-2018-2403-17.
[36]
SUN Y F. Efficient pricing and hedging under the double Heston stochastic volatility jump-diffusion model[J]. International Journal of Computer Mathematics, 2015, 92(12): 2551-2574. DOI: 10.1080/00207160.2015.1079311.
[37]
DAVIES G B, SATCHELL S E. The behavioural components of risk aversion[J]. Journal of Mathematical Psychology, 2007, 51(1): 1-13. DOI: 10.1016/j.jmp.2006.10.003.
[38]
PRELEC D. The probability weighting function[J]. Econometrica, 1998, 66(3): 497-527. DOI: 10.2307/2998573.
[39]
GIL-PELAEZ J. Note on the inversion theorem[J]. Biometrika, 1951, 38(3-4): 481-482. DOI: 10.1093/biomet/38.3-4.481.
[40]
WITKOVSKY V. Numerical inversion of a characteristic function: an alternative tool to form the probability distribution of output quantity in linear measurement models[J]. Acta IMEKO, 2016, 5(3): 32-44.