Agile Inexact Methods for Spectral Projector-Based Graph Fourier Transforms

We propose an inexact method for the graph Fourier transform of a graph signal, as defined by the signal decomposition over the Jordan subspaces of the graph adjacency matrix. This method projects the signal over the generalized eigenspaces of the adjacency matrix, which accelerates the transform computation over large, sparse, and directed adjacency matrices. The trade-off between execution time and fidelity to the original graph structure is discussed. In addition, properties such as a generalized Parseval's identity and total variation ordering of the generalized eigenspaces are discussed. The method is applied to 2010-2013 NYC taxi trip data to identify traffic hotspots on the Manhattan grid. Our results show that identical highly expressed geolocations can be identified with the inexact method and the method based on eigenvector projections, while reducing computation time by a factor of 26,000 and reducing energy dispersal among the spectral components corresponding to the multiple zero eigenvalue.

[1]  Axel Ruhe Properties of a matrix with a very ill-conditioned eigenproblem , 1970 .

[2]  José M. F. Moura,et al.  Graph Equivalence Classes for Spectral Projector-Based Graph Fourier Transforms , 2017, ArXiv.

[3]  Pascal Frossard,et al.  The emerging field of signal processing on graphs: Extending high-dimensional data analysis to networks and other irregular domains , 2012, IEEE Signal Processing Magazine.

[4]  Zhonggang Zeng Regularization and Matrix Computation in Numerical Polynomial Algebra , 2009 .

[5]  Vivek K Goyal,et al.  Foundations of Signal Processing , 2014 .

[6]  Michael G. Rabbat,et al.  Approximating signals supported on graphs , 2012, 2012 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[7]  F. R. Gantmakher The Theory of Matrices , 1984 .

[8]  J. H. Wilkinson The algebraic eigenvalue problem , 1966 .

[9]  Markus Püschel,et al.  Algebraic Signal Processing Theory: Foundation and 1-D Time , 2008, IEEE Transactions on Signal Processing.

[10]  José M. F. Moura,et al.  Discrete Signal Processing on Graphs , 2012, IEEE Transactions on Signal Processing.

[11]  Axel Ruhe An algorithm for numerical determination of the structure of a general matrix , 1970 .

[12]  P. P. Vaidyanathan,et al.  Extending Classical Multirate Signal Processing Theory to Graphs—Part I: Fundamentals , 2017, IEEE Transactions on Signal Processing.

[13]  Peter Lancaster,et al.  The theory of matrices , 1969 .

[14]  Santiago Segarra,et al.  Sampling of Graph Signals With Successive Local Aggregations , 2015, IEEE Transactions on Signal Processing.

[15]  Franz Franchetti,et al.  Discrete fourier transform on multicore , 2009, IEEE Signal Processing Magazine.

[16]  José M. F. Moura,et al.  Big Data Analysis with Signal Processing on Graphs: Representation and processing of massive data sets with irregular structure , 2014, IEEE Signal Processing Magazine.

[17]  L. Trefethen,et al.  Spectra and Pseudospectra , 2020 .

[18]  Jon Wilkening,et al.  An algorithm for computing Jordan chains and inverting analytic matrix functions , 2007 .

[19]  Louis L. Scharf,et al.  Signal processing applications of oblique projection operators , 1994, IEEE Trans. Signal Process..

[20]  José M. F. Moura,et al.  Algebraic Signal Processing Theory: 1-D Space , 2008, IEEE Transactions on Signal Processing.

[21]  José M. F. Moura,et al.  Spectral Projector-Based Graph Fourier Transforms , 2017, IEEE Journal of Selected Topics in Signal Processing.

[22]  Franz Franchetti,et al.  Big data computation of taxi movement in New York City , 2016, 2016 IEEE International Conference on Big Data (Big Data).

[23]  Santiago Segarra,et al.  Reconstruction of Graph Signals Through Percolation from Seeding Nodes , 2015, IEEE Transactions on Signal Processing.

[24]  José M. F. Moura,et al.  Discrete Signal Processing on Graphs: Frequency Analysis , 2013, IEEE Transactions on Signal Processing.

[25]  José M. F. Moura,et al.  Signal Recovery on Graphs: Variation Minimization , 2014, IEEE Transactions on Signal Processing.

[26]  Sunil K. Narang,et al.  Perfect Reconstruction Two-Channel Wavelet Filter Banks for Graph Structured Data , 2011, IEEE Transactions on Signal Processing.

[27]  Zhonggang Zeng,et al.  Sensitivity and Computation of a Defective Eigenvalue , 2016, SIAM J. Matrix Anal. Appl..

[28]  Charles R. Johnson,et al.  Matrix analysis , 1985, Statistical Inference for Engineers and Data Scientists.

[29]  Gene H. Golub,et al.  Ill-conditioned eigensystems and the computation of the Jordan canonical form , 1975, Milestones in Matrix Computation.

[30]  Michael L. Overton,et al.  An Efficient Algorithm for Computing the Generalized Null Space Decomposition , 2015, SIAM J. Matrix Anal. Appl..

[31]  Robert H. Halstead,et al.  Matrix Computations , 2011, Encyclopedia of Parallel Computing.

[32]  Lada A. Adamic,et al.  The political blogosphere and the 2004 U.S. election: divided they blog , 2005, LinkKDD '05.

[33]  W. Kahan Conserving Confluence Curbs Ill-Condition , 1972 .