BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

George Papanicolaou

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

ResearcherUnclaimed source record
math.NAmath.OCeess.SPmath.PRmath-phmath.APmath.MPInformation TheoryMachine Learningmath.ITphysics.comp-phcond-mat.stat-mecheess.IVnlin.AONumerical Analysisphysics.med-phphysics.opticsphysics.soc-phq-fin.MFq-fin.RMq-fin.STq-fin.TRSocial and Information Networks

Catalog footprint

What is connected

23works
23topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

23 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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Generalized correlation based Imaging for satellites

We consider imaging of fast moving small objects in space, such as low earth orbit satellites. The imaging system consists of ground based, asynchronous sources of radiation and several passive receivers above the dense atmosphere. We use the cross correlation of the received signals to reduce distortions from ambient medium fluctuations. Imaging with correlations also has the advantage of not requiring any knowledge about the probing pulse and depends weakly on the emitter positions. We account for the target's orbital velocity by introducing the necessary Doppler compensation. We show that over limited imaging regions, a constant Doppler factor can be used, resulting in an efficient data structure for the correlations of the recorded signals. We then investigate and analyze different imaging methods using the cross-correlation data structure. Specifically, we show that using a generalized two point migration of the cross correlation data, the top eigenvector of the migrated data matrix provides superior image resolution compared to the usual single-point migration scheme. We carry out a theoretical analysis that illustrates the role of the two point migration methods as well as that of the inverse aperture in improving resolution. Extensive numerical simulations support the theoretical results and assess the scope of the imaging methodology.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

PCA for Implied Volatility Surfaces

Principal component analysis (PCA) is a useful tool when trying to construct factor models from historical asset returns. For the implied volatilities of U.S. equities there is a PCA-based model with a principal eigenportfolio whose return time series lies close to that of an overarching market factor. The authors show that this market factor is the index resulting from the daily compounding of a weighted average of implied-volatility returns, with weights based on the options' open interest (OI) and Vega. The authors also analyze the singular vectors derived from the tensor structure of the implied volatilities of S&P500 constituents, and find evidence indicating that some type of OI and Vega-weighted index should be one of at least two significant factors in this market.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Synthetic aperture imaging and motion estimation using tensor methods

We consider a synthetic aperture imaging configuration, such as synthetic aperture radar (SAR), where we want to first separate reflections from moving targets from those coming from a stationary background, and then to image separately the moving and the stationary reflectors. For this purpose, we introduce a representation of the data as a third order tensor formed from data coming from partially overlapping sub-apertures. We then apply a tensor robust principal component analysis (TRPCA) to the tensor data which separates them into the parts coming from the stationary and moving reflectors. Images are formed with the separated data sets. Our analysis shows a distinctly improved performance of TRPCA, compared to the usual matrix case. In particular, the tensor decomposition can identify motion features that are undetectable when using the conventional motion estimation methods, including matrix RPCA. We illustrate the performance of the method with numerical simulations in the X-band radar regime.

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Consensus Convergence with Stochastic Effects

We consider a stochastic, continuous state and time opinion model where each agent's opinion locally interacts with other agents' opinions in the system, and there is also exogenous randomness. The interaction tends to create clusters of common opinion. By using linear stability analysis of the associated nonlinear Fokker-Planck equation that governs the empirical density of opinions in the limit of infinitely many agents, we can estimate the number of clusters, the time to cluster formation and the critical strength of randomness so as to have cluster formation. We also discuss the cluster dynamics after their formation, the width and the effective diffusivity of the clusters. Finally, the long term behavior of clusters is explored numerically. Extensive numerical simulations confirm our analytical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Limit order trading with a mean reverting reference price

Optimal control models for limit order trading often assume that the underlying asset price is a Brownian motion since they deal with relatively short time scales. The resulting optimal bid and ask limit order prices tend to track the underlying price as one might expect. This is indeed the case with the model of Avellaneda and Stoikov (2008), which has been studied extensively. We consider here this model under the condition when the underlying price is mean reverting. Our main result is that when time is far from the terminal, the optimal price for bid and ask limit orders is constant, which means that it does not track the underlying price. Numerical simulations confirm this behavior. When the underlying price is mean reverting, then for times sufficiently far from terminal, it is more advantageous to focus on the mean price and ignore fluctuations around it. Mean reversion suggests that limit orders will be executed with some regularity, and this is why they are optimal. We also explore intermediate time regimes where limit order prices are influenced by the inventory of outstanding orders. The duration of this intermediate regime depends on the liquidity of the market as measured by specific parameters in the model.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Synthetic Aperture Imaging of Direction and Frequency Dependent Reflectivities

We introduce a synthetic aperture imaging framework that takes into consideration directional dependence of the reflectivity that is to be imaged, as well as its frequency dependence. We use an $\ell_1$ minimization approach that is coordinated with data segmentation so as to fuse information from multiple sub-apertures and frequency sub-bands. We analyze this approach from first principles and assess its performance with numerical simulations in an X-band radar regime.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A risk analysis for a system stabilized by a central agent

We formulate and analyze a multi-agent model for the evolution of individual and systemic risk in which the local agents interact with each other through a central agent who, in turn, is influenced by the mean field of the local agents. The central agent is stabilized by a bistable potential, the only stabilizing force in the system. The local agents derive their stability only from the central agent. In the mean field limit of a large number of local agents we show that the systemic risk decreases when the strength of the interaction of the local agents with the central agent increases. This means that the probability of transition from one of the two stable quasi-equilibria to the other one decreases. We also show that the systemic risk increases when the strength of the interaction of the central agent with the mean field of the local agents increases. Following the financial interpretation of such models and their behavior given in our previous paper (Garnier, Papanicolaou and Yang, SIAM J. Fin. Math. 4, 2013, 151-184), we may interpret the results of this paper in the following way. From the point of view of systemic risk, and while keeping the perceived risk of the local agents approximately constant, it is better to strengthen the interaction of the local agents with the central agent than the other way around.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Array imaging of localized objects in homogeneous and heterogeneous media

We present a comprehensive study of the resolution and stability properties of sparse promoting optimization theories applied to narrow band array imaging of localized scatterers. We consider homogeneous and heterogeneous media, and multiple and single scattering situations. When the media is homogeneous with strong multiple scattering between scatterers, we give a non-iterative formulation to find the locations and reflectivities of the scatterers from a nonlinear inverse problem in two steps, using either single or multiple illuminations. We further introduce an approach that uses the top singular vectors of the response matrix as optimal illuminations, which improves the robustness of sparse promoting optimization with respect to additive noise. When multiple scattering is negligible, the optimization problem becomes linear and can be reduced to a hybrid-$\ell_1$ method when optimal illuminations are used. When the media is random, and the interaction with the unknown inhomogeneities can be primarily modeled by wavefront distortions, we address the statistical stability of these methods. We analyze the fluctuations of the images obtained with the hybrid-$\ell_1$ method, and we show that it is stable with respect to different realizations of the random medium provided the imaging array is large enough. We compare the performance of the hybrid-$\ell_1$ method in random media to the widely used Kirchhoff migration and the multiple signal classification methods.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Imaging strong localized scatterers with sparsity promoting optimization

We study active array imaging of small but strong scatterers in homogeneous media when multiple scattering between them is important. We use the Foldy-Lax equations to model wave propagation with multiple scattering when the scatterers are small relative to the wavelength. In active array imaging we seek to locate the positions and reflectivities of the scatterers, that is, to determine the support of the reflectivity vector and the values of its nonzero elements from echoes recorded on the array. This is a nonlinear inverse problem because of the multiple scattering. We show in this paper how to avoid the nonlinearity and form images non-iteratively through a two-step process which involves $\ell_1$ norm minimization. However, under certain illuminations imaging may be affected by screening, where some scatterers are obscured by multiple scattering. This problem can be mitigated by using multiple and diverse illuminations. In this case, we determine solution vectors that have a common support. The uniqueness and stability of the support of the reflectivity vector obtained with single or multiple illuminations are analyzed, showing that the errors are proportional to the amount of noise in the data with a proportionality factor dependent on the sparsity of the solution and the mutual coherence of the sensing matrix, which is determined by the geometry of the imaging array. Finally, to filter out noise and improve the resolution of the images, we propose an approach that combines optimal illuminations using the singular value decomposition of the response matrix together with sparsity promoting optimization jointly for all illuminations. This work is an extension of our previous paper [5] on imaging using optimization techniques where we now account for multiple scattering effects.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Role of scattering in virtual source array imaging

We consider imaging in a scattering medium where the illumination goes through this medium but there is also an auxiliary, passive receiver array that is near the object to be imaged. Instead of imaging with the source-receiver array on the far side of the object we image with the data of the passive array on the near side of the object. The imaging is done with travel time migration using the cross correlations of the passive array data. We showed in [J. Garnier and G. Papanicolaou, Inverse Problems {28} (2012), 075002] that if (i) the source array is infinite, (ii) the scattering medium is modeled by either an isotropic random medium in the paraxial regime or a randomly layered medium, and (iii) the medium between the auxiliary array and the object to be imaged is homogeneous, then imaging with cross correlations completely eliminates the effects of the random medium. It is as if we imaged with an active array, instead of a passive one, near the object. The purpose of this paper is to analyze the resolution of the image when both the source array and the passive receiver array are finite. We show with a detailed analysis that for isotropic random media in the paraxial regime, imaging not only is not affected by the inhomogeneities but the resolution can in fact be enhanced. This is because the random medium can increase the diversity of the illumination. We also show analytically that this will not happen in a randomly layered medium, and there may be some loss of resolution in this case.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Anomalous Shock Displacement Probabilities for a Perturbed Scalar Conservation Law

We consider an one-dimensional conservation law with random space-time forcing and calculate using large deviations the exponentially small probabilities of anomalous shock profile displacements. Under suitable hypotheses on the spatial support and structure of random forces, we analyze the scaling behavior of the rate function, which is the exponential decay rate of the displacement probabilities. For small displacements we show that the rate function is bounded above and below by the square of the displacement divided by time. For large displacements the corresponding bounds for the rate function are proportional to the displacement. We calculate numerically the rate function under different conditions and show that the theoretical analysis of scaling behavior is confirmed. We also apply a large-deviation-based importance sampling Monte Carlo strategy to estimate the displacement probabilities. We use a biased distribution centered on the forcing that gives the most probable transition path for the anomalous shock profile, which is the minimizer of the rate function. The numerical simulations indicate that this strategy is much more effective and robust than basic Monte Carlo.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Large deviations for a mean field model of systemic risk

We consider a system of diffusion processes that interact through their empirical mean and have a stabilizing force acting on each of them, corresponding to a bistable potential. There are three parameters that characterize the system: the strength of the intrinsic stabilization, the strength of the external random perturbations, and the degree of cooperation or interaction between them. The latter is the rate of mean reversion of each component to the empirical mean of the system. We interpret this model in the context of systemic risk and analyze in detail the effect of cooperation between the components, that is, the rate of mean reversion. We show that in a certain regime of parameters increasing cooperation tends to increase the stability of the individual agents but it also increases the overall or systemic risk. We use the theory of large deviations of diffusions interacting through their mean field.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Motion Estimation and Imaging of Complex Scenes with Synthetic Aperture Radar

We study synthetic aperture radar (SAR) imaging and motion estimation of complex scenes consisting of stationary and moving targets. We use the classic SAR setup with a single antenna emitting signals and receiving the echoes from the scene. The known motion estimation methods for such setups work only in simple cases, with one or a few targets in the same motion. We propose to extend the applicability of these methods to complex scenes, by complementing them with a data pre-processing step intended to separate the echoes from the stationary targets and the moving ones. We present two approaches. The first is an iteration designed to subtract the echoes from the stationary targets one by one. It estimates the location of each stationary target from a preliminary image, and then uses it to define a filter that removes its echo from the data. The second approach is based on the robust principle component analysis (PCA) method. The key observation is that with appropriate pre-processing and windowing, the discrete samples of the stationary target echoes form a low rank matrix, whereas the samples of a few moving target echoes form a high rank sparse matrix. The robust PCA method is designed to separate the low rank from the sparse part, and thus can be used for the SAR data separation. We present a brief analysis of the two methods and explain how they can be combined to improve the data separation for extended and complex imaging scenes. We also assess the performance of the methods with extensive numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Probability of Failure in Hypersonic Engines Using Large Deviations

We consider a reduced order model of an air-breathing hypersonic engine with a time-dependent stochastic inflow that may cause the failure of the engine. The probability of failure is analyzed by the Freidlin-Wentzell theory, the large deviation principle for finite dimensional stochastic differential equations. We compute the asymptotic failure probability by numerically solving the constrained optimization related to the large deviation problem. A large-deviation-based importance sampling suggested by the most probable inflow perturbation is also implemented to compute the probability of failure of the engine. The numerical simulations show that the importance sampling method is much more efficient than the basic Monte Carlo method.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Synthetic Aperture Radar Imaging and Motion Estimation via Robust Principle Component Analysis

We consider the problem of synthetic aperture radar (SAR) imaging and motion estimation of complex scenes. By complex we mean scenes with multiple targets, stationary and in motion. We use the usual setup with one moving antenna emitting and receiving signals. We address two challenges: (1) the detection of moving targets in the complex scene and (2) the separation of the echoes from the stationary targets and those from the moving targets. Such separation allows high resolution imaging of the stationary scene and motion estimation with the echoes from the moving targets alone. We show that the robust principal component analysis (PCA) method which decomposes a matrix in two parts, one low rank and one sparse, can be used for motion detection and data separation. The matrix that is decomposed is the pulse and range compressed SAR data indexed by two discrete time variables: the slow time, which parametrizes the location of the antenna, and the fast time, which parametrizes the echoes received between successive emissions from the antenna. We present an analysis of the rank of the data matrix to motivate the use of the robust PCA method. We also show with numerical simulations that successful data separation with robust PCA requires proper data windowing. Results of motion estimation and imaging with the separated data are presented, as well.

0 citations0 reviewsNo code linkedOpen access