相关论文: Rao-Blackwellised Interacting Markov Chain Monte C…
Monte Carlo simulations of systems of particles such as hard spheres or soft spheres with singular kernels can display around a phase transition prohibitively long convergence times when using traditional Hasting-Metropolis reversible…
A new method, based on the simulated annealing algorithm and aimed at the inverse problem in the analysis of intergalactic (interstellar) complex spectra of hydrogen and metal lines, is presented. We consider the process of line formation…
We consider the simulation of Bayesian statistical inverse problems governed by large-scale linear and nonlinear partial differential equations (PDEs). Markov chain Monte Carlo (MCMC) algorithms are standard techniques to solve such…
The inverse problem for a disordered system involves determining the interparticle interaction parameters consistent with a given set of experimental data. Recently, Rutledge has shown (Phys. Rev. E63, 021111 (2001)) that such problems can…
Classic inversion methods adjust a model with a predefined number of parameters to the observed data. With transdimensional inversion algorithms such as the reversible-jump Markov Chain Monte Carlo (rjMCMC), it is possible to vary this…
Typical geophysical inversion problems are ill-posed, non-linear and non-unique. Sometimes the problem is trans-dimensional, where the number of unknown parameters is one of the unknowns, which makes the inverse problem even more…
The electromigration (EM) of metallic lines is studied in terms of competition between two percolative processes taking place in a random resistor network. The effects associated with the transport of mass and with the consequent growth of…
An example of full solution of the inverse scattering problem on the half line is presented. For this purpose, a simple analytically solvable model system (Morse potential) is used, which is expected to be a reasonable approximation to a…
We develop a multilevel Monte Carlo (MLMC)-FEM algorithm for linear, elliptic diffusion problems in polytopal domain $\mathcal D\subset \mathbb R^d$, with Besov-tree random coefficients. This is to say that the logarithms of the diffusion…
Near-infrared spectroscopy (NIRS) including diffuse optical tomography is an imaging modality which makes use of diffuse light propagation in random media. When optical properties of a random medium is investigated from boundary…
Monte Carlo (MC) sampling methods are widely applied in Bayesian inference, system simulation and optimization problems. The Markov Chain Monte Carlo (MCMC) algorithms are a well-known class of MC methods which generate a Markov chain with…
This paper is concerned with inverse acoustic scattering problem of inferring the position and shape of a sound-soft obstacle from phaseless far-field data. We propose the Bayesian approach to recover sound-soft disks, line cracks and…
Monte Carlo methods are widely used to estimate observables in many-body quantum systems. However, conventional sampling schemes often require a large number of samples to achieve sufficient accuracy. In this work we propose the…
The use of sparse precision (inverse covariance) matrices has become popular because they allow for efficient algorithms for joint inference in high-dimensional models. Many applications require the computation of certain elements of the…
We apply a Bayesian data analysis scheme known as the Markov Chain Monte Carlo (MCMC) to the tomographic reconstruction of quantum states. This method yields a vector, known as the Markov chain, which contains the full statistical…
A new approach is proposed to the solution of the quantum mechanical inverse scattering problem at fixed energy. The method relates the fixed energy phase shifts to those arising in an auxiliary Sturm-Liouville problem via the interpolation…
Scattering theory has had a major roll in twentieth century mathematical physics. Mathematical modeling and algorithms of direct,- and inverse electromagnetic scattering formulation due to biological tissues are investigated. The algorithms…
This paper is concerned with uniqueness for reconstructing a periodic inhomogeneous medium covered on a perfectly conducting plate. We deal with the problem in the frame of time-harmonic Maxwell systems without TE or TM polarization. An…
Smoothing in state-space models amounts to computing the conditional distribution of the latent state trajectory, given observations, or expectations of functionals of the state trajectory with respect to this distributions. For models that…
The shell model Monte Carlo (SMMC) method is a powerful technique for calculating the statistical and collective properties of nuclei in the presence of correlations in model spaces that are many orders of magnitude larger than those that…