English

Graph Fourier Transform: A Stable Approximation

Signal Processing 2020-12-02 v2

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 AA that, in GSP, generalizes the Discrete Signal Processing (DSP) shift operator z1z^{-1}. The (right) eigenvectors of the shift AA (graph spectral components) diagonalize AA and lead to a graph Fourier basis FF that provides a graph spectral representation of the graph signal. The inverse of the (matrix of the) graph Fourier basis FF is the Graph Fourier transform (GFT), F1F^{-1}. Often, including in real world examples, this diagonalization is numerically unstable. This paper develops an approach to compute an accurate approximation to FF and F1F^{-1}, 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 FF per iteration. Likewise, we can apply SGFA to AHA^H and, hence, approximate the stable left eigenvectors for the graph shift AA and directly compute the GFT. We evaluate empirically the quality of SGFA by applying it to graph shifts AA 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.

Keywords

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