A high-precision regional gravity field model is significant in various geodesy applications.In the field of modelling regional gravity fields,the spherical radial basis functions(SRBFs)approach has recently gained wi...A high-precision regional gravity field model is significant in various geodesy applications.In the field of modelling regional gravity fields,the spherical radial basis functions(SRBFs)approach has recently gained widespread attention,while the modelling precision is primarily influenced by the base function network.In this study,we propose a method for constructing a data-adaptive network of SRBFs using a modified Hierarchical Density-Based Spatial Clustering of Applications with Noise(HDBSCAN)algorithm,and the performance of the algorithm is verified by the observed gravity data in the Auvergne area.Furthermore,the turning point method is used to optimize the bandwidth of the basis function spectrum,which satisfies the demand for both high-precision gravity field and quasi-geoid modelling simultaneously.Numerical experimental results indicate that our algorithm has an accuracy of about 1.58 mGal in constructing the gravity field model and about 0.03 m in the regional quasi-geoid model.Compared to the existing methods,the number of SRBFs used for modelling has been reduced by 15.8%,and the time cost to determine the centre positions of SRBFs has been saved by 12.5%.Hence,the modified HDBSCAN algorithm presented here is a suitable design method for constructing the SRBF data adaptive network.展开更多
Dam-break flows pose significant threats to urban areas due to their potential for causing rapid and extensive flooding. Traditional numerical methods for simulating these events struggle with complex urban landscapes...Dam-break flows pose significant threats to urban areas due to their potential for causing rapid and extensive flooding. Traditional numerical methods for simulating these events struggle with complex urban landscapes. This paper presents an alternative approach using Radial Basis Functions to simulate dam-break flows and their impact on urban flood inundation. The proposed method adapts a new strategy based on Particle Swarm Optimization for variable shape parameter selection on meshfree formulation to enhance the numerical stability and convergence of the simulation. The method’s accuracy and efficiency are demonstrated through numerical experiments, including well-known partial and circular dam-break problems and an idealized city with a single building, highlighting its potential as a valuable tool for urban flood risk management.展开更多
To improve the nonlinear approximating ability of cerebellar model articulation controller(CMAC), by introducing the Gauss basis functions and the similarity measure based addressing scheme, a new kind of fuzzy CMAC...To improve the nonlinear approximating ability of cerebellar model articulation controller(CMAC), by introducing the Gauss basis functions and the similarity measure based addressing scheme, a new kind of fuzzy CMAC with Gauss basis functions(GFCMAC) was presented. Moreover, based upon the improvement of the self organizing feature map algorithm of Kohonen, the structural self organizing algorithm for GFCMAC(SOGFCMAC) was proposed. Simulation results show that adopting the Gauss basis functions and fuzzy techniques can remarkably improve the nonlinear approximating capacity of CMAC. Compared with the traditional CMAC,CMAC with general basis functions and fuzzy CMAC(FCMAC), SOGFCMAC has the obvious advantages in the aspects of the convergent speed, approximating accuracy and structural self organizing.展开更多
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.展开更多
The application of Tikhonov regularization method dealing with the ill-conditioned problems in the regional gravity field modeling by Poisson wavelets is studied. In particular, the choices of the regularization matri...The application of Tikhonov regularization method dealing with the ill-conditioned problems in the regional gravity field modeling by Poisson wavelets is studied. In particular, the choices of the regularization matrices as well as the approaches for estimating the regularization parameters are investigated in details. The numerical results show that the regularized solutions derived from the first-order regularization are better than the ones obtained from zero-order regularization. For cross validation, the optimal regularization parameters are estimated from L-curve, variance component estimation（VCE） and minimum standard deviation（MSTD） approach, respectively, and the results show that the derived regularization parameters from different methods are consistent with each other. Together with the firstorder Tikhonov regularization and VCE method, the optimal network of Poisson wavelets is derived, based on which the local gravimetric geoid is computed. The accuracy of the corresponding gravimetric geoid reaches 1.1 cm in Netherlands, which validates the reliability of using Tikhonov regularization method in tackling the ill-conditioned problem for regional gravity field modeling.展开更多
In this paper, radial basis functions are used to obtain the solution of evolution equations which appear in variational level set method based image segmentation. In this method, radial basis functions are used to in...In this paper, radial basis functions are used to obtain the solution of evolution equations which appear in variational level set method based image segmentation. In this method, radial basis functions are used to interpolate the implicit level set function of the evolution equation with a high level of accuracy and smoothness. Then, the original initial value problem is discretized into an interpolation problem. Accordingly, the evolution equation is converted into a set of coupled ordinary differential equations, and a smooth evolution can be retained. Compared with finite difference scheme based level set approaches, the complex and costly re-initialization procedure is unnecessary. Numerical examples are also given to show the efficiency of the method.展开更多
This paper introduces the use of partition of unity method for the development of a high order finite volume discretization scheme on unstructured grids for solving diffusion models based on partial differential equat...This paper introduces the use of partition of unity method for the development of a high order finite volume discretization scheme on unstructured grids for solving diffusion models based on partial differential equations.The unknown function and its gradient can be accurately reconstructed using high order optimal recovery based on radial basis functions.The methodology proposed is applied to the noise removal problem in functional surfaces and images.Numerical results demonstrate the effectiveness of the new numerical approach and provide experimental order of convergence.展开更多
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.展开更多
Ensemble simulations, which use multiple short independent trajectories from dispersive initial conformations, rather than a single long trajectory as used in traditional simulations, are expected to sample complex sy...Ensemble simulations, which use multiple short independent trajectories from dispersive initial conformations, rather than a single long trajectory as used in traditional simulations, are expected to sample complex systems such as biomolecules much more efficiently. The re-weighted ensemble dynamics（RED） is designed to combine these short trajectories to reconstruct the global equilibrium distribution. In the RED, a number of conformational functions, named as basis functions,are applied to relate these trajectories to each other, then a detailed-balance-based linear equation is built, whose solution provides the weights of these trajectories in equilibrium distribution. Thus, the sufficient and efficient selection of basis functions is critical to the practical application of RED. Here, we review and present a few possible ways to generally construct basis functions for applying the RED in complex molecular systems. Especially, for systems with less priori knowledge, we could generally use the root mean squared deviation（RMSD） among conformations to split the whole conformational space into a set of cells, then use the RMSD-based-cell functions as basis functions. We demonstrate the application of the RED in typical systems, including a two-dimensional toy model, the lattice Potts model, and a short peptide system. The results indicate that the RED with the constructions of basis functions not only more efficiently sample the complex systems, but also provide a general way to understand the metastable structure of conformational space.展开更多
The purpose of this paper is to present the class of atomic basis functions(ABFs)which are of exponential type and are denoted by EFupn(x,ω).While ABFs of the algebraic type are already represented in the numerical m...The purpose of this paper is to present the class of atomic basis functions(ABFs)which are of exponential type and are denoted by EFupn(x,ω).While ABFs of the algebraic type are already represented in the numerical modeling of various problems inmathematical physics and computationalmechanics,ABFs of the exponential type have not yet been sufficiently researched.These functions,unlike the ABFs of the algebraic type Fupn(x),contain the tension parameterω,which gives them additional approximation properties.Exponential monomials up to the nth degree can be described exactly by the linear combination of the functions EFupn(x,ω).The function EFupn for n=0 is called the“mother”ABF of the exponential type,i.e.,EFup0(x,ω)≡Eup(x,ω).In other words,the functions EFupn(x,ω)are elements of the linear vector space EUPn and retain all the properties of their“mother”function Eup(x,ω).Thus,this paper,in terms of its content and purpose,can be understood as a sequel of the article by Brajcic Kurbasa et al.,which shows the basic properties and application of the basis function Eup(x,ω).This paper presents,in an analogous way,the development and application of the exponential basis functions EFupn(x,ω).Here,for the first time,expressions for calculating the values of the functions EFupn(x,ω)and their derivatives are given in a form suitable for application in numerical analyses,which is shown in the verification examples of the approximations of known functions.展开更多
It is important to reconstruct a continuous surface representation of the point cloud scanned from a human body. In this paper a new implicit surface method is proposed to reconstruct the human body surface from the p...It is important to reconstruct a continuous surface representation of the point cloud scanned from a human body. In this paper a new implicit surface method is proposed to reconstruct the human body surface from the points based on the combination of radial basis functions （RBFs） and adaptive partition of unity （PoU）. The whole 3D domain of the scanned human body is firstly subdivided into a set of overlapping subdomalns based on the improved octrees. The smooth local surfaces are then computed in the subdomalns based on RBFs. And finally the global human body surface is reconstructed by blending the local surfaces with the adaptive PoU functions. This method is robust for the surface reconstruction of the scanned human body even with large or non-uniform point cloud which has a sharp density variation.展开更多
The radial basis functions(RBFs)play an important role in the numerical simulation processes of partial differential equations.Since the radial basis functions are meshless algorithms,its approximation is easy to impl...The radial basis functions(RBFs)play an important role in the numerical simulation processes of partial differential equations.Since the radial basis functions are meshless algorithms,its approximation is easy to implement and mathematically simple.In this paper,the commonly⁃used multiquadric RBF,conical RBF,and Gaussian RBF were applied to solve boundary value problems which are governed by partial differential equations with variable coefficients.Numerical results were provided to show the good performance of the three RBFs as numerical tools for a wide range of problems.It is shown that the conical RBF numerical results were more stable than the other two radial basis functions.From the comparison of three commonly⁃used RBFs,one may obtain the best numerical solutions for boundary value problems.展开更多
The explicit expression of the G3 basis function is presented in this paper. It is derived by constructing the conversion matrix between G3 basis function and Brzier representation. After the matrix decomposition, equ...The explicit expression of the G3 basis function is presented in this paper. It is derived by constructing the conversion matrix between G3 basis function and Brzier representation. After the matrix decomposition, equations for constructing G3 splines can be presented independently of geometric shape parameters＇ values. It makes the equation＇s solving easier. It is also known that the general form of the G3spline basis function is given in the first time. Its geometric construction method is presented.展开更多
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.展开更多
Continuously differentiable radial basis functions (C<sup>∞</sup>-RBFs), while being theoretically exponentially convergent are considered impractical computationally because the coefficient matrices are ...Continuously differentiable radial basis functions (C<sup>∞</sup>-RBFs), while being theoretically exponentially convergent are considered impractical computationally because the coefficient matrices are full and can become very ill- conditioned. Similarly, the Hilbert and Vandermonde have full matrices and become ill-conditioned. The difference between a coefficient matrix generated by C<sup>∞</sup>-RBFs for partial differential or integral equations and Hilbert and Vandermonde systems is that C<sup>∞</sup>-RBFs are very sensitive to small changes in the adjustable parameters. These parameters affect the condition number and solution accuracy. The error terrain has many local and global maxima and minima. To find stable and accurate numerical solutions for full linear equation systems, this study proposes a hybrid combination of block Gaussian elimination (BGE) combined with arbitrary precision arithmetic (APA) to minimize the accumulation of rounding errors. In the future, this algorithm can execute faster using preconditioners and implemented on massively parallel computers.展开更多
The Radial Basis Function（RBF） method with data reduction is an effective way to perform mesh deformation. However, for large deformations on meshes of complex aerodynamic configurations, the efficiency of the RBF m...The Radial Basis Function（RBF） method with data reduction is an effective way to perform mesh deformation. However, for large deformations on meshes of complex aerodynamic configurations, the efficiency of the RBF mesh deformation method still needs to be further improved to fulfill the demand of practical application. To achieve this goal, a multistep RBF method based on a multilevel subspace RBF algorithm is presented to further improve the efficiency of the mesh deformation method in this research. A whole deformation is divided into a series of steps, and the supporting radius is adjusted in accordance with the maximal displacement error. Furthermore, parallel computing is applied to the interpolation to enhance the efficiency. Typical deformation problems of the NASA Common Research Model（CRM） configuration, the DLR-F6 wing-body-nacellepylon configuration, and the DLR-F11 high-lift configuration are tested to verify the feasibility of this method. Test results show that the presented multistep RBF mesh deformation method is efficient and robust in dealing with large deformation problems over complex geometries.展开更多
The present work describes the application of the method of fundamental solutions （MFS） along with the analog equation method （AEM） and radial basis function （RBF） approximation for solving the 2D isotropic and ...The present work describes the application of the method of fundamental solutions （MFS） along with the analog equation method （AEM） and radial basis function （RBF） approximation for solving the 2D isotropic and anisotropic Helmholtz problems with different wave numbers. The AEM is used to convert the original governing equation into the classical Poisson＇s equation, and the MFS and RBF approximations are used to derive the homogeneous and particular solutions, respectively. Finally, the satisfaction of the solution consisting of the homogeneous and particular parts to the related governing equation and boundary conditions can produce a system of linear equations, which can be solved with the singular value decomposition （SVD） technique. In the computation, such crucial factors related to the MFS-RBF as the location of the virtual boundary, the differential and integrating strategies, and the variation of shape parameters in multi-quadric （MQ） are fully analyzed to provide useful reference.展开更多
Based on our previous study,the accuracy of derivatives of interpolating functions are usually very poor near the boundary of domain when Compactly Supported Radial Basis Functions (CSRBFs)are used,so that it could re...Based on our previous study,the accuracy of derivatives of interpolating functions are usually very poor near the boundary of domain when Compactly Supported Radial Basis Functions (CSRBFs)are used,so that it could result in significant error in solving partial differential equations with Neumann boundary conditions.To overcome this drawback,the Consistent Compactly Supported Radial Basis Functions(CCSRBFs)are developed,which satisfy the predetermined consistency con- ditions.Meshless method based on point collocation with CCSRBFs is developed for solving partial differential equations.Numerical studies show that the proposed method improves the accuracy of approximation significantly.展开更多
The Radial Basis Functions Neural Network （RBFNN） is used to establish the model of a response system through the input and output data of the system. The synchronization between a drive system and the response syst...The Radial Basis Functions Neural Network （RBFNN） is used to establish the model of a response system through the input and output data of the system. The synchronization between a drive system and the response system can be implemented by employing the RBFNN model and state feedback control. In this case, the exact mathematical model, which is the precondition for the conventional method, is unnecessary for implementing synchronization. The effect of the model error is investigated and a corresponding theorem is developed. The effect of the parameter perturbations and the measurement noise is investigated through simulations. The simulation results under different conditions show the effectiveness of the method.展开更多
A boundary integral method with radial basis function approximation is proposed for numerically solving an important class of boundary value problems governed by a system of thermoelastostatic equations with variable ...A boundary integral method with radial basis function approximation is proposed for numerically solving an important class of boundary value problems governed by a system of thermoelastostatic equations with variable coe?cients. The equations describe the thermoelastic behaviors of nonhomogeneous anisotropic materials with properties that vary smoothly from point to point in space. No restriction is imposed on the spatial variations of the thermoelastic coe?cients as long as all the requirements of the laws of physics are satis?ed. To check the validity and accuracy of the proposed numerical method, some speci?c test problems with known solutions are solved.展开更多
基金funded by The Fundamental Research Funds for Chinese Academy of surveying and mapping(AR2402)Open Fund of Wuhan,Gravitation and Solid Earth Tides,National Observation and Research Station(No.WHYWZ202213)。
文摘A high-precision regional gravity field model is significant in various geodesy applications.In the field of modelling regional gravity fields,the spherical radial basis functions(SRBFs)approach has recently gained widespread attention,while the modelling precision is primarily influenced by the base function network.In this study,we propose a method for constructing a data-adaptive network of SRBFs using a modified Hierarchical Density-Based Spatial Clustering of Applications with Noise(HDBSCAN)algorithm,and the performance of the algorithm is verified by the observed gravity data in the Auvergne area.Furthermore,the turning point method is used to optimize the bandwidth of the basis function spectrum,which satisfies the demand for both high-precision gravity field and quasi-geoid modelling simultaneously.Numerical experimental results indicate that our algorithm has an accuracy of about 1.58 mGal in constructing the gravity field model and about 0.03 m in the regional quasi-geoid model.Compared to the existing methods,the number of SRBFs used for modelling has been reduced by 15.8%,and the time cost to determine the centre positions of SRBFs has been saved by 12.5%.Hence,the modified HDBSCAN algorithm presented here is a suitable design method for constructing the SRBF data adaptive network.
文摘Dam-break flows pose significant threats to urban areas due to their potential for causing rapid and extensive flooding. Traditional numerical methods for simulating these events struggle with complex urban landscapes. This paper presents an alternative approach using Radial Basis Functions to simulate dam-break flows and their impact on urban flood inundation. The proposed method adapts a new strategy based on Particle Swarm Optimization for variable shape parameter selection on meshfree formulation to enhance the numerical stability and convergence of the simulation. The method’s accuracy and efficiency are demonstrated through numerical experiments, including well-known partial and circular dam-break problems and an idealized city with a single building, highlighting its potential as a valuable tool for urban flood risk management.
文摘To improve the nonlinear approximating ability of cerebellar model articulation controller(CMAC), by introducing the Gauss basis functions and the similarity measure based addressing scheme, a new kind of fuzzy CMAC with Gauss basis functions(GFCMAC) was presented. Moreover, based upon the improvement of the self organizing feature map algorithm of Kohonen, the structural self organizing algorithm for GFCMAC(SOGFCMAC) was proposed. Simulation results show that adopting the Gauss basis functions and fuzzy techniques can remarkably improve the nonlinear approximating capacity of CMAC. Compared with the traditional CMAC,CMAC with general basis functions and fuzzy CMAC(FCMAC), SOGFCMAC has the obvious advantages in the aspects of the convergent speed, approximating accuracy and structural self organizing.
文摘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.
基金supported by the National Natural Science Foundation of China (Nos.41374023,41131067,41474019)the National 973 Project of China (No.2013CB733302)+2 种基金the China Postdoctoral Science Foundation (No.2016M602301)the Key Laboratory of Geospace Envi-ronment and Geodesy,Ministry of Education,Wuhan University (No.15-02-08)the State Scholarship Fund from Chinese Scholarship Council (No.201306270014)
文摘The application of Tikhonov regularization method dealing with the ill-conditioned problems in the regional gravity field modeling by Poisson wavelets is studied. In particular, the choices of the regularization matrices as well as the approaches for estimating the regularization parameters are investigated in details. The numerical results show that the regularized solutions derived from the first-order regularization are better than the ones obtained from zero-order regularization. For cross validation, the optimal regularization parameters are estimated from L-curve, variance component estimation（VCE） and minimum standard deviation（MSTD） approach, respectively, and the results show that the derived regularization parameters from different methods are consistent with each other. Together with the firstorder Tikhonov regularization and VCE method, the optimal network of Poisson wavelets is derived, based on which the local gravimetric geoid is computed. The accuracy of the corresponding gravimetric geoid reaches 1.1 cm in Netherlands, which validates the reliability of using Tikhonov regularization method in tackling the ill-conditioned problem for regional gravity field modeling.
基金Project supported by the National Natural Science Foundation of China (Grant No.11101454)the Educational Commission Foundation of Chongqing City,China (Grant No.KJ130626)the Program of Innovation Team Project in University of Chongqing City,China (Grant No.KJTD201308)
文摘In this paper, radial basis functions are used to obtain the solution of evolution equations which appear in variational level set method based image segmentation. In this method, radial basis functions are used to interpolate the implicit level set function of the evolution equation with a high level of accuracy and smoothness. Then, the original initial value problem is discretized into an interpolation problem. Accordingly, the evolution equation is converted into a set of coupled ordinary differential equations, and a smooth evolution can be retained. Compared with finite difference scheme based level set approaches, the complex and costly re-initialization procedure is unnecessary. Numerical examples are also given to show the efficiency of the method.
基金supported by PRIN-MIUR-Cofin 2006by University of Bologna"Funds for selected research topics"
文摘This paper introduces the use of partition of unity method for the development of a high order finite volume discretization scheme on unstructured grids for solving diffusion models based on partial differential equations.The unknown function and its gradient can be accurately reconstructed using high order optimal recovery based on radial basis functions.The methodology proposed is applied to the noise removal problem in functional surfaces and images.Numerical results demonstrate the effectiveness of the new numerical approach and provide experimental order of convergence.
基金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.
基金Project supported by the National Natural Science Foundation of China(Grant No.11175250)
文摘Ensemble simulations, which use multiple short independent trajectories from dispersive initial conformations, rather than a single long trajectory as used in traditional simulations, are expected to sample complex systems such as biomolecules much more efficiently. The re-weighted ensemble dynamics（RED） is designed to combine these short trajectories to reconstruct the global equilibrium distribution. In the RED, a number of conformational functions, named as basis functions,are applied to relate these trajectories to each other, then a detailed-balance-based linear equation is built, whose solution provides the weights of these trajectories in equilibrium distribution. Thus, the sufficient and efficient selection of basis functions is critical to the practical application of RED. Here, we review and present a few possible ways to generally construct basis functions for applying the RED in complex molecular systems. Especially, for systems with less priori knowledge, we could generally use the root mean squared deviation（RMSD） among conformations to split the whole conformational space into a set of cells, then use the RMSD-based-cell functions as basis functions. We demonstrate the application of the RED in typical systems, including a two-dimensional toy model, the lattice Potts model, and a short peptide system. The results indicate that the RED with the constructions of basis functions not only more efficiently sample the complex systems, but also provide a general way to understand the metastable structure of conformational space.
基金supported through Project KK.01.1.1.02.0027a project co-financed by the Croatian Government and the European Union through the European Regional Development Fund-the Competitiveness and Cohesion Operational Programme.
文摘The purpose of this paper is to present the class of atomic basis functions(ABFs)which are of exponential type and are denoted by EFupn(x,ω).While ABFs of the algebraic type are already represented in the numerical modeling of various problems inmathematical physics and computationalmechanics,ABFs of the exponential type have not yet been sufficiently researched.These functions,unlike the ABFs of the algebraic type Fupn(x),contain the tension parameterω,which gives them additional approximation properties.Exponential monomials up to the nth degree can be described exactly by the linear combination of the functions EFupn(x,ω).The function EFupn for n=0 is called the“mother”ABF of the exponential type,i.e.,EFup0(x,ω)≡Eup(x,ω).In other words,the functions EFupn(x,ω)are elements of the linear vector space EUPn and retain all the properties of their“mother”function Eup(x,ω).Thus,this paper,in terms of its content and purpose,can be understood as a sequel of the article by Brajcic Kurbasa et al.,which shows the basic properties and application of the basis function Eup(x,ω).This paper presents,in an analogous way,the development and application of the exponential basis functions EFupn(x,ω).Here,for the first time,expressions for calculating the values of the functions EFupn(x,ω)and their derivatives are given in a form suitable for application in numerical analyses,which is shown in the verification examples of the approximations of known functions.
基金the National Natural Science Foundation of China (No. 50575139)the Shanghai Special Fund of Informatization (No. 088)
文摘It is important to reconstruct a continuous surface representation of the point cloud scanned from a human body. In this paper a new implicit surface method is proposed to reconstruct the human body surface from the points based on the combination of radial basis functions （RBFs） and adaptive partition of unity （PoU）. The whole 3D domain of the scanned human body is firstly subdivided into a set of overlapping subdomalns based on the improved octrees. The smooth local surfaces are then computed in the subdomalns based on RBFs. And finally the global human body surface is reconstructed by blending the local surfaces with the adaptive PoU functions. This method is robust for the surface reconstruction of the scanned human body even with large or non-uniform point cloud which has a sharp density variation.
基金the Natural Science Foundation of Anhui Province(Grant No.1908085QA09)the University Natural Science Research Project of Anhui Province(KJ2019A0591).
文摘The radial basis functions(RBFs)play an important role in the numerical simulation processes of partial differential equations.Since the radial basis functions are meshless algorithms,its approximation is easy to implement and mathematically simple.In this paper,the commonly⁃used multiquadric RBF,conical RBF,and Gaussian RBF were applied to solve boundary value problems which are governed by partial differential equations with variable coefficients.Numerical results were provided to show the good performance of the three RBFs as numerical tools for a wide range of problems.It is shown that the conical RBF numerical results were more stable than the other two radial basis functions.From the comparison of three commonly⁃used RBFs,one may obtain the best numerical solutions for boundary value problems.
基金Supported by National Natural Science Foundation of China(Grants 61100129)Open Program of Key Laboratory of Intelligent Information Processing,Institute of Computing Technology,Chinese Academy of Sciences(IIP2014-7)
文摘The explicit expression of the G3 basis function is presented in this paper. It is derived by constructing the conversion matrix between G3 basis function and Brzier representation. After the matrix decomposition, equations for constructing G3 splines can be presented independently of geometric shape parameters＇ values. It makes the equation＇s solving easier. It is also known that the general form of the G3spline basis function is given in the first time. Its geometric construction method is presented.
文摘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.
文摘Continuously differentiable radial basis functions (C<sup>∞</sup>-RBFs), while being theoretically exponentially convergent are considered impractical computationally because the coefficient matrices are full and can become very ill- conditioned. Similarly, the Hilbert and Vandermonde have full matrices and become ill-conditioned. The difference between a coefficient matrix generated by C<sup>∞</sup>-RBFs for partial differential or integral equations and Hilbert and Vandermonde systems is that C<sup>∞</sup>-RBFs are very sensitive to small changes in the adjustable parameters. These parameters affect the condition number and solution accuracy. The error terrain has many local and global maxima and minima. To find stable and accurate numerical solutions for full linear equation systems, this study proposes a hybrid combination of block Gaussian elimination (BGE) combined with arbitrary precision arithmetic (APA) to minimize the accumulation of rounding errors. In the future, this algorithm can execute faster using preconditioners and implemented on massively parallel computers.
基金co-supported by the ‘‘111" Project of China (No. B17037)the National Natural Science Foundation of China (No. 11772265)
文摘The Radial Basis Function（RBF） method with data reduction is an effective way to perform mesh deformation. However, for large deformations on meshes of complex aerodynamic configurations, the efficiency of the RBF mesh deformation method still needs to be further improved to fulfill the demand of practical application. To achieve this goal, a multistep RBF method based on a multilevel subspace RBF algorithm is presented to further improve the efficiency of the mesh deformation method in this research. A whole deformation is divided into a series of steps, and the supporting radius is adjusted in accordance with the maximal displacement error. Furthermore, parallel computing is applied to the interpolation to enhance the efficiency. Typical deformation problems of the NASA Common Research Model（CRM） configuration, the DLR-F6 wing-body-nacellepylon configuration, and the DLR-F11 high-lift configuration are tested to verify the feasibility of this method. Test results show that the presented multistep RBF mesh deformation method is efficient and robust in dealing with large deformation problems over complex geometries.
文摘The present work describes the application of the method of fundamental solutions （MFS） along with the analog equation method （AEM） and radial basis function （RBF） approximation for solving the 2D isotropic and anisotropic Helmholtz problems with different wave numbers. The AEM is used to convert the original governing equation into the classical Poisson＇s equation, and the MFS and RBF approximations are used to derive the homogeneous and particular solutions, respectively. Finally, the satisfaction of the solution consisting of the homogeneous and particular parts to the related governing equation and boundary conditions can produce a system of linear equations, which can be solved with the singular value decomposition （SVD） technique. In the computation, such crucial factors related to the MFS-RBF as the location of the virtual boundary, the differential and integrating strategies, and the variation of shape parameters in multi-quadric （MQ） are fully analyzed to provide useful reference.
基金The project supported by the National Natural Science Foundation of China (10172052)
文摘Based on our previous study,the accuracy of derivatives of interpolating functions are usually very poor near the boundary of domain when Compactly Supported Radial Basis Functions (CSRBFs)are used,so that it could result in significant error in solving partial differential equations with Neumann boundary conditions.To overcome this drawback,the Consistent Compactly Supported Radial Basis Functions(CCSRBFs)are developed,which satisfy the predetermined consistency con- ditions.Meshless method based on point collocation with CCSRBFs is developed for solving partial differential equations.Numerical studies show that the proposed method improves the accuracy of approximation significantly.
基金This project was supported in part by the Science Foundation of Shanxi Province (2003F028)China Postdoctoral Science Foundation (20060390318).
文摘The Radial Basis Functions Neural Network （RBFNN） is used to establish the model of a response system through the input and output data of the system. The synchronization between a drive system and the response system can be implemented by employing the RBFNN model and state feedback control. In this case, the exact mathematical model, which is the precondition for the conventional method, is unnecessary for implementing synchronization. The effect of the model error is investigated and a corresponding theorem is developed. The effect of the parameter perturbations and the measurement noise is investigated through simulations. The simulation results under different conditions show the effectiveness of the method.
文摘A boundary integral method with radial basis function approximation is proposed for numerically solving an important class of boundary value problems governed by a system of thermoelastostatic equations with variable coe?cients. The equations describe the thermoelastic behaviors of nonhomogeneous anisotropic materials with properties that vary smoothly from point to point in space. No restriction is imposed on the spatial variations of the thermoelastic coe?cients as long as all the requirements of the laws of physics are satis?ed. To check the validity and accuracy of the proposed numerical method, some speci?c test problems with known solutions are solved.
