1,188 research outputs found
Adaptive transient solution of nonuniform multiconductor transmission lines using wavelets
AbstractâThis paper presents a highly adaptive algorithm for the transient simulation of nonuniform interconnects loaded with arbitrary nonlinear and dynamic terminations. The discretization of the governing equations is obtained through a weak formula-tion using biorthogonal wavelet bases as trial and test functions. It is shown how the multiresolution properties of wavelets lead to very sparse approximations of the voltages and currents in typical transient analyzes. A simple yet effective timeâspace adaptive al-gorithm capable of selecting the minimal number of unknowns at each time iteration is described. Numerical results show the high degree of adaptivity of the proposed scheme. Index TermsâElectromagnetic (EM) transient analysis, multi-conductor transmission lines (TLs), wavelet transforms. I
Wavelet-Fourier CORSING techniques for multi-dimensional advection-diffusion-reaction equations
We present and analyze a novel wavelet-Fourier technique for the numerical
treatment of multidimensional advection-diffusion-reaction equations based on
the CORSING (COmpRessed SolvING) paradigm. Combining the Petrov-Galerkin
technique with the compressed sensing approach, the proposed method is able to
approximate the largest coefficients of the solution with respect to a
biorthogonal wavelet basis. Namely, we assemble a compressed discretization
based on randomized subsampling of the Fourier test space and we employ sparse
recovery techniques to approximate the solution to the PDE. In this paper, we
provide the first rigorous recovery error bounds and effective recipes for the
implementation of the CORSING technique in the multi-dimensional setting. Our
theoretical analysis relies on new estimates for the local a-coherence, which
measures interferences between wavelet and Fourier basis functions with respect
to the metric induced by the PDE operator. The stability and robustness of the
proposed scheme is shown by numerical illustrations in the one-, two-, and
three-dimensional case
Recovering edges in ill-posed inverse problems: optimality of curvelet frames
We consider a model problem of recovering a function from noisy Radon data. The function to be recovered is assumed smooth apart from a discontinuity along a curve, that is, an edge. We use the continuum white-noise model, with noise level .
Traditional linear methods for solving such inverse problems behave poorly in the presence of edges. Qualitatively, the reconstructions are blurred near the edges; quantitatively, they give in our model mean squared errors (MSEs) that tend to zero with noise level only as as . A recent innovation--nonlinear shrinkage in the wavelet domain--visually improves edge sharpness and improves MSE convergence to . However, as we show here, this rate is not optimal.
In fact, essentially optimal performance is obtained by deploying the recently-introduced tight frames of curvelets in this setting. Curvelets are smooth, highly anisotropic elements ideally suited for detecting and synthesizing curved edges. To deploy them in the Radon setting, we construct a curvelet-based biorthogonal decomposition of the Radon operator and build "curvelet shrinkage" estimators based on thresholding of the noisy curvelet coefficients. In effect, the estimator detects edges at certain locations and orientations in the Radon domain and automatically synthesizes edges at corresponding locations and directions in the original domain.
We prove that the curvelet shrinkage can be tuned so that the estimator will attain, within logarithmic factors, the MSE as noise level . This rate of convergence holds uniformly over a class of functions which are except for discontinuities along curves, and (except for log terms) is the minimax rate for that class. Our approach is an instance of a general strategy which should apply in other inverse problems; we sketch a deconvolution example
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
Self-similar prior and wavelet bases for hidden incompressible turbulent motion
This work is concerned with the ill-posed inverse problem of estimating
turbulent flows from the observation of an image sequence. From a Bayesian
perspective, a divergence-free isotropic fractional Brownian motion (fBm) is
chosen as a prior model for instantaneous turbulent velocity fields. This
self-similar prior characterizes accurately second-order statistics of velocity
fields in incompressible isotropic turbulence. Nevertheless, the associated
maximum a posteriori involves a fractional Laplacian operator which is delicate
to implement in practice. To deal with this issue, we propose to decompose the
divergent-free fBm on well-chosen wavelet bases. As a first alternative, we
propose to design wavelets as whitening filters. We show that these filters are
fractional Laplacian wavelets composed with the Leray projector. As a second
alternative, we use a divergence-free wavelet basis, which takes implicitly
into account the incompressibility constraint arising from physics. Although
the latter decomposition involves correlated wavelet coefficients, we are able
to handle this dependence in practice. Based on these two wavelet
decompositions, we finally provide effective and efficient algorithms to
approach the maximum a posteriori. An intensive numerical evaluation proves the
relevance of the proposed wavelet-based self-similar priors.Comment: SIAM Journal on Imaging Sciences, 201
Image Compression by Wavelet Transform.
Digital images are widely used in computer applications. Uncompressed digital images require considerable storage capacity and transmission bandwidth. Efficient image compression solutions are becoming more critical with the recent growth of data intensive, multimedia-based web applications.
This thesis studies image compression with wavelet transforms. As a necessary background, the basic concepts of graphical image storage and currently used compression algorithms are discussed. The mathematical properties of several types of wavelets, including Haar, Daubechies, and biorthogonal spline wavelets are covered and the Enbedded Zerotree Wavelet (EZW) coding algorithm is introduced. The last part of the thesis analyzes the compression results to compare the wavelet types
- âŠ