9,866 research outputs found
On the Estimation of Nonrandom Signal Coefficients from Jittered Samples
This paper examines the problem of estimating the parameters of a bandlimited
signal from samples corrupted by random jitter (timing noise) and additive iid
Gaussian noise, where the signal lies in the span of a finite basis. For the
presented classical estimation problem, the Cramer-Rao lower bound (CRB) is
computed, and an Expectation-Maximization (EM) algorithm approximating the
maximum likelihood (ML) estimator is developed. Simulations are performed to
study the convergence properties of the EM algorithm and compare the
performance both against the CRB and a basic linear estimator. These
simulations demonstrate that by post-processing the jittered samples with the
proposed EM algorithm, greater jitter can be tolerated, potentially reducing
on-chip ADC power consumption substantially.Comment: 11 pages, 8 figure
Fast computation of effective diffusivities using a semi-analytical solution of the homogenization boundary value problem for block locally-isotropic heterogeneous media
Direct numerical simulation of diffusion through heterogeneous media can be
difficult due to the computational cost of resolving fine-scale
heterogeneities. One method to overcome this difficulty is to homogenize the
model by replacing the spatially-varying fine-scale diffusivity with an
effective diffusivity calculated from the solution of an appropriate boundary
value problem. In this paper, we present a new semi-analytical method for
solving this boundary value problem and computing the effective diffusivity for
pixellated, locally-isotropic, heterogeneous media. We compare our new solution
method to a standard finite volume method and show that equivalent accuracy can
be achieved in less computational time for several standard test cases. We also
demonstrate how the new solution method can be applied to complex heterogeneous
geometries represented by a grid of blocks. These results indicate that our new
semi-analytical method has the potential to significantly speed up simulations
of diffusion in heterogeneous media.Comment: 29 pages, 4 figures, 5 table
On the convergence rates of Gauss and Clenshaw-Curtis quadrature for functions of limited regularity
We study the optimal general rate of convergence of the n-point quadrature
rules of Gauss and Clenshaw-Curtis when applied to functions of limited
regularity: if the Chebyshev coefficients decay at a rate O(n^{-s-1}) for some
s > 0, Clenshaw-Curtis and Gauss quadrature inherit exactly this rate. The
proof (for Gauss, if 0 < s < 2, there is numerical evidence only) is based on
work of Curtis, Johnson, Riess, and Rabinowitz from the early 1970s and on a
refined estimate for Gauss quadrature applied to Chebyshev polynomials due to
Petras (1995). The convergence rate of both quadrature rules is up to one power
of n better than polynomial best approximation; hence, the classical proof
strategy that bounds the error of a quadrature rule with positive weights by
polynomial best approximation is doomed to fail in establishing the optimal
rate.Comment: 7 pages, the figure of the revision has an unsymmetric example, to
appear in SIAM J. Numer. Ana
On the numerical calculation of the roots of special functions satisfying second order ordinary differential equations
We describe a method for calculating the roots of special functions
satisfying second order linear ordinary differential equations. It exploits the
recent observation that the solutions of a large class of such equations can be
represented via nonoscillatory phase functions, even in the high-frequency
regime. Our algorithm achieves near machine precision accuracy and the time
required to compute one root of a solution is independent of the frequency of
oscillations of that solution. Moreover, despite its great generality, our
approach is competitive with specialized, state-of-the-art methods for the
construction of Gaussian quadrature rules of large orders when it used in such
a capacity. The performance of the scheme is illustrated with several numerical
experiments and a Fortran implementation of our algorithm is available at the
author's website
- …