2,534 research outputs found
Chebyshev Polynomial Approximation for Distributed Signal Processing
Unions of graph Fourier multipliers are an important class of linear
operators for processing signals defined on graphs. We present a novel method
to efficiently distribute the application of these operators to the
high-dimensional signals collected by sensor networks. The proposed method
features approximations of the graph Fourier multipliers by shifted Chebyshev
polynomials, whose recurrence relations make them readily amenable to
distributed computation. We demonstrate how the proposed method can be used in
a distributed denoising task, and show that the communication requirements of
the method scale gracefully with the size of the network.Comment: 8 pages, 5 figures, to appear in the Proceedings of the IEEE
International Conference on Distributed Computing in Sensor Systems (DCOSS),
June, 2011, Barcelona, Spai
Chebyshev expansion for Impurity Models using Matrix Product States
We improve a recently developed expansion technique for calculating real
frequency spectral functions of any one-dimensional model with short-range
interactions, by postprocessing computed Chebyshev moments with linear
prediction. This can be achieved at virtually no cost and, in sharp contrast to
existing methods based on the dampening of the moments, improves the spectral
resolution rather than lowering it. We validate the method for the exactly
solvable resonating level model and the single impurity Anderson model. It is
capable of resolving sharp Kondo resonances, as well as peaks within the
Hubbard bands when employed as an impurity solver for dynamical mean-field
theory (DMFT). Our method works at zero temperature and allows for arbitrary
discretization of the bath spectrum. It achieves similar precision as the
dynamical density matrix renormalization group (DDMRG), at lower cost. We also
propose an alternative expansion, of 1-exp(-tau H) instead of the usual H,
which opens the possibility of using established methods for the time evolution
of matrix product states to calculate spectral functions directly.Comment: 13 pages, 9 figure
CayleyNets: Graph Convolutional Neural Networks with Complex Rational Spectral Filters
The rise of graph-structured data such as social networks, regulatory
networks, citation graphs, and functional brain networks, in combination with
resounding success of deep learning in various applications, has brought the
interest in generalizing deep learning models to non-Euclidean domains. In this
paper, we introduce a new spectral domain convolutional architecture for deep
learning on graphs. The core ingredient of our model is a new class of
parametric rational complex functions (Cayley polynomials) allowing to
efficiently compute spectral filters on graphs that specialize on frequency
bands of interest. Our model generates rich spectral filters that are localized
in space, scales linearly with the size of the input data for
sparsely-connected graphs, and can handle different constructions of Laplacian
operators. Extensive experimental results show the superior performance of our
approach, in comparison to other spectral domain convolutional architectures,
on spectral image classification, community detection, vertex classification
and matrix completion tasks
Spectral functions and time evolution from the Chebyshev recursion
We link linear prediction of Chebyshev and Fourier expansions to analytic
continuation. We push the resolution in the Chebyshev-based computation of
many-body spectral functions to a much higher precision by deriving a
modified Chebyshev series expansion that allows to reduce the expansion order
by a factor . We show that in a certain limit the Chebyshev
technique becomes equivalent to computing spectral functions via time evolution
and subsequent Fourier transform. This introduces a novel recursive time
evolution algorithm that instead of the group operator only involves
the action of the generator . For quantum impurity problems, we introduce an
adapted discretization scheme for the bath spectral function. We discuss the
relevance of these results for matrix product state (MPS) based DMRG-type
algorithms, and their use within dynamical mean-field theory (DMFT). We present
strong evidence that the Chebyshev recursion extracts less spectral information
from than time evolution algorithms when fixing a given amount of created
entanglement.Comment: 12 pages + 6 pages appendix, 11 figure
Sidescan Sonar Image Enchancement Using a Decomposition Based on Orthogonal Functions. Applications with Chebyshev Polynomials
A method is presented to remove from sidescan sonar images of the seafloor, artifacts that are clearly unrelated to the backscattering properties of the seafloor. A spectral analysis performed on a ping by ping basis proved to be well suited to the problem. The technique relies on a decomposition using Chebyshev polynomials. This stochastic method does not require a priori knowledge of deterministic parameters. It deals with the low spatial frequency components of the image whose wavelengths are not very small compared to the swath width. Applications to sidescan sonar images obtained with the SeaMARC LI system are presented
On the resolution power of Fourier extensions for oscillatory functions
Functions that are smooth but non-periodic on a certain interval possess
Fourier series that lack uniform convergence and suffer from the Gibbs
phenomenon. However, they can be represented accurately by a Fourier series
that is periodic on a larger interval. This is commonly called a Fourier
extension. When constructed in a particular manner, Fourier extensions share
many of the same features of a standard Fourier series. In particular, one can
compute Fourier extensions which converge spectrally fast whenever the function
is smooth, and exponentially fast if the function is analytic, much the same as
the Fourier series of a smooth/analytic and periodic function.
With this in mind, the purpose of this paper is to describe, analyze and
explain the observation that Fourier extensions, much like classical Fourier
series, also have excellent resolution properties for representing oscillatory
functions. The resolution power, or required number of degrees of freedom per
wavelength, depends on a user-controlled parameter and, as we show, it varies
between 2 and \pi. The former value is optimal and is achieved by classical
Fourier series for periodic functions, for example. The latter value is the
resolution power of algebraic polynomial approximations. Thus, Fourier
extensions with an appropriate choice of parameter are eminently suitable for
problems with moderate to high degrees of oscillation.Comment: Revised versio
- …