1,001 research outputs found
Poisson inverse problems
In this paper we focus on nonparametric estimators in inverse problems for
Poisson processes involving the use of wavelet decompositions. Adopting an
adaptive wavelet Galerkin discretization, we find that our method combines the
well-known theoretical advantages of wavelet--vaguelette decompositions for
inverse problems in terms of optimally adapting to the unknown smoothness of
the solution, together with the remarkably simple closed-form expressions of
Galerkin inversion methods. Adapting the results of Barron and Sheu [Ann.
Statist. 19 (1991) 1347--1369] to the context of log-intensity functions
approximated by wavelet series with the use of the Kullback--Leibler distance
between two point processes, we also present an asymptotic analysis of
convergence rates that justifies our approach. In order to shed some light on
the theoretical results obtained and to examine the accuracy of our estimates
in finite samples, we illustrate our method by the analysis of some simulated
examples.Comment: Published at http://dx.doi.org/10.1214/009053606000000687 in the
Annals of Statistics (http://www.imstat.org/aos/) by the Institute of
Mathematical Statistics (http://www.imstat.org
Post-Processing Techniques and Wavelet Applications for Hammerstein Integral Equations
This dissertation is focused on the varieties of numerical solutions of nonlinear Hammerstein integral equations. In the first part of this dissertation, several acceleration techniques for post-processed solutions of the Hammerstein equation are discussed. The post-processing techniques are implemented based on interpolation and extrapolation. In this connection, we generalize the results in [29] and [28] to nonlinear integral equations of the Hammerstein type. Post-processed collocation solutions are shown to exhibit better accuracy. Moreover, an extrapolation technique for the Galerkin solution of Hammerstein equation is also obtained. This result appears new even in the setting of the linear Fredholm equation.
In the second half of this dissertation, the wavelet-collocation technique of solving nonlinear Hammerstein integral equation is discussed. The main objective is to establish a fast wavelet-collocation method for Hammerstein equation by using a \u27linearization\u27 technique. The sparsity in the Jacobian matrix takes place in the fast wavelet-collocation method for Hammerstein equation with smooth as well as weakly singular kernels. A fast algorithm is based upon the block truncation strategy which was recently proposed in [10]. A multilevel augmentation method for the linearized Hammerstein equation is subsequently proposed which further accelerates the solution process while maintaining the order of convergence. Numerical examples are given throughout this dissertation
On adaptive wavelet boundary element methods
The present article is concerned with the numerical solution of boundary integral equa- tions by an adaptive wavelet boundary element method. This method approximates the solution with a computational complexity that is proportional to the solutionâs best N -term approximation. The focus of this article is on algorithmic issues which includes the crucial building blocks and details about the efficient implementation. By numerical examples for the Laplace equation and the Helmholtz equation, solved for different geometries and right-hand sides, we validate the feasibility and efficiency of the adaptive wavelet boundary element method
Fast wavelet collocation methods for second kind integral equations on polygons
In this thesis we develop fast wavelet collocation methods for integral equations of the second kind with weakly singular kernels over polygons . For this purpose, we construct multiscale wavelet functions and collocation functionals having vanishing moments. Moreover, we propose several truncation strategies, which lead to fast algorithms, for the coefficient matrix of the corresponding discrete system. Critical issues for numerical implementation of such methods are considered, such as choices of practical truncation strategies, numerical integration of weakly singular integrals, error controls of numerical quadrature and numerical solutions of resulting compressed linear systems. Numerical experiments are given to demonstrate proposed ideas and methods. Finally, parallel computing using developed methods is investigated.;That this work received partial support from the US NSF grant EPSCoR-0132740
An optimal adaptive wavelet method for First Order System Least Squares
In this paper, it is shown that any well-posed 2nd order PDE can be
reformulated as a well-posed first order least squares system. This system will
be solved by an adaptive wavelet solver in optimal computational complexity.
The applications that are considered are second order elliptic PDEs with
general inhomogeneous boundary conditions, and the stationary Navier-Stokes
equations.Comment: 40 page
- âŠ