Papers › A Data-Driven McMillan Degree Lower Bound
A Data-Driven McMillan Degree Lower Bound
Jeffrey M. Hokanson
The archive published only this paper's code-link row. Authors, date and abstract are from arXiv's metadata (CC0), read from the Kaggle arXiv metadata snapshot of 2026-09-12 where its title matched the archive's; the title is the archive's.
Given measurements of a linear time-invariant system, the McMillan degree is the dimension of the smallest such system that reproduces these observed dynamics. Using impulse response measurements where the system has been started in some (unknown) state and then allowed to evolve freely, a classical result by Ho and Kalman reveals the McMillan degree as the rank of a Hankel matrix built from these measurements. However, if measurements are contaminated by noise, this Hankel matrix will almost surely be full rank. Hence practitioners often estimate the rank of this matrix---and thus the McMillan degree---by manually setting a threshold between the large singular values that correspond to the non-zero singular values of the noise-free Hankel matrix and the small singular values that are pertubations of the zero singular values. Here we introduce a probabilistic upper bound on the perturbation of the singular values of this Hankel matrix when measurements are corrupted by additive Gaussian noise, and hence provide guidance on setting the threshold to obtain a lower bound on the McMillan degree. This result is powered by a new, probabilistic bound on the 2-norm of a random Hankel matrix with normally distributed entries. Unlike existing results for random Hankel matrices, this bound features no unknown constants and, moreover, is within a small factor of the empirically observed bound when entries are independent and identically distributed. This bound on the McMillan degree provides an inexpensive alternative to more general model order selection techniques such as the Akaike Information Criteria (AIC).
Code
Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.
Code Syntology ran Syntology
Not run by Syntology. Nothing on this page verifies that the listed code works.
Results from the paper archive 2025-07-28
No leaderboard rows for this paper in the archive.
Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections