【FinE】债券久期和凸性

债券久期

久期(Duration)是由Macaulay在1938年提出，又被称为Macaulay久期，用来评估债券的平均还款期限，计算公式为
$$D=\sum _{t=1}^{T}t\frac{C_t}{(1+y{)}^t}/\sum _{t=1}^T\frac{C_t}{(1+y{)}^t}$$
方程符号表如下

符号	含义
$C_t$	各期现金流，当$t<T$时，$C_t=C$；当$t=T$时，$C_t=C+M$
$t$	现金流支付的时间
$y_t$	贴现率
$T$	债券的有效期，在债券发行之前，有效期等于债券期限；债券发行之后，有效期小于期限

案例

面值为$100$元，期限为$5$年，年利率为$10\%$，贴现率为$9\%$，每半年付息一次的债券，求该债券的久期

解析

根据债券久期公式
$$\begin{gathered} \{\begin{array}{rl}& M=100\\ & C=100\ast 10\%/2=5\\ & r=0.09\\ & T=10\end{array} \\ D=\sum _{t=1}^{T}t\frac{C_t}{(1+y{)}^t}/\sum _{t=1}^T\frac{C_t}{(1+y{)}^t} \end{gathered}$$

代码

#include<iostream>
#include<vector>
#include<cmath>
#include<cstring>

using namespace std;

double bond_duration(const vector<double> &cf, const double &y){
    double up=0.00;
    double down=0.00;
    for(int i=0; i<cf.size(); ++i){
        up+=(i+1)*cf[i]/pow(1+y, i+1);
        down+=cf[i]/pow(1+y, i+1);
    }
    return up/down; 
}

int main(){
    double M=100;
    double r=0.10;
    double y=0.09;
    int years=5;
    int tm=2;
    vector<double> cf;
    for(int i=0; i<tm*years; ++i){
        if(i<tm*years-1) cf.push_back(M*r/tm);
        else cf.push_back(M*r/tm+M);
    }
    cout<<"duration ="<<bond_duration(cf, y)<<endl;
    
    return 0;
}

修正久期

久期($D$)，到期收益率($y$)和债券价格($P$)之间的关系为
$$\frac{\Delta P}{P}=-\frac{D}{1+y}\Delta y$$
定义$D^{\ast }=\frac{D}{1+y}$为修正久期，则有
$$\frac{\Delta P}{P}=-D^{\ast }\Delta y$$
可以通过修正久期和YTM的变化快速求出债券价格的变化率.

债券凸性

只用$\frac{\Delta P}{P}$的一阶展开并不可以完全描述债券价格变化率，因此需要进行二阶展开. 定义债券的凸性$C$为债券价格$P$对收益率$y$的二阶导数与价格$P$的比值.
【连续型债券凸性】
$$C=\frac{1}{P}\cdot \frac{{\partial }^{2}P}{\partial y^2}$$
【离散型债券凸性】
$$C=\frac{1}{P(1+y{)}^2}\sum _{t=1}^{T}t(t+1)\frac{C_t}{(1+y{)}^t}$$

案例

债券面值是$100$，期限为$3$年，每年付息$1$次，年利率为$10\%$，到期收益率为$9\%$，计算债券的凸性.

解析

根据离散型债券凸性公式
$$\{\begin{array}{rl}& M=100\\ & T=3\\ & r=10\%\\ & y=9\%\\ & P=\sum _{t=1}^T\frac{C_t}{(1+y{)}^t}\\ & C=\frac{1}{P(1+y{)}^2}\sum _{t=1}^{T}t(t+1)\frac{C_t}{(1+y{)}^t}\end{array}$$

代码

#include<iostream>
#include<vector>
#include<cmath>
#include<cstring>

using namespace std;

double bond_convexity(const vector<double> &cf, const double &y){
    double P=0.00; // 债券价格
    for(int i=0; i<cf.size(); ++i){
        P+=cf[i]/pow(1+y, i+1);
    }
    double C=0.00; // 债券凸性
    for(int i=0; i<cf.size(); ++i){
        C+=(i+1)*(i+2)*cf[i]/pow(1+y, i+1);
    }
    return C/(P*pow(1+y, 2));
}

int main(){
    double M=100;
    double r=0.10;
    double y=0.09;
    int years=3;
    int tm=1;
    vector<double> cf;
    for(int i=0; i<tm*years; ++i){
        if(i<tm*years-1) cf.push_back(M*r/tm);
        else cf.push_back(M*r/tm+M);
    }
    cout<<"convexity ="<<bond_convexity(cf, y)<<endl;
    
    return 0;
}

近似关系

如果一年内支付了$m$次利息，则凸性的年度值为

凸性的年度值=每年$m$期内的凸性/$m^2$

债券价格变化与收益率变化之间的关系可以近似表示为
$$\frac{dP}{P}=D^{\ast }dy+\frac{1}{2}C(dy{)}^2$$
当收益率变化不大时，收益率变化与价格变化之间的关系可以近似表示为
$$\frac{\Delta P}{P}=-D^{\ast }\Delta y+\frac{1}{2}C(\Delta y{)}^2$$