QuantLib 金融計算——案例之普通利率互換分析（2）

来源:https://www.cnblogs.com/xuruilong100/archive/2020/05/22/12940396.html
-Advertisement-
Play Games

[TOC] QuantLib 金融計算——案例之普通利率互換分析（2） 概述 QuantLib 中涉及利率互換的功能大致分為兩大類： 對存續的利率互換合約估值； 根據利率互換合約的成交報價推算隱含的期限結構。 這兩類功能是緊密聯繫的，根據最新報價推算出的期限結構通常可以用來對存續合約進行估值。 本文 ...


目錄

QuantLib 金融計算——案例之普通利率互換分析（2）

概述

QuantLib 中涉及利率互換的功能大致分為兩大類：

  • 對存續的利率互換合約估值；
  • 根據利率互換合約的成交報價推算隱含的期限結構。

這兩類功能是緊密聯繫的，根據最新報價推算出的期限結構通常可以用來對存續合約進行估值。

本文接下來介紹如何根據合約報價推算出隱含的利率期限結構，並以建信金科的技術文檔 《利率互換貼現因數曲線的構造模型》 中圖 16 的結果作為比較基準。

圖 16 的結果：

合約條款

通過合約報價推算期限結構的過程稱為“bootstrap”，其思想和實踐非常類似於理論證明中用到的“數學歸納法”，大體過程如下：

  1. 首先將要用到的已知利率和金融工具根據期限升序排列；
  2. 假設已經求得期限結構上的第 \(n\) 個值——\(TS_n\)，對應於第 \(n\) 個報價；
  3. 令 \(TS_{n+1}\) 是個待定參數 \(x\)，並給定一個初值；
  4. 用已知的期限結構數據——\(TS_1,\dots,TS_n,x\)，對第 \(n+1\) 個金融工具進行估值；
  5. 調整 \(x\)，使得估值結果與報價達到一致，不存在套利空間；
  6. 此時，\(x\) 便是 \(TS_{n+1}\)；
  7. 以此類推。

其中 \(TS\) 可以是即期利率、遠期利率、貼現因數三者中的任意一個，而可用的插值方法更是五花八門，例如線性插值、樣條插值、對數線性插值和常數插值等等。兩個維度相互搭配可以產生非常多的組合，QuantLib 通過模板技術實現兩個維度的自由搭配，具體選擇哪種組合要視業務需要而定。

需要註意的是，利率互換的估值對合約條款比較敏感。

示例中的合約均是 Shibor 3M 的利率互換，條款細則如下：

  • 浮動利率：Shibor 3M
  • 利差：0.0%
  • 估值日期：2020-01-15
  • 結算延遲：1 天
  • 重置延遲：1 天
  • 浮動端支付頻率：季度
  • 浮動端天數計算規則：ACT/360
  • 固定端支付頻率：季度
  • 固定端天數計算規則：ACT/365
  • 日曆：中國銀行間市場
  • 工作日轉換規則：Modified Following（MFL）

實踐

完整的代碼請見 QuantLibEx 項目的 example.cpp 文件。

using namespace QuantLib;
using namespace std;

Calendar calendar = China(China::IB);
Date today(15, January, 2020);
Settings::instance().evaluationDate() = today;

Natural delayDays = 1;

Date settlementDate = calendar.advance(
    today, delayDays, Days);
// must be a business day
settlementDate = calendar.adjust(settlementDate);

cout << "Today: " << today << endl;
cout << "Settlement date: " << settlementDate << endl;

Today: January 15th, 2020
Settlement date: January 16th, 2020

設置 RateHelper

QuantLib 中 bootstrap 計算的核心是為 PiecewiseYieldCurve 模板類配置 RateHelper 對象，不同的金融工具要使用對應的派生類。對於已知利率通常用 DepositRateHelper 類，而普通互換則用 SwapRateHelper 類。

示例沒有使用 QuantLib 提供的 Shibor 類，而是自己根據合約重新配置了一個對象。（查看源代碼的話，這實際上正是 QuantLib 中 IborIndex 派生類的普遍構造方式）

另外，Actual365_25 是 QuantLib 中未提供的，要自己實現，幾乎就是 Actual365Fixed 的翻版。

DayCounter termStrcDayCounter = Actual365_25();

Period mn1(1, Months), mn3(3, Months), mn6(6, Months), mn9(9, Months),
    yr1(1, Years), yr2(2, Years), yr3(3, Years), yr4(4, Years),
    yr5(5, Years), yr7(7, Years), yr10(10, Years);

ext::shared_ptr<Quote>
    m1Rate(new SimpleQuote(2.7990 / 100.0)),
    m3Rate(new SimpleQuote(2.8650 / 100.0)),
    s6mRate(new SimpleQuote(2.8975 / 100.0)),
    s9mRate(new SimpleQuote(2.9125 / 100.0)),
    s1yRate(new SimpleQuote(2.9338 / 100.0)),
    s2yRate(new SimpleQuote(3.0438 / 100.0)),
    s3yRate(new SimpleQuote(3.1639 / 100.0)),
    s4yRate(new SimpleQuote(3.2805 / 100.0)),
    s5yRate(new SimpleQuote(3.3876 / 100.0)),
    s7yRate(new SimpleQuote(3.5575 / 100.0)),
    s10yRate(new SimpleQuote(3.7188 / 100.0));

Handle<Quote>
    m1RateHandle(m1Rate),
    m3RateHandle(m3Rate),
    s6mRateHandle(s6mRate),
    s9mRateHandle(s9mRate),
    s1yRateHandle(s1yRate),
    s2yRateHandle(s2yRate),
    s3yRateHandle(s3yRate),
    s4yRateHandle(s4yRate),
    s5yRateHandle(s5yRate),
    s7yRateHandle(s7yRate),
    s10yRateHandle(s10yRate);

DayCounter depositDayCounter = Actual360();

ext::shared_ptr<RateHelper>
    m1(new DepositRateHelper(
        m1RateHandle, mn1, delayDays, calendar,
        ModifiedFollowing, false, depositDayCounter)),
    m3(new DepositRateHelper(
        m3RateHandle, mn3, delayDays, calendar,
        ModifiedFollowing, false, depositDayCounter));

Frequency fixedLegFreq = Quarterly;
BusinessDayConvention fixedLegConv = ModifiedFollowing;
DayCounter fixedLegDayCounter = Actual365Fixed();

ext::shared_ptr<IborIndex> shiborIndex(
    new IborIndex(
        "Shibor", mn3, delayDays, CNYCurrency(),
        calendar, Unadjusted, false, Actual360()));

ext::shared_ptr<RateHelper>
    s6m(new SwapRateHelper(
        s6mRateHandle, mn6, calendar, fixedLegFreq, fixedLegConv,
        fixedLegDayCounter, shiborIndex)),
    s9m(new SwapRateHelper(
        s9mRateHandle, mn9, calendar, fixedLegFreq, fixedLegConv,
        fixedLegDayCounter, shiborIndex)),
    s1y(new SwapRateHelper(
        s1yRateHandle, yr1, calendar, fixedLegFreq, fixedLegConv,
        fixedLegDayCounter, shiborIndex)),
    s2y(new SwapRateHelper(
        s2yRateHandle, yr2, calendar, fixedLegFreq, fixedLegConv,
        fixedLegDayCounter, shiborIndex)),
    s3y(new SwapRateHelper(
        s3yRateHandle, yr3, calendar, fixedLegFreq, fixedLegConv,
        fixedLegDayCounter, shiborIndex)),
    s4y(new SwapRateHelper(
        s4yRateHandle, yr4, calendar, fixedLegFreq, fixedLegConv,
        fixedLegDayCounter, shiborIndex)),
    s5y(new SwapRateHelper(
        s5yRateHandle, yr5, calendar, fixedLegFreq, fixedLegConv,
        fixedLegDayCounter, shiborIndex)),
    s7y(new SwapRateHelper(
        s7yRateHandle, yr7, calendar, fixedLegFreq, fixedLegConv,
        fixedLegDayCounter, shiborIndex)),
    s10y(new SwapRateHelper(
        s10yRateHandle, yr10, calendar, fixedLegFreq, fixedLegConv,
        fixedLegDayCounter, shiborIndex));

vector<ext::shared_ptr<RateHelper>> instruments;

instruments.push_back(m1);
instruments.push_back(m3);
instruments.push_back(s6m);
instruments.push_back(s9m);
instruments.push_back(s1y);
instruments.push_back(s2y);
instruments.push_back(s3y);
instruments.push_back(s4y);
instruments.push_back(s5y);
instruments.push_back(s7y);
instruments.push_back(s10y);

Bootstrap

Bootstrap 的過程很簡單，這裡選用 PiecewiseYieldCurve<ForwardRate, BackwardFlat>，與示例一致，將得到一個“階梯狀”的遠期期限結構。

ext::shared_ptr<YieldTermStructure> termStrc(
    new PiecewiseYieldCurve<ForwardRate, BackwardFlat>(
        today,
        instruments,
        termStrcDayCounter));

驗證

Date curveNodeDate = calendar.adjust(settlementDate + mn1);

cout << setw(4) << curveNodeDate - today << ", "
        << termStrc->discount(curveNodeDate) << ", "
        << termStrc->zeroRate(curveNodeDate, termStrcDayCounter, Continuous).rate() * 100.0
        << endl;
curveNodeDate = calendar.adjust(settlementDate + mn3);
cout << setw(4) << curveNodeDate - today << ", "
        << termStrc->discount(curveNodeDate) << ", "
        << termStrc->zeroRate(curveNodeDate, termStrcDayCounter, Continuous).rate() * 100.0
        << endl;
curveNodeDate = calendar.adjust(settlementDate + mn6);
cout << setw(4) << curveNodeDate - today << ", "
        << termStrc->discount(curveNodeDate) << ", "
        << termStrc->zeroRate(curveNodeDate, termStrcDayCounter, Continuous).rate() * 100.0
        << endl;
curveNodeDate = calendar.adjust(settlementDate + mn9);
cout << setw(4) << curveNodeDate - today << ", "
        << termStrc->discount(curveNodeDate) << ", "
        << termStrc->zeroRate(curveNodeDate, termStrcDayCounter, Continuous).rate() * 100.0
        << endl;
curveNodeDate = calendar.adjust(settlementDate + yr1);
cout << setw(4) << curveNodeDate - today << ", "
        << termStrc->discount(curveNodeDate) << ", "
        << termStrc->zeroRate(curveNodeDate, termStrcDayCounter, Continuous).rate() * 100.0
        << endl;
curveNodeDate = calendar.adjust(settlementDate + yr2);
cout << setw(4) << curveNodeDate - today << ", "
        << termStrc->discount(curveNodeDate) << ", "
        << termStrc->zeroRate(curveNodeDate, termStrcDayCounter, Continuous).rate() * 100.0
        << endl;
curveNodeDate = calendar.adjust(settlementDate + yr3);
cout << setw(4) << curveNodeDate - today << ", "
        << termStrc->discount(curveNodeDate) << ", "
        << termStrc->zeroRate(curveNodeDate, termStrcDayCounter, Continuous).rate() * 100.0
        << endl;
curveNodeDate = calendar.adjust(settlementDate + yr4);
cout << setw(4) << curveNodeDate - today << ", "
        << termStrc->discount(curveNodeDate) << ", "
        << termStrc->zeroRate(curveNodeDate, termStrcDayCounter, Continuous).rate() * 100.0
        << endl;
curveNodeDate = calendar.adjust(settlementDate + yr5);
cout << setw(4) << curveNodeDate - today << ", "
        << termStrc->discount(curveNodeDate) << ", "
        << termStrc->zeroRate(curveNodeDate, termStrcDayCounter, Continuous).rate() * 100.0
        << endl;
curveNodeDate = calendar.adjust(settlementDate + yr7);
cout << setw(4) << curveNodeDate - today << ", "
        << termStrc->discount(curveNodeDate) << ", "
        << termStrc->zeroRate(curveNodeDate, termStrcDayCounter, Continuous).rate() * 100.0
        << endl;
curveNodeDate = calendar.adjust(settlementDate + yr10);
cout << setw(4) << curveNodeDate - today << ", "
        << termStrc->discount(curveNodeDate) << ", "
        << termStrc->zeroRate(curveNodeDate, termStrcDayCounter, Continuous).rate() * 100.0
        << endl;

  33, 0.997441, 2.83629
  92, 0.992733, 2.89565
 183, 0.985631, 2.88875
 275, 0.978375, 2.90377
 369, 0.970721, 2.94142
 733, 0.940822, 3.03969
1097, 0.909515, 3.15787
1462, 0.877015, 3.27853
1828, 0.843917, 3.39077
2560, 0.778373, 3.57473
3654, 0.687352, 3.74755

與建信金科專家們的模型結果非常接近了，只有一個日期出現了不一致。

差異可能的來源

由於工作日轉換規則是 MFL，對假期比較敏感，QuantLib 中包含中國假期的日曆類是 China，它所記錄的假期可能和建信金科系統的假期不一致。

下一步

  • 分析國內市場上掛鉤的 FR007 的利率互換。
  • 分析國內市場上掛鉤的 LPR 的利率互換。

您的分享是我們最大的動力!

-Advertisement-
Play Games
更多相關文章
  • Deno會在短期內取代Node嗎？
    轉載請註明出處：葡萄城官網，葡萄城為開發者提供專業的開發工具、解決方案和服務，賦能開發者。 原文出處：https://blog.bitsrc.io/what-is-deno-and-will-it-replace-nodejs-a13aa1734a74 Deno是什麼？ Deno v1.0.0已於5 ...
  • BUAA_OO_2020_UNIT3
    一、JML初探 ​ 作為一種形式化語言，可以約束 代碼中類和方法的狀態和行為形成規格，通過將一系列具體代碼實現抽象成明確的行為介面，可以形成一種契約式編程模式， 設計者無需考慮實際的數據結構與演算法，可以聚焦於程式的整體邏輯， 形式化語言的無二義性能讓實現者準確理解介面功能，根據問題需要選擇合適的實現 ...
  • 消息隊列RabbitMQ安裝
    首先安裝Erlang環境 因為 RabbitMQ 需要 erlang 環境的⽀持，所以必須先安裝 erlang 。 如果只是使用RabbitMQ，個人推薦使用RabbitMQ公司維護的 "erlang" 版本，該版本只保留了與RabbltMQ相關的功能， centOS6與7版本的都有，還有erlan ...
  • 手把手教系列之FIR濾波器設計
    【導讀】：前面的文章介紹了移動平均濾波器、IIR濾波器、梳狀濾波器，今天來談談FIR濾波器的設計實現。 本篇文章依然採用4W1H進行描述，從 What Why Where When How 幾個維度展開。為了便於理解4W1H，依然把5W1H的圖附上。 FIR濾波器之What？ LTI線性時不變系統沖 ...
  • 手把手教系列之梳狀濾波器設計實現
    [導讀]：前面一篇文章關於IIR/移動平均濾波器設計的文章。本文來聊一聊陷波濾波器，該濾波器在混入諧波干擾時非常有用，演算法簡單，實現代價低。本文來一探其在機理、應用場景。 註：儘量在每篇文章寫寫摘要，方便閱讀。信息時代，大家時間都很寶貴，如此亦可節約粉絲們的寶貴時間。 前文所說學習的倡導2W1H原則 ...
  • 《java併發編程藝術》學習筆記
    本文是學習筆記中的思維導圖，感覺思維導圖是個好東西，可以加深記憶、構建知識體系、還能為將來的複習提供便利。 完整筆記的地址在文末。說明下，併發編程的底層原理和JMM我打算背下來。 完整筆記地址: "java併發編程藝術學習筆記gitbook" "java併發編程藝術學習筆記github" 下麵這個是 ...
  • JAVA連載117-反射的應用與打破封裝新
    一、反射的概念 1.反射的概念是有Smith在1982年首次提出，主要是指程式可以訪問、檢測和修改它本身狀態或者行為的一種能力，並且根據自身行為的狀態和結果，調整或者修改應用所描述行為的狀態和相關的語義。 2.Java中，反射是一種強大的工具，它使您能夠創建靈活的代碼，這些代碼可以在運行時進行裝配， ...
  • 彙編實驗十（王爽）
    設計子程式 子程式一：在指定的位置，用指定的顏色，顯示一個用0結束的字元串 舉例：在屏幕的8行3列，用綠色顯示data段中的字元串 assume cs:code data segment db 'Welcome to masm!',0 data ends code segment start: mo ...
一周排行
    -Advertisement-
    Play Games
  • 移動開發（一）：使用.NET MAUI開發第一個安卓APP
    移動開發（一）：使用.NET MAUI開發第一個安卓APP 對於工作多年的C#程式員來說，近來想嘗試開發一款安卓APP，考慮了很久最終選擇使用.NET MAUI這個微軟官方的框架來嘗試體驗開發安卓APP，畢竟是使用Visual Studio開發工具，使用起來也比較的順手，結合微軟官方的教程進行了安卓 ...
  • wpf ToggleButton選中效果和一個登錄界面
    前言 QuestPDF 是一個開源 .NET 庫，用於生成 PDF 文檔。使用了C# Fluent API方式可簡化開發、減少錯誤並提高工作效率。利用它可以輕鬆生成 PDF 報告、發票、導出文件等。 項目介紹 QuestPDF 是一個革命性的開源 .NET 庫，它徹底改變了我們生成 PDF 文檔的方 ...
  • 彈幕樹洞項目功能新增篇
    項目地址 項目後端地址： https://github.com/ZyPLJ/ZYTteeHole 項目前端頁面地址： ZyPLJ/TreeHoleVue (github.com) https://github.com/ZyPLJ/TreeHoleVue 目前項目測試訪問地址： http://tree ...
  • 第27篇 sqlserver2022詳細安裝步驟
    話不多說，直接開乾 一.下載 1.官方鏈接下載： https://www.microsoft.com/zh-cn/sql-server/sql-server-downloads 2.在下載目錄中找到下麵這個小的安裝包 SQL2022-SSEI-Dev.exe，運行開始下載SQL server； 二. ...
  • .NET 開源高性能 MQTT 類庫
    前言 隨著物聯網（IoT）技術的迅猛發展，MQTT（消息隊列遙測傳輸）協議憑藉其輕量級和高效性，已成為眾多物聯網應用的首選通信標準。 MQTTnet 作為一個高性能的 .NET 開源庫，為 .NET 平臺上的 MQTT 客戶端與伺服器開發提供了強大的支持。 本文將全面介紹 MQTTnet 的核心功能 ...
  • Serilog文檔翻譯系列（六） - 可用的接收器、增強器、格式化輸出
    Serilog支持多種接收器用於日誌存儲，增強器用於添加屬性，LogContext管理動態屬性，支持多種輸出格式包括純文本、JSON及ExpressionTemplate。還提供了自定義格式化選項，適用於不同需求。 ...
  • 警惕 Visual Studio 屬性求值副作用導致邏輯不符合預期
    目錄簡介獲取 HTML 文檔解析 HTML 文檔測試參考文章 簡介 動態內容網站使用 JavaScript 腳本動態檢索和渲染數據，爬取信息時需要模擬瀏覽器行為，否則獲取到的源碼基本是空的。 本文使用的爬取步驟如下： 使用 Selenium 獲取渲染後的 HTML 文檔 使用 HtmlAgility ...
  • [使用目前最新版]HybridCLR6.9.0+YooAsset2.2.4實現純C# Unity熱更新方案 (一)
    1.前言 什麼是熱更新 游戲或者軟體更新時，無需重新下載客戶端進行安裝，而是在應用程式啟動的情況下，在內部進行資源或者代碼更新 Unity目前常用熱更新解決方案 HybridCLR，Xlua，ILRuntime等 Unity目前常用資源管理解決方案 AssetBundles，Addressable， ...
  • 在 ASP.NET Core Web API 中使用操作篩選器統一處理通用操作
    本文章主要是在C# ASP.NET Core Web API框架實現向手機發送驗證碼簡訊功能。這裡我選擇是一個互億無線簡訊驗證碼平臺，其實像阿裡雲，騰訊雲上面也可以。 首先我們先去 互億無線 https://www.ihuyi.com/api/sms.html 去註冊一個賬號 註冊完成賬號後，它會送 ...
  • 第28篇 如何.net中實現高效可靠數據同步api
    通過以下方式可以高效，並保證數據同步的可靠性 1.API設計 使用RESTful設計，確保API端點明確，並使用適當的HTTP方法（如POST用於創建，PUT用於更新）。 設計清晰的請求和響應模型，以確保客戶端能夠理解預期格式。 2.數據驗證 在伺服器端進行嚴格的數據驗證，確保接收到的數據符合預期格 ...
所有分類