Graph Fourier Transform: A Stable Approximation
Abstract
In Graph Signal Processing (GSP), data dependencies are represented by a graph whose nodes label the data and the edges capture dependencies among nodes. The graph is represented by a weighted adjacency matrix that, in GSP, generalizes the Discrete Signal Processing (DSP) shift operator . The (right) eigenvectors of the shift (graph spectral components) diagonalize and lead to a graph Fourier basis that provides a graph spectral representation of the graph signal. The inverse of the (matrix of the) graph Fourier basis is the Graph Fourier transform (GFT), . Often, including in real world examples, this diagonalization is numerically unstable. This paper develops an approach to compute an accurate approximation to and , while insuring their numerical stability, by means of solving a non convex optimization problem. To address the non-convexity, we propose an algorithm, the stable graph Fourier basis algorithm (SGFA) that we prove to exponentially increase the accuracy of the approximating per iteration. Likewise, we can apply SGFA to and, hence, approximate the stable left eigenvectors for the graph shift and directly compute the GFT. We evaluate empirically the quality of SGFA by applying it to graph shifts drawn from two real world problems, the 2004 US political blogs graph and the Manhattan road map, carrying out a comprehensive study on tradeoffs between different SGFA parameters. We also confirm our conclusions by applying SGFA on very sparse and very dense directed Erd\H os-R\'enyi graphs.
Cite
@article{arxiv.2001.05042,
title = {Graph Fourier Transform: A Stable Approximation},
author = {João Domingos and José M. F. Moura},
journal= {arXiv preprint arXiv:2001.05042},
year = {2020}
}
Comments
16 pages, 17 figures. Originally submitted in -IEEE Transactions on Signal Processing- on 01-Aug-2019. Resubmitted on 12-Jan-2020. Accept with mandatory minor revisions. Resubmitted again on 30-April-2020