A periodic map for linear barycentric rational trigonometric interpolation

被引:11
|
作者
Berrut, Jean-Paul [1 ]
Elefante, Giacomo [1 ]
机构
[1] Univ Fribourg, Dept Math, CH-1700 Fribourg, Switzerland
关键词
Barycentric rational interpolation; Trigonometric interpolation; Conformal maps; CHEBYSHEV PSEUDOSPECTRAL METHOD; POINT SHIFTS; FORMULAS;
D O I
10.1016/j.amc.2019.124924
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
Consider the set of equidistant nodes in [0, 2 pi), theta(k) :=k . 2 pi/n, k = 0, ..., n-1. For an arbitrary 2 pi-periodic function f(theta), the barycentric formula for the corresponding trigonometric interpolant between the theta(k)'s is T[f](theta) = Sigma(n-1)(k=0)(-1)(k) cst(theta-theta(k)/2) f(theta(k))/Sigma(n-1)(k=0)(-1)(k) cst(theta-theta(k)/2), where cst(.) := ctg(.) if the number of nodes n is even, and cst(.) := csc(.) if n is odd. Baltensperger [3] has shown that the corresponding barycentric rational trigonometric interpolant given by the right-hand side of the above equation for arbitrary nodes introduced in [9] converges exponentially toward f when the nodes are the images of the theta(k)'s under a periodic conformal map. In the present work, we introduce a simple periodic conformal map which accumulates nodes in the neighborhood of an arbitrarily located front, as well as its extension to several fronts. Despite its simplicity, this map allows for a very accurate approximation of smooth periodic functions with steep gradients. (C) 2019 Elsevier Inc. All rights reserved.
引用
收藏
页数:8
相关论文
共 50 条