---
title: Graph Fourier Transform (GFT) Overview
url: https://www.emergentmind.com/topics/graph-fourier-transform-gft
type: topic
---

# Graph Fourier Transform (GFT) Overview

The Graph Fourier Transform (GFT) generalizes the classical discrete Fourier transform (DFT) to signals defined on the nodes of a graph, capturing frequency analysis, filtering, and spectral decomposition in networks with complex topology. In the GFT framework, the underlying graph structure replaces traditional time or spatial regularity, leveraging the eigendecomposition of graph operators (typically Laplacians or adjacency matrices) to define spectra, Fourier modes, and graph frequencies. GFT constructions have been extensively developed for undirected graphs and, more recently, for directed graphs, often relying on various matrix representations, optimization-theoretic formulations, and spectral characterizations to address the intricacies of network directionality, asymmetry, and singularity.

## 1. Fundamental Principles and Classical GFT

The foundational classical GFT is defined for undirected, weighted graphs $\mathcal{G} = (\mathcal{V}, \mathcal{E}, W)$ with $N = |\mathcal{V}|$ nodes, adjacency matrix $W$, and degree matrix $D = \mathrm{diag}(d_1, ..., d_N)$ where $d_i = \sum_j W_{ij}$. The combinatorial Laplacian is $L = D - W$, a symmetric positive semidefinite matrix admitting orthonormal eigendecomposition $L = U \Lambda U^T$, with $U = [u_1, ..., u_N]$ and $\Lambda = \mathrm{diag}(\lambda_1, ..., \lambda_N)$, $0 = \lambda_1 \leq ... \leq \lambda_N$. The GFT of a signal $x \in \mathbb{R}^N$ is $\hat{x} = U^T x$ and the inverse is $x = U \hat{x}$; the eigenvectors $u_i$ are the "graph Fourier modes,” and the eigenvalues $\lambda_i$ are interpreted as "graph frequencies" governing the smoothness of the corresponding mode. Modes with smaller $\lambda_i$ represent slowly varying (low-frequency) patterns, while larger $\lambda_i$ reflect more oscillatory structure [2503.13164].

For directed graphs, the Laplacian becomes asymmetric, complicating the spectral theory and necessitating alternative approaches for defining the GFT.

## 2. GFTs on Directed Graphs: Laplacian, SVD, Jordan, and Polar Approaches

Directed graphs introduce nonsymmetric adjacency and Laplacian matrices. Several seminal GFT construction strategies for digraphs include:

### Directed Laplacian (Jordan-based) GFT

Given a digraph with adjacency $W$ (where $W_{ij}$ reflects the edge $j \to i$), in-degree matrix $D_{\mathrm{in}}$, the directed Laplacian is $L = D_{\mathrm{in}} - W$. As $L$ is typically nondiagonalizable (asymmetric), the graph harmonics are chosen as the (generalized) Jordan eigenvectors $v_i$ of $L$ with associated eigenvalues $\lambda_i$, where $L v_i = \lambda_i v_i$ [1601.03204]. The GFT is then defined as $\hat{f} = V^{-1}f$ for a signal $f$, with the inverse $f = V\hat{f}$, $V$ being the Jordan basis of $L$. Frequency ordering is obtained by considering the total variation $\mathrm{TV}_G(f) = \|L f\|_1$; for normalized harmonics $v_i$, $\mathrm{TV}_G(v_i) = |\lambda_i|$, yielding a smoothness-to-oscillation spectrum.

### SVD-based GFT

The SVD-based approach, applicable to general directed Laplacians or unified graph representation matrices (\textit{UGRM}s), decomposes the (potentially non-normal) matrix $L$ (or $P^{\alpha,k}$ in the UGRM framework) as $L = U \Sigma V^T$. Here, $\Sigma = \mathrm{diag}(\sigma_0, ..., \sigma_{N-1})$ contains nonnegative singular values ordered as graph frequencies, $V$ (resp.\ $U$) are right (left) singular vectors. The GFT is then defined by mixing these bases, e.g. $\mathcal{F}x = [(U^T+V^T)x/2; (U^T-V^T)x/2]$, yielding a unitary transform with Parseval equality even when $L$ is asymmetric [2205.06242, 2510.10532]. 

Such SVD constructions generalize classical GFTs, with convergence to known DFTs for circulant digraphs, and facilitate efficient implementation and parameterized family definitions (UGRM-GFT) which interpolate between degree, Laplacian, and adjacency representations [2510.10532].

### Polar Decomposition GFTs

Alternative digraph GFT designs leverage the polar decomposition of the adjacency matrix $A = QF = PU$, where $Q$, $U$ are unitary, $F$, $P$ are Hermitian positive semi-definite. Three GFTs are derived from the eigendecomposition of $P$, $F$, and $Q$, each associated with well-defined node-domain variation metrics—common-in-link, common-out-link, and shift-residual, respectively. This yields mutually orthogonal transforms, resolving non-orthogonality and numerical instability issues inherent in Jordan-based GFTs [2304.04350].

### Projector-based and Variation-minimizing GFTs

For defective or non-diagonalizable adjacency matrices, the projector-based GFT decomposes signals into spectral projectors onto Jordan subspaces of $A$, admitting generalized Parseval equality and unique, coordinate-free harmonic definitions [1701.02690]. Lovász extension and related convex relaxations are exploited to define orthonormal GFT bases for digraphs by minimizing piecewise-linear (l1-type) directed variations, with tractable optimization (e.g., SOC, PAMAL) yielding numerically stable, cluster-constant transforms [1601.05972].

## 3. GFTs for Product and Multi-dimensional Graphs

Large-scale graph signals often reside on cartesian product graphs (e.g., space × time), requiring multi-dimensional or joint GFTs. For graphs $G_1$, $G_2$ with Laplacians $L_1$, $L_2$, the joint Laplacian is $L_\boxtimes = L_1 \otimes I + I \otimes L_2$, leading to Kronecker-factorized GFTs with spectral indices as frequency tuples $(\lambda_i^{(1)}, \lambda_j^{(2)})$ [1712.07811, 2209.01336]. SVD-based methods allow either full product SVD (global basis, high expressiveness, high cost) or separable Kronecker-GFT (lower complexity, effective for locally coupled signals) [2510.10532]. These frameworks support multidimensional filtering, stationarity, and invariant analysis, and enable efficient denoising and subspace projections for spatiotemporal data [2209.01336].

## 4. Graph Frequencies, Frames, and Basis Design

One challenge in GFT analysis is the nonuniform distribution and possible localization of GFT frequencies/eigenvalues—contrasting with the uniform frequency grid of the DFT. To address biases in spectral sampling and sparse representation, overcomplete frames (denser-graph-frequency GFFs/DGFFs) have been proposed. These augment the usual GFT basis with additional oscillating vectors interpolating intermediate frequencies, constructed via linear or analytic interpolation, yielding improved dispersion, sparser representations, and boosting filtering performance [2503.13164].

Spread-frequency GFTs for directed graphs aim to construct orthonormal bases with directed variation (DV) frequencies as evenly spread as possible over the attainable frequency range. This is achieved via nonconvex dispersion minimization on the Stiefel manifold, or scalable greedy algorithms exploiting supermodularity and matroid constraints, with provable approximation guarantees [1705.10821, 1804.03000].

## 5. Computational and Numerical Considerations

The dense eigendecomposition required by standard GFTs incurs $O(N^2)$ time and memory complexity ($N$: number of nodes). Fast GFTs use symmetry, bipartition, and Haar/butterfly structures to reduce costs to $O(N \log N)$ on graphs with suitable structure (e.g., symmetric lines, cycles, grids) [1907.07875]. On arbitrary graphs, truncated Jacobi algorithms provide approximate fast GFTs controlled by the number of Givens rotations, enabling trade-offs between accuracy, storage, and run-time [1711.00386]. Numerical stability on digraphs is tackled via optimization-based frameworks (e.g., the Stable Graph Fourier Algorithm, SGFA), which guarantee well-conditioned spectral bases through alternating projections and contraction mapping, systematically balancing diagonalization error and minimum singular value [2001.05042].

## 6. GFT in Sampling, Recovery, and Signal Processing Applications

The GFT underpins sampling and recovery theory for graph signals, especially in settings where signals are bandlimited in the graph spectrum. In domains like water distribution networks, where dynamics exhibit low-rank structure rather than strict Laplacian bandlimiting, data-driven GFTs construct bases tailored to observed rank, guaranteeing perfect recovery from a minimal set of node measurements—a threshold often unattainable by classical or compressed sensing-based approaches [1904.03437]. Denser spectral frames and overcomplete GFFs further enhance reconstruction, denoising, and inpainting, while facilitating signal-adaptive filter design [2503.13164].

GFTs directly inform the design and spectral interpretation of graph filters—LSI filters are implemented as polynomials in the graph shift operator (Laplacian, adjacency, UGRM), with the spectral response specified by polynomial evaluation on (generalized) graph frequencies. This framework extends naturally to convolutional architectures in GNNs, enabling spectral kernel learning aligned to graph topology [1601.03204, 2510.10532].

## 7. Extensions and Theoretical Generalizations

Recent work connects finite-graph GFTs to continuum analogues—the graphon Fourier transform—by analyzing the convergence of GFT spectra and modes as graph size grows, illuminating the transferability of spectral analysis and filter design to infinite dense-graph families [1910.10195]. Generalized GFTs decouple the definitions of signal variation and energy, introducing flexible, application-driven orthonormal bases via generalized eigenproblems parameterized by positive-definite matrices reflecting sampling irregularity, energy metrics, or domain-specific structure [1802.10220]. For localized or task-aware analysis, sparse GFTs optimize orthogonality, energy capture, and vertex sparsity via lasso-regularized regression reformulations [1811.08609], providing interpretable, subgraph-focused spectral components.

---

In summary, the Graph Fourier Transform provides a unifying spectral framework for signal processing on graphs, with diverse formulations and extensions handling directed and undirected graphs, multidimensional domains, irregular topologies, and domain-specific requirements. Technical advances in basis design, spectral frame enrichment, fast algorithms, and robust optimization underpin practical applications in sampling, recovery, denoising, and GNNs, while also establishing connections to broader spectral and functional analytic theory. Key milestones include Laplacian and SVD-based GFTs for digraphs [1601.03204, 2205.06242, 2510.10532], spread-frequency frameworks [1705.10821, 1804.03000], numerically stable spectral optimization [2001.05042], and denser-frequency frames [2503.13164], consolidating the GFT’s foundational role in modern graph signal processing.

Source: https://www.emergentmind.com/topics/graph-fourier-transform-gft