使用Java与QuantLib 进行雪球结构产品的评价

最近业界朋友委托我，开发一个雪球结构产品的评价程序，因为是业界IT公司，因此偏好使用Java语言来撰写。刚好我想练习一下使用SWIG封装的QuantLib链接库，因此答应用Java来开发，另外我也试了C#版的QuantLib，当作一个比较，结果与Java计算的一样。我发现中金固收的文件(【中金固收·固收+】超额险中求：雪球结构产品介绍与历史回测)写得还蛮清楚，因此借用他的图表(图一)来简单提示产品规格。

由于这个产品的偿付条件是路径相关的，而且也相当复杂，因此只能使用蒙地卡罗模拟法来评价。考虑到Basel III即将施行，这些产品未来需要依照FRTB的要求，计算SBM法市场风险资本的计提，因此最好以面向对象的方法来开发。我就比照QuantLib的格式，写了一个SnowBall对象。这个版本主要是用作测试，因此不打算搞得太复杂，日历的推算不考虑假日的问题。

图二是一个简单的测试案例，契约发行日2022/3/1，评价日也是2022/3/1，到期日为2023/3/1，一年期产品，每月比价敲入，每日比价敲出。期初标的价格100.0，契约本金100.0。标的资产收益率2%，以Act/365计息，采用不考虑假日的日历。行352~390，建立一个即期利率期限结构对象。行396~419，建立一个波动性期限结构对象。这两个对象都是QuantLib内建的基础对象，只需依其规格输入便可使用，非常的方便好用。

期初价格设为100.0，敲出价位为期初价格的110%，敲入价位为期初价格的70%，或有票息15%，一倍杠杆。仿真1024条路径，每天走一步。计算得到的价格为97.2496。一共花了0.4393秒，不算快。下一个版本考虑直接使用GPU来模拟，根据以往的经验，应该有10倍以上的成效。

图四是整个类别的架构与SnowBall的建构函数，图五是平价的核心NPV()函数。行59~95，判断契约是一年或两年，用来决定比价期数与模拟的步数。行110与111叫用MT19937随机数生成器与InverseTransform函数，来产生常态分配随机数。MT19937是一个高质量的随机数生成器，QuantLib有提供这两个对象。

图六、七开始进行模拟与评价，考虑3种情况。行113~199是考虑期初评价的情况，行201~328是考虑发行后与到期之后的情况，这里只显示第一种情况。行127~136使用几何布朗运动过程仿真资产价格。行140~165逐条路径比较，是否敲出、敲入。行167~187根据比较的结果，计算每条路径的偿付。评价时都以单位面值计算，最后在行198，计算平均值时，再把契约面值调整回去，得到契约的价格MTM。

需注意的是，如果敲出，契约当月结束，因此需以该月比价时点的即期利率折现，其他情况是以契约到期日的即期利率折现。这个版本只是展示计算逻辑，并没有进行优化。有兴趣想要进一步了解使用QuantLib进行金融计算的朋友，可以跟我连络，我正在推广QuantLib的商业使用。我目前担任昀腾金融科技的技术长，台湾金融研训院的2022年菁英讲座，另外我是台湾大学财金所兼任的实践教授，WeChat：andydong3137。