3,114 research outputs found
On the convergence of the rescaled localized radial basis function method
The rescaled localized RBF method was introduced in Deparis, Forti, and Quarteroni (2014) for scattered data interpolation. It is a rational approximation method based on interpolation with compactly supported radial basis functions. It requires the solution of two linear systems with the same sparse matrix, which has a small condition number, due to the scaling of the basis function. Hence, it can be computed using an unpreconditioned conjugate gradient method in linear time. Numerical evidence provided in Deparis, Forti, and Quarteroni (2014) shows that the method produces good approximations for many examples but no theoretical results were provided. In this paper, we discuss the convergence of the rescaled localized RBF method in the case of quasi-uniform data and stationary scaling. As the method is not only interpolatory but also reproduces constants exactly, linear convergence is expected. We can show this linear convergence up to a certain conjecture
Shearlets and Optimally Sparse Approximations
Multivariate functions are typically governed by anisotropic features such as
edges in images or shock fronts in solutions of transport-dominated equations.
One major goal both for the purpose of compression as well as for an efficient
analysis is the provision of optimally sparse approximations of such functions.
Recently, cartoon-like images were introduced in 2D and 3D as a suitable model
class, and approximation properties were measured by considering the decay rate
of the error of the best -term approximation. Shearlet systems are to
date the only representation system, which provide optimally sparse
approximations of this model class in 2D as well as 3D. Even more, in contrast
to all other directional representation systems, a theory for compactly
supported shearlet frames was derived which moreover also satisfy this
optimality benchmark. This chapter shall serve as an introduction to and a
survey about sparse approximations of cartoon-like images by band-limited and
also compactly supported shearlet frames as well as a reference for the
state-of-the-art of this research field.Comment: in "Shearlets: Multiscale Analysis for Multivariate Data",
Birkh\"auser-Springe
A unified wavelet-based modelling framework for non-linear system identification: the WANARX model structure
A new unified modelling framework based on the superposition of additive submodels, functional components, and
wavelet decompositions is proposed for non-linear system identification. A non-linear model, which is often represented
using a multivariate non-linear function, is initially decomposed into a number of functional components via the wellknown
analysis of variance (ANOVA) expression, which can be viewed as a special form of the NARX (non-linear
autoregressive with exogenous inputs) model for representing dynamic input–output systems. By expanding each functional
component using wavelet decompositions including the regular lattice frame decomposition, wavelet series and
multiresolution wavelet decompositions, the multivariate non-linear model can then be converted into a linear-in-theparameters
problem, which can be solved using least-squares type methods. An efficient model structure determination
approach based upon a forward orthogonal least squares (OLS) algorithm, which involves a stepwise orthogonalization
of the regressors and a forward selection of the relevant model terms based on the error reduction ratio (ERR), is
employed to solve the linear-in-the-parameters problem in the present study. The new modelling structure is referred to
as a wavelet-based ANOVA decomposition of the NARX model or simply WANARX model, and can be applied to
represent high-order and high dimensional non-linear systems
Enhancing SPH using moving least-squares and radial basis functions
In this paper we consider two sources of enhancement for the meshfree
Lagrangian particle method smoothed particle hydrodynamics (SPH) by improving
the accuracy of the particle approximation. Namely, we will consider shape
functions constructed using: moving least-squares approximation (MLS); radial
basis functions (RBF). Using MLS approximation is appealing because polynomial
consistency of the particle approximation can be enforced. RBFs further appeal
as they allow one to dispense with the smoothing-length -- the parameter in the
SPH method which governs the number of particles within the support of the
shape function. Currently, only ad hoc methods for choosing the
smoothing-length exist. We ensure that any enhancement retains the conservative
and meshfree nature of SPH. In doing so, we derive a new set of
variationally-consistent hydrodynamic equations. Finally, we demonstrate the
performance of the new equations on the Sod shock tube problem.Comment: 10 pages, 3 figures, In Proc. A4A5, Chester UK, Jul. 18-22 200
A RBF partition of unity collocation method based on finite difference for initial-boundary value problems
Meshfree radial basis function (RBF) methods are popular tools used to
numerically solve partial differential equations (PDEs). They take advantage of
being flexible with respect to geometry, easy to implement in higher
dimensions, and can also provide high order convergence. Since one of the main
disadvantages of global RBF-based methods is generally the computational cost
associated with the solution of large linear systems, in this paper we focus on
a localizing RBF partition of unity method (RBF-PUM) based on a finite
difference (FD) scheme. Specifically, we propose a new RBF-PUM-FD collocation
method, which can successfully be applied to solve time-dependent PDEs. This
approach allows to significantly decrease ill-conditioning of traditional
RBF-based methods. Moreover, the RBF-PUM-FD scheme results in a sparse matrix
system, reducing the computational effort but maintaining at the same time a
high level of accuracy. Numerical experiments show performances of our
collocation scheme on two benchmark problems, involving unsteady
convection-diffusion and pseudo-parabolic equations
- …