(完整版)收益率曲线计算方法分析

收益率曲线计算方法浅析

目前债券绝大部分是银行间品种，只能在银行间交易，只有少量可以跨市场交易，考虑到我国债券市场的情况，直接以市价估值显然不合适，一个原因是，某些债券很少有交易，甚至一段时间都没有一笔发生，另外一个原因是，即使有交易，价格的真实性和代表性也不能保证，而报价也是如此，这样采用收益率曲线进行估值是比较合适的。收益率曲线的生成主要包括以下几种方法：1）存续期限法（Duration）

该方法由Macaulay提出，依照麦氏存续期限的定义，存续期限相同的债券，不论息票利率如何、不论是附息债还是零息债券，皆视为具有相同有效期限的债券，在其他因素如流动性等相同的情况下，市场对具有相同存续期限的附息债和零息债所要求的到期收益率必须是相等的，因此可以通过各债券的到期收益率和久期绘制出来，图形上相当于将到期收益率曲线向左移动（并非平移），如下图所示：

存续期限债券剩余期限

这种方法比较简单，但其对即期收益率的估计比较粗糙，计算久期的本身就蕴含了收益率曲线为水平形状的假设

2）一般计量方法

计量方法假定即期利率和时间因子存在着某种特定的函数关系，再以相应的计量方法对函数中的系数进行估计，从而得到一个适用于所有到期期限的即期利

率曲线。

在即期利率方式下，债券价格有如下的表达式：

∑=--+++=N i t N t i

N i t y M t y C PV 1))(1())(1(

其中，PV 为债券的现价（全价），C 为票息，M 为本金，N 为剩余附息次数，t i 为各期附息或还本的剩余时间，y(t i )为相应的即期利率。

令i

t i i t y t D -+=))(1()(为折现函数，则 )()(1N N

i i t D M t D C PV ⨯+⨯=∑=

我们可以假定函数D(t)可以近似表示为t 的多项式函数（任意一个连续可微函数可以用多项式逼近），这样，债券现值也就成为一个多元多项式函数，将市场上的各附息债券的现价、到期日期、票息等代入，就可以形成一组方程，其中包括待解的多项式的系数，利用上述的计量方法算出这些系数，就可以得到各期的即期收益率。

这种方法的优点是，除了能够在统计上检验曲线的拟合情况以及参数的显著性特征外，还可以得到即期利率的解析表达式，根据即期利率同远期利率的内在关系，能够方便的得到远期利率的表达式。但这种方法的缺点是模拟的误差较大，虽然可以通过提高多项式的次数来解决，但会造成结果的不稳定。

计量方法的具体算法如下：

设332210)(t a t a t a a t D +++=

由于D(0)=1,所以ɑ0=1,假设本金M=100元，票息不变（即为固定利率债券），

则有：

33221103

332312222211213

32213132121121))100(())100(())100(()100()

1)(100()1()

()100()()(a x a x a x x PV a t C Ct Ct a t C Ct Ct a t C Ct Ct NC t a t a t a C t a t a t a C t D C t CD t CD PV N N N N N N N +++=⇒+++++++++++++=+++++++++=++++= 对于每一个附息债券，我们可以计算出上面的x 0,x 1,x 2,x 3,而对于交易所市场

上的债券集，就可以得到一组含待定系数的等式，这样我们就可以通过最小二乘

估计或者极大似然估计得到这些待定系数的最优解。具体地，对于k 个债券的情形，我们得到k 个线性等式：

⎪⎪⎩

⎪⎪⎨⎧++=-++=-++=-33,22,11,0,33,222,211,20,2233,122,111,10,11a x a x a x x PV a x a x a x x PV a x a x a x x PV k k k k k

令⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛--=⎪⎪⎪⎭⎫ ⎝⎛=3210,0,113,2,1,3,12,11,1,,a a a x PV x PV P x x x x x x X k k

k k k β 在考虑进残差项，线性拟合方程为：

εβ+=X P ，其中⎪⎪⎪⎪⎪⎭

⎫ ⎝⎛=k εεεε 21,var(ε)=σ2I k ,I k 为k 阶单位矩阵。 这样就可以得出待定系数向量β的估计值为：

P X X X T T 1)(-=β （1） 实际中，各债券对应的误差项εi 的方差可能并不相同，而是随着债券到期

期限越长，变得越大，再用（1）来估计，拟合的误差就会较大。

现假设var(ε)=σ2Ω,其中⎪

⎪⎪⎪⎪⎭⎫ ⎝

⎛=Ω22221,0,,00,,,00,,0,k w w w ，则β的相应估计值为：

P X X X T T 111)(---ΩΩ=β （2） 上式中w i 可以去成如下的值：

222)()(i i i

i i PV D dy dPV w ⨯==，其中D i 为债券i 的修正久期。 根据计算得到的β，就可以计算各到期日的折现率，得出各个到期日的即期收益率，进而绘制出即期收益率曲线。而且对于参数的显著性，线性拟合的准确度可以做相应的统计推断。总体上，可以按照如下的两个步骤进行：

第一步：样本券的选择

样本按照各期限区间上个债券的发行规模、剩余期限、流动性水平、交易活跃程度等标准进行筛选。区间的制定可参考中央国债结算公司等，将指标年限划分为1年、2年、3年、5年、7年、10年、15年、20年等。以所有符合条件的债券的历史交易数据为依据，按照交易量、换手率等指标进行排序；同等条件下，优先选取交易量大而稳定、较近期发行的固定利率债券。并定期对筛选的样本债券进行调整。

第二步：根据上述的算法，计算各期的即期利率并绘制曲线图。

2）Spline 方法

这种方法相当于将上面的折现函数D(t)分成几段，然后对于每段采用上面的计量方法，同时保持整段曲线的连续，函数及其微分函数在各段交接点处连续。最后将折现函数表示成分段的多项式函数。采取这种方法，是为了提高拟合度，同时减少多项式的次数，但对于在何处分段，有不同的方法，可以按照市场债券到期的分布情况进行分段，或按照长、中、短期分段等。

现假设将贴现函数分成m 段，且在每段上表示为3次多项式的函数，这样在整个区间上的贴现函数可以表示为：

⎪⎪⎩⎪⎪⎨⎧∈∈∈=-]

,[),(],[),(],[),()(1212

101m m m T T t t D T T t t D T T t t D t D 相应的待定系数向量β也有m 个：β1，β2，…，βm,而贴现函数连续可微

的条件表示为：

⎪⎪⎪⎩

⎪⎪⎪⎨⎧-===-+-=-+-+-=-+-+-++++++++1,,1),,(0

)(3)(0)(3)(2)(0)()()(3,2,1,3,13,2,12,3,13,22,12,1,11,3,13,32,12,21,11,m j a a a a a T a a a a T a a T a a a a T a a T a a T T j j j j j j j j j j j j j j j j j j j j j j j j j j β 以上的限制条件共3m-3个方程，3m 个待定系数，有无数可行解。同前面一般计量方法，在k 个债券的情况下同样可以得到k 个线性等式，不同的是，每个等式中，至多将包含3m 个待定系数，这是因为每只债券不同期限的现金流将