文章快速检索    
  同济大学学报(自然科学版)  2018, Vol. 46 Issue (11): 1609-1614.  DOI: 10.11908/j.issn.0253-374x.2018.11.021
0

引用本文  

梁进, 邹宏春. 担保信用等级变换的利率互换衍生品定价[J]. 同济大学学报(自然科学版), 2018, 46(11): 1609-1614. DOI: 10.11908/j.issn.0253-374x.2018.11.021.
LIANG Jin, ZOU Hongchun. Pricing of Interest Rate Swap Derivatives for Assuring Credit Rating Migration[J]. Journal of Tongji University (Natural Science), 2018, 46(11): 1609-1614. DOI: 10.11908/j.issn.0253-374x.2018.11.021

基金项目

国家自然科学基金(11671301)

第一作者

梁进(1958—), 女, 教授, 博士生导师, 理学博士, 主要研究方向为金融数学及信用风险管理. E-mail:liang_jin@tongji.edu.cn

文章历史

收稿日期:2018-01-26
担保信用等级变换的利率互换衍生品定价
梁进 , 邹宏春     
同济大学 数学科学学院,上海 200092
摘要:考虑担保信用等级迁移风险的利率互换合约, 在结构化方法的框架下,建立了担保信用等级首次迁移所造成损失的保费定价模型,信用等级依赖于利率并具有高低2个等级.从一个新的角度将低等级下零息票利率互换的价值定义为保费价值的自变量,并利用对冲原理建立保费定价的偏微分方程模型.利用计价单位转换原理对方程进行降维,求出问题的半解析解,然后采用有限差分方法的显式格式对模型进行数值求解.最后, 讨论了保费价值关于利率参数的依赖关系.结果表明:保费价值和各参数间存在单调递减的关系.
关键词利率互换    信用等级变换    降维    有限差分方法    
Pricing of Interest Rate Swap Derivatives for Assuring Credit Rating Migration
LIANG Jin , ZOU Hongchun     
School of Mathematical Sciences, Tongji University, Shanghai 200092, China
Abstract: Considering the valuation of a protected swap on credit rating migration, under the framework of structural methods, a pricing model was established for protecting the loss caused by the first credit rating migration, where the credit ratings depend on the interest rate and have two grades. In the pricing model, an independent variable of the model was defined by the value of low-grade zero-coupon, which is a new perspective. A partial differential equation(PDE) pricing model was derived by the hedging method. A semi-closed solution was obtained by dimensionality reduction technique. The numerical solution was calculated by the explicit finite difference method. Finally, the dependency of parameters of the model was discussed and the results show that there is a monotonically decreasing relationship between premium value and each parameter.
Key words: interest rate swap    credit rating migration    dimensionality reduction    finite difference method    

自20世纪80年代初发生了第一笔互换合约后,互换市场发展迅猛,金融互换交易已成为目前国际金融市场最大的融资工具之一,其中利率互换和货币互换是目前市场上最核心的2种互换工具.虽然利率互换可以降低资金风险和利率风险,但是随着次贷危机和欧债危机的爆发,作为反映受评对象违约可能性大小的信用等级,越来越受到人们的关注.

信用等级确定的目的在于评估受评对象违约可能性的大小,一般由专门的信用评级机构进行评估.目前国际三大评级机构对信用评级的定义基本一致,均认为信用评级是对债务人偿债能力和偿债意愿的综合评价.三大评级机构认为:信用评级只是对受评对象信用风险的评价,不是对其资产价值的度量,不能单独用作投资操作的依据[1].

目前市场上对于信用风险的管理主要集中在违约风险上,对于信用等级变换风险管理的研究还比较少.在利率互换过程中,虽然利率作为一种公共因子存在,但是在利率互换这一特定的参考实体中,利率的高低会影响到互换一方违约的可能性(如利率下跌会使得固定利率支付方违约可能性变大),和公司的信用挂钩,因此定义利率信用等级.利率越低,固定利率支付方公司违约可能性越大,利率信用等级越低.

国内外学者对信用等级变换模型做了许多研究. Jarrow等[2]首次用Markov链模型描述信用等级迁移过程,并给出了在迁移强度为常数以及回收率为零的情况下公司债券价格的解析解.梁进等[3]研究了约化方法下信用等级迁移和违约风险的零息票债券定价,模型中信用等级迁移强度与利率公共因子有关,在假设参数为常数的情况下得出了每个信用等级债券价值的解析解.以上都是基于约化模型框架下对信用等级变换过程进行刻画,但信用等级变换与公司资产价值密切相关,所以采用结构化模型刻画信用等级变换比约化方法更加自然直接.目前为止,通过结构化模型刻画信用等级变换的研究比较少.曾楚琨[4]利用结构化方法讨论了具信用等级迁移的债券定价,分别研究了利率为常数和随机利率情况下的定价. Hu等[5]研究了具信用等级迁移债券定价的自由边界问题,通过结构化方法对债券定价并得到了自由边界问题的一些性质. Liang等[6]利用结构化方法研究了具信用等级迁移风险公司债券的效用无差异定价.

针对目前市场上存在的各种信用问题,投资者们自然想要采用某种手段来规避自己所面临的信用风险,于是信用违约互换(CDS)应运而生.它是一份为特定公司的违约风险提供保护的合约,特定公司称为参考公司.当投资者买入参考公司发行的参考债券时,为减少或避免参考公司违约对投资者造成的损失,投资者(即CDS买方)会向交易对手(即CDS卖方)买入一份CDS合约.合约规定,当参考公司违约时,交易对手须对投资者的损失进行赔偿,为此投资者须定期向交易对手支付保费,直至合约到期或信用事件的发生.作为目前市场上最流行的信用衍生产品,CDS研究受到广大学者的关注,研究主要集中在担保违约风险的情况[7-9].

对于利率互换中的违约风险问题,梁进等[10]利用约化方法讨论了信用攸关的利率互换的定价,并且对定价函数的性质和参数的依赖关系进行分析. Huang等[11]研究了无违约条款和违约条款的信用攸关的利率互换的定价,利用历史数据对参数进行估计并得出了合约价值的数值解.

在结构化方法框架下,本文给出在利率互换过程中担保信用等级的首次迁移所带来损失的合约定价模型.采用结构化方法刻画利率信用等级的首次迁移,并且在对合约价值的自变量的处理方法中,将低等级下零息票利率互换的价值作为合约价值函数的自变量.利用对冲技巧推导合约价值满足的偏微分方程,并且利用计价单位转换的方法得到合约价值的半解析解.最后, 通过对偏微分方程进行数值模拟得到合约价值随时间变化的关系图,并且分析了各个参数对合约价值的影响.

1 模型的建立

在利率互换中,浮动利率支付方支付浮动利率,固定利率支付方支付固定利率.随机利率时时刻刻都在变化,随机利率下跌使得固定利率支付方违约的可能性变大(相对应的随机利率上升导致浮动利率支付方的违约可能性变大),一旦随机利率低于某个特定值(信用等级边界),随机利率所遵循的随机过程就会发生变化.因此,对于浮动利率支付方而言,采用零息票定价法计算得到的利率互换的价值就会不一样,那么就有可能导致浮动利率支付方产生损失.

考虑这样一则合约,其目的是将利率互换过程中信用等级变换的风险转移给愿意承担这份风险的金融机构,在合约期限内,公司A(浮动利率支付方)与公司B(固定利率支付方)进行利率互换.公司A面临随机利率发生信用等级变换的风险,即若利率下跌,则公司B违约的可能性就会变大,公司B的信用等级下降.于是,公司A与另外一家金融机构C签订一张合约来规避利率信用等级变换带来的损失,其中公司A为利率保护买方,金融机构C为利率保护卖方.公司A一次性向金融机构C支付保费,在合约期间一旦随机利率发生信用等级变换,则合约终止,并且公司A利率互换的损失由金融机构C来支付.

本节主要从Vasicek利率模型出发,假设信用等级只有高低2种等级,并且只担保首次发生信用等级变换所带来的损失.利用结构化方法来定义利率信用等级迁移时刻,将低等级下利率互换的价值作为合约价值的自变量,以此来构建利率互换过程中利率等级发生变换的合约定价模型.

1.1 基本假设

作如下基本假设:

(1) 假设市场是完备的,不存在套利机会.

(2) 利率产品由一张面值为1、到期日为T的零息票债券构成.

(3) 担保的是固定利率支付方的信用等级,保费由利率信用等级保护买方(浮动利率支付方)在期初一次性支付给卖方(相应的担保浮动利率支付方也可以做类似考虑).

(4) 存在一个信用等级边界ra.当随机利率rra,利率处于高等级;当rra,利率处于低等级.初始时刻随机利率r0ra,即初始时刻利率处于高等级.

(5) 市场利率模型为Vasicek模型, 如下所示:

$ {\rm{d}}{r_t} = a\left( {\theta - {r_t}} \right){\rm{d}}t + \left( {{\sigma _1}{1_{\left\{ {{r_t} \ge {r_{\rm{a}}}} \right\}}} + {\sigma _2}{1_{\left\{ {{r_t} < {r_{\rm{a}}}} \right\}}}} \right){\rm{d}}{W_t} $

式中:aθ均为正常数;1A为事件A的示性函数;σ1σ2为常数,分别表示高、低等级下波动率, 有σ1σ2;{Wt: 0≤tT}为由概率空间(Ω, $\mathcal{F}$, Q)生成的标准Brown运动.

(6) 合约担保的是固定利率支付方信用等级的下降,并且只考虑担保信用等级的首次迁移.随机利率首次发生信用等级迁移的时刻定义为

$ \tau = \inf \left\{ {t\left| {{r_0} > {r_{\rm{a}}},{r_t} \le {r_{\rm{a}}}} \right.} \right\} $
1.2 现金流分析 1.2.1 零息票定价公式[12]

由于利率的不可交换性,因此利用利率的载体——零息票的价值来表示利率互换的价值.零息票的价值$ {P_t} = E\left( {{{\rm{e}}^{ - \int_t^T {{r_s}{\rm{d}}s} }}|r\left( t \right) = {r_t}} \right) $.

若不考虑信用等级变换,市场利率模型为以下Vasicek模型:

$ {\rm{d}}{r_t} = a\left( {\theta - {r_t}} \right){\rm{d}}t + \sigma {\rm{d}}{W_t} $

由Feynman-Kac公式可知,在鞅测度下函数P(r, t; T)适合以下偏微分方程Cauchy问题:

$ \left\{ \begin{array}{l} \frac{{\partial P}}{{\partial t}} + {{\cal L}_0}P = 0,\;\;\;r \in {\bf{R}},t \in \left[ {0,T} \right)\\ P\left( {r,T} \right) = 1,\;\;\;r \in {\bf{R}} \end{array} \right. $ (1)

其中,

$ {{\cal L}_0} = \frac{{{\sigma ^2}}}{2}\frac{{{\partial ^2}}}{{\partial {r^2}}} + a\left( {\theta - r} \right)\frac{\partial }{{\partial r}} - r $

上述Cauchy问题有如下形式的仿射结构解:

$ P\left( {r,t;T} \right) = A\left( t \right){{\rm{e}}^{ - rB\left( t \right)}} $ (2)

其中,

$ \begin{array}{l} A\left( t \right) = \exp \left( {\frac{1}{{{a^2}}}\left( {B{{\left( t \right)}^2} - \left( {T - t} \right)} \right)\left( {{a^2}\theta - \frac{{{\sigma ^2}}}{2}} \right) - } \right.\\ \;\;\;\;\;\;\;\;\;\;\;\left. {\frac{{{\sigma ^2}}}{{4a}}B{{\left( t \right)}^2}} \right) \end{array} $
$ B\left( t \right) = \frac{1}{a}\left( {1 - {{\rm{e}}^{ - a\left( {T - t} \right)}}} \right) $
1.2.2 现金流损失分析

由于在信用等级变换之前并没有产生损失,损失是由发生信用等级变换引起的,因此只需要分析信用等级变换以后直至到期日的时间段内2个不同信用等级下零息票的价值差.

高、低等级下利率分别用r1tr2t表示,相应地零息票利率互换的价值分别为P1tP2t.对于高等级而言,将低等级下零息票P2t=P(r2t, t; T)看作一种风险资产,由Itô公式可知[12]

$ \begin{array}{l} {\rm{d}}{P_{2t}} = \left( {\frac{{\partial {P_2}}}{{\partial t}} + \frac{{\sigma _2^2}}{2}\frac{{{\partial ^2}{P_2}}}{{\partial r_2^2}} + a\left( {\theta - {r_{2t}}} \right)\frac{{\partial {P_2}}}{{\partial {r_2}}}} \right){\rm{d}}t + \\ \;\;\;\;\;\;\;\;\;{\sigma _2}\frac{{\partial {P_2}}}{{\partial {r_2}}}{\rm{d}}{W_t} \end{array} $

由式(2)可知

$ \frac{{\partial {P_2}}}{{\partial {r_2}}} = - B\left( t \right){P_2} $

同时, 由于P2t满足式(1),立刻可得P2t满足的随机微分方程,如下所示:

$ \frac{{{\rm{d}}{P_{2t}}}}{{{P_{2t}}}} = {r_{2t}}{\rm{d}}t - B\left( t \right){\rm{d}}{W_t} $ (3)

若固定利率支付方因随机利率下降而在利率互换合约到期日T之前发生信用等级变换,变换时刻为τ,则利率保护买方(浮动利率支付方)将由于固定利率支付方信用等级变换而产生利差.若P1τP2τ,则利差为P1τ-P2τ;若P1τP2τ,则利差为零.因此,此时利率保护卖方应该向保护买方赔付的金额为(P1τ-P2τ)+.

由于赔付金额的表达式含低等级下零息票的价值P2,因此本文从一个新的角度对合约价值的自变量进行定义.信用等级一旦发生变化(从高等级到低等级),那么此时低等级下零息票对于高等级下利率可以看作是一种风险资产.式(3)为低等级零息票价值满足的随机微分方程.类似于债券和期权定价模型[13],将风险资产低等级下零息票价值P2作为合约价值的自变量,把合约的价值看作是时间t、高等级下利率r1和风险资产低等级下零息票价值P2的函数.因此,在t时刻合约的价值

$ V\left( {t,{r_1},{P_2}} \right) = E\left( {{{\rm{e}}^{ - \int_t^\tau {{r_{1s}}{\rm{d}}s} }}{{\left( {{P_{1\tau }} - {P_{2\tau }}} \right)}^ + }{1_{\left\{ {\tau < T} \right\}}}\left| {{{\cal F}_t}} \right.} \right) $
1.3 偏微分方程的推导

下面利用对冲技巧[13]来推导合约的价值Vt=V(t, r1, P2)所满足的偏微分方程.

构造投资组合

$ {\mathit{\Pi }_t} = {V_t} - {\Delta _{1t}}{P_{1t}} - {\Delta _{2t}}{P_{2t}} $

式中:Δ1t和Δ2t分别表示无等级变换风险高等级下零息票债券和低等级下风险资产的份额.通过选取适当的Δ1t和Δ2t,使得投资组合在高等级随机利率r1下是无风险的,如下所示:

$ {\rm{d}}{\mathit{\Pi }_t} = {\rm{d}}{V_t} - {\Delta _{1t}}{\rm{d}}{P_{1t}} - {\Delta _{2t}}{\rm{d}}{P_{2t}} = {r_{1t}}{\mathit{\Pi }_t}{\rm{d}}t $ (4)

由Itô公式,并结合式(1)和式(4),Vt=V(t, r1, P2)满足如下偏微分方程:

$ \frac{{\partial V}}{{\partial t}} + {{\cal L}_1}V + a\left( {\theta - {r_1}} \right)\frac{{\partial V}}{{\partial {r_1}}} + {r_1}{P_2}\frac{{\partial V}}{{\partial {P_2}}} - {r_1}V = 0 $ (5)

其中,

$ {{\cal L}_1} = \frac{{\sigma _1^2}}{2}\frac{{{\partial ^2}}}{{\partial r_1^2}} + \frac{{\sigma _2^2}}{2}B{\left( t \right)^2}P_2^2\frac{{{\partial ^2}}}{{\partial P_2^2}} - {\sigma _1}{\sigma _2}B\left( t \right){P_2}\frac{{{\partial ^2}}}{{\partial {r_1}\partial {P_2}}} $

因为r1的求解范围为r1ra,所以先将区域变换到坐标轴的右半部分,即作变换$\tilde r$1=r1-ra,则有

$ {\rm{d}}{{\tilde r}_1} = {\rm{d}}{r_1} = a\left( {\theta - {r_{\rm{a}}} - {{\tilde r}_{{\rm{1}}t}}} \right){\rm{d}}t + {\sigma _1}{\rm{d}}{W_t} $ (6)

将式(6)代入式(5),则在区域

$ {\mathit{\Sigma }_1}:\left\{ {{{\tilde r}_1} > 0,{P_2} > 0,0 \le t < T} \right\} $

$ \left\{ \begin{array}{l} \frac{{\partial V}}{{\partial t}} + {{\tilde {\cal L}}_1}V + a\left( {\theta - {r_{\rm{a}}} - {{\tilde r}_1}} \right)\frac{{\partial V}}{{\partial {{\tilde r}_1}}} + \\ \;\;\;\;\left( {{{\tilde r}_1} - {r_{\rm{a}}}} \right){P_2}\frac{{\partial V}}{{\partial {P_2}}} - \left( {{{\tilde r}_1} + {r_{\rm{a}}}} \right)V = 0\\ V\left( {t,0,{P_2}} \right) = f\left( t \right),\;\;\;\;{P_2} > 0,0 \le t < T\\ V\left( {T,{{\tilde r}_1},{P_2}} \right) = 0,\;\;\;\;\;{{\tilde r}_1} > 0,{P_2} > 0 \end{array} \right. $ (7)
$ {{\tilde {\cal L}}_1} = \frac{{\sigma _1^2}}{2}\frac{{{\partial ^2}}}{{\partial \tilde r_1^2}} + \frac{{\sigma _2^2}}{2}B{\left( t \right)^2}P_2^2\frac{{{\partial ^2}}}{{\partial P_2^2}} - {\sigma _1}{\sigma _2}B\left( t \right){P_2}\frac{{{\partial ^2}}}{{\partial {{\tilde r}_1}\partial {P_2}}} $
$ f\left( t \right) = {\left( {{P_{1t,{r_{\rm{a}}}}} - {P_{2t,{r_{\rm{a}}}}}} \right)^ + } $

式中:f(t)表示在信用等级变换边界r1=ra$\tilde r$1=0处,由高信用等级变为低信用等级利率互换的价值损失; P1t, raP2t, ra分别表示高、低信用等级在r=ra时的利率互换价值,分别有式(2)形式,即

$ {P_{1t,{r_{\rm{a}}}}} = {A_1}\left( t \right){{\rm{e}}^{ - {r_{\rm{a}}}B\left( t \right)}},\;\;\;\;{P_{2t,{r_{\rm{a}}}}} = {A_2}\left( t \right){{\rm{e}}^{ - {r_{\rm{a}}}B\left( t \right)}} $

式中:Ai(t)(i=1, 2)分别为A(t)中σσi(i=1, 2)的情形.

2 模型求解与数值分析 2.1 偏微分方程求解

接下来对偏微分方程(7)进行求解.偏微分方程(7)是一个二维偏微分方程.首先利用计价单位转化方法[12]将方程转化为一维问题,然后通过各种变量变换和函数变换将方程转化为对抛物型偏微分方程的第一类初边值问题的求解.

通过计价单位转化的方法将上述二维偏微分方程进行降维,作变换$ y = \frac{{{P_2}}}{{{{\tilde P}_1}}}, \tilde V = \frac{V}{{{{\tilde P}_1}}} $, 其中$ {{\tilde P}_1} $为以$ {{\tilde r}_1} $为利率的零息债券的价值函数.由变换可得到在区域Σ2: {yy0, 0≤tT}上的方程,如下所示:

$ \left\{ \begin{array}{l} \frac{{\partial \tilde V}}{{\partial t}} + \frac{1}{2}{\left( {{\sigma _2} - {\sigma _1}} \right)^2}B{\left( t \right)^2}{y^2}\frac{{{\partial ^2}\tilde V}}{{\partial {y^2}}} + \\ \;\;\;\;{r_{\rm{a}}}y\frac{{\partial \tilde V}}{{\partial y}} - {r_{\rm{a}}}\tilde V = 0\\ \tilde V\left( {t,{y_0}} \right) = g\left( t \right),\;\;\;\;\;0 \le t < T\\ \tilde V\left( {T,y} \right) = 0,\;\;\;\;\;y > {y_0} \end{array} \right. $ (8)

其中,

$ {y_0} = \frac{{{P_{2t,{r_{\rm{a}}}}}}}{{{{\tilde P}_{1t,0}}}} = \frac{{{P_{2t,{r_{\rm{a}}}}}}}{{{P_{1t,{r_{\rm{a}}}}}}} $
$ g\left( t \right) = \frac{{f\left( t \right)}}{{{{\tilde P}_{1t,0}}}} = \frac{{f\left( t \right)}}{{{P_{1t,0}}}} = {\left( {1 - \frac{{{P_{2t,{r_{\rm{a}}}}}}}{{{P_{1t,{r_{\rm{a}}}}}}}} \right)^ + } = {\left( {1 - {y_0}} \right)^ + } $

为了将上述方程(8)转化为一般的热传导方程的初边值问题,作如下变换[13]:

$ s = T - \int_0^t {\left( {{\sigma _2} - {\sigma _1}} \right)B\left( \omega \right){\rm{d}}\omega } ,\;\;\;x = \ln \left( {y{{\rm{e}}^{{r_{\rm{a}}}\left( {T - t} \right)}}} \right), $
$ u = \tilde V{{\rm{e}}^{\frac{1}{8}s - \frac{1}{2}x}}{{\rm{e}}^{{r_{\rm{a}}}\left( {T - t} \right)}} $

在上述变换下, 方程变为如下形式:

$ \left\{ \begin{array}{l} \frac{{\partial u}}{{\partial s}} - \frac{1}{2}\frac{{{\partial ^2}u}}{{\partial {x^2}}} = 0\\ u\left( {s,{x_0}} \right) = h\left( s \right),\;\;\;\;0 \le s < T\\ u\left( {0,x} \right) = 0,\;\;\;\;x \in {\bf{R}} \end{array} \right. $ (9)

其中,

$ {x_0} = \ln \left( {{y_0}{{\rm{e}}^{{r_{\rm{a}}}\left( {T - t} \right)}}} \right),\;\;\;\;h\left( s \right) = g\left( {T - s} \right){{\rm{e}}^{\frac{1}{8}s + \frac{1}{2}x}} $

方程(9)为抛物型方程第一初边值问题,作如下变换[14]:

$ z = x - {x_0},\;\;\;\;\;\varphi \left( {s,z} \right) = u\left( {s,z} \right) - h\left( s \right) $

代入方程(9)可以求解得到如下形式的解:

$ \begin{array}{*{20}{c}} {\varphi \left( {s,z} \right) = - \int_{ - \infty }^{ + \infty } {\mathit{\Gamma }\left( {z - \xi ,s} \right)h\left( 0 \right){\rm{d}}\xi } - }\\ {\int_0^s {h'\left( \omega \right){\rm{d}}\omega } \int_{ - \infty }^{ + \infty } {\mathit{\Gamma }\left( {z - \xi ,s - \omega } \right){\rm{d}}\xi } } \end{array} $

其中,

$ \mathit{\Gamma }\left( {z - \xi ,s - \omega } \right) = \left\{ \begin{array}{l} \frac{1}{{\sqrt {2{\rm{ \mathsf{ π} }}\left( {s - \omega } \right)} }}{{\rm{e}}^{ - \frac{{{{\left( {z - \xi } \right)}^2}}}{{2\left( {s - \omega } \right)}}}},\;\;\;\;\;s > \omega \\ 0,\;\;\;\;s \le \omega \end{array} \right. $

结合上述变换,有合约的价值

$ V = {{\tilde P}_1}{{\rm{e}}^{ - {r_{\rm{a}}}\left( {T - t} \right)}}{{\rm{e}}^{\frac{1}{8}s - \frac{1}{2}x}}\left( {\varphi \left( {s,x} \right) + h\left( s \right)} \right) $ (10)
2.2 数值结果

通过建立方程的有限差分的显式格式求出合约价值的数值解.首先, 为了对方程(7)进行数值计算,在方程(7)的基础上对t作变换η=T-t; 然后, 对方程建立有限差分的显式格式,直接计算V(η, ${\tilde r}$1, P2).

为了简单起见,采用等步长剖分,并且在${\tilde r}$1P2方向的步长相同,记Δ${\tilde r}$1P2=h时间步长为Δη=τ1.原问题的求解区域变为

$ \left( {\eta ,{{\tilde r}_1},{P_2}} \right) = \left[ {0,T} \right] \times \left[ {0,{R_1}} \right] \times \left[ {0,{R_2}} \right] $

式中:R1R2为足够大的常数.在网格点(n, i, j)上记V的值为Vn, i, j=V(1, ih, jh),其中$ 0 \le n \le \frac{T}{h}, 0 \le i \le \frac{{{R_1}}}{h}, 0 \le j \le \frac{{{R_2}}}{h} $, 可以得到定解问题(7)的九点显式差分格式,如下所示:

$ \begin{array}{l} {V_{n + 1,i,j}} = {\pi _1}{V_{n,i,j}} + {\pi _2}{V_{n,i + 1,j}} + {\pi _3}{V_{n,i - 1,j}} + \\ \;\;\;\;\;\;\;\;\;\;\;\;{\pi _4}{V_{n,i,j + 1}} + {\pi _5}{V_{n,i,j - 1}} + {\pi _6}\left( {{V_{n,i + 1,j + 1}} - } \right.\\ \;\;\;\;\;\;\;\;\;\;\;\;\left. {{V_{n,i - 1,j + 1}} - {V_{n,i + 1,j - 1}} + {V_{n,i - 1,j - 1}}} \right) \end{array} $

其中,

$ \left\{ \begin{array}{l} {\pi _1} = 1 - \left( {{r_{\rm{a}}} + ih} \right){\tau _1} - \frac{{\sigma _1^2{\tau _1}}}{{{h^2}}}\\ {\pi _2} = \frac{{\sigma _1^2{\tau _1}}}{{2{h^2}}} + a\left( {\theta - {r_{\rm{a}}} - ih} \right)\frac{{{\tau _1}}}{{2h}}\\ {\pi _3} = \frac{{\sigma _1^2{\tau _1}}}{{2{h^2}}} - a\left( {\theta - {r_{\rm{a}}} - ih} \right)\frac{{{\tau _1}}}{{2h}}\\ {\pi _4} = \frac{{\sigma _2^2B{{\left( {n{\tau _1}} \right)}^2}{j^2}{\tau _1}}}{2} + \left( {{r_{\rm{a}}} + ih} \right)j\frac{{{\tau _1}}}{2}\\ {\pi _5} = \frac{{\sigma _2^2B{{\left( {n{\tau _1}} \right)}^2}{j^2}{\tau _1}}}{2} - \left( {{r_{\rm{a}}} + ih} \right)j\frac{{{\tau _1}}}{2}\\ {\pi _6} = \frac{{{\sigma _1}{\sigma _2}B\left( {n{\tau _1}} \right)j{\tau _1}}}{{4h}} \end{array} \right. $

相应的定解条件和人工边界条件为

$ \begin{array}{*{20}{c}} {{V_{n,0,j}} = f\left( {n{\tau _1}} \right),\;\;\;\;{V_{0, - \frac{{{R_1}}}{h},j}} = 0,\;\;\;\;{V_{0,\frac{{{R_1}}}{h},j}} = 0,}\\ {{V_{0,i,j}} = 0,\;\;\;\;{V_{n,i, - \frac{{{R_2}}}{h}}} = {P_{1n{\tau _1},{r_{\rm{a}}}}},\;\;\;{V_{n,i,\frac{{{R_2}}}{h}}} = 0} \end{array} $

利用上面的九点显式差分格式,继而代回原变量,即得到利率保护买方的保费的预期支付额.

图 1显示了保费的预期支付额V与随机利率r和时间t的关系曲面.取参数a=0.5, θ=0.04, σ1=0.04, σ2=0.06, ra=0.04, T=5,时间步长τ1=0.05, 空间步长h=0.01.

图 1 保费预期支付额与随机利率和时间的关系 Fig.1 Variation of expected premiums with random interest rates and time

图 2显示了随机利率r取定后,保费的预期支付额V和时间t的关系.然而,由图 2并不能说明保费随时间的增长而变大.

图 2 保费预期支付额和时间的关系 Fig.2 Variation of expected premiums with time

图 3显示了利率回归速度a对保费的影响,a的大小反映了利率的变化速度.a越大,利率变化越快,保费的预期支付额越小.从图 3可以进一步看出,当a很小时,保费与时间的关系不再是单调的关系,这也就解释了之前的结论.

图 3 回归速度对保费预期支付额的影响 Fig.3 Effect of regression speed on expected premiums

图 4显示了利率回归均值θ对保费的影响.θ影响了利率最终趋于的值的大小,θ越大,说明随机利率最终趋于的值越大,也说明随机利率支付者的损失越少,所需支付的保费相应越少.

图 4 回归均值对保费预期支付额的影响 Fig.4 Effect of regression means on expected premiums

图 5显示了信用等级变换边界ra对保费的影响.由图 5可以看出,ra越大,保费的预期支付额越小.这是因为信用等级变换边界越大,在发生信用等级变换时对于浮动利率支付方来说损失就越少,故而保费的预期支付额就越小.

图 5 信用等级边界对保费预期支付额的影响 Fig.5 Effect of credit rating border on expected premiums
3 结语

在利率互换合约实施过程中,利率下跌造成固定利率支付方信用等级下降,因此浮动利率支付方与第三方机构签订一份担保合约.通过分析利率信用等级变换的首达时刻以及由于利率互换过程信用等级下降所产生的利差,给出了合约现金流的表达式,并且将低等级下零息票的价值作为合约价值的自变量,利用结构化方法推导出满足合约价值的二维偏微分方程.通过坐标变换和计价单位转换将问题化为一般抛物型偏微分方程的初边值问题,进而得到解析解的表达式.最后, 利用有限差分方法的显式格式得到问题的数值解,得出了合约价值随时间变化的关系图,并且分析了各个参数对合约价值的影响.结果表明, 合约价值与时间的单调关系与参数a有关,并且与参数aθra有着明显的单调递减关系.

参考文献
[1]
李丹, 伦杭, 聂逆, 等. 国际三大评级机构信用评级定义及方法研究[J]. 征信, 2013(8): 47
LI Dan, LUN Hang, NIE Ni, et al. Study on definitions and methods of credit rating for the three major international rating agencies[J]. Zhengxin, 2013(8): 47 DOI:10.3969/j.issn.1674-747X.2013.08.013
[2]
JARROW R A, LANDO D, TURNBULL S M. A Markov model for the term structure of credit risk spreads[J]. Review of Financial Studies, 1997, 10(2): 481 DOI:10.1093/rfs/10.2.481
[3]
梁进, 肖承志. 具有信用等级迁移风险的零息债券定价[J]. 同济大学学报(自然科学版), 2015, 43(8): 1284
LIANG Jin, XIAO Chengzhi. Valuation of zero-coupon bonds with credit rating migration risk[J]. Journal of Tongji University(Natural Science), 2015, 43(8): 1284
[4]
曾楚焜.基于结构化方法的信用等级迁移对债券及信用利差期权定价影响研究[D].上海: 同济大学, 2016.
ZENG Chukun. Impact analysis of credit rating migration on valuation of corporate bonds and credit spread option by structure approach[D]. Shanghai: Tongji University, 2016.
[5]
HU Bei, LIANG Jin, WU Yuan. A free boundary problem for corporate bond with credit rating migration[J]. Journal of Mathematical Analysis and Applications, 2015, 428: 896 DOI:10.1016/j.jmaa.2015.03.040
[6]
LIANG Jin, ZHANG Xudan, ZHAO Yuejuan. Utility indifference valuation of corporate bond with credit rating migration by structure approach[J]. Frontiers of Mathematics in China, 2015, 10(6): 1389 DOI:10.1007/s11464-015-0445-3
[7]
马俊美, 梁进. 一篮子信用违约互换定价的偏微分方程方法[J]. 高校应用数学学报:A辑, 2008, 23(4): 427
MA Junmei, LIANG Jin. Valuation of basket credit default swaps by partial differential equation method[J]. Applied Mathematics: A, 2008, 23(4): 427
[8]
周鹏, 梁进. 信用违约互换的定价方法[J]. 高校应用数学学报:A辑, 2007, 22(3): 311
ZHOU Peng, LIANG Jin. Valuation of credit default swap[J]. Applied Mathematics: A, 2007, 22(3): 311
[9]
吴森, 梁进. 抵押贷款信用违约互换的定价[J]. 高校应用数学学报:A辑, 2011, 26(3): 269
WU Sen, LIANG Jin. Valuation of mortgage loan CDS[J]. Applied Mathematics: A, 2011, 26(3): 269
[10]
梁进, 徐寅, 郭高月. 信用攸关的利率互换的定价[J]. 同济大学学报(自然科学版), 2011, 39(2): 299
LIANG Jin, XU Yin, GUO Gaoyue. Pricing for credit contingent interest rate swap[J]. Journal of Tongji University(Natural Science), 2011, 39(2): 299 DOI:10.3969/j.issn.0253-374x.2011.02.027
[11]
HUANG Haohan, HUANG Huaxiong, WANG E, et al. Credit contingent interest rate swap pricing[J]. Mathematics-in-Industry Case Studies, 2017, 8(1) DOI:10.1186/s40929-017-0015-x
[12]
姜礼尚, 徐承龙, 任学敏, 等. 金融衍生产品定价的数学模型与案例分析[M]. 2版. 北京: 高等教育出版社, 2013
JIANG Lishang, XU Chenglong, REN Xuemin, et al. Mathematical modelling and case analysis of financial derivatives pricing[M]. 2nd ed. Beijing: Higher Education Press, 2013
[13]
姜礼尚. 期权定价的数学模型和方法[M]. 2版. 北京: 高等教育出版社, 2008
JIANG Lishang. Mathematical modelling and methods of option pricing[M]. 2nd ed. Beijing: Higher Education Press, 2008
[14]
姜礼尚, 陈亚浙, 刘西垣, 等. 数学物理方程讲义[M]. 3版. 北京: 高等教育出版社, 2007
JIANG Lishang, CHEN Yazhe, LIU Xihuan, et al. Handouts of mathematical physics equations[M]. 3rd ed. Beijing: Higher Education Press, 2007