汪嘉骎

(廣西師范大學(xué)數(shù)學(xué)與統(tǒng)計學(xué)院，廣西桂林541004)

?

基于仿射跳擴散模型的利率衍生品定價

汪嘉骎

，鄧國和

(廣西師范大學(xué)數(shù)學(xué)與統(tǒng)計學(xué)院，廣西桂林541004)

本文考慮短期利率滿足一類仿射跳擴散期限結(jié)構(gòu)模型的利率衍生品定價。應(yīng)用Fourier變換方法和遠期測度技術(shù)，獲到了零息票債券及基于零息票債券為標(biāo)的資產(chǎn)的歐式債券期權(quán)價格的顯示解，并將此結(jié)果應(yīng)用于付息債券期權(quán)和利率期權(quán)的定價。最后，利用數(shù)值實例分別分析了債券和零息票債券期權(quán)價格受利率模型中各主要參數(shù)的影響，以及期權(quán)的隱含波動率問題。數(shù)值結(jié)果表明，跳躍風(fēng)險參數(shù)對利率衍生品價格和隱含波動率有顯著作用，這也驗證了仿射跳擴散的利率期限結(jié)構(gòu)模型具有較好擬合實際的能力。

仿射跳擴散模型；利率期限結(jié)構(gòu)；債券期權(quán)；Fourier變換；隱含波動率

0 引言

利率是金融市場中最基本和最重要的經(jīng)濟變量之一，許多金融衍生工具(例如股票型期權(quán)、信用衍生品等)的定價和設(shè)計從根本上來說都依賴于利率的變化過程，金融理論與實際應(yīng)用研究也離不開利率。利率模型(也稱利率期限結(jié)構(gòu)模型)是利率衍生品定價和風(fēng)險管理的重要工具，沒有利率理論模型，我們很難對許多實務(wù)問題做出合理的定價及尋找出理想的利率風(fēng)險對沖策略，更無法進行利率產(chǎn)品的金融創(chuàng)新設(shè)計。因此，對它的研究一直是金融數(shù)學(xué)與計量金融理論的熱點和難點。近年來，學(xué)者和業(yè)界投資者觀察到利率市場上債券的收益率曲線(yield curve)和隱含波動率曲面(implied volatility surface)出現(xiàn)了比以往更復(fù)雜的變化特征。這表明需要引入更復(fù)雜、彈性更好的數(shù)學(xué)模型以描述和研究利率的運動及其衍生產(chǎn)品價格的波動。

現(xiàn)有的許多經(jīng)典利率期限結(jié)構(gòu)模型，例如Vasicek[1]模型、CIR[2]模型、Brennan和Schwartz[3]模型以及其他擴展的模型[4-8]均假設(shè)狀態(tài)變量(例如短期利率、長期利率水平、波動率或其他)連續(xù)運動，且服從各種框架的連續(xù)形式擴散過程。盡管擴散過程在統(tǒng)計上有著良好的特性，并且計算方便，但越來越多的實證研究表明純擴散過程并不完全適合描述金融狀態(tài)變量(因為金融數(shù)據(jù)往往表現(xiàn)出尖峰厚尾特點，正態(tài)分布很難描述這類實際現(xiàn)象)。例如Hamilton[9]，Das[10]，Johannes[11]，Piazzesi[12]，Andersen、Benzoni和Lund[13]，Jarrow、Li和Zhao[14]等的研究發(fā)現(xiàn)利率動態(tài)過程具有明顯的跳躍特征存在。Lin和Yeh[15]研究了債券市場,發(fā)現(xiàn)跳躍擴散過程更加適合于描述利率過程。在定價和對沖金融衍生品時,跳躍擴散模型尤為重要,若忽視跳躍風(fēng)險部分,將會引起金融產(chǎn)品的定價和對沖風(fēng)險。因此，在描述利率的動態(tài)過程中引入跳躍風(fēng)險因素更為重要。

由于20世紀70年代西方各國放寬對利率的監(jiān)管，實行金融自由化政策，利率不確定性逐漸加劇,以致越來越多的金融機構(gòu)和投資者面臨利率風(fēng)險的沖擊。對沖和管理利率風(fēng)險的有效工具——利率衍生品(interest rate derivatives)開始逐漸產(chǎn)生，并在市場上受到了歡迎。利率衍生品是一種損益以某種方式依賴于利率水平變化的金融創(chuàng)新工具，市場上交易的產(chǎn)品類型很多，交易規(guī)模龐大，債券期權(quán)(bond options)就是一類重要的利率衍生工具，其買賣的標(biāo)的物是某一種債券(比如零息票債券、付息債券、違約債券等)。由于利率的動態(tài)變化，債券期權(quán)的定價比股票期權(quán)的定價更復(fù)雜。不同利率模型下的債券期權(quán)定價數(shù)學(xué)方法與計算也不盡相同。目前，關(guān)于債券期權(quán)的定價均是建立在利率滿足擴散過程驅(qū)動的期限結(jié)構(gòu)模型,研究成果豐富(如文獻[16])，但在利率的跳擴散模型下研究較少。Ahn和Thompson[17]首次將CIR[2]模型擴展為帶有跳躍部分的平方根利率過程(記為JCIR模型),并應(yīng)用線性化技術(shù)獲得了零息票債券的閉式解。Shirakawa[18]首次在利率跳擴散模型下開展債券期權(quán)定價研究，他們均假設(shè)短期利率具有有限個非隨機的跳躍振幅。后來，Baz和Das[19]在Vasicek[1]模型中加入跳躍部分擴展了純擴散利率過程(記為JVasicek模型)，并設(shè)定利率的跳躍振幅服從正態(tài)分布，應(yīng)用文獻[17]的近似技術(shù)獲得了零息票債券價格。Das和Foresi[20]也擴展了Vasicek[1]模型,并允許利率的跳躍振幅服從雙指數(shù)分布,獲得了零息票債券及債券期權(quán)的價格。Das[21-22]、Finnerty[23]推廣了Hull-White[4]模型,并設(shè)定利率的跳躍振幅服從正態(tài)分布或其他分布，同樣地獲得了零息票債券及債券期權(quán)的價格。最近，Beliaeva[24]、Nawalkha[25]和Deng[26]等在利率的跳擴散模型下研究了美式債券期權(quán)的定價方法。上述研究的一個共同點是假設(shè)利率動態(tài)過程滿足某個特定的跳擴散模型，在實際應(yīng)用過程中存在不足(比如，JVasicek模型可能導(dǎo)致利率出現(xiàn)負值，JCIR模型僅限于利率的正向跳躍風(fēng)險，等)。為了更加全面較準確地描述利率的動態(tài)行為特征，本文綜合上述利率的跳擴散模型，考慮一類更為廣泛的仿射跳擴散過程(此過程由Duffie等人[27]提出，在描述股價運動行為，及在股票型期權(quán)定價研究方面已有廣泛應(yīng)用)驅(qū)動利率的動態(tài)行為(記為JCIR-Ex模型)，該模型不僅融合了JVasicek和JCIR跳擴散模型的優(yōu)點，也克服了這兩個模型的不足，并且還能突現(xiàn)利率波動率的時變特征。本文的主要貢獻為: 首先，利用仿射跳擴散過程的特點建立短期利率行為的動力學(xué)模型，該模型中利率的波動率受內(nèi)部和外部因素的影響，且不同于JVasicek和JCIR跳擴散模型，并建立基于利率變量的任何歐式未定權(quán)益滿足的偏微分-積分方程，利用仿射結(jié)構(gòu)特點，獲得零息票債券價格的顯示解。其次，應(yīng)用Fourier反變換方法和遠期概率測度技術(shù)，研究歐式零息票債券期權(quán)的定價，并將結(jié)果應(yīng)用于付息債券期權(quán)和利率期權(quán)的定價，分別獲得這些期權(quán)價格計算的顯示公式，為金融風(fēng)險管理(如對沖策略的計算)提供便利。最后，給出數(shù)值計算實例，分析該利率模型中跳躍風(fēng)險因素對債券和債券期權(quán)價格的影響，以及債券期權(quán)的隱含波動率(即微笑現(xiàn)象)問題。

1 利率模型與債券價格

考慮一個無套利,無摩擦，且可連續(xù)進行交易的經(jīng)濟體，交易有限期限為τ=[0，T*]，T*<∞。該經(jīng)濟體中所有不確定性因素均定義在一個帶有信息濾波的給定完備概率空間(Ω，F(xiàn)，(Ft)t∈τ，Q)中，其中Q是風(fēng)險中性概率測度,信息濾波(Ft)t∈τ為滿足通常條件的完備參考族。假定在風(fēng)險中性概率Q下，瞬時利率rt的運動行為滿足下列仿射跳擴散模型：

(1)

為方便起見，記模型(1)為JCIR-Ex模型。該模型包括了現(xiàn)有一些應(yīng)用廣泛的瞬時利率模型：(ⅰ)如果模型中參數(shù)α=β=λ=0，則式(1)變?yōu)榇_定性利率變動模型。(ⅱ)如果模型中參數(shù)α=λ=0且β>0，則式(1)變?yōu)閂asicek[1]模型。(ⅲ)如果模型中參數(shù)β=λ=0且α>0，則式(1)變?yōu)镃IR[2]模型。(ⅳ)當(dāng)參數(shù)α=0且β>0，則式(1)變?yōu)镴Vasicek模型[19]；當(dāng)參數(shù)α>0且β=0，則式(1)變?yōu)镴CIR模型[17]。

按照風(fēng)險中性定價原理，該經(jīng)濟體中所有交易證券的折現(xiàn)過程是(Ω,Ft)-鞅過程。記該經(jīng)濟體中到期期限為T

(2)

F(T,r,T)=h(r)，

(3)

其中微分算子L定義為：

下面給出零息票債券的定價。記零息票債券t≤T時刻價格為P(t,r,T)。顯然,零息票債券在到期日T時的收益等于1，即P(T,r,T)=1，亦即H(rT)=1。于是由式(2)可得

(4)

為了計算債券及其衍生品的價格，先給出引理如下。

引理1 設(shè)短期利率rt滿足模型(1)。于是，對任意滿足exp{zr}有界的z∈C，復(fù)Fourier變換：

(5)

具有顯示表達式exp{A(t,T,z)+B(t,T,z)rt}，其中：

(6)

(7)

比較狀態(tài)變量前的系數(shù),進一步得到A(t,T,z)與B(t,T,z)分別滿足下列常微分方程：

(8)

(9)

現(xiàn)求解方程(9)。分兩種情形：α≠0和α≡0考慮。

(10)

利用Riccati方程理論，易得D具有表達式

(11)

將式(11)回代得B的表達式。

由于方程(9)中含有積分項，為了獲得系數(shù)A(t，T，z)的封閉式解，這里采用文獻[17]或文獻[19]中同樣的線性化技術(shù)近似方程(9)中的積分項，即應(yīng)用2階Taylor展開積分項中的指數(shù)函數(shù)：

于是

現(xiàn)令z≡0，并記τ=T-t，則零息票債券的價格為：

引理2T-期零息票債券t時刻的價格，P(t,r,t+τ),t

P(t,r,t+τ)=exp{a(τ)+b(τ)rt},

(12)

其中a(τ)=A(t,T,0)，b(τ)=B(t,T,0)。

2 債券期權(quán)定價

本節(jié)考慮標(biāo)的資產(chǎn)為零息票債券的歐式未定權(quán)益的定價。首先給出歐式零息票債券看漲期權(quán)的價格,再應(yīng)用歐式零息票債券看漲期權(quán)的價格研究付息債券期權(quán)及利率期權(quán)的定價。記歐式零息票債券看漲期權(quán)的到期日為T1(

(13)

并且有：

定理1 設(shè)短期利率rt滿足模型(1)，則

Ce(t,r,T1,T,K)=P(t,r,T)Π1-KP(t,r,T1)Π2,

(14)

φ1(u)=ea(T-T1)-a(τ)-b(τ)rΨ(t,T1,r,(iu+1)b(T-T1)),

φ2(u)=e-a(T1-t)-b(T1-t)rΨ(t,T1,r,iu,b(T-T1)),

R[·]表示[·]的實部，i是虛數(shù)單位。

證明 由式(13)可知

(15)

(16)

容易驗證 測度Q1和Q2都是概率測度Q的等價鞅概率測度。因此，在新的概率測度Q1和Q2下，式(15)可寫成：

Ce(t,r,T1,T,K)=P(t,r,T1)Q1(lnP(T1,r,T)≥lnK)-

KP(t,r,T)Q2(lnP(T1,r,T)≥lnK)=P(t,r,T1)Π1-KP(t,r,T)Π2,

(17)

其中Πj=Qj(lnP(T1,r,T)≥lnK)=Qj(b(T-T1)rT1≥K*)，j=1，2，且應(yīng)用Fourier反變換法，有

(18)

這里φj(u)，j=1，2是隨機變量b(T-T1)rT1分別在Q1和Q2下的特征函數(shù)。

下面計算在原概率測度Q下的特征函數(shù)φj(u)，j=1,2。由式(16)得：

ea(T1,T)-a(t,T)-b(t,T)rΨ(t,T1,r,(iu+1)b(T-T1))

(19)

以及

e-a(T1-t)-b(T1-t)rΨ(t,T1,r,iub(T-T1))。

(20)

類似容易證明，執(zhí)行價格為K，到期日為T1的歐式零息票債券看跌期權(quán)t時刻的價格為：

下面將歐式零息票期權(quán)定價的結(jié)果應(yīng)用于付息債券期權(quán)和利率期權(quán)的定價。

(21)

注2 式(21)表明付息債券看漲期權(quán)的價格是由零息票債券看漲期權(quán)價格的加權(quán)總和。

證明 因為

(22)

(23)

再由定理1得證。證畢。

(Ⅱ)利率期權(quán)(interest-rateoption)是一種與利率變化掛鉤的期權(quán)，到期時以現(xiàn)金或者與利率相關(guān)的合約(如利率期貨、遠期利率或者政府債券)進行結(jié)算。最早在場外市場交易的利率期權(quán)是1985年推出的利率上限期權(quán)(interest-rateCaps)。目前，利率期權(quán)除利率上限期權(quán)合約外，還有利率下限期權(quán)(interest-rateFloors)，利率雙限期權(quán)(interest-rateCollars)，利率互換期權(quán)(Swaps)等。這里僅給出利率下限期權(quán)的定價，其他產(chǎn)品的定價可類似證明。利率下限期權(quán)是由一系列的利率看跌期權(quán)(Floorlet)組合而成，而利率看跌期權(quán)的收益依賴于在時間t(

(24)

定理3 在短期利率rt滿足模型(1)下，則歐式利率看跌期權(quán)t(t

(25)

證明 由式(24)及條件期望的迭代性質(zhì)知：

(26)

(27)

注3 由定理3，很容易推導(dǎo)歐式利率上限期權(quán)(CapOption)和利率互換期權(quán)(Swaps)的定價公式分別為：

3 數(shù)值結(jié)果與分析

本節(jié)應(yīng)用數(shù)值計算實例考察仿射跳擴散利率模型(1)下債券及其標(biāo)準歐式看漲期權(quán)的價格影響因素。首先，分析債券價格分別在JVasicek與JCIR-Ex利率期限結(jié)構(gòu)模型下對跳躍風(fēng)險因素的性能表現(xiàn),并分別與純擴散過程的Vasicek和CIR利率期限結(jié)構(gòu)模型的結(jié)果進行比較。其次，分別在JVasicek與JCIR-Ex模型下研究歐式零息票債券看漲期權(quán)的價格受參數(shù)影響的變動情況。最后,討論債券期權(quán)的隱含波動率問題。數(shù)值計算選用Matlab 8.0軟件編程在Intel(R) Core(TM) i5-2500 CPU 3.30 GHz，4 GB RAM的聯(lián)想計算機上實現(xiàn)。模型(1)部分參數(shù)值來源于文獻[21](見表1)。

表1 模型(1)中基本參數(shù)值選擇

3.1 跳躍風(fēng)險參數(shù)對債券價格的影響

圖1 純擴散模型與跳擴散模型下零息票債券價格之比較Fig.1 Comparison of bond price from the Pure-diffusion and Jump-diffusion model

圖1考察了到期期限為10年的零息票債券價格分別在純擴散和跳擴散的利率期限結(jié)構(gòu)模型下的變動情況。圖1(a)選取參數(shù)α=0,β=0.002 5時比較了Vasicek利率模型與帶跳的Vasicek(即JVasicek)利率模型的債券價格。圖1(b)選取參數(shù)α=0.04,β=0或β≠0時比較了CIR利率模型與帶跳的CIR(即JCIR-Ex)利率模型的債券價格。計算結(jié)果表明：債券價格是期限值的減函數(shù)。另外，隨著期限的增加，JVasicek模型下的債券價格越來越高于純擴散過程驅(qū)動的Vasicek模型,這說明跳躍風(fēng)險對債券價格有正向影響，但JCIR或JCR-Ex模型下的債券價格卻越來越低于純擴散過程驅(qū)動的CIR模型,跳躍風(fēng)險對債券價格的作用是反向的。這是因為Vasicek利率模型服從正態(tài)分布，隨著期限的增加且跳躍風(fēng)險的存在,導(dǎo)致利率變小的機會增大,而CIR利率模型是服從卡方分布的,期限的增加與跳躍風(fēng)險的存在,促使利率值快速增加。

表2 債券價格對利率初始值及跳躍強度的依賴性

表2研究了利率初始值r0及跳躍強度λ對債券價格影響。顯然，在模型其他參數(shù)值不變的情形下，債券價格隨著利率初始值的增加而減少。此外,對于JVasicek利率模型債券價格隨著跳躍強度λ值減少而減少,且快速收斂到Vasicek模型時的債券價格,但對于JCIR-Ex利率模型債券價格則隨著λ的減少快速增大收斂到CIR模型的價格。這表明了跳躍強度參數(shù)λ對債券價格有顯著的影響,且對Vasicek與CIR利率期限結(jié)構(gòu)模型的債券收益有不同的影響效果。與此同時,在模型參數(shù)值相同情形下,JVasicek模型下的債券價格比JCIR-Ex模型下的債券價格要大得多，兩者相差近25倍，這是因為JVasicek模型下的利率值要小。

圖2 JVasicek模型與JCIR-Ex模型下參數(shù)μJ和σJ對零息票債券價格的影響Fig.2 Sensitivities of both JVasicek model and JCIR-Ex model for the bond price with μJ and σJ

圖2考查了債券價格依賴利率跳躍幅度參數(shù)μJ和σJ的敏感性。圖2(a)、(b)分別繪制了JVasicek模型與JCIR-Ex模型下債券價格隨參數(shù)μJ和σJ的變化情況。從圖示可以看出,跳擴散模型下的債券價格均是利率跳躍幅度參數(shù)σJ的減函數(shù)，但參數(shù)μJ在JVasicek模型與JCIR-Ex模型下對債券價格的影響是反向的。

3.2 期權(quán)價格的參數(shù)依賴性

表3分別考查了JVasicek模型和JCIR-Ex模型中的跳躍風(fēng)險參數(shù)λ、μJ和σJ對債券期權(quán)價格的影響,計算基本參數(shù)值見表1，這里選擇K=0.5,T1=1 a。從表中可以看出，歐式零息票債券看漲期權(quán)價格是利率跳躍振幅的期望μJ的減函數(shù)。對于JVasicek利率模型，債券期權(quán)的價格是利率跳躍振幅的方差σJ的增函數(shù)，但對于JCIR-Ex利率模型，結(jié)果卻相反。另外，一個有趣的發(fā)現(xiàn)：在JVasicek模型與JCIR-Ex模型中，利率跳躍振幅的期望μJ大小會影響跳躍強度參數(shù)λ對債券期權(quán)價格的作用效果。例如，在JVasicek模型中，當(dāng)μJ=0時，債券期權(quán)價格是跳躍強度參數(shù)λ的增函數(shù),但μJ≠0時，債券期權(quán)價格反而是跳躍強度參數(shù)λ的減函數(shù)。另一方面，在JCIR-Ex模型中，債券期權(quán)價格均是跳躍強度參數(shù)λ的減函數(shù),且不受μJ值的影響。

表3 跳躍風(fēng)險參數(shù)對零息債券看漲期權(quán)價格的影響

表4 JVasicek模型與JCIR-Ex模型下不同T1、r0和K值的零息票債券歐式看漲期權(quán)價格

表4考查了標(biāo)準歐式零息票債券看漲期權(quán)在JVasicek模型與JCIR-Ex模型下隨著T1、r0和K值變化時的價格。模型(1)計算的基本參數(shù)值見表1,選取利率初始值r0=0.04,0.06,0.08，期權(quán)到期日T1=1,3,5 a,以及期權(quán)執(zhí)行價K=0.5,0.7,0.9,1.0分別計算,結(jié)果見表4。從表中可以看出,零息票債券看漲期權(quán)價格隨著執(zhí)行價格K的增大(從價內(nèi)期權(quán)變化到價外期權(quán))而減少,這與股票型看漲期權(quán)的變化趨勢一致，并且期權(quán)價格是T1的增函數(shù)。另外,當(dāng)短期利率滿足JVasicek模型時,零息票債券看漲期權(quán)價格隨著r0的增加而變??；然而，當(dāng)短期利率滿足JCIR-Ex模型時,零息票債券看漲期權(quán)價格卻隨著r0的增加而變大。此外，JVasicek模型與JCIR-Ex模型下的標(biāo)準歐式零息票債券看漲期權(quán)價格相差很大,這表明JCIR-Ex模型比JVasicek模型在刻畫債券價格方面具有更小的波動或低債券價格(與表2、圖2的分析結(jié)果一致)。

3.3 債券期權(quán)的隱含波動率

金融市場通常以隱含波動率(implied volatility)的方式來引用看漲或看跌期權(quán)的價格變化。實際上，金融市場交易的歐式看漲期權(quán)的隱含波動率曲面依賴于期權(quán)敲定價格K和到期時間T。對相同的到期時間T，表現(xiàn)出波動率微笑/傾斜(smile/skew)特性，而且隨著時間演變，隱含波動率的曲面也在變動。隱含波動率問題已引起了業(yè)界和學(xué)術(shù)界對期權(quán)定價模型的推廣。下面分析仿射跳擴散利率模型下債券期權(quán)的隱含波動率問題，數(shù)值計算參數(shù)見表1。圖3(a)和(b)分別刻畫了JVasicek模型與JCIR-Ex模型下期權(quán)的隱含波動率變化。從隱含波動率的微笑曲面可以看出,仿射跳擴散的利率期限結(jié)構(gòu)模型具有較好擬合實際的能力，這與文獻[10，15]的實證結(jié)果一致。另外，從圖示還可以發(fā)現(xiàn)：①JVasicek模型下債券期權(quán)的隱含波動率要比JCIR-Ex模型下的小很多,且表現(xiàn)出較平坦；②期權(quán)到期期限短的隱含波動率的微笑現(xiàn)象比較顯著；③JVasicek模型下債券期權(quán)的隱含波動率是期權(quán)到期日T1的減函數(shù),而JCIR-Ex模型下債券期權(quán)的隱含波動率卻是期權(quán)到期日T1的增函數(shù)。

圖3 仿射跳擴散模型下債券期權(quán)的隱含波動率Fig.3 Implied volatilities of bond options in affine jump-diffusion model

圖3(c)和(d)在JCIR-Ex模型下分別考察跳躍強度參數(shù)λ=0及跳躍幅度的波動率σJ=0.02時的隱含波動率曲面。通過比較圖3(b)和(c)可以發(fā)現(xiàn)，跳躍強度參數(shù)對債券期權(quán)的隱含波動率影響很大，特別是在λ=0時(無跳躍風(fēng)險)債券期權(quán)的隱含波動率微笑現(xiàn)象顯著，這也表明引入跳躍風(fēng)險因素來捕捉市場利率的突發(fā)信息到達引發(fā)的利率劇烈波動是合理可行的。另外，圖3(b)和(d)表明JCIR-Ex模型下利率跳躍幅度的波動率對債券期權(quán)的隱含波動率的影響顯著，大波動率導(dǎo)致隱含波動率曲面更傾斜。

4 結(jié)論

本文在一類仿射跳擴散利率期限結(jié)構(gòu)模型下應(yīng)用Fourier變換方法和遠期測度技術(shù)獲得了債券及債券期權(quán)的價格顯示解,并應(yīng)用數(shù)值實例分析了模型中跳躍風(fēng)險參數(shù)對債券或期權(quán)價格的影響。數(shù)值計算結(jié)果表明，跳擴散過程驅(qū)動的利率期限結(jié)構(gòu)模型對利率衍生品的影響具有不同作用,也說明在連續(xù)的純擴散模型中引入不連續(xù)運動的跳躍風(fēng)險因素是必須的,且具有較好的擬合實際能力。本文僅在利率滿足單因素情形下研究利率衍生品的定價,模型還可以進一步擴展到多因素情形的跳擴散過程(如文獻[27])，以及研究其他利率衍生品的定價，這將是未來研究重點。

[1] VASICEK O. An equilibrium characterization of the term structure[J]. Journal of Financial Economics,1977，5(2): 177-188. DOI:10.1016/0304-405X(77)90016-2.

[2] COX J C，INGERSOLL J E，ROSS S A. A theory of the term structure of interest rates[J]. Econometrica，1985，53(2): 385-407. DOI:10.2307/1911242.

[3] BRENNAN M J，SCHWARTZ E S. A continuous time approach to the pricing of bonds[J]. Journal of Banking and Finance，1979，3(2): 133-155. DOI:10.1016/0378-4266(79)90011-6.

[4] HULL J，WHITE A. Pricing interest-rate-derivative securities[J]. The Review of Financial Studies，1990，3(4): 573-592. DOI:10.1093/rfs/3.4.573.

[5] HEATH D，JARROW R，MORTON A. Bond pricing and the term structure of interest rates: a new methodology for contingent claims valuation[J]. Econometrica，1992，60(1): 77-105. DOI:10.2307/2951677.

[6] CHEN R R，SCOTT L. Pricing interest rate options in a two-factor Cox-Ingersoll-Ross model of the term structure[J]. The Review of Financial Studies，1992，5(4): 613-636. DOI:10.1093/rfs/5.4.613.

[7] BANSAL R，ZHOU Hao. Term structure of interest rates with regime shifts[J]. The Journal of Finance，2002，57(5): 1997-2043. DOI:10.1111/0022-1082.00487.

[8] DUFFIE D，KAN R. A yield-factor model of interest rates[J]. Mathematical Finance，1996，6(4): 379-406. DOI:10.1111/j.1467-9965.1996.tb00123.x.

[9] HAMILTON J D. Rational-expectations econometric analysis of changes in regime: an investigation of the term structure of interest rates[J]. Journal of Economic Dynamics and Control，1988，12(2/3): 385-423. DOI:10.1016/0165-1889(88)90047-4.

[10] DAS S R. The surprise element: jumps in interest rates[J]. Journal of Econometrics，2002，106(1): 27-65. DOI:10.1016/S0304-4076(01)00085-9.

[11] JOHANNES M. The statistical and economic role of jumps in continuous-time interest rate models[J]. The Journal of Finance，2004，59(1): 227-260. DOI:10.1111/j.1540-6321.2004.00632.x.

[12] PIAZZESI M. Bond yields and the federal reserve[J]. Journal of Political Economy，2005，113(2): 311-344. DOI:10.1086/427466.

[13] ANDERSEN T，BENZONI L，LUND J. Stochastic volatility，mean drift，and jumps in the short-term interest rate[R]. Chicago: Northwestern University, 2004.

[14] JARROW R，LI Haitao，ZHAO Feng. Interest rate caps "smile" too! But can the LIBOR market models capture the smile?[J]. The Journal of Finance，2007，62(1): 345-382. DOI:10.1111/j.1540-6261.2007.01209.x.

[15] LIN B H，YEH S K. Jump-diffusion interest rate process: an empirical examination[J]. Journal of Business Finance and Accounting，1999，26(7/8): 967-995. DOI:10.1111/1468-5957.00282.

[16] 陳松男. 利率衍生品設(shè)計原理與應(yīng)用: 案例分析[M]. 北京: 機械工業(yè)出版社，2014.

[17] AHN C M，THOMPSON H E. Jump-diffusion processes and the term structure of interest rates[J]. The Journal of Finance，1988，43(1): 155-174. DOI:10.1111/j.1540-6261.1988.tb02595.x.

[18] SHIRAKAWA H. Interest rate option pricing with Poisson-Gaussian forward rate curve process[J]. Mathematical Finance，1991，1(4): 77-94. DOI:10.1111/j.1467-9965.1991.tb00020.x.

[19] BAZ J，DAS S R. Analytical approximations of the term structure for jump-diffusion process: a numerical analysis[J]. The Journal of Fixed Income，1996，6(1): 78-86. DOI:10.3905/jfi.1996.408164.

[20] DAS S R，F(xiàn)ORESI S. Exact solutions for bond and option prices with systematic jump risk[J]. Review of Derivatives Research，1996，1(1): 7-24. DOI:10.1007/BF01536393.

[21] DAS S R. Discrete-time bond and option pricing for jump-diffusion processes[J]. Review of Derivatives Research，1996，1(3): 211-243. DOI:10.1007/BF01531143.

[22] DAS S R. A direct discrete-time approach to Poisson-Gaussian bond option pricing in the Heath-Jarrow-Morton model[J]. J Econ Dyn Control，1998，23(3): 333-369. DOI:10.1016/S0165-1889(98)00031-1.

[23] FINNERTY J D. Exact formulas for pricing bonds and options when interest rate diffusions contain jumps[J]. The Journal of Financial Research，2005，28(3): 319-341. DOI:10.1111/j.1475-6803.2005.00127.x.

[24] BELIAEVA N A，NAWALKHA S K，SOTO G M. Pricing American interest rate options under the jump-extended Vasicek model[J]. Journal of Derivatives，2008，16(1): 29-43. DOI:10.3905/jod.2008.710896.

[25] BELIAEVA N，NAWALKHA S K. Pricing American interest rate options under the jump-extended constant-elasticity-of-variance short rate models[J]. Journal of Banking & Finance，2012，36(1): 151-163. DOI:10.1016/j.jbankfin.2011.06.012.

[26] DENG Guohe. Pricing American put option on zero-coupon bond in a jump-extended CIR model[J]. Commun Nonlinear Sci Numer Simulat，2015，22(1/2/3): 186-196. DOI:10.1016/j.cnsns.2014.10.003.

[27] DUFFIE D，PAN Jun，SINGLETON K. Transform analysis and asset pricing for affine jump-diffusions[J]. Econometrica，2000，68(6): 1343-1376. DOI:10.1111/1468-0262.00164.

[28] GEMAN H， KAROUI N E，ROCHET J C. Changes of numeraire，changes of probability measures and option pricing[J]. Journal of Applied Probability，1995，32(2): 443-458. DOI:10.2307/3215299.

(責(zé)任編輯 黃 勇)

Pricing of Interest Rate Derivatives Based on Affine Jump Diffusion Model

WANG Jiaqin, DENG Guohe

(School of Mathematics and Statistics，Guangxi Normal University, Guilin Guangxi 541004, China)

The pricing of interest rate derivatives is considered under an affine jump diffusion model. Using the Fourier transform method and the forward measure change technique，the closed explicit formulas for both the price of the default-free，zero-coupon bond and the value of the European options on the default-free，zero-coupon bond are obtained. Furthermore，pricing problems on both the European option on the coupon bond and the interest rate options are extended in this model by applying these explicit formulas above. Finally，the impacts of the key parameters in this model on prices for both the bond and bond option，and implied volatilities of bond options are analyzed by numerical examples，respectively. Numerical results show that the jump risks have more remarkable effects on the interest rate derivative prices and implied volatility，which show that the affine jump diffusion term structure model of the interest rate fits reality well.

affine jump diffusion model；term structure of interest rate；bond options；Fourier transform；implied volatility.

10.16088/j.issn.1001-6600.2016.03.011

2016-03-01

國家自然科學(xué)基金資助項目(11461008)；教育部人文社會科學(xué)研究規(guī)劃基金資助項目(13YJA910003)；廣西自然科學(xué)基金資助項目(2013GXNSFAA019005)；廣西高等學(xué)?？茖W(xué)技術(shù)研究重點項目(2013ZD010)

鄧國和(1969—)，男，湖南桂陽人，廣西師范大學(xué)教授，博士. E-mail:dengguohe@mailbox.gxnu.edu.cn

O211.9

A

1001-6600(2016)03-0074-12