全部版块 我的主页
论坛 经济学人 二区 外文文献专区
842 22
2022-05-31
英文标题:
《Incremental computation of block triangular matrix exponentials with
  application to option pricing》
---
作者:
Daniel Kressner, Robert Luce, Francesco Statti
---
最新提交年份:
2017
---
英文摘要:
  We study the problem of computing the matrix exponential of a block triangular matrix in a peculiar way: Block column by block column, from left to right. The need for such an evaluation scheme arises naturally in the context of option pricing in polynomial diffusion models. In this setting a discretization process produces a sequence of nested block triangular matrices, and their exponentials are to be computed at each stage, until a dynamically evaluated criterion allows to stop. Our algorithm is based on scaling and squaring. By carefully reusing certain intermediate quantities from one step to the next, we can efficiently compute such a sequence of matrix exponentials.
---
中文摘要:
我们以一种特殊的方式研究了块三角矩阵的矩阵指数的计算问题:从左到右逐块列计算。在多项式扩散模型中的期权定价背景下,这种评估方案的需要自然产生。在此设置中,离散化过程生成一系列嵌套块三角矩阵,并在每个阶段计算其指数,直到动态评估的标准允许停止。我们的算法基于缩放和平方。通过从一步到下一步仔细地重用某些中间量,我们可以有效地计算这样一系列矩阵指数。
---
分类信息:

一级分类:Mathematics        数学
二级分类:Numerical Analysis        数值分析
分类描述:Numerical algorithms for problems in analysis and algebra, scientific computation
分析和代数问题的数值算法,科学计算
--
一级分类:Quantitative Finance        数量金融学
二级分类:Pricing of Securities        证券定价
分类描述:Valuation and hedging of financial securities, their derivatives, and structured products
金融证券及其衍生产品和结构化产品的估值和套期保值
--

---
PDF下载:
-->
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2022-5-31 05:14:47
块三角矩阵指数的增量计算及其在期权定价中的应用*罗伯特·卢斯*弗朗西斯科·斯塔蒂*+2017年6月30日摘要我们以一种特殊的方式研究了块三角矩阵的矩阵指数计算问题:从左到右逐块列计算。在多项式微分模型中的期权定价背景下,这种评估方案的需要自然产生。在此设置中,离散化过程生成一系列嵌套块三角矩阵,并在每个阶段计算它们的分量,直到动态评估准则允许停止。我们的算法基于缩放和平方。通过在一个步骤到下一个步骤中小心地重用某些中间量,我们可以高效地计算这样一系列矩阵指数。关键词矩阵表达式、块三角矩阵、多项式微分模型、期权定价1引言我们研究了嵌套块三角矩阵序列的矩阵指数计算问题。为了给出预处理问题的公式,考虑块上三角r矩阵G,G,G。formGn的=G0,0G0,1··G0,nG1,1··G1,n。。。。。。Gn,n∈ Rdn×dn,(1)其中所有诊断块Gn,nare square。换句话说,矩阵giarisefromgi-1通过附加块列(并调整大小)。我们的目标是计算矩阵指数exp(G),exp(G),exp(G)。(2) 当然,可以使用标准技术单独计算每个指数(2)(参见【11】中的概述)。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2022-5-31 05:14:50
然而,矩阵指数序列(2)继承了矩阵Gnin(1)的嵌套结构,即。,*\'Ecole Polytechnique F'ed'erale de Lausane,Station 8,1015 Lausane,Switzerland({daniel.kressner,robert.luce,francesco.statti}@epfl.ch)+根据欧洲联盟第七框架计划(FP/2007-2013)/欧洲研究理事会第307465号拨款协议(ERC Grant Agreement n.307465-POLYTE)通过欧洲研究理事会支持的研究。exp(Gn-1) 是exp(Gn)的主要子矩阵。实际上,只需要计算exp(Gn)的last blockcolumn,本文的目标是解释如何以数字安全的方式实现这一点。在对角线块Gn、NAR的光谱分离的特殊情况下,Parlett方法[13]原则上产生了一个有效的计算方案:分别计算F0,0:=exp(G0,0)和F1,1:=exp(G1,1),然后exp(G)的缺失(1,2)块作为Sylvester方程G0,0X的唯一解X给出- XG1,1=F0,0G0,1- G0,1F1,1。以这种方式继续计算(2)所需的所有o f-对角线块可以从求解Sylvester方程中获得。然而,众所周知(见[8]第9章),只有当对角块的光谱被很好地分离时,Parlett方法才是数值安全的,因为所有涉及的Sylvester方程都是条件良好的。由于我们认为区块结构是固定的,因此施加此类条件将严重限制应用范围;我们下面讨论的应用程序肯定无法满足这一要求。指数增量计算的一类一般应用来自线性算子G:V的矩阵表示→ V受限于给定的有限维向量空间V的一系列嵌套的有限维子空间。更准确地说,我们从有限维子空间Vof V开始,以B为基础。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2022-5-31 05:14:53
接着,向量空间Vis扩展到V 五、 ···  V通过生成嵌套基序列B B B ···.假设GVn Vn对于所有n=0,1,考虑G相对于Bn的矩阵表示的序列。由于碱基的嵌套性,GNI由Gn构建-通过添加列来表示G toBn的操作-因此,我们得到了一个矩阵序列,其结构如(1)所示。上述塞纳里奥的一个具体例子出现在计算金融领域,即基于多项式差异模型的期权定价;参见【6】。正如我们在第3节中更详细地解释的那样,在这种情况下,G是随机微分方程(SDE)的生成器,V是多变量多项式的嵌套子空间。一些定价技术需要计算某些条件矩,这些条件矩可以从矩阵指数中提取出来(2)。虽然增加n可以更好地近似期权价格,但获得所需精度所需的n值通常是未知的。自适应选择n的算法依赖于整个序列的增量计算(2)。块三角矩阵的指数也在其他情况下进行了研究。对于two-by-two-block三角矩阵,Dieci和Papini在[4]中研究了条件问题,并在[3]中讨论了使用Pad'eapproximants到指数函数的缩放参数的选择。在ma矩阵也是oblock-Toeplitz的情况下,文[2]提出了一种快速求幂算法。本文的其余部分组织如下。在第2节中,我们详细描述了增量计算块三角矩阵指数的算法,如(1)所示。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2022-5-31 05:14:56
在第3节中,我们讨论了多项式微分模型,以及一些需要使用这种增量算法的定价技术。最后,数值结果在第4.2节增量缩放和平方中给出。由于共形分区块tr正则矩阵集形成一个代数,并且exp(Gn)是Gn中的多项式,因此矩阵exp(Gn)具有与Gn相同的块上三角结构,即exp(Gn)=exp(G0,0)* ··· *exp(G1,1)。。。。。。。。。*exp(Gn,n)∈ Rdn×dn。正如简介中所述,我们的目标是从左到右逐块列计算exp(Gn)block column。我们的算法基于缩放和平方方法,下面我们简要总结一下。2.1缩放和平方法概述缩放和平方法使用有理函数近似指数函数,通常包括三个步骤。用rk表示,m(z)=pk,m(z)qk,m(z)表示(k,m)-Pad'e对经验函数的逼近,这意味着分子是k次多项式,而表示符是m次多项式。这些Pad'e逼近非常接近原点,并且在第一步中,输入矩阵G将按2的幂进行缩放,因此k2-sGk足够小,可以保证精确的近似值rk,m(2-sG)≈ exp(2-sG)。第二步包括计算有理近似rk,m(2-最后,在第三步中,通过重复分析结果,即exp(G),获得exp(G)的近似值≈ rk,m(2-sG)s.标度参数s和近似度的不同选择,以及不同特征的产量方法。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2022-5-31 05:14:59
这些参数的选择对于近似质量和计算效率至关重要,请参见【8,第10章】。在下文中,我们描述了允许使用缩放和平方对块三角形矩阵ix(1)的矩阵分量进行增量评估的技术。这些技术可用于实际基础缩放和平方方法的任何选择,通过参数s、k和m.2.2工具定义,用于指数增量计算在解释算法之前,我们首先介绍一些贯穿始终的符号。矩阵Gnfrom(1)可以写成GN=G0,0···G0,n-1G0,n。。。。。。。。。Gn公司-1,n-1Gn-1,nGn,n=:Gn公司-1gn0 Gn,n(3) 其中Gn-1.∈ Rdn公司-1×dn-1,Gn,n∈ Rbn×bn,因此gn∈ Rdn公司-10亿。设s为标度参数,r=pq为有理函数,在近似中使用(对于隐式,我们通常会忽略指数k和m)。我们用▄Gn:=2表示缩放材料ix-我们按照(3)对其进行分区。该算法的出发点在于使用缩放和平方方法计算指数exp(G)=exp(G0,0)的Pad'e近似值。然后,通过重用先前获得的每一步量,递增地计算matr ix指数(2)的序列。所以更一般地说,假设exp(Gn-1) 已通过使用平方和平方法进行近似。获得exp(Gn)的Pad'e近似值的三个主要计算步骤是(i)计算多项式p(▄Gn),q(▄Gn),(ii)计算p(▄Gn)-1q(~Gn)和(iii)重复对其进行方形化。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群