Signal recovery from wavelet transform maxima

Cataloged from PDF version of article.This paper presents an iterative algorithm for signal recovery from discrete-time wavelet transform maxima. The signal recovery algorithm is developed by using the method of projections onto convex sets. Convergence of the algorithm is assured

Sparsity and Incoherence in Compressive Sampling

We consider the problem of reconstructing a sparse signal $x^0\in\R^n$ from a limited number of linear measurements. Given $m$ randomly selected samples of $U x^0$, where $U$ is an orthonormal matrix, we show that $\ell_1$ minimization recovers $x^0$ exactly when the number of measurements exceeds $$m\geq \mathrm{Const}\cdot\mu^2(U)\cdot S\cdot\log n,$$ where $S$ is the number of nonzero components in $x^0$, and $\mu$ is the largest entry in $U$ properly normalized: $\mu(U) = \sqrt{n} \cdot \max_{k,j} |U_{k,j}|$. The smaller $\mu$, the fewer samples needed. The result holds for ``most'' sparse signals $x^0$ supported on a fixed (but arbitrary) set $T$. Given $T$, if the sign of $x^0$ for each nonzero entry on $T$ and the observed values of $Ux^0$ are drawn at random, the signal is recovered with overwhelming probability. Moreover, there is a sense in which this is nearly optimal since any method succeeding with the same probability would require just about this many samples

Multi-Fractal Spectral Analysis of the 1987 Stock Market Crash

The multifractal model of asset returns captures the volatility persistence of many financial time series. Its multifractal spectrum computed from wavelet modulus maxima lines provides the spectrum of irregularities in the distribution of market returns over time and thereby of the kind of uncertainty or randomness in a particular market. Changes in this multifractal spectrum display distinctive patterns around substantial market crashes or drawdowns. In other words, the kinds of singularities and the kinds of irregularity change in a distinct fashion in the periods immediately preceding and following major market drawdowns. This paper focuses on these identifiable multifractal spectral patterns surrounding the stock market crash of 1987. Although we are not able to find a uniquely identifiable irregularity pattern within the same market preceding different crashes at different times, we do find the same uniquely identifiable pattern in various stock markets experiencing the same crash at the same time. Moreover, our results suggest that all such crashes are preceded by a gradual increase in the weighted average of the values of the Lipschitz regularity exponents, under low dispersion of the multifractal spectrum. At a crash, this weighted average irregularity value drops to a much lower value, while the dispersion of the spectrum of Lipschitz exponents jumps up to a much higher level after the crash. Our most striking result, therefore, is that the multifractal spectra of stock market returns are not stationary. Also, while the stock market returns show a global Hurst exponent of slight persistence 0.5Financial Markets, Persistence, Multi-Fractal Spectral Analysis, Wavelets

Extreme Value Analysis of Empirical Frame Coefficients and Implications for Denoising by Soft-Thresholding

Denoising by frame thresholding is one of the most basic and efficient methods for recovering a discrete signal or image from data that are corrupted by additive Gaussian white noise. The basic idea is to select a frame of analyzing elements that separates the data in few large coefficients due to the signal and many small coefficients mainly due to the noise \epsilon_n. Removing all data coefficients being in magnitude below a certain threshold yields a reconstruction of the original signal. In order to properly balance the amount of noise to be removed and the relevant signal features to be kept, a precise understanding of the statistical properties of thresholding is important. For that purpose we derive the asymptotic distribution of max_{\omega \in \Omega_n} || for a wide class of redundant frames (\phi_\omega^n: \omega \in \Omega_n}. Based on our theoretical results we give a rationale for universal extreme value thresholding techniques yielding asymptotically sharp confidence regions and smoothness estimates corresponding to prescribed significance levels. The results cover many frames used in imaging and signal recovery applications, such as redundant wavelet systems, curvelet frames, or unions of bases. We show that `generically' a standard Gumbel law results as it is known from the case of orthonormal wavelet bases. However, for specific highly redundant frames other limiting laws may occur. We indeed verify that the translation invariant wavelet transform shows a different asymptotic behaviour.Comment: [Content: 39 pages, 4 figures] Note that in this version 4 we have slightely changed the title of the paper and we have rewritten parts of the introduction. Except for corrected typos the other parts of the paper are the same as the original versions

2-D Prony-Huang Transform: A New Tool for 2-D Spectral Analysis

This work proposes an extension of the 1-D Hilbert Huang transform for the analysis of images. The proposed method consists in (i) adaptively decomposing an image into oscillating parts called intrinsic mode functions (IMFs) using a mode decomposition procedure, and (ii) providing a local spectral analysis of the obtained IMFs in order to get the local amplitudes, frequencies, and orientations. For the decomposition step, we propose two robust 2-D mode decompositions based on non-smooth convex optimization: a "Genuine 2-D" approach, that constrains the local extrema of the IMFs, and a "Pseudo 2-D" approach, which constrains separately the extrema of lines, columns, and diagonals. The spectral analysis step is based on Prony annihilation property that is applied on small square patches of the IMFs. The resulting 2-D Prony-Huang transform is validated on simulated and real data.Comment: 24 pages, 7 figure