Distributed Algorithm to Solve a System of Linear Equations With Unique or Multiple Solutions From Arbitrary Initializations

Abstract

A discrete-time distributed algorithm to solve a system of linear equations Ax = b is proposed with M-Fejer mappings. The algorithm can find a solution of Ax = b from arbitrary initializations at a geometric rate when Ax = b has either unique or multiple solutions. When Ax = b has a unique solution, the geometric convergence rate of the algorithm is proved by analyzing the mixed norm of homogeneous M-Fejer mappings. Then, when Ax = b has multiple solutions, the geometric convergence rate is proved through orthogonal decompositions of the agents' estimates onto the row space and null space of A, and the relationship between the initializations and the final convergence point is also specified. Quantitative upper bounds of the convergence rates for two special cases are given. Finally, some simulation examples are adopted to illustrate the effectiveness of the proposed algorithm.