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.展开更多
The shallow-water temperature profile is typically parameterized using a few empirical orthogonal function(EOF)coefficients.However,when the experimental area is poorly known or highly variable,the adaptability of the...The shallow-water temperature profile is typically parameterized using a few empirical orthogonal function(EOF)coefficients.However,when the experimental area is poorly known or highly variable,the adaptability of the EOFs will be significantly reduced.In this study,a new set of basis functions,generated by combining the internal-wave eigenmodes with the average temperature gradient,is developed for characterizing the temperature perturbations.Temperature profiles recorded by a thermistor chain in the South China Sea in 2015 are processed and analyzed.Compared to the EOFs,the new set of basis functions has higher reconstruction accuracy and adaptability;it is also more stable in ocean regions that have internal waves.展开更多
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.展开更多
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 crack fault is one of the most common faults in the rotor system,and researchers have paid close attention to its fault diagnosis.However,most studies focus on discussing the dynamic response characteristics cause...The crack fault is one of the most common faults in the rotor system,and researchers have paid close attention to its fault diagnosis.However,most studies focus on discussing the dynamic response characteristics caused by the crack rather than estimating the crack depth and position based on the obtained vibration signals.In this paper,a novel crack fault diagnosis and location method for a dual-disk hollow shaft rotor system based on the Radial basis function(RBF)network and Pattern recognition neural network(PRNN)is presented.Firstly,a rotor system model with a breathing crack suitable for a short-thick hollow shaft rotor is established based on the finite element method,where the crack's periodic opening and closing pattern and different degrees of crack depth are considered.Then,the dynamic response is obtained by the harmonic balance method.By adjusting the crack parameters,the dynamic characteristics related to the crack depth and position are analyzed through the amplitude-frequency responses and waterfall plots.The analysis results show that the first critical speed,first subcritical speed,first critical speed amplitude,and super-harmonic resonance peak at the first subcritical speed can be utilized for the crack fault diagnosis.Based on this,the RBF network and PRNN are adopted to determine the depth and approximate location of the crack respectively by taking the above dynamic characteristics as input.Test results show that the proposed method has high fault diagnosis accuracy.This research proposes a crack detection method adequate for the hollow shaft rotor system,where the crack depth and position are both unknown.展开更多
ions in the ECG signal.The cardiologist and medical specialistfind numerous difficulties in the process of traditional approaches.The specified restrictions are eliminated in the proposed classifier.The fundamental ai...ions in the ECG signal.The cardiologist and medical specialistfind numerous difficulties in the process of traditional approaches.The specified restrictions are eliminated in the proposed classifier.The fundamental aim of this work is tofind the R-R interval.To analyze the blockage,different approaches are implemented,which make the computation as facile with high accuracy.The information are recovered from the MIT-BIH dataset.The retrieved data contain normal and pathological ECG signals.To obtain a noiseless signal,Gaborfilter is employed and to compute the amplitude of the signal,DCT-DOST(Discrete cosine based Discrete orthogonal stock well transform)is implemented.The amplitude is computed to detect the cardiac abnormality.The R peak of the underlying ECG signal is noted and the segment length of the ECG cycle is identified.The Genetic algorithm(GA)retrieves the primary highlights and the classifier integrates the data with the chosen attributes to optimize the identification.In addition,the GA helps in performing hereditary calculations to reduce the problem of multi-target enhancement.Finally,the RBFNN(Radial basis function neural network)is applied,which diminishes the local minima present in the signal.It shows enhancement in characterizing the ordinary and anomalous ECG signals.展开更多
Radial Basis Function Neural Network(RBFNN)ensembles have long suffered from non-efficient training,where incorrect parameter settings can be computationally disastrous.This paper examines different evolutionary algor...Radial Basis Function Neural Network(RBFNN)ensembles have long suffered from non-efficient training,where incorrect parameter settings can be computationally disastrous.This paper examines different evolutionary algorithms for training the Symbolic Radial Basis Function Neural Network(SRBFNN)through the behavior’s integration of satisfiability programming.Inspired by evolutionary algorithms,which can iteratively find the nearoptimal solution,different Evolutionary Algorithms(EAs)were designed to optimize the producer output weight of the SRBFNN that corresponds to the embedded logic programming 2Satisfiability representation(SRBFNN-2SAT).The SRBFNN’s objective function that corresponds to Satisfiability logic programming can be minimized by different algorithms,including Genetic Algorithm(GA),Evolution Strategy Algorithm(ES),Differential Evolution Algorithm(DE),and Evolutionary Programming Algorithm(EP).Each of these methods is presented in the steps in the flowchart form which can be used for its straightforward implementation in any programming language.With the use of SRBFNN-2SAT,a training method based on these algorithms has been presented,then training has been compared among algorithms,which were applied in Microsoft Visual C++software using multiple metrics of performance,including Mean Absolute Relative Error(MARE),Root Mean Square Error(RMSE),Mean Absolute Percentage Error(MAPE),Mean Bias Error(MBE),Systematic Error(SD),Schwarz Bayesian Criterion(SBC),and Central Process Unit time(CPU time).Based on the results,the EP algorithm achieved a higher training rate and simple structure compared with the rest of the algorithms.It has been confirmed that the EP algorithm is quite effective in training and obtaining the best output weight,accompanied by the slightest iteration error,which minimizes the objective function of SRBFNN-2SAT.展开更多
Radial Basis Function methods for scattered data interpolation and for the numerical solution of PDEs were originally implemented in a global manner. Subsequently, it was realized that the methods could be implemented...Radial Basis Function methods for scattered data interpolation and for the numerical solution of PDEs were originally implemented in a global manner. Subsequently, it was realized that the methods could be implemented more efficiently in a local manner and that the local approaches could match or even surpass the accuracy of the global implementations. In this work, three localization approaches are compared: a local RBF method, a partition of unity method, and a recently introduced modified partition of unity method. A simple shape parameter selection method is introduced and the application of artificial viscosity to stabilize each of the local methods when approximating time-dependent PDEs is reviewed. Additionally, a new type of quasi-random center is introduced which may be better choices than other quasi-random points that are commonly used with RBF methods. All the results within the manuscript are reproducible as they are included as examples in the freely available Python Radial Basis Function Toolbox.展开更多
This paper concerns the implementation of the orthogonal polynomials using the Galerkin method for solving Volterra integro-differential and Fredholm integro-differential equations. The constructed orthogonal polynomi...This paper concerns the implementation of the orthogonal polynomials using the Galerkin method for solving Volterra integro-differential and Fredholm integro-differential equations. The constructed orthogonal polynomials are used as basis functions in the assumed solution employed. Numerical examples for some selected problems are provided and the results obtained show that the Galerkin method with orthogonal polynomials as basis functions performed creditably well in terms of absolute errors obtained.展开更多
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.展开更多
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.展开更多
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.展开更多
We use Radial Basis Functions (RBFs) to reconstruct smooth surfaces from 3D scattered data. An object's surface is defined implicitly as the zero set of an RBF fitted to the given surface data. We propose improveme...We use Radial Basis Functions (RBFs) to reconstruct smooth surfaces from 3D scattered data. An object's surface is defined implicitly as the zero set of an RBF fitted to the given surface data. We propose improvements on the methods of surface reconstruction with radial basis functions. A sparse approximation set of scattered data is constructed by reducing the number of interpolating points on the surface. We present an adaptive method for finding the off-surface normal points. The order of the equation decreases greatly as the number of the off-surface constraints reduces gradually. Experimental results are provided to illustrate that the proposed method is robust and may draw beautiful graphics.展开更多
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.展开更多
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.展开更多
A new method based on angular momentum theory was proposed to construct the basis functions of the irreducible representations(IRs) of point groups. The transformation coefficients, i. e., coefficients S, are the com...A new method based on angular momentum theory was proposed to construct the basis functions of the irreducible representations(IRs) of point groups. The transformation coefficients, i. e., coefficients S, are the components of the eigenvectors of some Hermitian matrices, and can be made as real numbers for all pure rotation point groups. The general formula for coefficient S was deduced, and applied to constructing the basis functions of single-valued irreducible representations of icosahedral group from the spherical harmonics with angular momentum j≤7.展开更多
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.展开更多
基金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.
基金The Natural Science Foundation of Shandong Province of China under contract Nos ZR2022MA051 and ZR2020MA090the Fund of China Postdoctoral Science Foundation under contract No.2020M670891+1 种基金the Shandong University of Science and Technology Research Fund under contract No.2019TDJH103the Talent Introduction Plan for Youth Innovation Team in Universities of Shandong Province(Innovation Team of Satellite Positioning and Navigation).
文摘The shallow-water temperature profile is typically parameterized using a few empirical orthogonal function(EOF)coefficients.However,when the experimental area is poorly known or highly variable,the adaptability of the EOFs will be significantly reduced.In this study,a new set of basis functions,generated by combining the internal-wave eigenmodes with the average temperature gradient,is developed for characterizing the temperature perturbations.Temperature profiles recorded by a thermistor chain in the South China Sea in 2015 are processed and analyzed.Compared to the EOFs,the new set of basis functions has higher reconstruction accuracy and adaptability;it is also more stable in ocean regions that have internal waves.
基金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.
文摘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.
基金Supported by National Natural Science Foundation of China (Grant No.11972129)National Science and Technology Major Project of China (Grant No.2017-IV-0008-0045)+1 种基金Heilongjiang Provincial Natural Science Foundation (Grant No.YQ2022A008)the Fundamental Research Funds for the Central Universities。
文摘The crack fault is one of the most common faults in the rotor system,and researchers have paid close attention to its fault diagnosis.However,most studies focus on discussing the dynamic response characteristics caused by the crack rather than estimating the crack depth and position based on the obtained vibration signals.In this paper,a novel crack fault diagnosis and location method for a dual-disk hollow shaft rotor system based on the Radial basis function(RBF)network and Pattern recognition neural network(PRNN)is presented.Firstly,a rotor system model with a breathing crack suitable for a short-thick hollow shaft rotor is established based on the finite element method,where the crack's periodic opening and closing pattern and different degrees of crack depth are considered.Then,the dynamic response is obtained by the harmonic balance method.By adjusting the crack parameters,the dynamic characteristics related to the crack depth and position are analyzed through the amplitude-frequency responses and waterfall plots.The analysis results show that the first critical speed,first subcritical speed,first critical speed amplitude,and super-harmonic resonance peak at the first subcritical speed can be utilized for the crack fault diagnosis.Based on this,the RBF network and PRNN are adopted to determine the depth and approximate location of the crack respectively by taking the above dynamic characteristics as input.Test results show that the proposed method has high fault diagnosis accuracy.This research proposes a crack detection method adequate for the hollow shaft rotor system,where the crack depth and position are both unknown.
文摘ions in the ECG signal.The cardiologist and medical specialistfind numerous difficulties in the process of traditional approaches.The specified restrictions are eliminated in the proposed classifier.The fundamental aim of this work is tofind the R-R interval.To analyze the blockage,different approaches are implemented,which make the computation as facile with high accuracy.The information are recovered from the MIT-BIH dataset.The retrieved data contain normal and pathological ECG signals.To obtain a noiseless signal,Gaborfilter is employed and to compute the amplitude of the signal,DCT-DOST(Discrete cosine based Discrete orthogonal stock well transform)is implemented.The amplitude is computed to detect the cardiac abnormality.The R peak of the underlying ECG signal is noted and the segment length of the ECG cycle is identified.The Genetic algorithm(GA)retrieves the primary highlights and the classifier integrates the data with the chosen attributes to optimize the identification.In addition,the GA helps in performing hereditary calculations to reduce the problem of multi-target enhancement.Finally,the RBFNN(Radial basis function neural network)is applied,which diminishes the local minima present in the signal.It shows enhancement in characterizing the ordinary and anomalous ECG signals.
基金This work is supported by Ministry of Higher Education(MOHE)through Fundamental Research Grant Scheme(FRGS)(FRGS/1/2020/STG06/UTHM/03/7).
文摘Radial Basis Function Neural Network(RBFNN)ensembles have long suffered from non-efficient training,where incorrect parameter settings can be computationally disastrous.This paper examines different evolutionary algorithms for training the Symbolic Radial Basis Function Neural Network(SRBFNN)through the behavior’s integration of satisfiability programming.Inspired by evolutionary algorithms,which can iteratively find the nearoptimal solution,different Evolutionary Algorithms(EAs)were designed to optimize the producer output weight of the SRBFNN that corresponds to the embedded logic programming 2Satisfiability representation(SRBFNN-2SAT).The SRBFNN’s objective function that corresponds to Satisfiability logic programming can be minimized by different algorithms,including Genetic Algorithm(GA),Evolution Strategy Algorithm(ES),Differential Evolution Algorithm(DE),and Evolutionary Programming Algorithm(EP).Each of these methods is presented in the steps in the flowchart form which can be used for its straightforward implementation in any programming language.With the use of SRBFNN-2SAT,a training method based on these algorithms has been presented,then training has been compared among algorithms,which were applied in Microsoft Visual C++software using multiple metrics of performance,including Mean Absolute Relative Error(MARE),Root Mean Square Error(RMSE),Mean Absolute Percentage Error(MAPE),Mean Bias Error(MBE),Systematic Error(SD),Schwarz Bayesian Criterion(SBC),and Central Process Unit time(CPU time).Based on the results,the EP algorithm achieved a higher training rate and simple structure compared with the rest of the algorithms.It has been confirmed that the EP algorithm is quite effective in training and obtaining the best output weight,accompanied by the slightest iteration error,which minimizes the objective function of SRBFNN-2SAT.
文摘Radial Basis Function methods for scattered data interpolation and for the numerical solution of PDEs were originally implemented in a global manner. Subsequently, it was realized that the methods could be implemented more efficiently in a local manner and that the local approaches could match or even surpass the accuracy of the global implementations. In this work, three localization approaches are compared: a local RBF method, a partition of unity method, and a recently introduced modified partition of unity method. A simple shape parameter selection method is introduced and the application of artificial viscosity to stabilize each of the local methods when approximating time-dependent PDEs is reviewed. Additionally, a new type of quasi-random center is introduced which may be better choices than other quasi-random points that are commonly used with RBF methods. All the results within the manuscript are reproducible as they are included as examples in the freely available Python Radial Basis Function Toolbox.
文摘This paper concerns the implementation of the orthogonal polynomials using the Galerkin method for solving Volterra integro-differential and Fredholm integro-differential equations. The constructed orthogonal polynomials are used as basis functions in the assumed solution employed. Numerical examples for some selected problems are provided and the results obtained show that the Galerkin method with orthogonal polynomials as basis functions performed creditably well in terms of absolute errors obtained.
文摘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.
基金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.
文摘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.
文摘We use Radial Basis Functions (RBFs) to reconstruct smooth surfaces from 3D scattered data. An object's surface is defined implicitly as the zero set of an RBF fitted to the given surface data. We propose improvements on the methods of surface reconstruction with radial basis functions. A sparse approximation set of scattered data is constructed by reducing the number of interpolating points on the surface. We present an adaptive method for finding the off-surface normal points. The order of the equation decreases greatly as the number of the off-surface constraints reduces gradually. Experimental results are provided to illustrate that the proposed method is robust and may draw beautiful graphics.
文摘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.
基金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.
文摘A new method based on angular momentum theory was proposed to construct the basis functions of the irreducible representations(IRs) of point groups. The transformation coefficients, i. e., coefficients S, are the components of the eigenvectors of some Hermitian matrices, and can be made as real numbers for all pure rotation point groups. The general formula for coefficient S was deduced, and applied to constructing the basis functions of single-valued irreducible representations of icosahedral group from the spherical harmonics with angular momentum j≤7.
文摘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.