收益率曲线计算方法简介

1、收益率曲线计算方法浅析杨华目前债券绝大部分是银行间品种，只能在银行间交易，只有少量可以跨市场 交易，考虑到我国债券市场的情况，直接以市价估值显然不合适，一个原因是， 某些债券很少有交易，甚至一段时间都没有一笔发生，另外一个原因是，即使有 交易，价格的真实性和代表性也不能保证， 而报价也是如此，这样采用收益率曲 线进行估值是比较合适的。收益率曲线的生成主要包括以下几种方法：1）存续期限法（Duration）该方法由Macaulay提出，依照麦氏存续期限的定义，存续期限相同的债券， 不论息票利率如何、不论是附息债还是零息债券，皆视为具有相同有效期限的债 券，在其他因素如流动性等相同的情况下， 市场

2、对具有相同存续期限的附息债和 零息债所要求的到期收益率必须是相等的， 因此可以通过各债券的到期收益率和 久期绘制出来，图形上相当于将到期收益率曲线向左移动（并非平移），如下图所示：这种方法比较简单，但其对即期收益率的估计比较粗糙， 计算久期的本身就蕴含 了收益率曲线为水平形状的假设2） 般计量方法计量方法假定即期利率和时间因子存在着某种特定的函数关系，再以相应的计量方法对函数中的系数进行估计，从而得到一个适用于所有到期期限的即期利 率曲线。在即期利率方式下，债券价格有如下的表达式：NPV 二' C (1 y(ti)M (1 y(tN ) _tNi =1其中，PV为债券的现价(全价)，C

3、为票息，M为本金，N为剩余附息次数， ti为各期附息或还本的剩余时间，y(ti)为相应的即期利率。令D (tj ) = (1 y(ti )丄为折现函数，则NPV = ' C D(t) M D(tN)i 4我们可以假定函数D(t)可以近似表示为t的多项式函数(任意一个连续可微 函数可以用多项式逼近)，这样，债券现值也就成为一个多元多项式函数，将市 场上的各附息债券的现价、到期日期、票息等代入，就可以形成一组方程，其中 包括待解的多项式的系数，利用上述的计量方法算出这些系数，就可以得到各期 的即期收益率。这种方法的优点是，除了能够在统计上检验曲线的拟合情况以及参数的显著 性特征外，还可以得

4、到即期利率的解析表达式，根据即期利率同远期利率的内在 关系，能够方便的得到远期利率的表达式。但这种方法的缺点是模拟的误差较大， 虽然可以通过提高多项式的次数来解决，但会造成结果的不稳定。计量方法的具体算法如下：设 D(t)二a0 a1t a2t2 a3t3由于D(0)=1,所以ao=1,假设本金M=10Q元，票息不变(即为固定利率债券)， 则有：PV 二 CD(tJ CD(t2)(100 C)D(tN)2 323二 C(1 a1t1 a2t1 ash )(100 C)(1 a azt”ast” )2 2 2= (NC 100) (6 Ct?(C 100)tN)a(CbCt?(C 100”n )

5、a23 33(Ct1Ct2(C 100)tN )a3=PV = x0 x1a1 x2a2 x3a3对于每一个附息债券，我们可以计算出上面的Xo,X 1,X2,X 3,而对于交易所市场 上的债券集，就可以得到一组含待定系数的等式， 这样我们就可以通过最小二乘 估计或者极大似然估计得到这些待定系数的最优解。具体地，对于k个债券的情 形，我们得到k个线性等式：PV1Xi,。=Xi,iai +Xi,2a2 +Xi,3a3PV2 X2,0 = X2,iai + X2,2a2+ X2,3a3，ZXi,iXi,2Xi,3',P =PVi Xi,0、,p =a2(Xk,iXk,2Xk,3>PVk

6、 - Xk,0令X -PVk - Xk,0 = xk,iai' Xk,2a2 Xk,3a3在考虑进残差项，线性拟合方程为:P =X；，其中；,var( £ )=/lk ,lk为k阶单位矩阵。这样就可以得出待定系数向量B的估计值为:(1)1 =(XTX)XTP实际中，各债券对应的误差项£ i的方差可能并不相同，而是随着债券到期 期限越长，变得越大，再用(1)来估计，拟合的误差就会较大。现假设var( £ )=(T 2Q ,其中-2W1 ,0,020,W2 ,0,则B的相应估计值为:匕=(XTi 2x)_xT，p(2)上式中w可以去成如下的值：Wj2 =(dP

7、Vi)2 =(Dj PVJ2，其中D为债券i的修正久期。dyi根据计算得到的B,就可以计算各到期日的折现率，得出各个到期日的即期收益率，进而绘制出即期收益率曲线。而且对于参数的显著性，线性拟合的准确 度可以做相应的统计推断。总体上，可以按照如下的两个步骤进行：第一步：样本券的选择样本按照各期限区间上个债券的发行规模、剩余期限、流动性水平、交易活 跃程度等标准进行筛选。区间的制定可参考中央国债结算公司等， 将指标年限划 分为1年、2年、3年、5年、7年、10年、15年、20年等。以所有符合条件的 债券的历史交易数据为依据，按照交易量、换手率等指标进行排序；同等条件下， 优先选取交易量大而稳定、较

8、近期发行的固定利率债券。并定期对筛选的样本债 券进行调整。第二步：根据上述的算法，计算各期的即期利率并绘制曲线图。2)Spline 方法这种方法相当于将上面的折现函数 D(t)分成几段，然后对于每段采用上面的 计量方法，同时保持整段曲线的连续，函数及其微分函数在各段交接点处连续。 最后将折现函数表示成分段的多项式函数。采取这种方法，是为了提高拟合度， 同时减少多项式的次数，但对于在何处分段，有不同的方法，可以按照市场债券 到期的分布情况进行分段，或按照长、中、短期分段等。现假设将贴现函数分成 m段，且在每段上表示为3次多项式的函数，这样 在整个区间上的贴现函数可以表示为：D1(t),Tc,T1

9、D2(t),"T1,T2D(t)二Dm(t),r Tm4,Tm相应的待定系数向量B也有m个：B 1，B 2，B m,而贴现函数连续可微的条件表示为:'.23Tj(aj,1 -aj 十,1)+Tj (aj,2 aj*2) +Tj (aj,3 a芦,3)= 02(aj,1 -aj卅J +2&,2 -aj+J +3©,3 -ajM =0*(aj,2 a*) +旺佝,3 -aw) =0Bj =(aj,1，aj,2aj,3)Tj =1,，m 1以上的限制条件共3m-3个方程，3m个待定系数，有无数可行解。同前面一般计量方法，在k个债券的情况下同样可以得到k个线性等式，

10、不同的是，每个等式中，至多将包含3m个待定系数，这是因为每只债券不同期限的现金流将根据到期期限的不同选择相应的贴现函数， 即选择相应的待定系数，特别地，对跨越整个区间的长期债券来说，等式中将取满所有 3m个待定系数假设:一GT2T,，；)T,XiPVi = Xi,0 (Xi,1a1,1 Xi,2a1,2-(xi ,1 , Xi ,3m )'冷3&1,3)(X ,3m _2am,1 ' Xj,3mam,2 ' Xj,3mam,3)-xi,0 Xi再设：P =PVX1,0 AlPVk -Xk,0 JX 二并将限制条件写成:R =0的形式，其中:,0*2323T1,T

11、1 ,T1 , _T1, _T1 7 -Tl ,0,23230 0 T T T_T_T_Tm 4 m J m Jmd m 4mV,02 21 , 2T1,3T1 ,-1,-2，-3，0,为(3m-1) x 3m矩阵2 20, ,0,1,2Tmd,3Tm，-1,-2匚小 -3匚二0,1,3，0,-1,-3£,0, ,00,0,0,1,3Tm 4,0, T, -3Tm/我们同样得到拟合方程:P根据误差协方差矩阵的不同，我们可以得到在上述限制条件下 的估计为：PT4TT4TT4T4TXT=(X X) X P -(X X) R R(X X) R R(X X) X P或者:=(X HX)4X

12、14(X 14X)4RTR(X 14X)4RT4R(XX)4X HP这样就可以估计出所有3m个待定系数，从而得到整个区间的贴现函数和即 期收益率函数在得到即期收益率的多项式表达后， 我们就可以绘制即期收益率曲线， 并且 对市场上的债券进行定价， 从而发现被高估或者低估的债券， 对于市场定价偏离 一定程度的债券， 可以给出买入或卖出的投资建议， 当然这种投资建议没有考虑 流动性等因素， 而是仅从即期利率结构出发给出的。 除了上述的方法外， 还可以 通过更复杂的利率模型来生成，这里省略。上述不同的方法生成的收益率曲线会有所不同，同时收益率曲线的生成结 果也和样本的选择有关， 从样本范围的选取来看， 应按照经验选择成交相对比较 活跃的品种，同时期限尽可能平均分布并跨越较大的范围， 这里有几个问题， 首 先是选取的时间长度问题， 如果仅仅以某一天的收盘价为样本数据， 可能样本量 严重不足， 因此时间长度可以扩展为一周或者两周， 既可以使样本量增大， 期限 涵盖加宽，同时也使得交易活跃的券种的影响力加大（权数增大） 。其次 对