The radial basis function （RBF） interpolation approach proposed by Freedman is used to solve inverse problems encountered in well-logging and other petrophysical issues. The approach is to predict petrophysical prop...The radial basis function （RBF） interpolation approach proposed by Freedman is used to solve inverse problems encountered in well-logging and other petrophysical issues. The approach is to predict petrophysical properties in the laboratory on the basis of physical rock datasets, which include the formation factor, viscosity, permeability, and molecular composition. However, this approach does not consider the effect of spatial distribution of the calibration data on the interpolation result. This study proposes a new RBF interpolation approach based on the Freedman＇s RBF interpolation approach, by which the unit basis functions are uniformly populated in the space domain. The inverse results of the two approaches are comparatively analyzed by using our datasets. We determine that although the interpolation effects of the two approaches are equivalent, the new approach is more flexible and beneficial for reducing the number of basis functions when the database is large, resulting in simplification of the interpolation function expression. However, the predicted results of the central data are not sufficiently satisfied when the data clusters are far apart.展开更多
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.Unlik...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.展开更多
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'...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.展开更多
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...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.展开更多
In this paper,a new quasi-interpolation with radial basis functions which satis- fies quadratic polynomial reproduction is constructed on the infinite set of equally spaced data.A new basis function is constructed by ...In this paper,a new quasi-interpolation with radial basis functions which satis- fies quadratic polynomial reproduction is constructed on the infinite set of equally spaced data.A new basis function is constructed by making convolution integral with a constructed spline and a given radial basis function.In particular,for twicely differ- entiable function the proposed method provides better approximation and also takes care of derivatives approximation.展开更多
Compactly supported radial basis function can enable the coefficient matrix of solving weigh linear system to have a sparse banded structure, thereby reducing the complexity of the algorithm. Firstly, based on the com...Compactly supported radial basis function can enable the coefficient matrix of solving weigh linear system to have a sparse banded structure, thereby reducing the complexity of the algorithm. Firstly, based on the compactly supported radial basis function, the paper makes the complex quadratic function （Multiquadric, MQ for short） to be transformed and proposes a class of compactly supported MQ function. Secondly, the paper describes a method that interpolates discrete motion capture data to solve the motion vectors of the interpolation points and they are used in facial expression reconstruction. Finally, according to this characteris- tic of the uneven distribution of the face markers, the markers are numbered and grouped in accordance with the density level, and then be interpolated in line with each group. The approach not only ensures the accuracy of the deformation of face local area and smoothness, but also reduces the time complexity of computing.展开更多
In many deformation analyses,the partial derivatives at the interpolated scattered data points are required.In this paper,the Gaussian Radial Basis Functions(GRBF)is proposed for the interpolation and differentiation ...In many deformation analyses,the partial derivatives at the interpolated scattered data points are required.In this paper,the Gaussian Radial Basis Functions(GRBF)is proposed for the interpolation and differentiation of the scattered data in the vertical deformation analysis.For the optimal selection of the shape parameter,which is crucial in the GRBF interpolation,two methods are used:the Power Gaussian Radial Basis Functions(PGRBF)and Leave One Out Cross Validation(LOOCV)(LGRBF).We compared the PGRBF and LGRBF to the traditional interpolation methods such as the Finite Element Method(FEM),polynomials,Moving Least Squares(MLS),and the usual GRBF in both the simulated and actual Interferometric Synthetic Aperture Radar(InSAR)data.The estimated results showed that the surface interpolation accuracy was greatly improved by LGRBF and PGRBF methods in comparison withFEM,polynomial,and MLS methods.Finally,LGRBF and PGRBF interpolation methods are used to compute invariant vertical deformation parameters,i.e.,changes in Gaussian and mean Curvatures in the Groningen area in the North of Netherlands.展开更多
This study evaluates the effectiveness of a new technique that transforms doma in integrals into boundary integrals that is applicable to the boundary element method.Si mulations were conducted in which two-dimensiona...This study evaluates the effectiveness of a new technique that transforms doma in integrals into boundary integrals that is applicable to the boundary element method.Si mulations were conducted in which two-dimensional surfaces were approximated by inter polation using radial basis functions with full and compact supports.Examples involving Poisson’s equation are presented using the boundary element method and the proposed te chnique with compact radial basis functions.The advantages and the disadvantages are e xamined through simulations.The effects of internal poles,the boundary mesh refinemen t and the value for the support of the radial basis functions on performance are assessed.展开更多
Ocean remote sensing satellites provide observations with high spatiotemporal resolution.However,the influence of clouds,fog,and haze frequently leads to significant data gaps.Accurate and effective estimation of thes...Ocean remote sensing satellites provide observations with high spatiotemporal resolution.However,the influence of clouds,fog,and haze frequently leads to significant data gaps.Accurate and effective estimation of these missing data is highly valuable for engineering and scientific research.In this study,the radial basis function(RBF)method is used to estimate the spatial distribution of total suspended matter(TSM)concentration in Hangzhou Bay using remote sensing data with severe data gaps.The estimation precision is validated by comparing the results with those of other commonly used interpolation methods,such as the Kriging method and the basic spline(B-spline)method.In addition,the applicability of the RBF method is explored.Results show that the estimation of the RBF method is significantly close to the observation in Hangzhou Bay.The average of the mean absolute error,mean relative error,and root mean square error in all the experiments is evidently smaller than those of the Kriging and B-spline interpolations,indicating that the proposed method is more appropriate for estimating the spatial distribution of the TSM in Hangzhou Bay.Finally,the TSM distribution in the blank observational area is predicted.This study can provide some reference values for handling watercolor remote sensing data.展开更多
In this paper,we present an approach for smooth surface reconstructions interpolating triangular meshes with ar-bitrary topology and geometry.The approach is based on the well-known radial basis functions(RBFs)and the...In this paper,we present an approach for smooth surface reconstructions interpolating triangular meshes with ar-bitrary topology and geometry.The approach is based on the well-known radial basis functions(RBFs)and the constructed surfaces are generalized thin-plate spline surfaces.Our algorithm first defines a pair of offset points for each vertex of a given mesh to en-hance the controUability of local geometry and to assure stability of the construction.A linear system is then solved by LU decomposi-tion and the implicit governing equation of interpolating surface is obtained.The constructed surfaces finally are visualized by a Marching Cubes based polygonizer.The approach provides a robust and efficient solution for smooth surface reconstruction from various 3 D meshes.展开更多
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 difficult...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.展开更多
In this paper,we introduce new non-polynomial basis functions for spectral approximation of time-fractional partial differential equations (PDEs). Different from many other approaches,the nonstandard singular basis fu...In this paper,we introduce new non-polynomial basis functions for spectral approximation of time-fractional partial differential equations (PDEs). Different from many other approaches,the nonstandard singular basis functions are defined from some generalised Birkhoff interpolation problems through explicit inversion of some prototypical fractional initial value problem (FIVP) with a smooth source term. As such,the singularity of the new basis can be tailored to that of the singular solutions to a class of time-fractional PDEs,leading to spectrally accurate approximation. It also provides the acceptable solution to more general singular problems.展开更多
A novel spatial interpolation method based on integrated radial basis function artificial neural networks （IRBFANNs） is proposed to provide accurate and stable predictions of heavy metals concentrations in soil at u...A novel spatial interpolation method based on integrated radial basis function artificial neural networks （IRBFANNs） is proposed to provide accurate and stable predictions of heavy metals concentrations in soil at un- sampled sites in a mountain region. The IRBFANNs hybridize the advantages of the artificial neural networks and the neural networks integration approach. Three experimental projects under different sampling densities are carried out to study the performance of the proposed IRBFANNs-based interpolation method. This novel method is compared with six peer spatial interpolation methods based on the root mean square error and visual evaluation of the distribution maps of Mn elements. The experimental results show that the proposed method performs better in accuracy and stability. Moreover, the proposed method can provide more details in the spatial distribution maps than the compared interpolation methods in the cases of sparse sampling density.展开更多
In the paper, firstly, based on new non-tensor-product-typed partially inverse divided differences algorithms in a recursive form, scattered data interpolating schemes are constructed via bivariate continued fractions...In the paper, firstly, based on new non-tensor-product-typed partially inverse divided differences algorithms in a recursive form, scattered data interpolating schemes are constructed via bivariate continued fractions with odd and even nodes, respectively. And equivalent identities are also obtained between interpolated functions and bivariate continued fractions. Secondly, by means of three-term recurrence relations for continued fractions, the characterization theorem is presented to study on the degrees of the numerators and denominators of the interpolating continued fractions. Thirdly, some numerical examples show it feasible for the novel recursive schemes. Meanwhile, compared with the degrees of the numera- tors and denominators of bivariate Thiele-typed interpolating continued fractions, those of the new bivariate interpolating continued fractions are much low, respectively, due to the reduc- tion of redundant interpolating nodes. Finally, the operation count for the rational function interpolation is smaller than that for radial basis function interpolation.展开更多
It＇s well-known that there is a very powerful error bound for Gaussians put forward by Madych and Nelson in 1992. It＇s of the form｜f（x） - s（x）｜≤（Cd）c/d｜｜f｜｜h where C, c are constants, h is the Gaussian ...It＇s well-known that there is a very powerful error bound for Gaussians put forward by Madych and Nelson in 1992. It＇s of the form｜f（x） - s（x）｜≤（Cd）c/d｜｜f｜｜h where C, c are constants, h is the Gaussian function, s is the interpolating function, and d is called fill distance which, roughly speaking, measures the spacing of the points at which interpolation occurs. This error bound gets small very fast as d → 0. The constants C and c are very sensitive. A slight change of them will result in a huge change of the error bound. The number c can be calculated as shown in [9]. However, C cannot be calculated, or even approximated. This is a famous question in the theory of radial basis functions. The purpose of this paper is to answer this question.展开更多
Accurate simulations of ultra-wideband （UWB） electromagnetic radiation from an antenna were developed based on a time-domain finite element method （TDFEM） based on p-step Lagrange interpolation for the temporal ex...Accurate simulations of ultra-wideband （UWB） electromagnetic radiation from an antenna were developed based on a time-domain finite element method （TDFEM） based on p-step Lagrange interpolation for the temporal expansion. The motivation was to utilize the good interpolation features and straightforward computations for UWB antenna simulations. Numerical results were obtained from the cases of the cavity resonance problem, a bowtie and a Sierpinski bowtie antenna. Comparisons with an existing TDFEM approach employed linear temporal basis functions show good agreement to demonstrate the validity of the present schemes. The TDFEM with 2-step Lagrange interpolation as the temporal basis functions achieves better numerical results with only a small increase to run time and memory use in terms of the relative errors of the resonant frequency in the cavity for the transverse electric mode and the radiation patterns of the bowtie antenna.展开更多
An improved interpolating complex variable element-frees Galerkin(IICVEFG)method for the two-dimensional elastic problems is developed.This method is based on the improved interpolating complex variable moving least-s...An improved interpolating complex variable element-frees Galerkin(IICVEFG)method for the two-dimensional elastic problems is developed.This method is based on the improved interpolating complex variable moving least-squares(IICVMLS)method and the integral form of the elastic problems.In the IICVEFG method,the proposed shape function has the interpolating feature.Therefore,the essential boundary conditions can be exerted directly.Additionally,the unnecessary t erms in the discrete mat rices are removed,which resul ts in a set of concise formulas.This method is verified by analyzing three elastic examples under different constraints and loads.The numerical results show that the IICVEFG method is superior in precision and efficiency to other non-interpolating meshless methods.展开更多
In this paper,general interpolating isogeometric boundary node method(IIBNM)and isogeometric boundary element method(IBEM)based on parameter space are proposed for 2D elasticity problems.In both methods,the integral c...In this paper,general interpolating isogeometric boundary node method(IIBNM)and isogeometric boundary element method(IBEM)based on parameter space are proposed for 2D elasticity problems.In both methods,the integral cells and elements are defined in parameter space,which can reproduce the geometry exactly at all the stages.In IIBNM,the improved interpolating moving leastsquare method(IIMLS)is applied for field approximation and the shape functions have the delta function property.The Lagrangian basis functions are used for field approximation in IBEM.Thus,the boundary conditions can be imposed directly in both methods.The shape functions are defined in 1D parameter space and no curve length needs to be computed.Besides,most methods for the treatment of the singular integrals in the boundary element method can be applied in IIBNM and IBEM directly.Numerical examples have demonstrated the accuracy of the proposed methods.展开更多
基金supported by the National Science and Technology Major Projects(No.2011ZX05020-008)Well Logging Advanced Technique and Application Basis Research Project of Petrochina Company(No.2011A-3901)
文摘The radial basis function （RBF） interpolation approach proposed by Freedman is used to solve inverse problems encountered in well-logging and other petrophysical issues. The approach is to predict petrophysical properties in the laboratory on the basis of physical rock datasets, which include the formation factor, viscosity, permeability, and molecular composition. However, this approach does not consider the effect of spatial distribution of the calibration data on the interpolation result. This study proposes a new RBF interpolation approach based on the Freedman＇s RBF interpolation approach, by which the unit basis functions are uniformly populated in the space domain. The inverse results of the two approaches are comparatively analyzed by using our datasets. We determine that although the interpolation effects of the two approaches are equivalent, the new approach is more flexible and beneficial for reducing the number of basis functions when the database is large, resulting in simplification of the interpolation function expression. However, the predicted results of the central data are not sufficiently satisfied when the data clusters are far apart.
基金supported by the National Key Research and Development Program of China (No.2016YFB0200701)the National Natural Science Foundation of China (Nos. 11532016 and 91530325)
文摘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.
文摘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.
文摘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.
文摘In this paper,a new quasi-interpolation with radial basis functions which satis- fies quadratic polynomial reproduction is constructed on the infinite set of equally spaced data.A new basis function is constructed by making convolution integral with a constructed spline and a given radial basis function.In particular,for twicely differ- entiable function the proposed method provides better approximation and also takes care of derivatives approximation.
基金Supported by the National Natural Science Foundation of China (No.60875046)by Program for Changjiang Scholars and Innovative Research Team in University(No.IRT1109)+5 种基金the Key Project of Chinese Ministry of Education (No.209029)the Program for Liaoning Excellent Talents in University(No.LR201003)the Program for Liaoning Science and Technology Research in University (No.LS2010008,2009S008,2009S009,LS2010179)the Program for Liaoning Innovative Research Team in University(Nos.2009T005,LT2010005,LT2011018)Natural Science Foundation of Liaoning Province (201102008)by "Liaoning BaiQianWan Talents Program(2010921010,2011921009)"
文摘Compactly supported radial basis function can enable the coefficient matrix of solving weigh linear system to have a sparse banded structure, thereby reducing the complexity of the algorithm. Firstly, based on the compactly supported radial basis function, the paper makes the complex quadratic function （Multiquadric, MQ for short） to be transformed and proposes a class of compactly supported MQ function. Secondly, the paper describes a method that interpolates discrete motion capture data to solve the motion vectors of the interpolation points and they are used in facial expression reconstruction. Finally, according to this characteris- tic of the uneven distribution of the face markers, the markers are numbered and grouped in accordance with the density level, and then be interpolated in line with each group. The approach not only ensures the accuracy of the deformation of face local area and smoothness, but also reduces the time complexity of computing.
文摘In many deformation analyses,the partial derivatives at the interpolated scattered data points are required.In this paper,the Gaussian Radial Basis Functions(GRBF)is proposed for the interpolation and differentiation of the scattered data in the vertical deformation analysis.For the optimal selection of the shape parameter,which is crucial in the GRBF interpolation,two methods are used:the Power Gaussian Radial Basis Functions(PGRBF)and Leave One Out Cross Validation(LOOCV)(LGRBF).We compared the PGRBF and LGRBF to the traditional interpolation methods such as the Finite Element Method(FEM),polynomials,Moving Least Squares(MLS),and the usual GRBF in both the simulated and actual Interferometric Synthetic Aperture Radar(InSAR)data.The estimated results showed that the surface interpolation accuracy was greatly improved by LGRBF and PGRBF methods in comparison withFEM,polynomial,and MLS methods.Finally,LGRBF and PGRBF interpolation methods are used to compute invariant vertical deformation parameters,i.e.,changes in Gaussian and mean Curvatures in the Groningen area in the North of Netherlands.
文摘This study evaluates the effectiveness of a new technique that transforms doma in integrals into boundary integrals that is applicable to the boundary element method.Si mulations were conducted in which two-dimensional surfaces were approximated by inter polation using radial basis functions with full and compact supports.Examples involving Poisson’s equation are presented using the boundary element method and the proposed te chnique with compact radial basis functions.The advantages and the disadvantages are e xamined through simulations.The effects of internal poles,the boundary mesh refinemen t and the value for the support of the radial basis functions on performance are assessed.
基金supported by the Open Funds for Hubei Key Laboratory of Marine Geological Resources,China University of Geosciences(No.MGR202308)the Natural Science Foundation of Shandong Province(No.ZR2020MD085)+3 种基金the National Natural Science Foundation of China(No.41821004)the Taishan Scholar Program(No.tstp2022114)the Shandong Provincial Natural Science Foundation(No.DKXZZ202206)the National Key Research and Development Program of China(No.2016YFC1402404).
文摘Ocean remote sensing satellites provide observations with high spatiotemporal resolution.However,the influence of clouds,fog,and haze frequently leads to significant data gaps.Accurate and effective estimation of these missing data is highly valuable for engineering and scientific research.In this study,the radial basis function(RBF)method is used to estimate the spatial distribution of total suspended matter(TSM)concentration in Hangzhou Bay using remote sensing data with severe data gaps.The estimation precision is validated by comparing the results with those of other commonly used interpolation methods,such as the Kriging method and the basic spline(B-spline)method.In addition,the applicability of the RBF method is explored.Results show that the estimation of the RBF method is significantly close to the observation in Hangzhou Bay.The average of the mean absolute error,mean relative error,and root mean square error in all the experiments is evidently smaller than those of the Kriging and B-spline interpolations,indicating that the proposed method is more appropriate for estimating the spatial distribution of the TSM in Hangzhou Bay.Finally,the TSM distribution in the blank observational area is predicted.This study can provide some reference values for handling watercolor remote sensing data.
基金This work presented in this paper is supported by research grant of National Natural Science Foundation ofChina(No.60503058).
文摘In this paper,we present an approach for smooth surface reconstructions interpolating triangular meshes with ar-bitrary topology and geometry.The approach is based on the well-known radial basis functions(RBFs)and the constructed surfaces are generalized thin-plate spline surfaces.Our algorithm first defines a pair of offset points for each vertex of a given mesh to en-hance the controUability of local geometry and to assure stability of the construction.A linear system is then solved by LU decomposi-tion and the implicit governing equation of interpolating surface is obtained.The constructed surfaces finally are visualized by a Marching Cubes based polygonizer.The approach provides a robust and efficient solution for smooth surface reconstruction from various 3 D meshes.
基金Supported by National Natural Science Youth Foundation (10401021).
文摘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.
基金the China Postdoctoral Science Foundation Funded Project (No.2017M620113)the National Natural Science Foundation of China (Nos.11801120,71773024 and 11771107)+4 种基金the Fundamental Research Funds for the Central Universities (Grant No.HIT.NSRIF.2019058)the Natural Science Foundation of Heilongjiang Province of China (No.G2018006)Singapore MOE AcRF Tier 2 Grants (MOE2017-T2-2-014 and MOE2018-T2-1-059)National Science Foundation of China (No.11371376)the Innovation-Driven Project and Mathematics.
文摘In this paper,we introduce new non-polynomial basis functions for spectral approximation of time-fractional partial differential equations (PDEs). Different from many other approaches,the nonstandard singular basis functions are defined from some generalised Birkhoff interpolation problems through explicit inversion of some prototypical fractional initial value problem (FIVP) with a smooth source term. As such,the singularity of the new basis can be tailored to that of the singular solutions to a class of time-fractional PDEs,leading to spectrally accurate approximation. It also provides the acceptable solution to more general singular problems.
基金The National Natural Science Foundation of China(No.61261007,61062005)the Key Program of Yunnan Natural Science Foundation(No.2013FA008)
文摘A novel spatial interpolation method based on integrated radial basis function artificial neural networks （IRBFANNs） is proposed to provide accurate and stable predictions of heavy metals concentrations in soil at un- sampled sites in a mountain region. The IRBFANNs hybridize the advantages of the artificial neural networks and the neural networks integration approach. Three experimental projects under different sampling densities are carried out to study the performance of the proposed IRBFANNs-based interpolation method. This novel method is compared with six peer spatial interpolation methods based on the root mean square error and visual evaluation of the distribution maps of Mn elements. The experimental results show that the proposed method performs better in accuracy and stability. Moreover, the proposed method can provide more details in the spatial distribution maps than the compared interpolation methods in the cases of sparse sampling density.
基金Supported by the Special Funds Tianyuan for the National Natural Science Foundation of China(Grant No.11426086)the Fundamental Research Funds for the Central Universities(Grant No.2016B08714)the Natural Science Foundation of Jiangsu Province for the Youth(Grant No.BK20160853)
文摘In the paper, firstly, based on new non-tensor-product-typed partially inverse divided differences algorithms in a recursive form, scattered data interpolating schemes are constructed via bivariate continued fractions with odd and even nodes, respectively. And equivalent identities are also obtained between interpolated functions and bivariate continued fractions. Secondly, by means of three-term recurrence relations for continued fractions, the characterization theorem is presented to study on the degrees of the numerators and denominators of the interpolating continued fractions. Thirdly, some numerical examples show it feasible for the novel recursive schemes. Meanwhile, compared with the degrees of the numera- tors and denominators of bivariate Thiele-typed interpolating continued fractions, those of the new bivariate interpolating continued fractions are much low, respectively, due to the reduc- tion of redundant interpolating nodes. Finally, the operation count for the rational function interpolation is smaller than that for radial basis function interpolation.
文摘It＇s well-known that there is a very powerful error bound for Gaussians put forward by Madych and Nelson in 1992. It＇s of the form｜f（x） - s（x）｜≤（Cd）c/d｜｜f｜｜h where C, c are constants, h is the Gaussian function, s is the interpolating function, and d is called fill distance which, roughly speaking, measures the spacing of the points at which interpolation occurs. This error bound gets small very fast as d → 0. The constants C and c are very sensitive. A slight change of them will result in a huge change of the error bound. The number c can be calculated as shown in [9]. However, C cannot be calculated, or even approximated. This is a famous question in the theory of radial basis functions. The purpose of this paper is to answer this question.
文摘Accurate simulations of ultra-wideband （UWB） electromagnetic radiation from an antenna were developed based on a time-domain finite element method （TDFEM） based on p-step Lagrange interpolation for the temporal expansion. The motivation was to utilize the good interpolation features and straightforward computations for UWB antenna simulations. Numerical results were obtained from the cases of the cavity resonance problem, a bowtie and a Sierpinski bowtie antenna. Comparisons with an existing TDFEM approach employed linear temporal basis functions show good agreement to demonstrate the validity of the present schemes. The TDFEM with 2-step Lagrange interpolation as the temporal basis functions achieves better numerical results with only a small increase to run time and memory use in terms of the relative errors of the resonant frequency in the cavity for the transverse electric mode and the radiation patterns of the bowtie antenna.
基金The authors sincerely acknowledge the financial support from the National Science Foundation of China(No.12002240)the National Science and Technology Major Project(No.2017-IV-0003-0040).
文摘An improved interpolating complex variable element-frees Galerkin(IICVEFG)method for the two-dimensional elastic problems is developed.This method is based on the improved interpolating complex variable moving least-squares(IICVMLS)method and the integral form of the elastic problems.In the IICVEFG method,the proposed shape function has the interpolating feature.Therefore,the essential boundary conditions can be exerted directly.Additionally,the unnecessary t erms in the discrete mat rices are removed,which resul ts in a set of concise formulas.This method is verified by analyzing three elastic examples under different constraints and loads.The numerical results show that the IICVEFG method is superior in precision and efficiency to other non-interpolating meshless methods.
基金The research for this paper was supported by(1)the National Natural Science Foundation of China(Grants Nos.51708429,51708428)the Open Projects Foundation(Grant No.2017-04-GF)of State Key Laboratory for Health and Safety of Bridge Structures+1 种基金Wuhan Institute of Technology Science Found(Grant No.K201734)the science and technology projects of Wuhan Urban and Rural Construction Bureau(Grants Nos.201831,201919).
文摘In this paper,general interpolating isogeometric boundary node method(IIBNM)and isogeometric boundary element method(IBEM)based on parameter space are proposed for 2D elasticity problems.In both methods,the integral cells and elements are defined in parameter space,which can reproduce the geometry exactly at all the stages.In IIBNM,the improved interpolating moving leastsquare method(IIMLS)is applied for field approximation and the shape functions have the delta function property.The Lagrangian basis functions are used for field approximation in IBEM.Thus,the boundary conditions can be imposed directly in both methods.The shape functions are defined in 1D parameter space and no curve length needs to be computed.Besides,most methods for the treatment of the singular integrals in the boundary element method can be applied in IIBNM and IBEM directly.Numerical examples have demonstrated the accuracy of the proposed methods.
