Assigning Stationary Distributions to Sparse Stochastic Matrices
Abstract
The target stationary distribution problem (TSDP) is the following: given an irreducible stochastic matrix and a target stationary distribution , construct a minimum norm perturbation, , such that is also stochastic and has the prescribed target stationary distribution, . In this paper, we revisit the TSDP under a constraint on the support of , that is, on the set of non-zero entries of . This is particularly meaningful in practice since one cannot typically modify all entries of . We first show how to construct a feasible solution that has essentially the same support as the matrix . Then we show how to compute globally optimal and sparse solutions using the component-wise norm and linear optimization. We propose an efficient implementation that relies on a column-generation approach which allows us to solve sparse problems of size up to in a few minutes. We illustrate the proposed algorithms with several numerical experiments.
Cite
@article{arxiv.2312.16011,
title = {Assigning Stationary Distributions to Sparse Stochastic Matrices},
author = {Nicolas Gillis and Paul Van Dooren},
journal= {arXiv preprint arXiv:2312.16011},
year = {2025}
}
Comments
29 pages, code available from https://gitlab.com/ngillis/TSDP. In this third version, we have added clarifications, corrections and remarks suggested to us by anonymous reviewers