Distributed State Estimation for Discrete-Time Systems With Unknown Inputs: An Optimization Approach
Pith reviewed 2026-05-10 15:15 UTC · model grok-4.3
The pith
A distributed framework lets sensor nodes estimate the full state of systems with unknown inputs by locally reconstructing reachable parts and optimizing their fusion.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The DUIO framework achieves bounded estimation error by having each node compute its maximal locally reconstructible state from local inputs and outputs, then fusing the partial estimates through a distributed optimization algorithm whose error bound depends on the number of communication iterations per time step and the strong-convexity constant determined by the system parameters; a normalization step mitigates curvature anisotropy from ill-conditioned geometry.
What carries the argument
The Distributed Unknown Input Observer (DUIO) that combines local maximal-state reconstruction with a distributed optimization step that fuses the partial estimates while preserving strong convexity.
If this is right
- The explicit error bound decreases monotonically with more communication iterations per time step.
- The normalization step inside the optimizer reduces the effect of poor system conditioning on convergence rate.
- The framework applies directly to large-scale discrete-time networks whose sensors and unknown inputs are spatially dispersed.
Where Pith is reading between the lines
- The same local-reconstruction-plus-optimization pattern could be tested on continuous-time or nonlinear systems without changing the core convexity argument.
- If the strong-convexity constant can be computed from local data alone, the method could run with minimal global knowledge.
- The bounded-error guarantee supplies a concrete starting point for designing distributed controllers that use the produced state estimates.
Load-bearing premise
Each node can independently compute the largest state segment reconstructible from its own local inputs and outputs, and the global fusion problem remains strongly convex with parameters set by the system.
What would settle it
Run the algorithm on a system where the local reconstructibility condition fails for at least one node and observe whether the global estimation error grows unbounded even as communication iterations increase.
Figures
read the original abstract
This paper proposes a novel Distributed Unknown Input Observer (DUIO) framework for state estimation in large-scale systems subject to local unknown inputs. We consider systems where outputs are measured by a network of spatially distributed sensors and inputs are introduced through multiple dispersed channels. In this framework, each local node utilizes only its local input and output measurements to estimate the maximal locally reconstructible state. Subsequently, nodes collaboratively reconstruct the whole system state via a distributed optimization algorithm that fuses these partial estimates. We provide a rigorous analysis showing that the estimation error is bounded, with the error bound explicitly dependent on the number of communication iterations per time step and strongly convexity constant determined by the system parameters. Furthermore, to counteract curvature anisotropy induced by poor conditioned system geometry, we embed a normalization step into the distributed optimization procedure. Simulation results demonstrate the effectiveness of the proposed framework and the performance improvements yielded by the normalization procedure.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a Distributed Unknown Input Observer (DUIO) framework for discrete-time large-scale systems with unknown inputs. Each local node uses only its local input and output measurements to compute the maximal locally reconstructible state subspace. Nodes then collaboratively reconstruct the full state by solving a distributed quadratic program that fuses the partial estimates. The manuscript claims a rigorous bounded-error analysis in which the estimation error is bounded by an expression that depends explicitly on the number of communication iterations per time step and on a strong-convexity constant μ determined solely by the system matrices (A,B,C,D). A normalization step is embedded in the distributed optimizer to mitigate curvature anisotropy, and simulation results are provided to illustrate effectiveness.
Significance. If the claimed error bound and its explicit dependence on iteration count and system-determined μ can be rigorously established, the work would supply a useful optimization-based route to distributed state estimation under unknown inputs, particularly for sensor networks where global observability is absent but local reconstructibility is exploitable. The normalization procedure for handling ill-conditioned geometry is a practical contribution. The manuscript does not supply machine-checked proofs or fully reproducible code, but the simulation validation and the attempt at an iteration-dependent bound are positive features.
major comments (2)
- [Abstract / error-bound derivation] Abstract and the section deriving the error bound: the headline claim requires that the fused quadratic program remain strongly convex with parameter μ > 0 fixed by the system matrices alone and independent of the communication graph. No rank, observability, or complementarity condition on the local reconstructible subspaces is stated that would guarantee the global Hessian is positive definite with such a μ. Without this, the derivation of the iteration-dependent error bound collapses.
- [Distributed optimization procedure] The normalization step (introduced to counteract curvature anisotropy): it is unclear whether this rescaling preserves the claimed strong-convexity constant μ or merely alters the effective conditioning while leaving the error-bound derivation intact. The manuscript must show that the normalized Hessian still satisfies the same μ lower bound used in the convergence analysis.
minor comments (1)
- [Abstract] Abstract: 'strongly convexity constant' is a grammatical error and should read 'strong-convexity constant'.
Simulated Author's Rebuttal
We thank the referee for the constructive and detailed comments. We address each major point below and will incorporate revisions to strengthen the theoretical claims.
read point-by-point responses
-
Referee: [Abstract / error-bound derivation] Abstract and the section deriving the error bound: the headline claim requires that the fused quadratic program remain strongly convex with parameter μ > 0 fixed by the system matrices alone and independent of the communication graph. No rank, observability, or complementarity condition on the local reconstructible subspaces is stated that would guarantee the global Hessian is positive definite with such a μ. Without this, the derivation of the iteration-dependent error bound collapses.
Authors: We agree that the conditions guaranteeing strong convexity of the fused quadratic program with μ determined solely by the system matrices (and independent of the communication graph) were not stated explicitly. In the revised manuscript we will introduce a new assumption on the complementarity of the local reconstructible subspaces (ensuring their direct sum spans the full state space) together with a supporting lemma that derives the explicit lower bound μ from the system matrices A, B, C, D alone. The iteration-dependent error bound derivation will then be restated under this assumption, restoring its validity. revision: yes
-
Referee: [Distributed optimization procedure] The normalization step (introduced to counteract curvature anisotropy): it is unclear whether this rescaling preserves the claimed strong-convexity constant μ or merely alters the effective conditioning while leaving the error-bound derivation intact. The manuscript must show that the normalized Hessian still satisfies the same μ lower bound used in the convergence analysis.
Authors: The normalization is a positive-definite diagonal scaling chosen to improve conditioning. We will add a short lemma in the revision proving that this scaling leaves the strong-convexity parameter μ unchanged: the minimal eigenvalue of the normalized Hessian is bounded below by the same μ derived from the original system matrices. Consequently the convergence analysis and the iteration-dependent error bound remain exactly as stated. revision: yes
Circularity Check
No circularity: error bound derived from stated assumptions on local subspaces and strong convexity
full rationale
The abstract and derivation outline a sequential process: local nodes compute maximal reconstructible states from their own (u,y) data, then solve a distributed QP whose Hessian is asserted strongly convex with μ fixed by (A,B,C,D). The error bound is then obtained explicitly in terms of iteration count and this μ. No equation reduces the bound to a fitted parameter, self-defined quantity, or prior self-citation chain; the convexity parameter is an input assumption rather than an output of the target result. The normalization step is an algorithmic addition, not a redefinition of the bound. The chain is therefore self-contained and non-circular.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Discrete-time linear system dynamics with unknown inputs and distributed sensor measurements
- domain assumption Distributed optimization converges with error bound governed by iteration count and strong-convexity constant
Reference graph
Works this paper leans on
-
[1]
, " * write output.state after.block = add.period write newline
ENTRY address author booktitle chapter doi edition editor eid howpublished institution journal key month note number organization pages publisher school series title type url volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 'mid.sent...
-
[2]
" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in capitalize " " * FUNCT...
-
[3]
Aybat, N.S., Wang, Z., Lin, T., and Ma, S. (2017). Distributed linearized alternating direction method of multipliers for composite convex consensus optimization. IEEE Transactions on Automatic Control, 63(1), 5--20
work page 2017
-
[4]
Bastianello, N., Carli, R., Schenato, L., and Todescato, M. (2020). Asynchronous distributed optimization over lossy networks via relaxed admm: Stability and linear convergence. IEEE Transactions on Automatic Control, 66(6), 2620--2635
work page 2020
-
[5]
Bastianello, N., Deplano, D., Franceschelli, M., and Johansson, K.H. (2024). Robust online learning over networks. IEEE Transactions on Automatic Control, 70(2), 933--946
work page 2024
-
[6]
Bhattacharyya, S. (2003). Observer design for linear systems with unknown inputs. IEEE transactions on Automatic Control, 23(3), 483--484
work page 2003
-
[7]
Cao, G. and Wang, J. (2023). Distributed unknown input observer. IEEE Transactions on Automatic Control, 68(12), 8244--8251
work page 2023
-
[8]
Disar \`o , G., Fattore, G., and Valcher, M.E. (2025). Distributed state estimation for discrete-time lti systems in the presence of unknown inputs. IEEE Transactions on Automatic Control
work page 2025
-
[9]
Farina, M., Ferrari-Trecate, G., and Scattolini, R. (2010). Distributed moving horizon estimation for linear constrained systems. IEEE Transactions on Automatic Control, 55(11), 2462--2475
work page 2010
-
[10]
Fioravanti, C., Makridis, E., Oliva, G., Vrakopoulou, M., and Charalambous, T. (2024). Distributed estimation and control for lti systems under finite-time agreement. IEEE Transactions on Automatic Control, 69(11), 7909--7916
work page 2024
-
[11]
Horn, R.A. and Johnson, C.R. (2012). Matrix analysis. Cambridge university press
work page 2012
-
[12]
Kim, T., Lee, D., and Shim, H. (2023). Decentralized design and plug-and-play distributed control for linear multichannel systems. IEEE Transactions on Automatic Control, 69(5), 2807--2822
work page 2023
-
[13]
Li, P. (2012). Geometric analysis, volume 134. Cambridge University Press
work page 2012
-
[14]
Liu, F., Wang, L., Fullmer, D., and Morse, A.S. (2023). Distributed feedback control of multichannel linear systems. IEEE Transactions on Automatic Control, 68(12), 7164--7178
work page 2023
-
[15]
Rego, F.F., Pascoal, A.M., Aguiar, A.P., and Jones, C.N. (2019). Distributed state estimation for discrete-time linear time invariant systems: A survey. Annual Reviews in Control, 48, 36--56
work page 2019
-
[16]
Wang, L., Sundaram, S., and LeGrand, K. (2024). Distributed optimization-based state estimation of nonlinear dynamical systems. In 2024 American Control Conference (ACC), 2816--2821. IEEE
work page 2024
-
[17]
Wu, Y., Su, H., and Zeng, Z. (2025). Almost sure distributed robust observer with unknown inputs. IEEE Transactions on Control of Network Systems
work page 2025
-
[18]
Yang, G., Barboni, A., Rezaee, H., and Parisini, T. (2022). State estimation using a network of distributed observers with unknown inputs. Automatica, 146, 110631
work page 2022
-
[19]
Yang, G., Ren, X., Bastianello, N., and Parisini, T. (2025). State estimation using a network of observers: A distributed optimization approach. In 2025 European Control Conference (ECC), 3151--3157. IEEE
work page 2025
-
[20]
Zhang, J., Zhao, X., Zheng, G., Zhu, F., and Dinh, T.N. (2025). On distributed prescribed-time unknown input observers. IEEE Transactions on Automatic Control
work page 2025
- [21]
- [22]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.