×

基于重心有理插值的最佳有理逼近算法。 (英文) Zbl 1473.65015号

摘要:我们提出了一种新的算法,用于计算零缺陷情况下实标量函数的最佳一致有理逼近。该方法被称为BRASIL(通过连续间隔长度调整获得的最佳有理逼近),其依据是对函数(f)的最佳有理有理逼近必须在一定数量的插值节点((x_j)处进行插值。此外,每个区间((x_{j-1},x_j))的局部最大误差序列必须相等振荡。该算法以均衡局部误差为目标,迭代地重新调整区间长度。使用重心有理公式以稳定的方式计算所需的有理插值。BRASIL算法可以被视为插值节点的定点迭代,并且线性收敛。我们证明了适当设计的重缩放和重启安德森加速(RAA)方法显著提高了其收敛速度。新算法具有良好的数值稳定性,仅使用标准IEEE双精度算法即可在几秒钟内计算出许多函数的高度最佳有理逼近。提供了一个免费的Python开源实现。我们通过与文献中的结果进行比较来验证算法。我们还证明了它在当前最先进的方法,即实现Remez算法重心变体的Chebfun包中的minimax函数失败的情况下快速收敛。

MSC公司:

65日第15天 函数逼近算法
41A20型 有理函数逼近
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 阿基瑟,N.I.:近似理论。多佛高等数学书籍。多佛出版社,ISBN 9780486671291(1992)
[2] Anderson,DG,非线性积分方程的迭代程序,J.ACM(JACM),12,4,547-560(1965)·Zbl 0149.11503号 ·doi:10.1145/321296.321305
[3] 贝鲁特,J-P;Trefethen,LN,重心拉格朗日插值,SIAM Rev.,46,3,501-517(2004)·Zbl 1061.65006号 ·doi:10.1137/s0036144502417715
[4] Berrut,J.-P.,Baltensperger,R.,Mittelmann,H.D.:重心有理插值的最新发展。在:构造逼近的趋势和应用。doi:10.1007/3-7643-7356-33,第27-51页。Birkhäuser,巴塞尔(2005年)·Zbl 1077.65009号
[5] Braess,D.,非线性近似理论(1986),柏林:Springer,柏林·Zbl 0656.41001号 ·doi:10.1007/978-3-642-61609-9
[6] Brezinski,C。;Redivo-Zaglia,M.,Padé-型有理插值和重心插值,数值。数学。,125, 1, 89-113 (2013) ·Zbl 1276.65005号 ·doi:10.1007/s00211-013-0535-7
[7] Brezinski,C。;Redivo Zaglia,M.,padépadé-type的新表示法和部分padé的近似,J.Compute。申请。数学。,284, 69-77 (2015) ·Zbl 1323.30044号 ·doi:10.1016/j.cam.2014.07.007
[8] Brezinski,C。;Redivo Zaglia,M。;Saad,Y.,Shanks序列变换和Anderson加速,SIAM Rev.,60,3,646-669(2018)·兹比尔1395.65001 ·数字对象标识代码:10.1137/17m1120725
[9] Carpenter,A.J.,Ruttan,A.,Varga,R.S.:有理逼近理论中1/9猜想的扩展数值计算。In:有理逼近和插值。doi:10.1007/bfb0072427,第383-411页。柏林施普林格(1984)·Zbl 0553.41022号
[10] 切尼,EW;Loeb,HL,有理逼近的两种新算法,Numer。数学。,3, 1, 72-75 (1961) ·Zbl 0114.32804号 ·doi:10.1007/bf01386002
[11] 德梅尔,J。;顾先生。;艾森斯塔特,S。;斯拉普尼查尔,I。;韦塞利奇,K。;Drmać,Z.,计算高相对精度的奇异值分解,线性代数应用。,299, 1-3, 21-80 (1999) ·Zbl 0952.65032号 ·doi:10.1016/s0024-3795(99)00134-2
[12] Driscoll,T.A.,Hale,N.,Trefethen,法律公告:Chebfun Guide。Pafnoty出版物。http://www.chebfun.org/docs/guide/ (2014)
[13] 方,H。;Saad,Y.,非线性加速度的两类多斜面方法,Numer。线性代数应用。,16, 3, 197-221 (2009) ·Zbl 1224.65134号 ·doi:10.1002/nla.617
[14] 菲利普,S-I;Y.Nakatsukasa。;特里芬,LN;Beckermann,B.,通过自适应重心表示的有理极小极大近似,SIAM J.Sci。计算。,40、4、A2427-A2455(2018)·Zbl 1396.41011号 ·数字对象标识代码:10.1137/17m1132409
[15] Golub,G.H.,Van Loan,C.F.:矩阵计算。约翰·霍普金斯大学出版社,第4版。ISBN 9781421408590(2012)·Zbl 0559.65011号
[16] Harizanov,S.,Lazarov,R.,Margenov,S..,Marinov,P.:最佳一致有理逼近:用于求解涉及椭圆算子分数幂的方程。arXiv:1910.13865(2019)·Zbl 1429.65064号
[17] Harizanov,S.,Lazarov,R.,Margenov,S..,Marinov,P.,Pasciak,J.:基于最佳一致有理逼近的谱分数阶椭圆方程数值方法分析。J.计算。物理学。doi:10.1016/j.jcp.2020.109285。在线提供(2020)·Zbl 1429.65064号
[18] Higham,NJ,QR因式分解,具有完整的枢轴和SVD的精确计算,线性代数应用。,309、1-3、153-174(2000年)·Zbl 0953.65025号 ·doi:10.1016/s0024-3795(99)00230-x
[19] 新泽西州海姆,《数值算法的准确性和稳定性》(2002),费城:工业和应用数学学会,费城·Zbl 1011.65010号 ·数字对象标识代码:10.1137/1.9780898718027
[20] 海厄姆,新泽西州,重心拉格朗日插值的数值稳定性,IMA J.Numer。分析。,24, 4, 547-556 (2004) ·Zbl 1067.65016号 ·doi:10.1093/imanum/24.547
[21] 新泽西州海姆;Strabic̀,N.,Anderson交替投影加速法,用于计算最近相关矩阵,Numer。算法,72,4,1021-1042(2015)·Zbl 1347.65074号 ·doi:10.1007/s11075-015-0078-3
[22] Hofreither,C.,分数扩散的一些数值方法的统一观点,计算数学应用,80,2,332-350(2020)·Zbl 1446.65153号 ·doi:10.1016/j.camwa.2019.07.025
[23] 爱奥尼亚,A.C.:拉格朗日有理插值及其在大规模动力系统近似中的应用。莱斯大学博士论文,休斯顿,TY(2013)
[24] Knockaert,L.,重心有理插值的一种简单而精确的算法,IEEE信号处理。莱特。,15, 154-157 (2008) ·doi:10.1109/lsp.2007.913583
[25] 压力,WH;Teukolsky,SA;韦特林,WT;英国石油公司弗兰纳里,《C中的数字配方》(1992),剑桥:剑桥大学出版社,剑桥·Zbl 0778.65003号
[26] 施耐德,C。;沃纳,W.,有理插值的一些新方面,数学。计算。,47, 175, 285-285 (1986) ·Zbl 0612.65005号 ·doi:10.1090/s0025-5718-1986-0842136-8
[27] Stahl,HR,[0,1]上xα的最佳一致有理逼近,数学学报,190,2241-306(2003)·Zbl 1077.41012号 ·doi:10.1007/bf02392691
[28] 托斯,A。;Kelley,CT,Anderson加速度收敛分析,SIAM J.Numer。分析。,53, 2, 805-819 (2015) ·Zbl 1312.65083号 ·doi:10.1137/130919398
[29] Trefethen,法律公告:近似理论和近似实践。应用数学的其他标题。SIAM,2013年。印尼盾9781611972405·Zbl 1264.41001号
[30] 瓦尔加,RS;Carpenter,AJ,关于[0,1]上xα最佳一致有理逼近的一些数值结果,Numer。算法,2,171-185(1992)·兹比尔0763.41025 ·doi:10.1007/bf02145384
[31] 瓦尔加,RS;Ruttan,A。;Carpenter,AD,[-1,1]上|x|最佳一致有理逼近的数值结果,苏联斯博尼克数学,74,2,271-290(1993)·Zbl 0774.65008号 ·doi:10.1070/SM1993v074n02ABEH003347
[32] Walker,HF;Ni,P.,Anderson定点迭代加速度,SIAM J.Numer。分析。,49, 4, 1715-1735 (2011) ·兹比尔1254.65067 ·数字对象标识码:10.1137/10078356x
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。