自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.当随机利率r≥ra，利率处于高等级；当r＜ra，利率处于低等级.初始时刻随机利率r0≥ra，即初始时刻利率处于高等级. (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≤t≤T}为由概率空间(Ω, $\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个不同信用等级下零息票的价值差. 高、低等级下利率分别用r1t和r2t表示，相应地零息票利率互换的价值分别为P1t和P2t.对于高等级而言，将低等级下零息票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的求解范围为r1＞ra，所以先将区域变换到坐标轴的右半部分，即作变换$\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, ra和P2t, 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: {y＞y0, 0≤t＜T}上的方程，如下所示: $ \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}$1、P2方向的步长相同，记Δ${\tilde r}$1=ΔP2=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] $ 式中:R1、R2为足够大的常数.在网格点(n, i, j)上记V的值为Vn, i, j=V(nτ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有着明显的单调递减关系.