DIMENSION AND DIFFERENTIABILIIY OF A CLASS OF FRACTAL INTERPOLATION FUNCTIONS

In this paper, a new iterated function system consisting of non-linear affinemaps is constructed. We investigate the fractal interpolation functions generated bysuch a system and get its differentiability, its box dimension, its packing dimension,and a lower bound of its Hansdorff dimension.

HOLDER PROPERTY OF FRACTAL INTERPOLATION FUNCTION

The purpose of this paper is to prove a Holder property about the fractal interpolation function L(x), ω(L,δ)=O(δ~α), and an approximate estimate |f-L|≤2{α(h)+||f||/1-h^(2-D)·h^(2-D)}, where D is a fractal dimension of L(x).

On cubic Hermite coalescence hidden variable fractal interpolation functions

Hermite interpolation is a very important tool in approximation theory and nu- merical analysis, and provides a popular method for modeling in the area of computer aided geometric design. However, the classical Hermite interpolant is unique for a prescribed data set, and hence lacks freedom for the choice of an interpolating curve, which is a crucial requirement in design environment. Even though there is a rather well developed fractal theory for Hermite interpolation that offers a large flexibility in the choice of interpolants, it also has the short- coming that the functions that can be well approximated are highly restricted to the class of self-affine functions. The primary objective of this paper is to suggest a gl-cubic Hermite in- terpolation scheme using a fractal methodology, namely, the coalescence hidden variable fractal interpolation, which works equally well for the approximation of a self-affine and non-self-affine data generating functions. The uniform error bound for the proposed fractal interpolant is established to demonstrate that the convergence properties are similar to that of the classical Hermite interpolant. For the Hermite interpolation problem, if the derivative values are not actually prescribed at the knots, then we assign these values so that the interpolant gains global G2-continuity. Consequently, the procedure culminates with the construction of cubic spline coalescence hidden variable fractal interpolants. Thus, the present article also provides an al- ternative to the construction of cubic spline coalescence hidden variable fractal interpolation functions through moments proposed by Chand and Kapoor [Fractals, 15（1） （2007）, pp. 41-53].

HAAR EXPANSIONS OF A CLASS OF FRACTAL INTERPOLATION FUNCTIONS AND THEIR LOGICAL DERIVATIVES

In this paper,we study a special class of fractal interpolation functions,and give their Haar-wavelet expansions.On the basis of the expansions,we investigate the H(o|¨)lder smoothness of such functions and their logical derivatives of order α.

Energy and Laplacian of fractal interpolation functions

Abstract. In this paper, we first characterize the finiteness of fractal interpolation functions （FIFs） on post critical finite self-similar sets. Then we study the Laplacian of FIFs with uniform vertical scaling factors on the Sierpinski gasket （SG）. As an application, we prove that the solution of the following Dirichlet problem on SG is a FIF with uniform vertical scaling factor 1/5 ：△=0 on SG ／ {q1, q2, q3}, and u（qi）=ai, i = 1, 2, 3, where qi, i=1, 2, 3, are boundary points of SG.

ON THE CONTINUITY AND DIFFERENTIABILITY OF AKIND OF FRACTAL INTERPOLATION FUNCTION

The sufficient conditions of Holder continuity of two kinds of fractal interpolation functions defined by IFS (Iterated Function System) were obtained. The sufficient and necessary condition for its differentiability was proved. Its derivative was a fractal interpolation function generated by the associated IFS, if it is differentiable.

Graph-Directed Coalescence Hidden Variable Fractal Interpolation Functions

Fractal interpolation function (FIF) is a special type of continuous function which interpolates certain data set and the attractor of the Iterated Function System (IFS) corresponding to a data set is the graph of the FIF. Coalescence Hidden-variable Fractal Interpolation Function (CHFIF) is both self-affine and non self-affine in nature depending on the free variables and constrained free variables for a generalized IFS. In this article, graph directed iterated function system for a finite number of generalized data sets is considered and it is shown that the projection of the attractors on is the graph of the CHFIFs interpolating the corresponding data sets.

Radial Basis Function Interpolation and Galerkin Projection for Direct Trajectory Optimization and Costate Estimation

This work presents a novel approach combining radial basis function(RBF)interpolation with Galerkin projection to efficiently solve general optimal control problems.The goal is to develop a highly flexible solution to optimal control problems,especially nonsmooth problems involving discontinuities,while accounting for trajectory accuracy and computational efficiency simultaneously.The proposed solution,called the RBF-Galerkin method,offers a highly flexible framework for direct transcription by using any interpolant functions from the broad class of global RBFs and any arbitrary discretization points that do not necessarily need to be on a mesh of points.The RBF-Galerkin costate mapping theorem is developed that describes an exact equivalency between the Karush-Kuhn-Tucker(KKT)conditions of the nonlinear programming problem resulted from the RBF-Galerkin method and the discretized form of the first-order necessary conditions of the optimal control problem,if a set of discrete conditions holds.The efficacy of the proposed method along with the accuracy of the RBF-Galerkin costate mapping theorem is confirmed against an analytical solution for a bang-bang optimal control problem.In addition,the proposed approach is compared against both local and global polynomial methods for a robot motion planning problem to verify its accuracy and computational efficiency.

HERMITE—BIRKHOFF INTERPOLATION OF SCATTERED DATA BY RADIAL BASIS FUNCTIONS

For Hermite-Birkhoff interpolation of scattered multidumensional data by radial basis function (?),existence and characterization theorems and a variational principle are proved. Examples include (?)(r)=r^b,Duchon's thin-plate splines,Hardy's multiquadrics,and inverse multiquadrics.

ON THE POINTWISE ESTIMATIONS OF APPROXIMATION OF FUNCTIONS AND THEIR DERIVATIVES BY HERMITE INTERPOLATION

The object of this paper is to establish the pointwise estimations of approximation of functions in C^1 and their derivatives by Hermite interpolation polynomials. The given orders have been proved to be exact in gen- eral.

OPTIMAL ERROR BOUNDS FOR THE CUBIC SPLINE INTERPOLATION OF LOWER SMOOTH FUNCTIONS(1)

In this paper,the kernel of the cubic spline interpolation is given.An optimal error bound for the cu- bic spline interpolation of lower smooth functions is obtained.

ON SIMULTANEOUS APPROXIMATION TO A DIFFERENTIABLE FUNCTION AND ITS DERIVATIVE BY INVERSE PAL-TYPE INTERPOLATION POLYNOMIALS

Let ξn-1<ξn-2 <ξn-2 <… < ξ1 be the zeros of the the (n -1)-th Legendre polynomial Pn-1(x) and - 1 = xn < xn-1 <… < x1 = 1 the zeros of the polynomial W n(x) =- n(n - 1) Pn-1(t)dt = (1 -x2)P'n-1(x). By the theory of the inverse Pal-Type interpolation, for a function f(x) ∈ C[-1 1], there exists a unique polynomial Rn(x) of degree 2n - 2 (if n is even) satisfying conditions Rn(f,ξk) = f(∈ek)(1≤ k≤ n - 1) ;R'n(f,xk) = f'(xk)(1≤ k≤ n). This paper discusses the simultaneous approximation to a differentiable function f by inverse Pal-Type interpolation polynomial {Rn(f,x)} (n is even) and the main result of this paper is that if f ∈ C'[1,1], r≥2, n≥ + 2> and n is even thenholds uniformly for all x ∈ [- 1,1], where h(x) = 1 +

Spatio-temporal variability of surface chlorophyll a in the Yellow Sea and the East China Sea based on reconstructions of satellite data of 2001-2020

Chlorophyll-a(Chl-a)concentration is a primary indicator for marine environmental monitoring.The spatio-temporal variations of sea surface Chl-a concentration in the Yellow Sea(YS)and the East China Sea(ECS)in 2001-2020 were investigated by reconstructing the MODIS Level 3 products with the data interpolation empirical orthogonal function(DINEOF)method.The reconstructed results by interpolating the combined MODIS daily+8-day datasets were found better than those merely by interpolating daily or 8-day data.Chl-a concentration in the YS and the ECS reached its maximum in spring,with blooms occurring,decreased in summer and autumn,and increased in late autumn and early winter.By performing empirical orthogonal function(EOF)decomposition of the reconstructed data fields and correlation analysis with several potential environmental factors,we found that the sea surface temperature(SST)plays a significant role in the seasonal variation of Chl a,especially during spring and summer.The increase of SST in spring and the upper-layer nutrients mixed up during the last winter might favor the occurrence of spring blooms.The high sea surface temperature(SST)throughout the summer would strengthen the vertical stratification and prevent nutrients supply from deep water,resulting in low surface Chl-a concentrations.The sea surface Chl-a concentration in the YS was found decreased significantly from 2012 to 2020,which was possibly related to the Pacific Decadal Oscillation(PDO).

A meshless model for transient heat conduction analyses of 3D axisymmetric functionally graded solids

A meshless numerical model is developed for analyzing transient heat conductions in three-dimensional （3D） axisymmetric continuously nonhomogeneous functionally graded materials （FGMs）. Axial symmetry of geometry and boundary conditions reduces the original 3D initial-boundary value problem into a two-dimensional （2D） problem. Local weak forms are derived for small polygonal sub-domains which surround nodal points distributed over the cross section. In order to simplify the treatment of the essential boundary conditions, spatial variations of the temperature and heat flux at discrete time instants are interpolated by the natural neighbor interpolation. Moreover, the using of three-node triangular finite element method （FEM） shape functions as test functions reduces the orders of integrands involved in domain integrals. The semi-discrete heat conduction equation is solved numerically with the traditional two-point difference technique in the time domain. Two numerical examples are investigated and excellent results are obtained, demonstrating the potential application of the proposed approach.

Large Scattered Data Fitting Based on Radial Basis Functions

Solving large radial basis function （RBF） interpolation problem with non-customized methods is computationally expensive and the matrices that occur are typically badly conditioned. In order to avoid these difficulties, we present a fitting based on radial basis functions satisfying side conditions by least squares, although compared with interpolation the method loses some accuracy, it reduces the computational cost largely. Since the fitting accuracy and the non-singularity of coefficient matrix in normal equation are relevant to the uniformity of chosen centers of the fitted RBE we present a choice method of uniform centers. Numerical results confirm the fitting efficiency.

SIMULTANEOUSE APPROXIMATION TO A DIFFERENTIABLE FUNCTION AND ITS DERIVATIVES BY LAGRANGE INTERPOLATING POLYNOMIALS

This paper establishes the following pointwise result for simultancous Lagrange imterpolating approxima- tion:,then |f^(k)(x)-P_n^(k)(f,x)|=O(1)△_n^(q-k)(x)ω where P_n(f,x)is the Lagrange interpolating potynomial of deereeon the nodes X_nUY_n(see the definition of the next).

COMPACTLY SUPPORTED NON-TENSOR PRODUCT FORM TWO-DIMENSION WAVELET FINITE ELEMENT

Some theorems of compactly supported non-tensor product form two-dimension Daubechies wavelet were analysed carefully. Compactly supported non-tensor product form two-dimension wavelet was constructed, then non-tensor product form two dimension wavelet finite element was used to solve the deflection problem of elastic thin plate. The error order was researched. A numerical example was given at last.

Enhancing PIV image and fractal descriptor for velocity and shear stresses propagation around a circular pier

In this study,the fractal dimensions of velocity fluctuations and the Reynolds shear stresses propagation for flow around a circular bridge pier are presented.In the study reported herein,the fractal dimension of velocity fluctuations（u′,v′,w′） and the Reynolds shear stresses（u′v′ and u′w′） of flow around a bridge pier were computed using a Fractal Interpolation Function（FIF） algorithm.The velocity fluctuations of flow along a horizontal plane above the bed were measured using Acoustic Doppler Velocity meter（ADV）and Particle Image Velocimetry（P1V）.The PIV is a powerful technique which enables us to attain high resolution spatial and temporal information of turbulent flow using instantaneous time snapshots.In this study,PIV was used for detection of high resolution fractal scaling around a bridge pier.The results showed that the fractal dimension of flow fluctuated significantly in the longitudinal and transverse directions in the vicinity of the pier.It was also found that the fractal dimension of velocity fluctuations and shear stresses increased rapidly at vicinity of pier at downstream whereas it remained approximately unchanged far downstream of the pier.The higher value of fractal dimension was found at a distance equal to one times of the pier diameter in the back of the pier.Furthermore,the average fractal dimension for the streamwise and transverse velocity fluctuations decreased from the centreline to the side wall of the flume.Finally,the results from ADV measurement were consistent with the result from PIV,therefore,the ADV enables to detect turbulent characteristics of flow around a circular bridge pier.

An efficient large-scale mesh deformation method based on MPI/Open MP hybrid parallel radial basis function interpolation

An efficient MPI/OpenMP hybrid parallel Radial Basis Function (RBF) strategy for both continuous and discontinuous large-scale mesh deformation is proposed to reduce the computational cost and memory consumption.Unlike the conventional parallel methods in which all processors use the same surface displacement and implement the same operation,the present method employs different surface points sets and influence radius for each volume point movement,accompanied with efficient geometry searching strategy.The deformed surface points,also called Control Points (CPs),are stored in each processor.The displacement of spatial points is interpolated by using only 20-50 nearest control points,and the local influence radius is set to 5-20 times the maximum displacement of control points.To shorten the searching time for the nearest control point clouds,an Alternating Digital Tree (ADT) algorithm for 3D complex geometry is designed based on an iterative bisection technique.Besides,an MPI/OpenMP hybrid parallel approach is developed to reduce the memory cost in each High-Performance Computing (HPC) node for large-scale applications.Three 3D cases,including the ONERA-M6 wing and a commercial transport airplane standard model with up to 2.5 billion hybrid elements,are used to test the present mesh deformation method.The robustness and high parallel efficiency are demonstrated by a wing deflection case with a maximum bending angle of 450 and more than 80% parallel efficiency with 1024 MPI processors.In addition,the availability for both continuous and discontinuous surface deformation is verified by interpolating the projecting displacement with opposite directions surface points to the spatial points.

A Note on a Theorem of J. Szabados

In this note, we establish a companion result to the theorem of J. Szabados on the maximum of fundamental functions of Lagrange interpolation based on Chebyshev nodes.