Source author record

Alexei Novikov

Alexei Novikov appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
14topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

Quantitative phase and absorption contrast imaging

We present an algorithm for coherent diffractive imaging with phaseless measurements. It treats the forward model as a combination of coherent and incoherent waves. The algorithm reconstructs absorption and phase contrast that quantifies the attenuation and the refraction of the waves propagating through an object. It requires coherent or partially coherent illuminations, and several detectors to record the intensity of the distorted wave that passes through the object under inspection. The diversity of illuminations, obtained by putting masks between the source and the object, provides enough information for imaging. The computational cost of our algorithm is linear in the number of pixels of the image. Therefore, it is efficient for high-resolution imaging. Our algorithm guarantees exact recovery if the image is sparse for a given basis. Numerical experiments in the setting of phaseless diffraction imaging of sparse objects validate the efficiency and the precision of the suggested algorithm.

preprint2019arXiv

Imaging with highly incomplete and corrupted data

We consider the problem of imaging sparse scenes from a few noisy data using an $l_1$-minimization approach. This problem can be cast as a linear system of the form $A \, ρ=b$, where $A$ is an $N\times K$ measurement matrix. We assume that the dimension of the unknown sparse vector $ρ\in {\mathbb{C}}^K$ is much larger than the dimension of the data vector $b \in {\mathbb{C}}^N$, i.e, $K \gg N$. We provide a theoretical framework that allows us to examine under what conditions the $\ell_1$-minimization problem admits a solution that is close to the exact one in the presence of noise. Our analysis shows that $l_1$-minimization is not robust for imaging with noisy data when high resolution is required. To improve the performance of $l_1$-minimization we propose to solve instead the augmented linear system $ [A \, | \, C] ρ=b$, where the $N \times Σ$ matrix $C$ is a noise collector. It is constructed so as its column vectors provide a frame on which the noise of the data, a vector of dimension $N$, can be well approximated. Theoretically, the dimension $Σ$ of the noise collector should be $e^N$ which would make its use not practical. However, our numerical results illustrate that robust results in the presence of noise can be obtained with a large enough number of columns $Σ\approx 10 K$.

preprint2019arXiv

The Noise Collector for sparse recovery in high dimensions

The ability to detect sparse signals from noisy high-dimensional data is a top priority in modern science and engineering. A sparse solution of the linear system $A ρ= b_0$ can be found efficiently with an $l_1$-norm minimization approach if the data is noiseless. Detection of the signal's support from data corrupted by noise is still a challenging problem, especially if the level of noise must be estimated. We propose a new efficient approach that does not require any parameter estimation. We introduce the Noise Collector (NC) matrix $C$ and solve an augmented system $A ρ+ C η= b_0 + e$, where $ e$ is the noise. We show that the $l_1$-norm minimal solution of the augmented system has zero false discovery rate for any level of noise and with probability that tends to one as the dimension of $ b_0$ increases to infinity. We also obtain exact support recovery if the noise is not too large, and develop a Fast Noise Collector Algorithm which makes the computational cost of solving the augmented system comparable to that of the original one. Finally, we demonstrate the effectiveness of the method in applications to passive array imaging.

preprint2016arXiv

A fractional kinetic process describing the intermediate time behaviour of cellular flows

This paper studies the intermediate time behaviour of a small random perturbation of a periodic cellular flow. Our main result shows that on time scales shorter than the diffusive time scale, the limiting behaviour of trajectories that start close enough to cell boundaries is a fractional kinetic process: A Brownian motion time changed by the local time of an independent Brownian motion. Our proof uses the Freidlin-Wentzell framework, and the key step is to establish an analogous averaging principle on shorter time scales. As a consequence of our main theorem, we obtain a homogenization result for the associated advection-diffusion equation. We show that on intermediate time scales the effective equation is a fractional time PDE that arises in modelling anomalous diffusion.

preprint2016arXiv

Coherent imaging without phases

In this paper we consider narrow band, active array imaging of weak localized scatterers when only the intensities are recorded at an array with N transducers. We consider that the medium is homogeneous and, hence, wave propagation is fully coherent. This work is an extension of our previous paper, where we showed that using linear combinations of intensity-only measurements imaging of localized scatterers can be carried out efficiently using MUSIC or sparsity promoting optimization. Here we show the same strategy can be accomplished with only 3N-2 illuminations, therefore reducing enormously the data acquisition process. Furthermore, we show that in the paraxial regime one can form the images by using six illuminations only. In particular, this paraxial regime includes Fresnel and Fraunhofer diffraction. The key point of this work is that if one controls the illuminations, imaging with intensity-only can be easily reduced to a imaging with phases and, therefore, one can apply standard imaging techniques. Detailed numerical simulations illustrate the performance of the proposed imaging strategy with and without data noise.

preprint2016arXiv

Multifrequency interferometric imaging with intensity-only measurements

We propose an illumination strategy for interferometric imaging that allows for robust depth recovery from intensity-only measurements. For an array with colocated sources and receivers, we show that all the possible interferometric data for multiple sources, receivers and frequencies can be recovered from intensity-only measurements provided that we have sufficient source location and frequency illumination diversity. There is no need for phase reconstruction in this approach. Using interferometric imaging methods we show that in homogeneous media there is no loss of resolution when imaging with intensities-only. If in these imaging methods we reduce incoherence by restricting the multifrequency interferometric data to nearby array elements and nearby frequencies we obtain robust images in weakly inhomogeneous background media with a somewhat reduced resolution.

preprint2015arXiv

Anomalous diffusion in fast cellular flows at intermediate time scales

It is well known that on long time scales the behaviour of tracer particles diffusing in a cellular flow is effectively that of a Brownian motion. This paper studies the behaviour on "intermediate" time scales before diffusion sets in. Various heuristics suggest that an anomalous diffusive behaviour should be observed. We prove that the variance on intermediate time scales grows like $O(\sqrt{t})$. Hence, on these time scales the effective behaviour can not be purely diffusive, and is consistent with an anomalous diffusive behaviour.

preprint2014arXiv

Illumination strategies for intensity-only imaging

We propose a new strategy for narrow band, active array imaging of localized scat- terers when only the intensities are recorded and measured at the array. We consider a homogeneous medium so that wave propagation is fully coherent. We show that imaging with intensity-only measurements can be carried out using the time reversal operator of the imaging system, which can be obtained from intensity measurements using an appropriate illumination strategy and the polarization identity. Once the time reversal operator has been obtained, we show that the images can be formed using its singular value decomposition (SVD). We use two SVD-based methods to image the scatterers. The proposed approach is simple and efficient. It does not need prior information about the sought image, and guarantees exact recovery in the noise-free case. Furthermore, it is robust with respect to additive noise. Detailed numerical simulations illustrate the performance of the proposed imaging strategy when only the intensities are captured.

preprint2012arXiv

A differential equations approach to $l_1$-minimization with applications to array imaging

We present an ordinary differential equations approach to the analysis of algorithms for constructing $l_1$ minimizing solutions to underdetermined linear systems of full rank. It involves a relaxed minimization problem whose minimum is independent of the relaxation parameter. An advantage of using the ordinary differential equations is that energy methods can be used to prove convergence. The connection to the discrete algorithms is provided by the Crandall-Liggett theory of monotone nonlinear semigroups. We illustrate the effectiveness of the discrete optimization algorithm in some sparse array imaging problems.

preprint2012arXiv

The regularizing effects of resetting in a particle system for the Burgers equation

We study the dissipation mechanism of a stochastic particle system for the Burgers equation. The velocity field of the viscous Burgers and Navier-Stokes equations can be expressed as an expected value of a stochastic process based on noisy particle trajectories [Constantin and Iyer Comm. Pure Appl. Math. 3 (2008) 330-345]. In this paper we study a particle system for the viscous Burgers equations using a Monte-Carlo version of the above; we consider N copies of the above stochastic flow, each driven by independent Wiener processes, and replace the expected value with $\frac{1}{N}$ times the sum over these copies. A similar construction for the Navier-Stokes equations was studied by Mattingly and the first author of this paper [Iyer and Mattingly Nonlinearity 21 (2008) 2537-2553]. Surprisingly, for any finite N, the particle system for the Burgers equations shocks almost surely in finite time. In contrast to the full expected value, the empirical mean $\frac{1}{N}\sum_1^N$ does not regularize the system enough to ensure a time global solution. To avoid these shocks, we consider a resetting procedure, which at first sight should have no regularizing effect at all. However, we prove that this procedure prevents the formation of shocks for any $N\geq2$, and consequently as $N\to\infty$ we get convergence to the solution of the viscous Burgers equation on long time intervals.

preprint2011arXiv

From homogenization to averaging in cellular flows

We consider an elliptic eigenvalue problem with a fast cellular flow of amplitude $A$, in a two-dimensional domain with $L^2$ cells. For fixed $A$, and $L \to \infty$, the problem homogenizes, and has been well studied. Also well studied is the limit when $L$ is fixed, and $A \to \infty$. In this case the solution equilibrates along stream lines. In this paper, we show that if \textit{both} $A \to \infty$ and $L \to \infty$, then a transition between the homogenization and averaging regimes occurs at $A \approx L^4$. When $A\gg L^4$, the principal Dirichlet eigenvalue is approximately constant. On the other hand, when $A\ll L^4$, the principal eigenvalue behaves like ${\bar σ(A)}/L^2$, where $\bar σ(A) \approx \sqrt{A} I$ is the effective diffusion matrix. A similar transition is observed for the solution of the exit time problem. The proof in the homogenization regime involves bounds on the second correctors. Miraculously, if the slow profile is quadratic, these estimates can be obtained using drift independent $L^p \to L^\infty$ estimates for elliptic equations with an incompressible drift. This provides effective sub and super-solutions for our problem.

preprint2010arXiv

Homogenization of the G-equation with incompressible random drift in two dimensions

We study the homogenization limit of solutions to the G-equation with random drift. This Hamilton-Jacobi equation is a model for flame propagation in a turbulent fluid in the regime of thin flames. For a fluid velocity field that is statistically stationary and ergodic, we prove sufficient conditions for homogenization to hold with probability one. These conditions are expressed in terms of travel times for the associated control problem. When the spatial dimension is equal to two and the fluid velocity is divergence-free, we verify that these conditions hold under suitable assumptions about the growth of the random stream function.