非高斯单因子短期利率模型正则化参数估计.pdfVIP

第29卷第6期 计 算 物 理 Vol29．NO．6 CHINESE OF 2012年11月 JOURNALCOMPUTATIONALPHYSICS Nov．，2012 文章编号：1001—246X(2012)06-0837-08 非高斯单因子短期利率模型正则化参数估计 江 良1’2， 徐承龙1 (1．同济大学数学系，上海200092；2莆田学院数学系，福建351100) 摘要：应用正则化方法估计带约束条件非高斯单因子短期利率模型中含时间变量参数估计．基于变分原理，证 明该正则化方法具有稳定性和收敛性．由于该问题是一个带有约束条件的反问题，本文应用罚方法把原问题转化 为非约束条件问题以方便计算．最后，通过数值模拟确认该算法的有效性并给出实证结果． 关键词：非高斯短期利率模型；正则化方法；罚方法；光滑算子 中图分类号：0241．83；029 文献标识码：A O 引言 众所周知，对于单因子的短期利率模型可以分为高斯模型和非高斯模型两大类．前者具有很好的数学 模型口1．不过CIR模型虽然保证了瞬时利率非负，但是实证结果表明该模型很难较好地拟合市场数据，如 零息债券，这归咎于该模型中的参数是常数(详细参考本文的实证部分)．为了更好地拟合市场数据，Hull 方面该模型只允许长期回归均值是时间变量的函数，而其他所有的参数是常数．另一方面，在实际应用中， 基于约化模型定价风险资产价格时，需要考虑违约强度和短期瞬时利率之间关系，此时瞬时利率被假设基 于ECIR模型(或者CIR模型)¨1，这归咎于该模型保证瞬时利率是正的，这样就使得计算违约概率有意义． 因此估计ECIR模型的参数变得尤为重要． 使用(零息)债券的数据来估计模型中的参数．假设市场数据和理论价格完全匹配i这可导出带有约束 smoothing 定性第二类Volterra积分问题．虽然通过对积分方程求二阶导数可以推出带有约束条件具有稳定性的第二 类Votterra积分问题，但是过程需要对输入数据求二阶导数，即要求市场数据具有一定光滑性．然而输人数 据求导数是一个不稳定的过程哺1．众所周知，市场数据都含有噪声．Rogers∽1也明确地指出任何涉及到对于 市场数据直接求导数的数值结果是不可行的．因此求解上述问题可分为两类方法：1)通过光滑算子将输入 数据光滑化¨…；2)直接使用正则化方法求解．本文考虑后者的方法，这主要基于以下的原因：①正则化方 法具有一般性；②很容易使用一些数学技巧证明其稳定性和收敛性；③正则化方法能够通过正则化参数来 控制数值结果的稳定性及其精度．关于正则方法求解不适定问题可参考文献[1l一12]，其基本思想是通过 引入正则项及正则化参数控制原问题的不稳定性． 本文使用一阶微分算子作为正则项来控制数值结果的稳定性．证明该正则化方法的稳定性和收敛性． 由于研究的问题是带有约束条件的正则化问题，因此通过引入罚函数，把原问题转化为无约束条件的非线 性问题．这类技巧也适用于美式期权数值计算¨””3．其主要的思想是：通过增大罚函数的参数使得近似解 落在可行域中，即近似值解满足约束条件．在第三部分，数值模拟结果确认对于一个足够大的罚函数参数， 其近似数值解很好地逼近原问题的数值解．由于转化后的非线性问题在一些点上是不可导的，导致一些经 收稿日期：2012—03—11；修回Ft期：2012—06—28 作者葱介：江良(t978一)，男，博士生，主要研究计算数学及金融数学，E-rnail：ptjiangliang@yahoo，com．cn 万方数据 计 算 物 理 第29卷 典的数值方法无法直接使用(如牛顿迭代算法)．本文引入光滑算子，把非线性问题转化为具有光滑性的问 题，从而使得应用一些经典的算法(见本文第三部分)具有可行性．而数值模拟结果确认上述变化是可 行的． 1利率与债券定价模型 假设短期利率r在风险中性测度下满足如下的随机微分方程∞1 dr=(0(t)一ar)dt+Or4rdW(t)， (1) 2 到期日T支付为1单位的债券价格有如