目前许多代码和资源在进行债券定价时，大多以债券发行日作为贴现值的时间点，但在实际应用中我们常常需要对早就发行的债券进行定价，这就需要计算准确的现金流、现金流日距离当前贴现时间点的时间距离和不同时间距离的零息利率，这一过程会较多使用插值法，还涉及到时间的转换。此外，本文还进行了麦考利久期、修正久期和凸度的计算，并使用久期凸度分析了价格的敏感性。

目录

1 主要步骤概述

2 使用插值法计算特殊点零息利率，并画出收益率曲线

3  对债券进行精确定价

4 计算债券的麦考利久期、修正久期和凸度，并进行价格敏感性分析

国家开发银行发行了200亿元的“22国开05”债券（代码220205.IB），在“中国货币网”查找该债券的关键信息：面值、期限、票面利率、付息频次、起息日等。并请用python程序完成了下列任务：

1、登录“中国货币网”，查找“基准利率”中的“实时收益率曲线”，下载“政策性金融债”的实时收益率数据到excel，python读入这些数据，并利用数据作图将其可视化。

2、计算“22国开05”各个未来现金流的剩余期限，利用插值法得到这些时刻对应的即期利率（即零息债券的收益率）。然后，再利用前面得到的零息利率对“22国开05”债券进行精确定价（净价）。

3、利用前面得到的数据计算该债券的修正久期和凸性值；假定即期利率曲线（零息利率曲线）发生了平移，即每期都增加100个基点，请用久期和凸性方法重新估计债券价格。

以下数据为2022年5月10日国开行的不同剩余期限的基准债券的数据。我们以一下数据作为初始数据，并使用插值法，计算0~30年中特殊关键时点的零息利率，并画出图。

数据导入：

使用插值法将收益率补齐：

 画收益率曲线图：（这里用到了我最近发现的一个更好地数据可视化工具plotly，十分好用，后面另写文章介绍）

 输出：

债券价值=未来各期利息收入的现值合计+未来到期本金或售价的现值，债券价值是指进行债券投资时投资者预期可获得的现金流入的现值。

债券的现金流入主要包括利息和到期收回的本金或出售时获得的现金两部分。

因此对债券进行精确定价，需要知道该债券的每一现金流，以及可以获得的每一现金流距离贴现点的时间长度，以及这些时间长度所对应的零息利率。如上图中的公式，C(票息)和M（债券面值）为每一现金流，实际上每一期的C可能不同，ri对应不同时间长度的零息利率，也就是每一现金流有不同的r。

计算债券的现金流列表，每一现金流对应的零息利率，每一现金流距离指定时间点间的时间距离的函数：

 债券精确定价函数：

配置22国开05债券的信息并调用函数：

输出：

麦考利久期(Macaulay duration) 。久期的概念最早是麦考利 （Frederick Robertson Macaulay (1882.8.12–1970.3) )在1938年提出来的，所以又称麦考利久期（简记为D）。麦考利久期是使用加权平均数的形式计算债券的平均到期时间。它是债券在未来产生现金流的时间的加权平均,其权重是各期现值在债券价格中所占的比重。

修正久期是对于给定的到期收益率的微小变动，债券价格的相对变动与其麦考利久期的比例。这种比例关系是一种近似的比例关系，以债券的到期收益率很小为前提。是在考虑了收益率的基础上对麦考利久期进行的修正，是债券价格对于利率变动灵敏性的更加精确的度量。

凸性（ Convexity ）是收益率变化 1 %所引起的久期的变化。用来衡量债券价格收益率曲线的曲度。凸性越大，债券价格曲线弯曲程度越大，用修正久期度量的利率风险所产生的误差越大。

 麦考利久期代码实现：

 计算修正久期（涉及到到期收益率ytm的计算）：

计算凸性：

计算利率变化时的债券价格变动：

调用函数进行计算：

输出：