1,674 research outputs found
Wavelet-Fourier CORSING techniques for multi-dimensional advection-diffusion-reaction equations
We present and analyze a novel wavelet-Fourier technique for the numerical
treatment of multidimensional advection-diffusion-reaction equations based on
the CORSING (COmpRessed SolvING) paradigm. Combining the Petrov-Galerkin
technique with the compressed sensing approach, the proposed method is able to
approximate the largest coefficients of the solution with respect to a
biorthogonal wavelet basis. Namely, we assemble a compressed discretization
based on randomized subsampling of the Fourier test space and we employ sparse
recovery techniques to approximate the solution to the PDE. In this paper, we
provide the first rigorous recovery error bounds and effective recipes for the
implementation of the CORSING technique in the multi-dimensional setting. Our
theoretical analysis relies on new estimates for the local a-coherence, which
measures interferences between wavelet and Fourier basis functions with respect
to the metric induced by the PDE operator. The stability and robustness of the
proposed scheme is shown by numerical illustrations in the one-, two-, and
three-dimensional case
Non-equispaced B-spline wavelets
This paper has three main contributions. The first is the construction of
wavelet transforms from B-spline scaling functions defined on a grid of
non-equispaced knots. The new construction extends the equispaced,
biorthogonal, compactly supported Cohen-Daubechies-Feauveau wavelets. The new
construction is based on the factorisation of wavelet transforms into lifting
steps. The second and third contributions are new insights on how to use these
and other wavelets in statistical applications. The second contribution is
related to the bias of a wavelet representation. It is investigated how the
fine scaling coefficients should be derived from the observations. In the
context of equispaced data, it is common practice to simply take the
observations as fine scale coefficients. It is argued in this paper that this
is not acceptable for non-interpolating wavelets on non-equidistant data.
Finally, the third contribution is the study of the variance in a
non-orthogonal wavelet transform in a new framework, replacing the numerical
condition as a measure for non-orthogonality. By controlling the variances of
the reconstruction from the wavelet coefficients, the new framework allows us
to design wavelet transforms on irregular point sets with a focus on their use
for smoothing or other applications in statistics.Comment: 42 pages, 2 figure
Near-best univariate spline discrete quasi-interpolants on non-uniform partitions
International audienceUnivariate spline discrete quasi-interpolants (abbr. dQIs) are approximation operators using B-spline expansions with coefficients which are linear combinations of discrete values of the function to be approximated. When working with nonuniform partitions, the main challenge is to find dQIs which have both good approximation orders and bounded uniform norms independent of the given partition. Near-best dQIs are obtained by minimizing an upper bound of the infinite norm of dQIs depending on a certain number of free parameters, thus reducing this norm. This paper is devoted to the study of some families of near-best dQIs of approximation order 2
Adaptive transient solution of nonuniform multiconductor transmission lines using wavelets
Abstract—This paper presents a highly adaptive algorithm for the transient simulation of nonuniform interconnects loaded with arbitrary nonlinear and dynamic terminations. The discretization of the governing equations is obtained through a weak formula-tion using biorthogonal wavelet bases as trial and test functions. It is shown how the multiresolution properties of wavelets lead to very sparse approximations of the voltages and currents in typical transient analyzes. A simple yet effective time–space adaptive al-gorithm capable of selecting the minimal number of unknowns at each time iteration is described. Numerical results show the high degree of adaptivity of the proposed scheme. Index Terms—Electromagnetic (EM) transient analysis, multi-conductor transmission lines (TLs), wavelet transforms. I
Fast multi-dimensional scattered data approximation with Neumann boundary conditions
An important problem in applications is the approximation of a function
from a finite set of randomly scattered data . A common and powerful
approach is to construct a trigonometric least squares approximation based on
the set of exponentials . This leads to fast numerical
algorithms, but suffers from disturbing boundary effects due to the underlying
periodicity assumption on the data, an assumption that is rarely satisfied in
practice. To overcome this drawback we impose Neumann boundary conditions on
the data. This implies the use of cosine polynomials as basis
functions. We show that scattered data approximation using cosine polynomials
leads to a least squares problem involving certain Toeplitz+Hankel matrices. We
derive estimates on the condition number of these matrices. Unlike other
Toeplitz+Hankel matrices, the Toeplitz+Hankel matrices arising in our context
cannot be diagonalized by the discrete cosine transform, but they still allow a
fast matrix-vector multiplication via DCT which gives rise to fast conjugate
gradient type algorithms. We show how the results can be generalized to higher
dimensions. Finally we demonstrate the performance of the proposed method by
applying it to a two-dimensional geophysical scattered data problem
- …