A sequential least squares algorithm for ARMAX dynamic network identification

Identification of dynamic networks in prediction error setting often requires the solution of a non-convex optimization problem, which can be difficult to solve especially for large-scale systems. Focusing on ARMAX models of dynamic networks, we instead employ a method based on a sequence of least-squares steps. For single-input single-output models, we show that the method is equivalent to the recently developed Weighted Null Space Fitting, and, drawing from the analysis of that method, we conjecture that the proposed method is both consistent as well as asymptotically efficient under suitable assumptions. Simulations indicate that the sequential least squares estimates can be of high quality even for short data sets.