董艷

（陜西鐵路工程職業(yè)技術學院基礎部，陜西渭南714000）

?

非線性Black-Scholes模型下Bala期權定價

董艷

（陜西鐵路工程職業(yè)技術學院基礎部，陜西渭南714000）

在非線性Black-Scholes模型下，研究了Bala期權定價問題.首先利用雙參數(shù)攝動方法，將Bala期權適合的偏微分方程分解成一系列常系數(shù)拋物方程.其次通過計算這些常系數(shù)拋物型方程的解，給出了Bala期權的近似定價公式.最后利用Green函數(shù)分析了近似結論的誤差估計.

Bala期權；非線性Black-Scholes模型；Green函數(shù)；誤差估計

§1　引　言

近些年來有關障礙期權的研究工作已有一些文獻，文獻［1］研究了分數(shù)布朗運動環(huán)境下歐式障礙期權定價問題，利用偏微分方程方法給出了歐式障礙期權的顯示解，并給出了看張看跌之間的平價關系.文獻［2］和文獻［3］則是從數(shù)值模擬的角度出發(fā)，利用蒙特卡羅模擬研究了障礙期權定價問題.胡文偉和李湛在文獻［4］中采用二叉樹方法研究了障礙期權定價問題.文獻［5］則放棄了原生資產(chǎn)波動率和收益率的常數(shù)假設，在較為一般的模型基礎之上，研究了障礙期權定價問題，將障礙期權滿足的拋物初值問題做近似化處理，給出了障礙期權的一個近似解析表達式，并給出了相應的誤差估計.

遺憾的是上述文獻均是研究普通障礙期權所得結論，近些年來，越來越多的學者發(fā)現(xiàn)了障礙期權存在的某些不足:對于敲出障礙期權的持有者而言，當原生資產(chǎn)價格觸及到障礙值時會導致期權立刻失效，他就失去了全部的投資.為了克服這一不足，金融實踐家們衍生出了Bala期權.

下降敲出類型的Bala期權設置了另一種緩沖措施，該期權合約設置了一個“期權失效計時器”.當原生資產(chǎn)價格St在障礙值B的下方停留總時間達到規(guī)定的時長D后，期權全部失效，否則期權轉化為歐式看跌期權.因此期權在到期日T時的收益為

（K - ST）+I｛τT＜D｝，

其中τt的定義見公式（2）.易得當D→0時，下降敲出類型的Bala期權就轉化為下降敲出看跌障礙期權；當時，該Bala期權則轉化為標準的歐式看跌期權.除此Bala期權還有很多類型，本文只研究其中的一種情況，至于其他Bala期權的相關結論可類推本文證明.

目前近些年來有關Bala期權定價還沒有專門的文獻，姜禮尚在文獻［6］中給出了Bala期權適合的偏微分方程模型，并給出了相應的優(yōu)先差分格式.由于Bala期權適合的偏微分方程結構復雜（它甚至都不是拋物方程），以目前的研究手段尋找Bala期權的精確解析解是困難的.

既然對于Bala期權不能得到其精確解析定價公式，那么找到一個數(shù)學上一致性和嚴密的近似解也不失為一個很好的研究方向.基于此，本文在混合分數(shù)布朗運動驅動的非線性Black-Scholes模型下研究階梯期權定價問題，其中原生資產(chǎn)價格St滿足非線性隨機微分方程

無風險資產(chǎn)價格Mt滿足

其中S0＞0為原生資產(chǎn)的初始價格，wt為標準布朗運動，μ（t，St）表示風險資產(chǎn)的期望收益率，σ0（t，St；ε）表示風險資產(chǎn)的波動率，ε（0＜ε＜1）為攝動參數(shù)，μ（t，St）和σ0（t，St；ε）為時間t和股票價格St的一般函數(shù).

如果用常數(shù)r表示無風險資產(chǎn)的利率，P（t，τ，S）表示下降敲出類型的看跌Bala期權的價格，用C（t，τ，S）表示相應的看漲Bala期權的價格，由文獻［6］易得P（t，τ，S）和C（t，τ，S）為如下兩個偏微分方程初值問題的解

其中算子L為

接下來，將采用多尺度方法給出兩種階梯期權的定價公式，并采用Feymann-Kac公式給出結論的誤差估計.本文提供的方法同樣適合其他期權，例如歐式期權、亞式期權以及回望期權.

§2　多尺度方法

為方便敘述，首先定義如下的符號

其中

容易看出公式（5）提供了一個由P（t，τ，S）到V（t，τ，x）的變換，并且在該變換下偏微分方程問題（1）變?yōu)?/p>

假定當ε= 0時f（t，x；ε）為與t和x無關的常數(shù)，即f（t，x；0）= f0（0，0），并且f（t，x；ε）可以根據(jù)泰勒公式在ε= 0附近展成如下形式的冪級數(shù)

fn（t，x）表示f（t，x；ε）對ε的n階導數(shù)在ε= 0處的值.由于f（t，x；ε）=σ2（t，ex；ε）＞0，所以這里假設f0（0，0）＞0.

下面以公式（6）為研究對象采用多尺度方法（雙參數(shù)攝動展開方法）去處理Bala期權定價問題.根據(jù)多尺度方法的中心思想［7-8］，將V（t，τ，x）在（ε，η）=（0，0）附近展開成冪級數(shù)

并將上述級數(shù)代入偏微分方程初值問題（6）的主方程，有

比較ε和η同次冪的系數(shù)，可以得到

由x和t的任意性，可令ε和η的相同次冪項的系數(shù)為零，從而

其中

易見，當（n，m）=（0，0）時，gn，m（t，τ，x）= 0，同時經(jīng)過整理gn，m（t，τ，x）也可以寫成如下形式

最后根據(jù)公式（8）將偏微分方程初值問題（6）的初值條件進行分解，則拋物型方程（6）的求解問題轉化為一系列常系數(shù)拋物方程初值問題，并且這些常系數(shù)拋物初值問題是可解的，其中V0，0（t，τ，x）為拋物方程初值問題

的解，Vn，m（t，τ，x）是如下拋物方程初值問題的解

下面依次計算Vn，m，n = 0，1，2，···，m = 0，1，2，···.另一方面，觀察拋物初值問題（12）和（13）可以看出，主方程中不含變量τ，它只作為參量出現(xiàn)在初值條件中，所以在以后的證明過程中將其視為常數(shù).此外，公式（11）中gn，m（t，τ，x）不包含Vi，j（i＜n，j＜m）這樣的項，這是最后一節(jié)能進行誤差分析的關鍵.

§3　期權定價公式

引理1若V0，0（t，τ，x）是常系數(shù)拋物方程初值問題（12）的解，則它可表示為

其中N（x）表示標準正態(tài)分布的概率累積函數(shù)，

證對于首項V0，0（t，τ，x），首先考慮如下的變換

其中

若將τ視為常數(shù)，則常系數(shù)拋物初值問題（12）在變換（14）下可以寫為

從而利用Poisson公式可知，拋物初值問題（15）的解為

其中

將初值條件帶入并進行積分運算，可得

利用變換（14）回到可得定理結論.

接下來求解Vn，m（t，τ，x）（n，m = 0，1，2，···）.觀察公式（11）可知，gn，m（t，τ，x）僅僅與

有關，而與?xxVn，m-?xVn，m，0＜η＜1無關，這啟發(fā)采用遞推法完成如下結論的證明.

引理2線性拋物方程（13）的解可表示為

其中

證注意到計算Vn，m（t，τ，x）時，gn，m（t，τ，x）是已知的，因此將τ視為常數(shù)，并作變換

則拋物方程初值問題（30）可以轉化為如下形式

注意到gn，m（t，τ，x）是解析的以及?gn，m（s，τ，x）= exp｛-αx-βs｝gn，m（T -s，τ，x），根據(jù)文獻［9-10］，拋物初值問題（18）的解可表示為

對變換（17）進行逆變換可完成定理證明.

觀察引理1，引理2以及公式（11），易得

這里δ（·）為狄利克雷函數(shù)，將上式代入gn，m（t，τ，x）的表達式有如下結論成立.

定理3若P（t，S）表示下降敲出類型的Bala看跌期權的價格，則它可以表示為

其中

Vn，m（t，τ，x）見引理2.

定理4若C（t，τ，S）表示下降敲出類型的Bala看漲期權的價格，則它可以表示為

其中

Vn，m（t，τ，x）的表示方法與定理3相同.

觀察定理3和定理4容易驗證，當τ＞D時，

從而當τ＞D時，有C（t，τ，S）= P（t，τ，S）= 0，這恰與文獻［6］描述的事實相符.

§4　誤差估計

本節(jié)，將在如下假設的基礎之上進一步考慮近似結論的誤差估計問題:

（2）假設fn（t，x）在［0，T］×R+上一致有界，即對任意的（t，x）∈［0，T］×R+以及任意正整數(shù)n，存在正常數(shù)M0，使得

并假定對任意的（t，x）∈［0，T］×R+有

上面的假設也意味著f（t，x；ε）和f（t，x；ε）均滿足Lipschtiz條件和線性增長條件.事實上，由公式（7）和公式（23）可知對任意的（t，x）∈［0，T］×R有

顯然f（t，x；ε）和f（t，x；ε）在（t，x）∈［0，T］×R上關于x滿足線性增長條件.由公式（43）可得

因此f（t，x；ε）也關于變量x滿足一致Lipschtiz條件.

為了方便證明，還需要兩個有關Rn，m（t，τ，x）=?xxVn，m-?xVn，m的引理，在以后的證明中經(jīng)常用到它們.

引理5對任意的（t，τ，x）∈［0，T］×［0，T］×R，存在不依賴t，τ和x的正常數(shù)M1使得

證由引理1直接計算容易得到

進一步，利用公式（27）和公式（28）以及引理1有

將公式（29）和公式（30）代入R0，0（t，τ，x）可得

容易驗證對任意的（t，τ，x）∈［0，T］×［0，T］×R都有

因此exp｛x｝N′（a2）在［0，T］×R上有界，即存在正常數(shù)M3使得

進一步，由文獻［10］第四章第十六節(jié)可知，存在正常數(shù)M4和M5，使得

將公式（32）代入公式（31）得到

因此公式（26）中第二個不等式成立.

引理6當T足夠小時，對任意的非負整數(shù)n和m（n2+ m2＞0），存在不依賴x和t的非負常數(shù)M6和M7，使得

證明可類推文獻［11］之引理2.6，只需用引理3代替文獻［11］之引理2.4.

值得注意的是，上述證明得以進行是因為構造的多尺度方法使得gn，m（t，τ，x）的表達式可以寫成Vn，j（t，τ，x）和Vi，m（t，τ，x）的組合，其表達式并不包含Vi，j（t，τ，x）這樣的項，0≤i＜n，0≤j＜m.

引理7對任意的（t，τ，x）∈［0，T］×［0，T］×R有

證由引理6可知

進一步聯(lián)立公式（39）可以得到

因此利用文獻［12］中定理11.1，可得

引理7表明: V（t，τ，x）對變量τ的偏導數(shù)與V（t，τ，x）成線性關系，這意味著Lε，ηV（t，τ，x）具有拋物算子性質，即

這說明算子Lε，η為拋物算子，因此本文可以根據(jù)拋物方程的Feyman-Kac公式研究V（t，τ，x）的近似解的誤差估計問題，下面給出本節(jié)的兩個主要結論.

定理8當T足夠小時，存在不依賴τ，x和t的正常數(shù)，使得

證考慮C（t，τ，S）和V（t，τ，x）的對應關系，這里只需證明

成立.令x = ln S，E0，0（t，τ，x）= V0，0（t，τ，x）- V（t，τ，x），顯然

又因為Lε，ηV（t，τ，x）= 0，L0V0，0（t，τ，x）= 0，所以

故聯(lián)立公式（34）和公式（35），E0（t，x）為偏微分方程初值問題

的解，其中

利用公式（33），上述偏微分方程初值問題可以被重新寫為

一方面，由引理5可得

又由公式（22），對任意的（t，x）∈［0，T］×R+，

注意到0＜η＜1，由中值定理可知，存在θ∈（0，η）使得

因此將公式（26）和公式（40）代入公式（37），對任意的（t，x）∈［0，T］×R+，h0，0（t，x）滿足

另一方面，運用Feymann-Kac公式可知［13］，拋物初值問題（38）的解可以表示為如下形式的條件概率問題

其中Xt適合隨機微分方程

t≤s≤T，Xt= x.因為公式（23）-公式（25）已經(jīng)證明f（t，x；ε）和f（t，x；ε）關于變量x滿足線性增長條件和一致Lipschtiz條件，所以由文獻［13］可知，隨機微分方程（44）的解是存在且唯一的.最后，將公式（42）代入公式（43）可得E（t，x）的一個估計

定理9假定T足夠小并且存在正整數(shù)ω0使得

證類似定理7，這里只需證明

其中

顯然En（t，τ，x）滿足如下的初值條件

又因為Lε，ηV（t，τ，x）= 0，

所以有

其中

故聯(lián)立公式（46）和公式（47），En（t，x）為如下拋物初值問題的解

為估計hn（t，x），利用公式（22），可得對任意的（t，τ，x）∈［0，T］×［0，T］×R有

又因為0＜η＜1，利用帶有Lagrange余項的泰勒公式

將公式（50）和公式（51）代入公式（48），hn（x，t）滿足如下形式的估計

進一步將引理5和引理6中的結論代入公式（52）可得

由Feynman-Kac公式可知，拋物初值問題（49）的解可表示為如下形式的條件概率問題

其中Xτ滿足隨機微分方程（66）.將公式（53）代入公式（54）可知，En（t，x）滿足估計式

［1］霍海峰，溫鮮，鄧國和.分數(shù)次布朗運動的歐式障礙期權定價［J］.經(jīng)濟數(shù)學，2009，26（4）: 97-103.

［2］徐騰飛，曹小龍，胡云姣.離散障礙期權定價的蒙特卡羅模擬［J］.北京化工大學學報（自然科學版），2013，40（3）: 123-127.

［3］張利花，張衛(wèi)國，許文坤.美式障礙期權定價的總體最小二乘擬蒙特卡羅模擬方法［J］.數(shù)理統(tǒng)計與管理，2013，32（5）: 923-930.

［4］胡文偉，李湛.基于二叉樹方法的障礙期權與標準期權差價分析模型［J］.上海交通大學學報，2012，46（5）: 825-831.

［5］孫玉東，師義民，吳敏.參數(shù)依賴股票價格情形下的障礙期權定價［J］.數(shù)學物理學報，2013，33（5）: 912-925.

［6］姜禮尚.期權定價的數(shù)學模型和方法（第二版）［M］.北京:高等教育出版社，2008.

［7］瓦西里耶娃，布圖索夫.奇異攝動方程解的漸近展開［M］.北京：高等教育出版社，2008.

［8］鄭連存，張欣欣.非線性偏微分方程近代分析方法［M］.北京：科學出版社，2011.

［9］Chen Yazhe. The second order parabolic partial differential equations［M］. Beijing: Peking University Press，2002.

［10］Ladyzenskaja O A，Solonikov V A，Uralceva N N. Linear and Quasilinear Equations of Parabolic Type［M］. New York: American Mathematical Society，1964.

［11］孫玉東.非線性Black-Scholes期權定價模型研究［D］.西安:西北工業(yè)大學，2014.

［12］歐陽光中，朱學炎，金福臨，等.數(shù)學分析第三版上冊［M］.北京:高等教育出版社，2015.

［13］金志明.隨機分析基礎［M］.北京:國防工業(yè)出版社，2001.

MR Subject Classification: 60H10；90A06

The pricing of Bala options under the nonlinear Black-Scholes model

DONG Yan
（Department of Basic，Shaanxi Railway Institute，Weinan 714000，China）

In this paper，the pricing problems of Bala options are studied under the nonlinear Black-Scholes model. Firstly，the partial differential equations for the Bala options are transformed into a series of parabolic equations with constant coefficients by the perturbation method of doubleparameter. Secondly，the approximate pricing formulae of the Bala options are given by solving those parabolic equations with constant coefficients. Finally，the error estimates of the approximate solutions are given by using Green function.

Bala options；nonlinear Black-Scholes model；Green function；error estimates

O211.6；F830.9

A

1000-4424（2016）01-0009-12

2015-03-22

2015-10-25

陜西鐵路工程職業(yè)技術學院基金（2015-08）