Table of Contents
Fetching ...

Learning Time-Varying Graphs from Incomplete Graph Signals

Chuansen Peng, Xiaojing Shen

TL;DR

This work tackles the problem of learning time-varying graph topologies from incomplete graph signals by formulating a single non-convex optimization that jointly recovers a sequence of Laplacians and imputes missing signals. A structural fusion regularizer enforces shared topology patterns across graphs, while an ADMM algorithm with closed-form updates yields scalable, provably convergent solutions to the joint problem. The authors establish convergence to stationary points and derive non-asymptotic error bounds that quantify how sampling, sparsity, signal smoothness, and cross-graph similarity affect estimation accuracy. Empirical results on synthetic and real spatio-temporal data demonstrate faster convergence and superior joint graph learning and signal recovery compared to state-of-the-art baselines, highlighting strong practical impact for dynamic networks under partial observability.

Abstract

This paper tackles the challenging problem of jointly inferring time-varying network topologies and imputing missing data from partially observed graph signals. We propose a unified non-convex optimization framework to simultaneously recover a sequence of graph Laplacian matrices while reconstructing the unobserved signal entries. Unlike conventional decoupled methods, our integrated approach facilitates a bidirectional flow of information between the graph and signal domains, yielding superior robustness, particularly in high missing-data regimes. To capture realistic network dynamics, we introduce a fused-lasso type regularizer on the sequence of Laplacians. This penalty promotes temporal smoothness by penalizing large successive changes, thereby preventing spurious variations induced by noise while still permitting gradual topological evolution. For solving the joint optimization problem, we develop an efficient Alternating Direction Method of Multipliers (ADMM) algorithm, which leverages the problem's structure to yield closed-form solutions for both the graph and signal subproblems. This design ensures scalability to large-scale networks and long time horizons. On the theoretical front, despite the inherent non-convexity, we establish a convergence guarantee, proving that the proposed ADMM scheme converges to a stationary point. Furthermore, we derive non-asymptotic statistical guarantees, providing high-probability error bounds for the graph estimator as a function of sample size, signal smoothness, and the intrinsic temporal variability of the graph. Extensive numerical experiments validate the approach, demonstrating that it significantly outperforms state-of-the-art baselines in both convergence speed and the joint accuracy of graph learning and signal recovery.

Learning Time-Varying Graphs from Incomplete Graph Signals

TL;DR

This work tackles the problem of learning time-varying graph topologies from incomplete graph signals by formulating a single non-convex optimization that jointly recovers a sequence of Laplacians and imputes missing signals. A structural fusion regularizer enforces shared topology patterns across graphs, while an ADMM algorithm with closed-form updates yields scalable, provably convergent solutions to the joint problem. The authors establish convergence to stationary points and derive non-asymptotic error bounds that quantify how sampling, sparsity, signal smoothness, and cross-graph similarity affect estimation accuracy. Empirical results on synthetic and real spatio-temporal data demonstrate faster convergence and superior joint graph learning and signal recovery compared to state-of-the-art baselines, highlighting strong practical impact for dynamic networks under partial observability.

Abstract

This paper tackles the challenging problem of jointly inferring time-varying network topologies and imputing missing data from partially observed graph signals. We propose a unified non-convex optimization framework to simultaneously recover a sequence of graph Laplacian matrices while reconstructing the unobserved signal entries. Unlike conventional decoupled methods, our integrated approach facilitates a bidirectional flow of information between the graph and signal domains, yielding superior robustness, particularly in high missing-data regimes. To capture realistic network dynamics, we introduce a fused-lasso type regularizer on the sequence of Laplacians. This penalty promotes temporal smoothness by penalizing large successive changes, thereby preventing spurious variations induced by noise while still permitting gradual topological evolution. For solving the joint optimization problem, we develop an efficient Alternating Direction Method of Multipliers (ADMM) algorithm, which leverages the problem's structure to yield closed-form solutions for both the graph and signal subproblems. This design ensures scalability to large-scale networks and long time horizons. On the theoretical front, despite the inherent non-convexity, we establish a convergence guarantee, proving that the proposed ADMM scheme converges to a stationary point. Furthermore, we derive non-asymptotic statistical guarantees, providing high-probability error bounds for the graph estimator as a function of sample size, signal smoothness, and the intrinsic temporal variability of the graph. Extensive numerical experiments validate the approach, demonstrating that it significantly outperforms state-of-the-art baselines in both convergence speed and the joint accuracy of graph learning and signal recovery.
Paper Structure (14 sections, 11 theorems, 135 equations, 12 figures, 2 tables, 1 algorithm)

This paper contains 14 sections, 11 theorems, 135 equations, 12 figures, 2 tables, 1 algorithm.

Key Result

Proposition 1

The objective function with respect to the matrix variable $\mathbf{X}$ is given by $f(\mathbf{X})=\frac{1}{n}\cdot\frac{1}{\sigma^2}\mathrm{tr}((\mathbf{Y}_M-\mathbf{M}\odot\mathbf{X})(\mathbf{Y}_M-\mathbf{M}\odot\mathbf{X})^T)+\frac{n_k}{n}\cdot\mathrm{tr}(\mathbf{F}^T\mathbf{X}\mathbf{X}^T\mathbf where the constituent matrix $\mathbf{T}$ and vector $\mathbf{e}$ are defined as

Figures (12)

  • Figure 1: The F-score and the relative error of the estimated Laplacian matrix $\mathbf{L}$ are reported for a synthetic data model across varying noise levels $\sigma$, with the sampling ratio fixed at $\mathrm{SR}=0.8$.
  • Figure 2: The F-score and relative error of the estimated Laplacian matrix $\mathbf{L}$ are reported for a synthetic data model across varying sampling ratio $\mathrm{SR}$, with fixed noise level $\sigma=0.1$.
  • Figure 3: NMSE and SNR curves illustrating reconstruction performance of the data matrix $\mathbf{X}$ in the synthetic model, shown for various noise levels $\sigma$ with the sampling ratio fixed at $\mathrm{SR}=0.8$.
  • Figure 4: NMSE and SNR curves illustrating reconstruction performance of the data matrix $\mathbf{X}$ in the synthetic model, shown for various sampling ratio $\mathrm{SR}$ with the noise level fixed at $\sigma=0.1$.
  • Figure 5: The comparison of computation time for different $N$.
  • ...and 7 more figures

Theorems & Definitions (16)

  • Remark 1
  • Proposition 1
  • Theorem 1
  • Corollary 1
  • Theorem 2
  • Remark 2
  • Remark 3
  • Remark 4
  • Remark 5
  • Lemma 1
  • ...and 6 more