Model-Free $H_\infty$ Nash Equilibrium Seeking in Network Games: Double-Homotopy Policy Iteration and Output Regulation

Khanh Thai1,*

1School of Electrical & Electronic Engineering, Hanoi University of Science and Technology (HUST), Vietnam

*Corresponding author

Paper Coming soon Code Coming soon

Working draft. Not yet public and not yet submitted for peer review — the paper and code will be posted here once ready.

No initial stabilizing policy needed $\hat\gamma$ within $\sim\!1\%$ of $\gamma^\star$, from data $30\%$ smaller disturbance gain than a data-driven LMI baseline Small-gain certified for graphs with cycles
Convergence of P, K, L across six continuation levels of the mu-homotopy for player 1, plus the final level compared to the model-based optimum
Player 1's whole learning run: six continuation levels in $\mu=\gamma^{-2}$, each converging quadratically (blue: successive difference; orange: distance to that level's saddle point). Panel (d) shows the final level closing in on the model-based optimum — using data only, with no initial stabilizing gain and no system identification.

Overview

Nash equilibrium seeking, learned from data, robust to what you didn't model.

Distributed Nash equilibrium (NE) seeking for players with dynamics has recently been cast as a cooperative output regulation problem, where the pseudo-gradient of the game plays the role of the regulated error. We make this fully data-driven and add what that reformulation was missing: an $H_\infty$-optimal attenuation of disturbances the model does not capture, a network-level guarantee that survives communication graphs with cycles, and a plain account of what happens when the data themselves are not exact.

Model-Free $H_\infty$ NE Seeking

Each player learns, from its own local input–state data, the saddle-point policy of a zero-sum game on its internal-model-augmented dynamics — exact rejection of exosystem disturbances plus a certified $\mathcal L_2$ attenuation of unmodelled ones, with no system identification, no initial stabilizing gain, and no communication during learning.

Double-Homotopy Policy Iteration

A state-shift homotopy builds a stabilizing gain from scratch; a second continuation in $\mu=\gamma^{-2}$ then carries the policy from the LQR solution to the $H_\infty$ solution, with every iterate certified stabilizing directly from data. As a by-product, it returns a data-based estimate of the optimal attenuation level $\gamma^\star$.

Beyond Acyclic Graphs

An explicit disturbance-to-pseudo-gradient $\mathcal L_2$ gain, automatic for acyclic graphs and holding under a small-gain condition for graphs with cycles — both certified from each player's own local data, with no model and no central computation.

An Honest Account of Inexact Data

An exact perturbation bound for the data solve, with an amplification factor computable from the same data. The estimator is errors-in-variables biased, so collecting more data does not remove the error — verified with $16\times$ more data. Averaging the learning identity over shifted start times buys about a factor $30$ in tolerable noise.

Method

Two homotopies: one for stability, one for optimal attenuation.

Policy iteration is Newton's method on a Riccati equation and converges quadratically — but it must start from a stabilizing policy, which is itself model information. For zero-sum games a second difficulty appears: the simultaneous update is Newton's method on a Riccati equation with an indefinite quadratic term, so a stabilizing pair alone is not enough for convergence. We remove both obstacles with two nested homotopies, using only input–state data collected once, off-policy.

Phase I: the stabilizer homotopy driving the state shift to zero while the certified stability bound holds. Phase II: the attenuation homotopy taking gamma from infinity (LQR) down to the target level in five accepted continuation steps
Double-homotopy policy iteration. (a) Phase I (state-shift homotopy): the shift $\bar\beta-\sum\alpha$ is driven to zero while a data-based certified bound on the closed-loop spectral abscissa holds, producing a stabilizing gain with no model. (b) Phase II (attenuation homotopy): a continuation in $\mu=\gamma^{-2}$ carries the policy from $\gamma=\infty$ (the LQR solution) down to the target level in a handful of accepted steps, each certified stabilizing from data before it is accepted. Click to open full size.

Theorem (convergence of the attenuation homotopy)

Under a standard excitation condition on the collected data, the continuation in $\mu=\gamma^{-2}$ accepts only finitely many steps, every policy pair it uses is certified stabilizing, and the PI at the target level converges quadratically to the true saddle point — for every $\gamma>\gamma^\star$, with no initial stabilizing policy required. As a corollary, running the continuation past the true optimum and taking the last accepted level gives a data-based estimate $\hat\gamma\to\gamma^\star$.

A network-level counterpart holds too: if every player's learned gain is $\gamma_i>\gamma_i^\star$, the closed-loop network is exponentially stable and $\|e\|_{\mathcal L_2}\le\Gamma\|d\|_{\mathcal L_2}$ for an explicit $\Gamma$ — automatically on acyclic graphs, and under a small-gain condition, certified from local data alone, on graphs with cycles.

Why the second homotopy is needed

Simultaneous zero-sum policy iteration started from a genuinely stabilizing pair need not stay stabilizing. On our benchmark, starting from the LQR gain at $\gamma=1.1\gamma^\star$ or $1.3\gamma^\star$ drives the closed loop unstable already at the first update, and the iteration diverges. This is a known initialization difficulty for the game Riccati equation's indefinite quadratic term; the attenuation homotopy sidesteps it by following the solution path in $\gamma$ instead of jumping straight to the target level.

Simultaneous zero-sum policy iteration started directly at the target level from the stabilizing LQR pair, without the mu-homotopy: for gamma at most 1.3 times gamma star the iterates leave the stabilizing set
Started directly at the target level, without the continuation: for $\gamma\le1.3\gamma^\star$ the iterates leave the stabilizing set within a few updates.

Results

Five heterogeneous robots, learned $H_\infty$ gains, and a fair fight with a robust baseline.

Five mobile robots with different friction coefficients play a network game on a communication graph, each with an unknown linear model, an unmodelled $\mathcal L_2$ disturbance channel, and a known exosystem to reject. Every player learns its gain from $600$ data intervals with no initial stabilizing policy, and the resulting network is compared against the data-driven LMI design of Wang et al. (2026), which stabilizes but is not optimal and needs an acyclic graph for its stability guarantee.

Learning results per player (model-based $\gamma^\star$, $K^\star$ shown only for validation — never used by the algorithm):

Player $\gamma^\star$ $\hat\gamma$ (data) $\mu$-steps PI iters $\|K-K^\star\|/\|K^\star\|$
12.1772.197529$9.5\times10^{-4}$
22.4062.436530$2.7\times10^{-3}$
32.3662.391529$1.7\times10^{-3}$
42.4332.462427$2.9\times10^{-3}$
52.3892.411530$3.3\times10^{-3}$

The data-based estimate $\hat\gamma$ is within about $1\%$ of the model-based $\gamma^\star$ for every player, and the learned gain $K$ matches the true GARE solution $K^\star$ to $3\times10^{-3}$ relative error — all without ever solving a model-based Riccati equation.

Network comparison with the data-driven LMI design (acyclic graph):

Proposed LMI baseline
$\|T_{d\to e}\|_\infty$3.625.15
$\|T_{d\to u}\|_\infty$3.053.86
empirical $\|e\|_2/\|d\|_2$ (rest)1.772.47
$\int_0^{60}\|e\|^2dt$422703
$\int_0^{60}\|u\|^2dt$20801380
final $\|y-y^\star\|_\infty$$1.3\times10^{-7}$$2.0\times10^{-7}$

The learned $H_\infty$ gains give a $30\%$ smaller disturbance-to-error gain and a $30\%$ smaller disturbance-to-input gain than the LMI baseline. The one line where the baseline wins is total control energy over the full scenario — the learned gains are larger, which costs more effort during the initial transient.

Citation

Cite this work

@misc{thai2026hinfne,
  title={Model-Free $H_\infty$ Nash Equilibrium Seeking in Network Games via
         Double-Homotopy Policy Iteration and Output Regulation},
  author={Thai, Khanh},
  year={2026},
  note={Working draft}
}