Next Article in Journal
On Weak Limiting Distributions for Random Walks on a Spider
Next Article in Special Issue
Legendre-Gould Hopper-Based Sheffer Polynomials and Operational Methods
Previous Article in Journal
An Algebraic Decision Support Model for Inventory Coordination in the Generalized n-Stage Non-Serial Supply Chain with Fixed and Linear Backorders Costs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rational Approximation on Exponential Meshes

by
Umberto Amato
1 and
Biancamaria Della Vecchia
2,*
1
Istituto di Scienze Applicate e Sistemi Intelligenti, Consiglio Nazionale delle Ricerche, 80131 Napoli, Italy
2
Dipartimento di Matematica ‘Guido Castelnuovo’, Università degli Studi ‘La Sapienza’, 00185 Roma, Italy
*
Author to whom correspondence should be addressed.
Submission received: 7 November 2020 / Revised: 1 December 2020 / Accepted: 2 December 2020 / Published: 4 December 2020
(This article belongs to the Special Issue Special Functions and Polynomials)

Abstract

:
Error estimates of pointwise approximation, that are not possible by polynomials, are obtained by simple rational operators based on exponential-type meshes, improving previous results. Rational curves deduced from such operators are analyzed by Discrete Fourier Transform and a CAGD modeling technique for Shepard-type curves by truncated DFT and the PIA algorithm is developed.

1. Introduction

In [1], simple rational operators r n and q n were constructed, reaching the following pointwise approximation error estimates x [ 0 , 1 ] and f C ( [ 0 , 1 ] )
| f ( x ) r n ( x ) | C ω f ; ϕ ( x ) n , ϕ ( x ) = x ( 1 x ) α , 0 < α < 1 ,
and
| f ( x ) q n ( x ) | C ω f ; ψ ( x ) n , ψ ( x ) = | 2 x 1 | β , 0 < β < 1 .
Here, C denotes a positive constant independent of n and ω ( f ) is the usual modulus of continuity of f .
We remark that if 1 2 < α < 1 in (1) or 0 < β < 1 in (2), such error estimates are not achievable by polynomials, as proved by Gopengauz (see [1,2]).
Moreover, in [1] the authors proved that the exponents α and β in (1) and (2), respectively, cannot be equal to 1 (see remark to Corollary 2.5 proposed by Totik [1]). A crucial key for such estimates is the careful choice of node meshes of algebraic-type.
Zeros of orthogonal polynomials as nodes for rational approximation were examined in [3], but corresponding results are weaker than (1) and (2).
It was an open problem to construct simple rational operators, reaching pointwise approximation error estimates in terms of functions vanishing at any fixed point in [ 0 , 1 ] faster than in (1) and (2), like
χ ( x ) = x ( 1 x ) 1 log x ( 1 x ) λ , 0 < λ ,
or
χ ( x ) = | 2 x 1 | 1 log | 2 x 1 | μ , 0 < μ .
On the other hand, it was an unsolved question to get results of convergence and pointwise estimates of the approximation error by the above rational operators, based on exponential-type meshes. In such cases, the main difficulty is to balance the rational nature of above operators with the exponential mesh behavior.
The present paper aims at giving a first successful answer to the above problems, considering easy rational operators based on exponential-type meshes.
In Section 2, convergence statements and pointwise approximation error estimates, involving functions of type (3) or (4), are given in Theorems 1–4.
Then, in Section 3, we use above operators to construct rational curves by Discrete Fourier Transform (DFT) and design a truncated progressive iterative approximation technique, useful in CAGD modeling.
Recently, parametric Shepard-type curves were studied and properties and algorithms, interesting in CAGD and image processing, were deduced in [4,5]. In Section 3, we study how to transform Shepard-type curves into a different domain using Discrete Fourier Transform through the well-known Fast Fourier Transform algorithm. Theorem 6 presents the inner structure of Shepard-type curves in the form of other Shepard-type curves, known as base Shepard-type curves in the transform form. These base Shepard-type curves have as control points twiddle factors of DFT; hence, the geometry of these curves is determined by the geometry of the corresponding star polygons. The modeling power of such basic Shepard-type curves overcomes the analogous Bézier and B-spline ones, as in original Shepard-type curves case (cf. [4]). Then, in Section 3.1, by truncated DFT and the progressive iterative approximation (PIA in short) algorithm introduced in [4], we deduce a method to construct new Shepard-type curves having good shaping behavior. Indeed, the number of base Fourier functions and the iteration level of PIA algorithm handle modeling the outline of Shepard-type curve, getting as a limit case the Shepard-type interpolating curve.
The proofs of the above results are included in Section 4. They are based on smart choices of nodes meshes, careful estimates for our operators and suitable matrix formulations of DFT and the PIA algorithm. Finally, Section 5 shows numerical experiments, illustrating achieved results.

2. Rational Approximation on Exponential-Type Meshes

For n N , we introduce the nodes matrix X ¯ = x ¯ n , k = x ¯ k , k = 0 , , n [ a , b ] , a , b R . Then, for any function f C ( [ a , b ] ) , we consider the Shepard operator S n as
S n X ¯ ; f ; x = S n ( f ; x ) = k = 0 n f ( x ¯ k ) ( x x ¯ k ) s k = 0 n 1 ( x x ¯ k ) s ,
with x [ a , b ] and s > 2 even. From the definition, we can deduce that S n is a positive, linear operator, preserving constants; S n ( f ) is a rational function having degree ( s n , s n ) , interpolating f at x ¯ k , k = 0 , , n and is stable in Fejér sense, i.e., min a x b | f ( x ) | | S n ( f ; x ) | max a x b | f ( x ) | . If the mesh is equispaced on [ a , b ] , and we assume [ a , b ] = [ 0 , 1 ] for example, then S n is a symmetric operator, i.e.,
k = 0 n f ( x ¯ k ) ( x x ¯ k ) s k = 0 n 1 ( x x ¯ k ) s = k = 0 n f ( x ¯ n k ) ( 1 x x ¯ k ) s k = 0 n 1 ( 1 x x ¯ k ) s .
This operator is widely used in classical approximation theory and problems of scattered data interpolation (see, e.g., [6,7,8,9]).
When the nodes’ mesh is equispaced, direct and converse results for the operator S n exist (e.g., [7,10]). For nodes mesh of algebraic type, pointwise estimates were obtained in [1].
The case of nodes exponentially thicker near one fixed point was an open problem. Now we prove that, by nodes exponentially dense, S n operator achieves results of convergence and estimates of pointwise approximation error, involving functions of type (3) and (4).
To this aim, we assume [ a , b ] = [ 0 , 1 ] and denote by | | f | | the usual supremum norm on [ 0 , 1 ] of f C ( [ 0 , 1 ] ) . Moreover, C denotes a positive constant that can assume different values, even when it appears more times in the same formula. Then let
X = x n , k = x k : = exp n k η + 1 , k = 1 , , n , x 0 = 0 , n N [ 0 , 1 ] , η > 0 ,
be the nodes’ matrix. This mesh is exponentially thick near 0 . Then, we consider operator S n ( X ) defined by (5) based on the nodes matrix X in (7). We have
Theorem 1.
Let f L i p 1 ( [ 0 , 1 ] ) and S n ( X ) be the operator defined by (5)–(7). Then
lim n | | f S n ( X , f ) | | = 0 .
Moreover, x [ 0 , 1 ] , x x k , k = 0 , , n ,
| f ( x ) S n ( X ; f ; x ) | C χ 1 ( x ) n , x = C , χ 1 ( x ) n + χ ¯ 1 ( x ) , x = o ( 1 ) ,
where χ 1 ( x ) = x ( 1 log x ) 1 + 1 / η , χ ¯ 1 ( x ) = x ( 1 log x ) and x = o ( 1 ) means x vanishing, as n .
Remark 1.
Equation (8) allows us to deduce the uniform convergence of S n ( X ; f ) to f, as n , f L i p 1 ( [ 0 , 1 ] ) . Moreover, estimate (9) gives an answer to the problems presented in Introduction. The thickness of the exponential-type mesh near 0 affects estimate of the error (due to the presence of functions χ 1 and χ ¯ 1 at the r.h.s. in (9)).
As remarked in [1], the left-end point 0 is not a special point and we can get analogous results for the right-endpoint 1 case, both endpoints case and any interior point case. For example, let
y k = 1 x k , k = 0 , , n ,
with x k given in (7). Now consider the matrix Y = ( y k , k = 0 , , n , n N ) . This mesh is exponentially finer near 1 . Moreover, for any f C ( [ 0 , 1 ] ) , consider the operator S n ( Y ) based on the matrix Y . We can state
Theorem 2.
Let f L i p 1 ( [ 0 , 1 ] ) and S n ( Y ) be the operator defined by (5) and (10). Then
lim n | | f S n ( Y ; f ) | | = 0 .
Furthermore, x [ 0 , 1 ] , x y k , k = 0 , , n ,
| f ( x ) S n ( Y ; f ; x ) | C χ 2 ( x ) n , 1 x = C , χ 2 ( x ) n + χ ¯ 2 ( x ) , 1 x = o ( 1 ) ,
with χ 2 ( x ) = ( 1 x ) 1 log ( 1 x ) 1 + 1 / η and χ ¯ 2 ( x ) = ( 1 x ) ( 1 log ( 1 x ) ) .
Remark 2.
Obviously, from (11) we deduce the uniform convergence of S n ( Y ; f ) to f, as n , f L i p 1 ( [ 0 , 1 ] ) . Moreover, error estimate is affected by the mesh thickness near 1, because of functions χ 2 ( x ) and χ ¯ 2 at the r.h.s. of (12).
Combining the above results, we can consider the matrix Z = z k , k = 0 , , n , n N , with n even,
z k = 1 2 x n / 2 , k , k = 0 , , n 2 , 1 z n k , k = n 2 + 1 , , n ,
x n / 2 , k as in (7). Then, denote by S n ( Z ) the operator S n based on the matrix Z . When n is odd, we can replace S n ( Z ) by S n + 1 ( Z ) . We have
Theorem 3.
Let f L i p 1 ( [ 0 , 1 ] ) and S n ( Z ) be the operator defined by (5) and (13). Then
lim n | | f S n ( Z ; f ) | | = 0 .
In addition, x [ 0 , 1 ] , x z k , k = 0 , , n ,
| f ( x ) S n ( Z ; f ; x ) | C χ 3 ( x ) n , x ( 1 x ) = C , χ 3 ( x ) n + χ ¯ 3 ( x ) , x ( 1 x ) = o ( 1 ) ,
with χ 3 ( x ) = x ( 1 x ) 1 log ( x ( 1 x ) ) 1 + 1 / η and χ ¯ 3 ( x ) = x ( 1 x ) ( 1 log ( x ( 1 x ) ) ) .
Remark 3.
From (14), we have the uniform convergence of S n ( Z ; f ) to f , as n , f L i p 1 ( [ 0 , 1 ] ) . Moreover, the mesh thickness near 0 and 1 affects estimate of the error (due to functions χ 3 ( x ) and χ ¯ 3 at the r.h.s. of Equation (15)). Hence, Theorem 3 successfully answers the problems posed in the introduction. For example, if x ( 1 x ) exp ( n ) and f L i p 1 ( [ 0 , 1 ] ) , then by (15) | f ( x ) S n ( Z ; f ; x ) | C exp ( n ) n , which is faster than O ( exp ( α n ) n ) coming from (1).
Additionally for n even, let
t k = 1 / 2 1 exp n 2 k η + 1 , k = 0 , , n 2 , 1 t n k , k = n 2 + 1 , , n .
This mesh is thicker near the inner point 1 / 2 . Denote by S n ( T ) the operator S n based on the matrix T = ( t k , k = 0 , , n , n N ) [ 0 , 1 ] . If n is odd, we can replace S n ( T ) by S n + 1 ( T ) .
We prove
Theorem 4.
Let f L i p 1 ( [ 0 , 1 ] ) and S n ( T ) be the operator defined by (5) and (16). Then
lim n | | f S n ( T ; f ) | | = 0 .
Moreover, x [ 0 , 1 ] , x t k , k = 0 , , n ,
| f ( x ) S n ( T ; f ; x ) | C χ 4 ( x ) n , | 2 x 1 | = C , χ 4 ( x ) n + χ ¯ 4 ( x ) , | 2 x 1 | = o ( 1 ) ,
with χ 4 ( x ) = | 2 x 1 | 1 log | 2 x 1 | 1 + 1 / η and χ ¯ 4 ( x ) = | 2 x 1 | 1 log | 2 x 1 | .
Remark 4.
Equation (17) proves the uniform convergence of S n ( T ; f ) to f, as n , f L i p 1 ( [ 0 , 1 ] ) . The error estimate is affected by the mesh thickness near 1 / 2 , because of functions χ 4 ( x ) and χ ¯ 4 at the r.h.s. of (18). Theorem 4 answers to the problems posed in the introduction. For example if | 2 x 1 | exp ( n ) and f L i p 1 ( [ 0 , 1 ] ) , then by (18) | f ( x ) S n ( T ; f ; x ) | C exp ( n ) n , which is better than O exp β n n descending from (2).

3. Shepard-Type Curves and Discrete Fourier Transform

Let us first recall some properties of parametric curves of Shepard-type (see [4]). Let A n ( t ) = [ A n , 0 ( t ) , A n , 1 ( t ) , , A n , n ( t ) ] , where
A n , i ( t ) = 1 ( t t i ) s + λ k = 0 n 1 ( t t i ) s + λ ,
for 0 i n , n N , t [ 0 , 1 ] , t i = i n , i = 0 , 1 , , n , 0 < n s λ C , where C is any fixed positive constant, and s > 2 is even.
Let P = [ P 0 , P 1 , , P n ] T , P i R d , i = 0 , 1 , , n , d 2 , be the control vector. Then the parametric curve of Shepard-type S n [ P , t ] is defined as
S n [ P , t ] = i = 0 n A n , i ( t ) P i = A n ( t ) · P .
In [4], some properties of such curves were studied, that are interesting in CAGD. In particular, S n [ P ] is a rational curve that reproduces points and lies within the convex hull of the control polygon P. Such a curve satisfies the pseudo-local control property; that is, each function A n , j ( t ) , 0 j n , reaches its maximum value close to 1 at t = t j ; therefore the point P j strongly affects the shape of the curve close to t = t j . From (6), Shepard-type curves are symmetric, i.e., S n [ P , t ] = S n [ P ˜ , 1 t ] , with P ˜ = [ P n , P n 1 , , P 0 ] T . The choice 0 < n s λ C makes S n [ P ] a near-interpolating curve of the original control polygon. This fixes the flat spots artifact that affects the original Shepard operator (cf. [4]).
In [4], we introduced and investigated a PIA technique for curves of Shepard-type. The PIA process starts with an initial Shepard-type curve; then through iterations, it adjusts the control points at each iteration, resulting in a sequence of fitting curves. The limiting of the curves at different iterations is the Shepard-type curve interpolating the data. In details, given the control vector P and the basis A n , i ( t ) , i = 0 , , n , defined by (19), the initial curve is generated as
γ 0 ( t ) = i = 0 n A n , i ( t ) P i 0 = S n [ P , t ] ,
with P i 0 = P i , i = 0 , , n . Then the remaining curves of the sequence γ k + 1 ( t ) , for k 0 , are computed as
γ k + 1 ( t ) = i = 0 n P i k + 1 A n , i ( t ) ,
with
P i k + 1 = P i k + Δ i k ,
and Δ i k the adjusting vectors given by
Δ i k = P i γ k ( t i ) , i = 0 , 1 , , n .
The iterative process can be also expressed in matrix form as
Δ 0 k , Δ 1 k , , Δ n k T = ( I B ) Δ 0 k 1 , Δ 1 k 1 , , Δ n k 1 T = ( I B ) k Δ 0 0 , Δ 1 0 , , Δ n 0 T ,
with B the collocation matrix of A n ( t ) basis, i.e.,
B = A n , 0 ( t 0 ) A n , 1 ( t 0 ) A n , n ( t 0 ) A n , 0 ( t 1 ) A n , 1 ( t 1 ) A n , n ( t 1 ) A n , 0 ( t n ) A n , 1 ( t n ) A n , n ( t n ) .
We remark that B is a symmetric, stochastic, diagonally dominant matrix (see [4]).
We say that the γ 0 curve satisfies the PIA property iff lim k γ k ( t i ) = P i , i = 0 , , n . We can state (cf. [4])
Theorem 5.
Curve γ 0 satisfies the PIA property.
Remark 5.
The proof of Theorem 5 is based on the well-known result
M 1 = I + ( I M ) + ( I M ) 2 + ( I M ) 3 + ,
where M is any nonsingular n × n matrix such that ρ ( I M ) < 1 , with ρ ( I M ) being the spectral radius of I M , I the identity matrix, ( I M ) i the i-th power of matrix ( I M ) and M 1 the inverse matrix, i.e., M 1 M = I .
Thanks to the PIA property, we can construct a sequence of control polygons that converge to the control polygon generating a Shepard-type interpolating curve. In addition, we can use k as a shape parameter so to model several outlines; the extreme cases are the curve of Shepard-type and the global interpolating curve of Shepard-type. As remarked in [4], the rate of convergence of such a procedure is faster than the Bézier case.
Then we recall some preliminaries on Fourier analysis, namely the resulting Fourier Transform. The Fourier Transform is the decomposition of a function into components at higher and higher frequencies. Discrete Fourier Transform (DFT) is the discrete counterpart of the Fourier Transform, yielding an estimate starting from a finite sample. DFT maps an ordered set of n + 1 complex number to a different one. Let x ( k ) x 0 , x 1 , x 2 , , x n be a series of n + 1 complex numbers. We assume a periodicity condition that outside the range 0 , n the series is extended n + 1 periodic, i.e., x L = x L + n + 1 , L . We denote DFT of this series x ( L ) . The forward transform is defined as
x ( L ) = 1 n + 1 i = 0 n x ( i ) w L i , L = 0 , 1 , , n ,
with w = exp ( j 2 π n + 1 ) known as twiddle factor and j the solution of j 2 = 1 (complex unity).
In practical applications, DFT is computed by the well-known Fast Fourier Transform algorithm.
Now we are ready to analyze Shepard-type curves by DFT.
Theorem 6.
Let S n [ P ] be the Shepard-type curve defined in (20). We can write
S n [ P , t ] = L = 0 n P ¯ L A ˜ n , L ( t ) ,
where
A ˜ n , L ( t ) = 1 n + 1 i = 0 n A n , i ( t ) w L i , L = 0 , 1 , , n ,
is called the “base Shepard-type curve” and
P ¯ L = 1 n + 1 i = 0 n P i w L i .
In matrix notation, we can write
S n [ P , t ] = A ˜ · P ¯ , A ˜ = A W H , P ¯ = W P ,
where A ˜ = A ˜ n , 0 ( t ) , A ˜ n , 1 ( t ) , , A ˜ n , n ( t ) , P ¯ = P ¯ 0 , P ¯ 1 , , P ¯ n T , A = A n ( t ) and W is the Fourier matrix
W L , i = 1 n + 1 w L i , i , L = 0 , , n ,
with H denoting the Hermitian.
Remark 6.
From Theorem 6, S n [ P , t ] curve can be expressed as weighted average of base Shepard-type curves, by DFT of the original control points.
It can be easily proved that A ˜ behaves as a basis for Shepard-type curve by polygon points in transformed form. Thus, the original Shepard-type curves are generated by these base Shepard-type curves.
The control points that are obtained for them are w L , w 2 L , w 3 L , , w ( n + 1 ) L , so the corresponding polygon is a regular or star one. The relative base curves of Shepard-type lie in the convex-hull of these regular polygons or star polygons, (cf. [11]).
Since w ( ( n + 1 ) L ) i = w L i , obtained Shepard-type curves A ˜ n , L ( t ) and A ˜ n , n + 1 L ( t ) are mirror images of each other about x axis.
Analogous Fourier analysis for Bézier, B-splines and rational B-splines curves was presented in [12,13], but at lower modeling power than our curves (see Section 5).

3.1. Shepard-Type Curves via DFT and PIA

Consider the global interpolating Shepard-type curve (cf. [4])
G n [ P , t ] = k = 0 n A n , k ( t ) Q k ,
with Q k R d , d 2 , such that
G n [ P , t i ] = P i , i = 0 , , n ,
i.e., by (21)
B Q = P , Q = [ Q 0 , Q 1 , , Q n ] T ,
or equivalently
Q = B 1 P = V P , V = B 1 .
From (25), (24) in Theorem 6 and (26), we can write
G n [ P , t ] = A ˜ · Q ¯ = A W H W Q = A W H W V P .
Now let us consider only the first k, k n , Fourier basis functions in W, or in matrix notation W ( k ) R k × n + 1 , with index ( k ) denoting the truncation procedure. This choice is made in analogy with truncation occurring in some statistical contexts involving FT. Then the truncated interpolating curve of Shepard-type is introduced:
G n ( k ) [ P , t ] = A n ( t ) W ( k ) H W ( k ) V P .
Obviously, if k = n , W ( n ) H W ( n ) = I and we get back (25). Hence, varying 0 k n in (27), we get different curves approaching the interpolating Shepard-type curve.
Representation (27) suggests to use the above PIA technique (see Theorem 5, (22) and consequent remark) to construct a method giving a pencil of curves of Shepard-type, modeling the original data points.
We sketch the procedure. In (27), we replace V by
V r = k = 0 r ( I B ) k ,
with r being the iteration level. From [4], lim r V r = V and convergence rate is fast, so a few iterations are enough to go close to V.
Then playing on two shape parameters, the number k of basis Fourier functions in (27) and the iteration level of the PIA algorithm—the index r in (28), the designer can get intermediate contours not reachable by original PIA format (see Section 5).

4. Proofs

Proof of Theorem 1. 
Due to the interpolatory behavior of S n at x k , k = 0 , , n , we can assume that x x k , k = 0 , , n . Let x j , 0 j n , be the closest knot to x, x j < x . The case when x j + 1 is the closest knot to x, x < x j + 1 , can be treated similarly. First, we prove that for 0 j n 1
| x x j | C χ 1 ( x ) n ,
with χ 1 ( x ) = x ( 1 log x ) 1 + 1 / η . Indeed, letting g ( δ ) = exp δ η + 1 , we put x = g ( ϑ ) and x j = g j n , with j n < ϑ < j + 1 n . Hence, by Lagrange’s theorem, for ξ j n , ϑ
| x x j | = | g ( ϑ ) g j n = | g ( ξ ) | ϑ j n C n exp ξ η + 1 ξ η 1 C n exp θ η + 1 θ η 1 C n x ( 1 log x ) 1 + 1 / η = C n χ 1 ( x ) ,
that is (29).
Now,
| f ( x ) S n ( X ; f ; x ) | k = 0 n | f ( x ) f ( x k ) | ( x x k ) s k = 0 n 1 ( x x k ) s | f ( x ) f ( x j ) | ( x x j ) s k = 0 n 1 ( x x k ) s + | f ( x ) f ( x j + 1 ) | ( x x j + 1 ) s k = 0 n 1 ( x x k ) s + k j , j + 1 | f ( x ) f ( x k ) | ( x x k ) s k = 0 n 1 ( x x k ) s : = Σ 1 + Σ 2 + Σ 3 .
Since
1 k = 0 n 1 ( x x k ) s ( x x j ) s ,
it follows from (30)
Σ 1 ( x x j ) s | f ( x ) f ( x j ) | ( x x j ) s ω f ; χ 1 ( x ) n .
By | x x j + 1 | > | x x j | and a well known property of modulus of continuity, from (30) and (29)
Σ 2 ( x x j ) s | f ( x ) f ( x j + 1 ) | ( x x j + 1 ) s ( x x j ) s | x x j + 1 | s 1 ω ( f ; | x x j + 1 ) | | x x j + 1 | C ( x x j ) s | x x j + 1 | s 1 ω ( f ; | x x j | ) | x x j | C ( x x j ) s | x x j | s 1 ω ( f ; | x x j | ) | x x j | C ω ( f ; | x x j | ) C ω f ; χ 1 ( x ) n .
Moreover, from (30) it follows
Σ 3 ( x x j ) s k = 0 j 1 + k = j + 2 n ω f ; | x x k | ( x x k ) s : = Σ 4 + Σ 5 .
Now in Σ 5
x k x C χ 1 ( x ) k j n .
Hence, by usual calculations (see e.g., [1])
Σ 5 C ω f ; χ 1 ( x ) n χ 1 ( x ) s n s k = j + 2 n ( k j ) n s ( k j ) s χ 1 ( x ) s C ω f ; χ 1 ( x ) n .
Now we estimate Σ 4 . First let x = C . If x k = C , then
| x x k | C j k n ,
so
| x x j | s ω f ; | x x k | | x x k | s C n s n ω f ; 1 n n s 1 ( j k ) s 1 C ω f ; 1 n ( j k ) s 1 .
If x k = o ( 1 ) , then | x x k | > C , consequently
| x x j | s ω f ; | x x k | | x x k | s C n s ω f ; 1 n n .
Collecting the above estimations, from (31) we get
Σ 4 C ω f ; 1 n C ω f ; χ 1 ( x ) n .
Now we assume x = o ( 1 ) , which implies j = o ( n ) . Then by | x x j | | x x k | , k = 0 , , j 1 , and f L i p 1 ( [ 0 , 1 ] ) , it follows that
Σ 4 | x x j | s k = 0 j 1 1 ( x x k ) s 1 C j | x x j | C χ 1 ( x ) ( 1 log x ) 1 / η .
Collecting above estimations, the theorem is proved. ☐
Proof of Theorem 2. 
Analogously to Theorem 1, we have
| x y j | | x y j + 1 | C χ 2 ( x ) n ,
with χ 2 ( x ) = ( 1 x ) ( 1 log ( 1 x ) ) 1 + 1 / η .
Then, as in the proof of Theorem 1,
| x y k | C χ 2 ( x ) ( j k ) n , k = 0 , , j 1 .
Hence,
( x y j ) s k = 0 j 1 ω f ; | x y k | ( x y k ) s C ω f ; χ 2 ( x ) n .
Then we assume 1 x = o ( 1 ) . By | x y j | | x y k | , k = j + 2 , , n and f L i p 1 ( [ 0 , 1 ] ) ,
( x y j ) s k = j + 2 n ω f ; | x y k | ( x y k ) s C ( n j 2 ) | x y j | C χ 2 ( x ) ( 1 log ( 1 x ) ) 1 / η .
Other cases can be handled as in Theorem 1 and we deduce the statement. ☐
Proof of Theorem 3. 
Proceeding as in Theorems 1 and 2, we obtain
| x z j | C χ 3 ( x ) n .
Working as before (cf. [1]), the assertion is proved. ☐
Proof of Theorem 4. 
Working as in the proof of Theorems 1 and 2, we get
| x t j | C χ 4 ( x ) n .
Following the above demonstrations (cf. [1]), we derive the statement. ☐
Proof of Theorem 6. 
Let S n [ P , t ] be the Shepard-type curve defined by (20). The control points P i , i = 0 , 1 , , n , can be considered as points on complex plane. Therefore, each P i , i = 0 , 1 , , n , can be equivalently seen as the inverse transform of points P ¯ i , with P ¯ i being the Fourier Transform of points P i , i.e.,
P i = 1 n + 1 L = 0 n P ¯ L w L i ,
with w = exp j 2 π n + 1 the twiddle factor. Thus, Equation (20) can be written as
S n [ P , t ] = 1 n + 1 i = 0 n A n , i ( t ) L = 0 n P ¯ L w L i .
As P ¯ L ’s are independent of i, we get
S n [ P , t ] = 1 n + 1 L = 0 n P ¯ L i = 0 n A n , i ( t ) w L i
and by (23) the assertion follows. ☐

5. Numerical Experiments

In this section, we show some examples of Shepard-type curves in the transformed domain.

5.1. Example 1

Let n = 3 . Then the control points of basic Shepard-type curves are 1 , w 1 , w 1 2 , w 1 3 , with w 1 = exp ( 2 π 4 ) , or 1 , i , 1 , i . In this case, A ˜ n , 0 ( t ) degenerates to the single point ( 1 , 0 ) in the plane, while the basic curve A ˜ n , 2 ( t ) describes the line segment from ( 1 , 0 ) to ( 1 , 0 ) and back (cf. [12] for the corresponding degenerated Bézier case).

5.2. Example 2

Here, we present basic Shepard-type curves defined by (20) for s = 4 corresponding to cases
(a)
n = 9 , j = 7 , λ = 10 5 (Figure 1);
(b)
n = 11 , j = 3 , λ = 10 5 (Figure 2, left);
(c)
n = 11 , j = 10 , λ = 10 5 (Figure 2, right);
(d)
n = 13 , j = 3 , λ = 10 6 (Figure 3, left);
(e)
n = 13 , j = 5 , λ = 10 6 (Figure 3, right).
From such figures, the modeling outperformance of basic Shepard-type curves with respect to Bézier case is clear (cf. [12]).

5.3. Example 3

Let n = 6 , s = 4 and λ = 5 × 10 5 . Figure 4 shows the corresponding basic Shepard-type curves satisfying properties mentioned in Section 3. Comparing Figure 4 with Figure 1 in [13] for the corresponding base B-spline curves, we can note that basic Shepard-type curves show superior shaping capabilities.

5.4. Example 4

Consider Archimedes spiral given by
( x ( t ) , y ( t ) ) = ( t cos 6 π t , t sin 6 π t ) , t [ 0 , 1 ] .
We extract from the spiral a sample of 11 ( n = 10 ) control points as
( x ( s i ) , y ( s i ) ) , s i = i n , i = 0 , , n .
They form the control polygon in Figure 5.
We start from these control points and fit the spiral by a sequence of three curves that are defined by the truncated iterative procedure introduced in Section 3.1 with k = 9 , 10, 11, respectively, in (27), r = 1 , s = 4 and λ = 10 5 (see Figure 5, left). From Figure 5 (left), we can see the modeling behavior of our procedure; indeed for k = 11 and r = 1 , we find back the curve obtained by one iteration of the PIA algorithm, while the choices k = 9, 10 and r = 1 give intermediate curves approximating the given data points. Analogous curves for the choice r = 2 , s = 4 and λ = 10 5 are shown in Figure 5 (right). In this way, the designer has different choices for outlines modeling the spiral.

5.5. Example 5

We sample the spiral of Example 4 with 19 ( n = 18 ) control points, given by
( x ( s i ) , y ( s i ) ) , s i = i n , i = 0 , , n .
Then we start from these control points and fit the spiral by a sequence of three curves defined by our truncated iterative procedure with s = 4 , λ = 3 × 10 6 , k = 17 , 18 , 19 in (27) and r = 1 (Figure 6, left) and r = 3 (Figure 6, right), respectively. Comparison of Figure 6 (left) with Figure 6 (right) shows the influence of shape parameters in slight changes of our curves, giving new contours fitting the spiral.

5.6. Example 6

Consider the Golden spiral curve in a sequence of 15 ( n = 14 ) points defined by
( x ( t i ) , y ( t i ) ) = a φ 2 τ i / π cos τ i , a φ 2 τ i / π sin τ i , φ = 1 2 , τ i = 5 + 10 i n , i = 0 , , n ,
forming the control polygon in Figure 7. Actually, since the curve lacks periodicity, we virtually extend the control points at the right and/or left with mirror points. As is well known, the enlarged function is now even, then the Fourier Transform reverts to the Cosine Fourier Transform [14], with the favorable consequence that the coefficients and the inverse transform keep real. From the computational point of view this is obtained by simply computing the Discrete Cosine Fourier Transform of the original 15 points.
We fit these points by 6 sequences of curves defined by the above procedure with s = 4 , λ = 10 5 , k = 10 , 11, …, 15 and r = 1 (see Figure 7), from which one can see several contours approximating Golden Spiral.

Author Contributions

Conceptualization, B.D.V.; methodology, U.A. and B.D.V.; software, U.A.; validation, U.A.; investigation, U.A. and B.D.V. ; writing—original draft preparation, B.D.V.; writing—review and editing, U.A. and B.D.V.; visualization, U.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors are grateful to the anonymous reviewers for their interesting comments that improved the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
DFTDiscrete Fourier Transform
CAGDComputer Aided Geometric Design
FTFourier Transform
PIAProgressive Iterative Approximation

References

  1. Della Vecchia, B.; Mastroianni, G. Pointwise simultaneous approximation by rational operators. J. Approx. Theory 1991, 65, 140–150. [Google Scholar] [CrossRef] [Green Version]
  2. Gopengauz, I. A theorem of A.F. Timan on the approximation of functions by polynomials on a finite segment. Math. Notes 1967, 1, 110–116. [Google Scholar] [CrossRef]
  3. Criscuolo, G.; Mastroianni, G. Estimates of the Shepard interpolatory procedure. Acta Math. Hungar. 1993, 6, 79–91. [Google Scholar] [CrossRef]
  4. Amato, U.; Della Vecchia, B. Modelling by Shepard-type curves and surfaces. J. Comp. Anal. Applic. 2016, 20, 611–634. [Google Scholar]
  5. Amato, U.; Della Vecchia, B. On Shepard-Gupta-type operators. J. Ineq. Appl. 2018, 232. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Allasia, G. A class of interpolatory positive linear operators: Theoretical and computational aspects. In Approximation Theory, Wavelets and Applications; NATO ASI Series C; Springer: Dordrecht, The Netherlands, 1995; Volume 454, pp. 1–36. [Google Scholar]
  7. Szabados, J. On a problem of R. DeVore. Acta Math. Acad. Sci. Hungar. 1976, 27, 219–223. [Google Scholar] [CrossRef]
  8. Tian, Y.; Mou, J.; Shi, T.; Xia, Q. Shepard Interpolation based on geodesic distance for optimization of fiber reinforced composite structures with non-convex shape. Appl. Compos. Mater. 2019, 26, 575–590. [Google Scholar] [CrossRef]
  9. Wang, H.; Bettens, R. Modelling potential energy surfaces for small clusters using Shepard interpolation with Gaussian-form nodal functions. Phys. Chem. Chem. Phys. 2019, 8, 4513–4522. [Google Scholar] [CrossRef] [PubMed]
  10. Somorjai, G. On a saturation problem. Acta Math. Acad. Sci. Hungar. 1978, 32, 377–381. [Google Scholar] [CrossRef]
  11. Coxeter, H.S.M. Introduction to Geometry, 2nd ed.; Wiley: New York, NY, USA, 1969. [Google Scholar]
  12. Barry, P. The Fourier Analysis of Bézier Curves. J. Visual Mathem. 2003, 5. Available online: http://www.mi.sanu.ac.rs/vismath/barry/index.html (accessed on 25 November 2020).
  13. Ganguly, A.; Arondekar, P. Analysis of B-spline curve using Discrete Fourier Transform. Math. Comp. Appl. 2010, 15, 127–136. [Google Scholar] [CrossRef] [Green Version]
  14. Rao, K.R.; Yip, P. Discrete Cosine Transform: Algorithms, Advantages, Applications; Academic Press: Boston, MA, USA, 1990. [Google Scholar]
Figure 1. Basic Shepard curve for n = 9 , λ = 10 5 and j = 7 .
Figure 1. Basic Shepard curve for n = 9 , λ = 10 5 and j = 7 .
Symmetry 12 01999 g001
Figure 2. Basic Shepard curve for n = 11 , λ = 10 5 and j = 3 (left), j = 10 (right).
Figure 2. Basic Shepard curve for n = 11 , λ = 10 5 and j = 3 (left), j = 10 (right).
Symmetry 12 01999 g002
Figure 3. Basic Shepard curve for n = 13 , λ = 10 6 and j = 3 (left), j = 5 (right).
Figure 3. Basic Shepard curve for n = 13 , λ = 10 6 and j = 3 (left), j = 5 (right).
Symmetry 12 01999 g003
Figure 4. Basic Shepard curves for n + 1 = 7 , λ = 5 × 10 5 and j = 1 , , 6 .
Figure 4. Basic Shepard curves for n + 1 = 7 , λ = 5 × 10 5 and j = 1 , , 6 .
Symmetry 12 01999 g004
Figure 5. Archimedes spiral, n = 10 , λ = 10 5 , k = 9 , 10, 11 and r = 1 (left), r = 2 (right).
Figure 5. Archimedes spiral, n = 10 , λ = 10 5 , k = 9 , 10, 11 and r = 1 (left), r = 2 (right).
Symmetry 12 01999 g005
Figure 6. Archimedes spiral, n = 18 , λ = 3 × 10 6 , k = 17, 18, 19 and r = 1 (left), r = 3 (right).
Figure 6. Archimedes spiral, n = 18 , λ = 3 × 10 6 , k = 17, 18, 19 and r = 1 (left), r = 3 (right).
Symmetry 12 01999 g006
Figure 7. Golden spiral n = 14 , λ = 10 5 , r = 1 and k = 10, 11, …, 15.
Figure 7. Golden spiral n = 14 , λ = 10 5 , r = 1 and k = 10, 11, …, 15.
Symmetry 12 01999 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Amato, U.; Della Vecchia, B. Rational Approximation on Exponential Meshes. Symmetry 2020, 12, 1999. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12121999

AMA Style

Amato U, Della Vecchia B. Rational Approximation on Exponential Meshes. Symmetry. 2020; 12(12):1999. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12121999

Chicago/Turabian Style

Amato, Umberto, and Biancamaria Della Vecchia. 2020. "Rational Approximation on Exponential Meshes" Symmetry 12, no. 12: 1999. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12121999

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop