21-利率期限结构模型[金融计算与建模]

第21章利率期限结构模型清华大学经管学院朱世武[e-mail]样本数据：论坛：利率期限结构模型简介利率期限结构相关符号表：(,):BtT在未来时间T到期的零息票债券在时间t的价格，即在未来时间T支付单位1的债券在时间t的价格。ˆ(,)Rt起息日为时间t，剩余到期期限为年的零息票债券利率。有：1(,)ˆ[1(,)]BttRt(,)Rt起息日为时间t，剩余到期期限为年的连续复合利率。有：(,)exp[(,)]BttRt(,,)FtsTs在时间t计算的，起息日为时间s，剩余到期期限为T-s的远期利率。有：(,)(,,)exp[()(,,)](,)fBtTBtsTTsFtsTsBts(,)fts在时间t计算的，起息日为时间s的瞬时远期利率。有：0(,)lim(,,)TsftsFtsTsln(,)(,)Btsftsstr即期利率，时间t计算的，剩余到期期限无限小时的零息票债券的连续符合内部收益率。有：0lim(,)trRtln(,)|tTtBtTrTtR起息日为时间t，剩余到期期限为年的连续复合利率。有：t(,)tRRtt(,)tT贴现债券价格在时间t的预期瞬间收益。(,)BtT(,)tT贴现债券价格在时间t的瞬时波动。(,)BtTˆ,WW标准布朗运动。(,)tT瞬间远期利率的波动。有：(,)ftT(,)(,)tTtTT(,)Vt贴现债券利率的波动。(,)Rt(,)iBnT重组树中，在第i种状态下，剩余到期期限为T的贴现债券在时间n的均衡价格。注意，与的定义不同，此处T表示的是剩余到期期限，而非到期日。(,)BtT利率期限结构的概念利率（interestrate）是经济和金融领域的一个核心变量，它实质上是资金的价格，反映了资金的供求关系。利率期限结构（termstructureofinterestrates），又称收益率曲线（yieldcurve），是指在相同风险水平下，利率与到期期限之间的关系，或者说是理论上的零息债券利率曲线。常见的利率期限结构有以下四种：贴现因子曲线（discountfactorcurve）：；零息票收益曲线（zero-couponyieldcurve），（常用）：或；远期利率曲线（forwardratescurve）：瞬时远期利率期限结构（instantaneousforwardtermstructure），（常用）：。(,)Btt(,)Rtˆ(,)Rt(,,)TFtsTs(,)sfts静态模型动态模型样条函数模型节约型模型指数样条法（Vasicek&Fong，1982）均衡模型套利模型Vasicek模型（Vasicek，1977）CIR模型（Cox、Ingersoll&Ross，1985）Ho-Lee模型（Ho&Lee，1986）Hull-White模型（Hull&White，1990）HJM模型（Heath,Jarrow&Morton，1992）Nelson-Siegel模型（Nelsen&Siegel，1987）Svensson扩展模型（Svensson，1994）B样条法,（Steeley，1991）多项式样条法（McCulloch，1971，1975）利率期限结构模型静态利率期限结构模型静态利率期限结构模型概述静态利率期限结构模型以当天市场的债券价格信息为基础，构造利率曲线函数，利用所构造的利率曲线得到理论价格来逼近债券的市场价格，从而得出符合当天价格信息的利率期限结构。静态利率期限结构模型最为常见的有样条函数模型和节约型模型，样条函数模型主要包括多项式样条法、指数样条法和B样条法，节约型模型的主要代表是Nelson-Siegel模型及其扩展模型。通常，使用静态模型拟合利率期限结构的具体过程如下：首先，从市场上选出一组无违约风险的附息债券。设该组附息债券在时间t的市场价格为，在时间s的现金流入为，其中，，j表示该组的第j支债券。jtP()jsFst由于期限结构指的是零息债券的收益率与其到期日间之关系，因此必须先调整“息票效应”（CouponEffect）。息票效应是指：对于剩余到期期限相同的债券来说，它们的到期收益率不仅与当前的利率期限结构有关，还与它们的票面利率水平有关。对于相同的即期利率期限结构而言，到期收益率是这些即期利率的加权平均，而权重是各个现金流的现值。于是，假想出贴现函数或零息票债券利率的具体形式，其中和为参数向量。然后利用假想出的具体形式，来推导附息债券的理论价格，当推导出的理论价格与给定的市场价格最为接近时，就可以估计出由和构成的参数向量，即：1(,)(;)Btsfst2(,)(;)Rtstgst1212*21ˆˆargmin()njjttjPP其中，是从模型ˆjtP()ˆ(;)jjtssPFfst或模型()ˆ(;)jjtssPFfst()ˆexp[()(;)]jjtssPFstgst推导出的附息债券理论价格。显然，债券样本中长期品种的价格波动性应大于短期品种，而由此带来的结果是：以上述方法中表示长期债券的定价误差往往大于短期债券。这就是在进行收益率曲线拟合时无法避免的样本异方差特征，导致的结果往往是收益率曲线在远端出现“过度拟合”（Overfitting）的情况，而在近端则无法很好地表现短期债的实际情况。为了解决这一问题，应该对短期债券赋予较高的权重，而对长期债券赋予较低的权重，从而允许长期债券存在较高的误差。在Bolder和Streliski(1999)的论文中，设定了如下的权重系数：1/1/jjjDurwDur而将参数的估计过程定义为：2*21ˆ()ˆargminjjnttjjPPwn多项式样条法多项式样条函数假设折现因子是到期期限s的多项式分段连续函数。在运用此函数时，仔细选择多项式的阶数是至关重要的。阶数的多少决定了利率曲线的平滑程度和拟合程度，同时也影响到待估参数的数量。本书将多项式样条函数的阶数定为3。这是因为，当多项式样条函数为二阶时，的二阶导数是离散的；当阶数过高（四阶或五阶）时，验证三阶或四阶导数是否连续的难度将增大，待估参数的数量也将增大。()Ds()Ds(2)()Ds一般选用如下形式的多项式样条函数：231111123222222333333()()()()...........DsabscsdsDsabscsdsDsDsabscsds11223[0,][,][,]sTsTTsTT注意，对于即期贴现率函数来说，显然有。另外，为了保证分段函数的平滑性以及在分段点的平滑过渡，必须保证贴现函数在整个定义域内连续且一、二阶可导，还需要满足如下约束条件：()Ds(0)1D1111221()()()()()()iiiiiiiiiiiiDTDTDTDTDTDT其中的一阶导数和二阶导数。(1)(2)(),()()iiiDTDTDT分别为例如，考虑30年期的贴现率函数，可以用三次多项式分段拟合如下：2301111235222223103333()[0,5]()()[5,10]()[10,30]DsabscsdssDsDsabscsdssDsabscsdss其中，函数必须满足以下的7个约束条件：()()055100(5)(5)(10)(10)0,1,2(0)1iiiiDDDDiD从而，我们可以将互相独立的参数缩减到5个：23011133235111232310111332()10,5()15-55,10()()15-510DsbscsdssDsbscsdssdssDsDsbscsdssdss33-1010,30dss指数样条法指数样条函数是VasicekandFong(1982)提出的。与在多项式样条函数部分所述的原因相同，也采用三阶指数形式样条函数，其形式为：231111123222222333333()()()()...........usususususususususDsabecedeDsabecedeDsDsabecede11223[0,][,][,]sTsTTsTT模型中，除了u也是一个参数，并且有明显的经济含义。VasicekandFong(1982)证明了如下等式：,,,iiiiabcd外，lim(0,)sufs即，u可以被认为是当前的起息日为未来无限远时的瞬间远期利率。同样，指数样条法也必须满足如下约束条件：1111221(0)1()()()()()()iiiiiiiiiiiiDDTDTDTDTDTDT(1)(2)(),()()iiiDTDTDT分别为其中，的一阶导数和二阶导数。选择样条函数的分段数量和取样条分界点在指数样条法中也同样十分重要，其方法可以参见多项式样条法。并且，指数样条模型也容易导致远期利率曲线不稳定。不同于多项式样条法的是，其参数估计必须采用非线性最优化。Nelson-Siegel模型及其扩展形式Nelson-Siegel模型可以由一个公式来说明，该公式的形式与那些描述动态利率的普通微分方程的解的表达式十分类似。该公式为：012111(0,)exp()exp()f其中，表示即期计算的，在未来时间时发生的瞬间远期利率。均为待估参数。利用(0,)f0012,,1以及01(0,)(0,)Rfsds可以得到：110121111exp()1exp()(0,)expR这就是Nelson-Siegel模型的基本表达形式。当固定时，通过的不同组合，利用这个模型，可以推出四种不同形状的零息票债券收益曲线：递增、递减、水平和倒置。但是，这个模型无法推导出形状更为复杂的收益曲线，例如V形收益曲线和驼峰收益曲线。012和为了克服上述缺点，Svensson(1994)将上述模型扩展如下：012311122(0,)exp()exp()exp()f于是，可以得到：1101211123221exp()1exp()(0,)exp1exp()expR-4.0-3.0-2.0-1.00.01.02.03.04.05.06.001234567891011121314151617181920%B0B1B2Yield动态利率期限结构模型动态利率期限模型包括均衡模型和套利模型。均衡模型是一种由均衡分析方法得出的模型，它从假设一些经济变量开始，推出短期无风险利率的一个过程，然后寻找该过程对债券价格和期权价格的含义。均衡模型利用以下三步来为利率或有要求权定价：利用已建立好的因子模型来推导出理论零息票债券收益率曲线。利用参考债券的市场价格来校准模型并推出模型的参数值。最后，利用已确定的参数来为金融衍生品定价。套利模型由无套利分析方法得出，它是利用市场上的价格信息来推导出利率随机微分方程的形式的。均衡模型根据状态变量集中随机变量的个数，可以将利率期限结构模型区分为单因素和两(多)因素模型两大类。一般单因素模型(,)(,)ttttdrtrdttrdW对取不同