×

在笛卡尔网格上使用样条函数对一般有界区域进行有效的函数逼近。 (英语) Zbl 1504.65041号

摘要:在科学计算中,有界域上的函数通常在适应几何形状的网格上使用分段多项式近似。我们研究了在边界框上定义的规则但过采样的网格上使用样条函数逼近的问题。这种方法允许使用高阶和高度结构化的样条函数作为分段多项式的基础。该方法类似于傅里叶扩展,在边界框上使用傅里叶级数,从而获得平滑函数的光谱精度。然而,傅里叶扩展近似涉及求解一个高度病态的线性系统,这是一个昂贵的步骤。最新算法的计算复杂度在一维中是(mathcal{O}左(N\log^2\左(N \right)\右)),在二维中是(mathcal{O}左\)在1-D中,在2-D中是(mathcal{O}左(N^{3/2}右),在(D)-D中是更一般的(mathcal{O}\左(N_{3(D-1)/D}\右),带有\(D>1\)。通过对相关线性系统使用直接稀疏QR解算器,我们还发现在实际中计算复杂度可以近似线性。这是以仅实现代数收敛速度为代价的。我们的陈述得到了数值实验的证实,并且可以使用Julia代码。

MSC公司:

65日第15天 函数逼近算法
65D07年 使用样条曲线进行数值计算
65层20 超定系统伪逆的数值解
65T40型 三角逼近和插值的数值方法
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] 阿德科克,B。;Huybrechs,D.,《框架与数值近似》,SIAM Rev.,61,3,443-473(2019)·Zbl 1421.42015年 ·doi:10.1137/17M1114697
[2] Adcock,B.,Huybrechs,D.:框架和数值近似II:广义采样。J.傅里叶分析。申请。,26(87) (2020) ·Zbl 1462.42048号
[3] 阿德科克,B。;Huybrechs,D。;Martín-Vaquero,J.,《关于傅里叶扩张的数值稳定性》,Found。计算。数学。,14, 4, 635-687 (2014) ·Zbl 1298.65198号 ·doi:10.1007/s10208-013-9158-8
[4] Adcock,B.,Shadrin,A.:关于通过多项式框架从等距样本快速稳定逼近分析函数的可能性。技术报告arXiv:2110.03755(2021)
[5] 阿尔德鲁比,A。;Unser,M。;Eden,M.,基数样条滤波器:理想正弦插值的稳定性和收敛性,信号处理。,28, 2, 127-138 (1992) ·兹比尔0755.41007 ·doi:10.1016/0165-1684(92)90030-Z
[6] Averbuch,A.Z.,Neitaanmaki,P.:样条和样条小波方法及其在信号和图像处理中的应用:第一卷:周期样条,第一卷Springer(2014)·Zbl 1316.65019号
[7] Averbuch,A.Z.,Neitaanmȧki,P.:V.A Zheludev。样条和样条小波方法及其在信号和图像处理中的应用:第二卷:非周期样条,第二卷Springer(2016)·Zbl 1338.65032号
[8] Boyd,JP,Fourier嵌入式域方法:将定义在不规则区域上的函数扩展为矩形,以便扩展在空间上具有周期性,C,Appl。数学。计算。,161, 2, 591-597 (2005) ·Zbl 1061.65127号
[9] 布鲁诺,OP;韩,Y。;Pohlman,MM,《通过傅里叶延拓分析实现复杂三维曲面的精确高阶表示》,J.Compute。物理。,227, 2, 1094-1125 (2007) ·Zbl 1128.65017号 ·doi:10.1016/j.jcp.2007.08.029
[10] Chen,M.,关于循环线性系统的解,SIAM J.Numer。分析。,24, 3, 668-683 (1987) ·Zbl 0622.65019号 ·数字对象标识代码:10.1137/0724044
[11] Coppe,V.:BSPlineextension.jl v0.1.0https://github.com/framefunVC/BSplineextension.jl (2019)
[12] Coppé,V.,Huybrechs,D.:在笛卡尔网格上使用小波对一般有界域进行有效函数逼近。(arXiv:2004.03537),2020年,提交
[13] Coppé,V.,Huybrechs,D.,Matthysen,R.,Webb,M.:具有已知不完全广义逆的最小二乘系统的AZ算法。SIAM矩阵分析与应用杂志。即将出现(2020年)·Zbl 1461.65058号
[14] Davis,PJ,循环矩阵。《纯粹与应用数学》(1979),纽约:威利出版社,纽约·Zbl 0418.15017号
[15] Davis,T.:多前沿多线程等级稀疏QR分解。摘自:Dagstuhl研讨会论文集。Dagstuhl-Leibniz-Zentrum Für Informatik修道院(2009)
[16] De Boor,C.,样条配置矩阵的全正性,印第安纳大学数学杂志,25,6,541-551(1976)·Zbl 0311.41008号 ·doi:10.1512/iumj.1976.25.25043
[17] Demko,S.,带矩阵的逆和样条投影的局部收敛,SIAM J.Numer。分析。,14, 4, 616-619 (1977) ·Zbl 0367.65024号 ·doi:10.1137/0714041
[18] DeVore,R.A.,Lorentz,G.G.:《建构近似》,第303卷,施普林格科学与商业媒体(1993)·兹伯利0797.41016
[19] Edelman,A。;McCorquodale,P。;托莱多,S.,《未来快速傅里叶变换》,SIAM J.Sci。计算。,20, 3, 1094-1114 (1999) ·Zbl 0926.65145号 ·doi:10.1137/S106482759731666
[20] GMES RDA项目(EU-DEM)。欧洲数字高程模型(eu dem),2014年。从欧洲环境署检索的数据,https://www.eea.europa.eu/data-and-maps/data/copernicus-land-monitoring-service-eu-dem
[21] 指南:花键实用指南。申请。数学。科学。,27 (1978) ·Zbl 0406.41003号
[22] Höllig,K.:B样条有限元方法。SIAM(2003)·Zbl 1020.65085号
[23] Höllig,K。;赖夫,美国。;Wipper,J.,Neumann问题的B样条逼近,印前,2001年2月(2001年)·Zbl 0996.65119号
[24] Huybrechs,D.,关于非周期函数的Fourier扩展,SIAM J.Numer。分析。,47, 6, 4326-4355 (2010) ·Zbl 1209.65153号 ·doi:10.1137/090752456
[25] JuliaMatrices公司。Lowrank约为jl v0.4https://github.com/JuliaMatrices/LowRankApprox.jl (2019)
[26] Lyon,M.,傅里叶延拓的快速算法,SIAM J.Sci。计算。,33, 6, 3241-3260 (2011) ·Zbl 1255.65253号 ·数字对象标识码:10.1137/1082436X
[27] 马提森,R。;Huybrechs,D.,计算任意长度傅里叶扩展的快速算法,SIAM科学计算杂志,38,2,A899-A922(2016)·Zbl 1337.65181号 ·doi:10.1137/15M1030923
[28] 马提森,R。;Huybrechs,D.,《使用傅里叶扩展框架在任意域上的函数逼近》,SIAM J.Numer。分析。,56, 3, 1360-1385 (2018) ·兹比尔1404.33019 ·doi:10.137/17M1134809
[29] Parvizian,J。;杜斯特,A。;Rank,E.,有限单元法,计算力学,41,1,121-133(2007)·Zbl 1162.74506号 ·doi:10.1007/s00466-007-0173-y
[30] 席林格,D。;Ruess,M.,《有限单元法:CAD和基于图像的几何模型的高阶结构分析综述》,工程计算方法档案,22,3391-455(2015)·Zbl 1348.65056号 ·doi:10.1007/s11831-014-9115-y
[31] 勋伯格,IJ,基数插值和样条函数,近似理论杂志,2167-206(1969)·Zbl 0202.34803号 ·doi:10.1016/0021-9045(69)90040-9
[32] Schumaker,L.:样条函数:剑桥数学图书馆基础理论(1981)·Zbl 0449.41004号
[33] Slepian,D.,《Prolate椭球波函数、傅里叶分析和不确定性V:离散情况》,《贝尔系统技术期刊》,571371-1430(1978)·Zbl 0378.33006号 ·doi:10.1002/j.1538-7305.1978.tb02104.x
[34] Unser,M。;阿尔德鲁比,A。;Eden,M.,《用于连续图像表示和插值的快速B样条变换》,IEEE Trans。模式分析。机器。智力。,13, 3, 277-285 (1991) ·数字对象标识代码:10.1109/34.75515
[35] Unser,M。;阿尔德鲁比,A。;Eden,M.,《B样条信号处理:第一部分-理论》,IEEE Trans。信号处理。,41, 2, 821-833 (1993) ·Zbl 0825.94086号 ·数字对象标识代码:10.1109/78.193220
[36] Unser,M。;阿尔德鲁比,A。;Eden,M.,《B样条信号处理:第二部分-高效设计和应用》,IEEE Trans。信号处理。,41, 2, 843-848 (1993) ·Zbl 0825.94087号 ·数字对象标识代码:10.1109/78.193221
[37] Yannakakis,M.,《计算最小填充是NP-complete》,SIAM J.Alg。光盘。方法。,2, 77-79 (1981) ·Zbl 0496.68033号 ·doi:10.1137/0602010
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。