Analysis and extraction of strongly frequency modulated signals have been a challenging problem for Adaptive Data Analysis methods; e.g., Empirical Mode Decomposition (EMD) [13]. In fact, many of the Newtonian dynamical systems, including conservative mechanical systems, are sources of signals with low to strong levels of frequency modulation. Analysis of such signals is an important issue in system identification problems. In this paper, we present a novel method to accurately extract Intrawave Signals. This method is a descendant of Sparse Time-Frequency Representation (STFR) methods [8, 7]. We will present numerical examples to show the performance of this new algorithm. Theoretical analysis of convergence of the algorithm is also presented as a support for the method. We will show that the algorithm is stable to noise perturbation, as well.

We propose an extension of the computational fluid mechanics approach to the Monge-Kantorovich mass transfer problem, which was developed by Benamou-Brenier. Our extension allows optimal transfer of unnormalized and unequal masses. We obtain a one-parameter family of simple modifications of the formulation in [4]. This leads us to a new Monge-Ampere type equation and a new Kantorovich duality formula. These can be solved efficiently by, for example, the Chambolle-Pock primal-dual algorithm. This solution to the extended mass transfer problem gives us a simple metric for computing the distance between two unnormalized densities. The L1 version of this metric was shown in [23] (which is a precursor of our work here) to have desirable properties.

No abstract uploaded!

We propose an efficient eigensolver for computing densely distributed spectra of the two-dimensional transmission eigenvalue problem (TEP), which is derived from Maxwell's equations with Tellegen media and the transverse magnetic mode. The governing equations, when discretized by the standard piecewise linear finite element method, give rise to a large-scale quadratic eigenvalue problem (QEP). Our numerical simulation shows that half of the positive eigenvalues of the QEP are densely distributed in some interval near the origin. The quadratic JacobiDavidson method with a so-called non-equivalence deflation technique is proposed to compute the dense spectrum of the QEP. Extensive numerical simulations show that our proposed method processes the convergence efficiently, even when it needs to compute more than 5000 desired eigenpairs. Numerical results also illustrate that the computed

A class of upwind flux splitting methods in the Euler equations of compressible flow is considered in this paper. Using the property that Euler flux F(U) is a homogeneous function of degree one in U, we reformulate the splitting fluxes with F^{+}=A^{+} U, F^{-}=A^{-} U, and the corresponding matrices are either symmetric or symmetrizable and keep only non-negative and non-positive eigenvalues. That leads to the conclusion that the first order schemes are positive in the sense of Lax-Liu [18], which implies that it is L^2-stable in some suitable sense. Moreover, the second order scheme is a stable perturbation of the first order scheme, so that the positivity of the second order schemes is also established, under a CFL-like condition. In addition, these splitting methods preserve the positivity of density and energy.