Hankel and Toeplitz Rank-1 Decomposition of Arbitrary Matrices with Applications to Signal Direction-of-Arrival Estimation
Pith reviewed 2026-05-07 13:26 UTC · model grok-4.3
The pith
Optimal rank-1 Hankel and Toeplitz approximations under L2 and L1 norms produce maximum-likelihood DoA estimators for few-shot scenarios.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We consider the problems of computing the optimal rank-1 Hankel and Toeplitz-structured approximation of arbitrary matrices under L2 and L1-norm error. Such problems arise naturally in engineered systems, including the basic few-shot signal Direction-of-Arrival (DoA) estimation problem that is of importance to modern autonomous systems applications. We develop accurate and computationally efficient structured matrix decomposition algorithms for both formulations and then derive analytically grounded small-sample-support DoA estimators for practical sensing system deployments. The resulting estimators under the L2 and L1 norms are formally shown to be maximum-likelihood optimal under white G
What carries the argument
Rank-1 Hankel and Toeplitz structured matrix decomposition under L2 and L1 norms, which extracts the signal component while respecting the structure induced by the array geometry or time series.
If this is right
- DoA estimators become maximum-likelihood optimal under matching noise distributions.
- Performance improves in few-shot regimes with limited samples.
- Applicable to both Gaussian and heavy-tailed Laplace noise models.
- Computationally efficient algorithms allow real-time use in autonomous systems.
- Validated on both simulated and real-world data for reliability.
Where Pith is reading between the lines
- The decomposition technique could be adapted for other structured approximations in array signal processing.
- It separates the signal structure from the noise statistics, allowing flexible noise modeling.
- In practice, this might lead to better integration with adaptive sensing systems that adjust to observed noise.
- Extensions to higher-rank structures or multi-dimensional arrays could follow similar principles.
Load-bearing premise
The observed matrix is a noisy version of an exact rank-1 matrix with Hankel or Toeplitz structure generated by the signal model.
What would settle it
Running the proposed estimators on data from a known DoA scenario with added white Gaussian noise and comparing their accuracy to the theoretical maximum-likelihood bound or to unstructured estimators would show if they achieve the claimed optimality and superiority in small-sample cases.
Figures
read the original abstract
We consider the problems of computing the optimal rank-$1$ Hankel and Toeplitz-structured approximation of arbitrary matrices under $L_2$ and $L_1$-norm error. Such problems arise naturally in engineered systems, including the basic few-shot signal Direction-of-Arrival (DoA) estimation problem that is of importance to modern autonomous systems applications. We develop accurate and computationally efficient structured matrix decomposition algorithms for both formulations and then derive analytically grounded small-sample-support DoA estimators for practical sensing system deployments. The resulting estimators under the $L_2$ and $L_1$ norms are formally shown to be maximum-likelihood optimal under white Gaussian and Laplace noise, respectively. The estimators are further validated through extensive simulation studies and real-world data experiments in few-shot DoA inference.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops algorithms for computing the optimal rank-1 Hankel and Toeplitz structured approximations to arbitrary matrices under L2 and L1 norms. It applies these to derive small-sample-support DoA estimators and analytically shows that the L2-based estimator is maximum-likelihood optimal under white Gaussian noise while the L1-based estimator is optimal under Laplace noise. The claims are supported by algorithmic descriptions, theoretical derivations, simulations, and real-world experiments.
Significance. If the ML optimality claims hold under the physical sensing model, the work would provide a significant contribution by delivering parameter-free, theoretically grounded DoA estimators suitable for few-shot regimes in autonomous systems. The structured decomposition algorithms may also have broader utility in signal processing applications involving Hankel or Toeplitz matrices. The paper merits credit for combining analytical derivations with extensive empirical validation on both simulated and real data.
major comments (1)
- [DoA application and ML optimality derivation] § on DoA signal model and ML derivation (the section deriving optimality from the likelihood): The claim that the L2 (resp. L1) rank-1 Hankel/Toeplitz approximant is ML optimal under white Gaussian (Laplace) noise assumes iid noise across the observed matrix entries. However, the motivating DoA model forms the matrix via Hankel/Toeplitz embedding of the short vector y = a(θ)s + n where n is white vector noise; overlapping entries therefore contain identical noise samples and the effective noise on the matrix is correlated. This mismatch means the formal optimality proof applies to an observation model that does not match the physical sensing model used to motivate the estimators. This is load-bearing for the central claim in the abstract.
minor comments (2)
- [Abstract] The abstract asserts that the estimators are 'formally shown' to be ML optimal but does not cite the specific theorem, equation, or subsection containing the derivation, which reduces readability.
- [Problem formulation] Notation for the embedding operators that construct the Hankel and Toeplitz matrices from the vector y should be introduced with an explicit definition and example in the problem formulation section to aid clarity.
Simulated Author's Rebuttal
We thank the referee for the careful reading and for highlighting this important distinction between the matrix-level noise model used in the optimality derivation and the vector-level white noise model that motivates the DoA application. We address the comment below.
read point-by-point responses
-
Referee: [DoA application and ML optimality derivation] § on DoA signal model and ML derivation (the section deriving optimality from the likelihood): The claim that the L2 (resp. L1) rank-1 Hankel/Toeplitz approximant is ML optimal under white Gaussian (Laplace) noise assumes iid noise across the observed matrix entries. However, the motivating DoA model forms the matrix via Hankel/Toeplitz embedding of the short vector y = a(θ)s + n where n is white vector noise; overlapping entries therefore contain identical noise samples and the effective noise on the matrix is correlated. This mismatch means the formal optimality proof applies to an observation model that does not match the physical sensing model used to motivate the estimators. This is load-bearing for the central claim in the abstract.
Authors: We agree that the formal ML optimality derivation in the paper is performed under the assumption of i.i.d. noise across the entries of the observed matrix. In the DoA section, the observed matrix is obtained by Hankel/Toeplitz embedding of a short vector y whose noise component is white (i.i.d. across vector entries). Because of the repeated entries in the embedding, the induced noise on the matrix entries is correlated and the Frobenius (or element-wise L1) distance on the matrix is not equivalent to the standard Euclidean (or L1) distance on the original vector. Consequently, the rank-1 structured approximant is exactly maximum-likelihood only for the matrix observation model with i.i.d. entries; it is not strictly ML for the vector white-noise model that physically motivates the DoA problem. We will revise the manuscript as follows: (i) explicitly state the noise model under which the optimality holds, (ii) qualify the abstract claim to indicate that the estimators are ML-optimal under the matrix i.i.d. model and are motivated by (but not exactly optimal for) the embedded vector model, and (iii) add a short discussion of the relationship between the two models together with a note that the matrix-norm formulation remains a principled and parameter-free estimator for few-shot DoA. These changes directly address the load-bearing concern for the central claim. revision: yes
Circularity Check
No circularity: ML optimality derived from standard likelihood under explicit iid matrix noise model
full rationale
The paper derives the L2 and L1 rank-1 Hankel/Toeplitz approximants analytically and states that they are ML optimal under white Gaussian and Laplace noise, respectively. This follows directly from the standard log-likelihood for iid entries (Gaussian or Laplace) plus a rank-1 structured signal; the derivation does not reduce to a fitted parameter renamed as prediction, nor to any self-citation chain, nor to a self-definitional loop. The DoA application simply instantiates the same structured-matrix model; any mismatch between vector noise and induced matrix-entry correlations is an external modeling assumption, not a circular reduction inside the paper's equations.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The underlying signal model results in a rank-1 Hankel or Toeplitz matrix structure.
- domain assumption Noise is white (independent and identically distributed) Gaussian for L2 and Laplace for L1.
Reference graph
Works this paper leans on
-
[1]
Markovsky, Low-Rank Approximation: Algorithms, Imple- mentation, Applications
I. Markovsky, Low-Rank Approximation: Algorithms, Imple- mentation, Applications. London, U.K.: Springer, 2011
work page 2011
-
[2]
The approximation of one matrix by another of lower rank,
C. Eckart and G. Young, “The approximation of one matrix by another of lower rank,” Psychometrika, vol. 1, no. 3, pp. 211–218, Sep. 1936
work page 1936
-
[3]
G. H. Golub and C. F. V. Loan, Matrix Computations, 3rd ed. Baltimore, MD, USA: The Johns Hopkins University Press, 1996
work page 1996
-
[4]
Liii. on lines and planes of closest fit to systems of points in space,
K. Pearson, “Liii. on lines and planes of closest fit to systems of points in space,” The London, Edinburgh, and Dublin philosophical magazine and journal of science, vol. 2, no. 11, pp. 559–572, Nov. 1901
work page 1901
-
[5]
A global geometric framework for nonlinear dimensionality reduction,
J. B. Tenenbaum, V. d. Silva, and J. C. Langford, “A global geometric framework for nonlinear dimensionality reduction,” Science, vol. 290, no. 5500, pp. 2319–2323, Dec. 2000
work page 2000
-
[6]
Laplacian eigenmaps for dimension- ality reduction and data representation,
M. Belkin and P. Niyogi, “Laplacian eigenmaps for dimension- ality reduction and data representation,” Neural computation, vol. 15, no. 6, pp. 1373–1396, Jun. 2003
work page 2003
-
[7]
H. L. Van Trees, Optimum array processing: Part IV of detec- tion, estimation, and modulation theory. John Wiley & Sons, Apr. 2002
work page 2002
-
[8]
On the L1-norm approximation of a matrix by another of lower rank,
N. Tsagkarakis, P. P. Markopoulos, and D. A. Pados, “On the L1-norm approximation of a matrix by another of lower rank,” in Proc. 15th IEEE Int. Conf. Mach. Learn. Appl. (ICMLA), Anaheim, CA, USA, Dec. 2016, pp. 768–773
work page 2016
-
[9]
Q. Ke and T. Kanade, “Robust L1-norm factorization in the presence of outliers and missing data by alternative convex programming,” in Proc. IEEE Comput. Soc. Conf. Comput. Vis. Pattern Recognit. (CVPR), San Diego, CA, USA, Jun. 2005, pp. 739–746
work page 2005
-
[10]
A pure L1-norm principal component analysis,
J. P. Brooks, J. H. Dulá, and E. L. Boone, “A pure L1-norm principal component analysis,” Comput. Statist. Data Anal., vol. 61, pp. 83–98, May 2013
work page 2013
-
[11]
Parameter estimation techniques: A tutorial with application to conic fitting,
Z. Zhang, “Parameter estimation techniques: A tutorial with application to conic fitting,” Image Vis. Comput., vol. 15, no. 1, pp. 59–76, Jan. 1997
work page 1997
-
[12]
Y. Ma, S. Soatto, J. Košecká, and S. Sastry, An Invitation to 3-D Vision: From Images to Geometric Models. New York, NY, USA: Springer, Nov. 2003
work page 2003
-
[13]
Shape and motion from image streams: A factorization method,
C. Tomasi and T. Kanade, “Shape and motion from image streams: A factorization method,” Proc. Natl. Acad. Sci., vol. 90, no. 21, pp. 9795–9802, Nov. 1993
work page 1993
-
[14]
Hankel matrix rank minimization with applications to system identification and realization,
M. Fazel, T. K. Pong, D. Sun, and P. Tseng, “Hankel matrix rank minimization with applications to system identification and realization,” SIAM J. Matrix Anal. Appl., vol. 34, no. 3, pp. 946–977, Jul. 2013
work page 2013
-
[15]
M. R. Laskar and A. K. Dutta, “A proposed quantum frame- work for low-complexity quantum simulation and spectrum es- timation of Hankel-patterned systems,” IEEE Trans. Quantum Eng., vol. 4, pp. 1–18, Nov. 2023
work page 2023
-
[16]
W. Sun, Y. Zhou, J. Xiang, B. Chen, and W. Feng, “Hankel matrix-based condition monitoring of rolling element bearings: An enhanced framework for time-series analysis,” IEEE Trans. Instrum. Meas., vol. 70, pp. 1–10, Feb. 2021
work page 2021
-
[17]
System approxi- mation via restructured hankel matrix,
R. S. Sengar, K. Chatterjee, and J. Singh, “System approxi- mation via restructured hankel matrix,” Circuits Syst. Signal Process., vol. 40, no. 12, pp. 6354–6370, Jun. 2021. 11
work page 2021
-
[18]
Time-series analysis with small and faulty data: L1-norm decompositions of hankel matrices,
G. I. Orfanidis, D. A. Pados, and G. Sklivanitis, “Time-series analysis with small and faulty data: L1-norm decompositions of hankel matrices,” in Proc. SPIE Big Data IV: Learning, Analytics, and Applications, vol. 12097, Orlando, FL, USA, May 2022, pp. 97–104
work page 2022
-
[19]
Single-sample direction-of- arrival estimation by Hankel-matrix decompositions,
G. I. Orfanidis, D. A. Pados, G. Sklivanitis, E. S. Bentley, J. Suprenant, and M. J. Medley, “Single-sample direction-of- arrival estimation by Hankel-matrix decompositions,” in Proc. 56th Asilomar Conf. Signals, Syst., Comput., Pacific Grove, CA, USA, Nov. 2022, pp. 1026–1030
work page 2022
-
[20]
S. Mazokha, S. Naderi, G. I. Orfanidis, G. Sklivanitis, D. A. Pados, and J. O. Hallstrom, “Single-sample direction-of-arrival estimation for fast and robust 3d localization with real measure- ments from a massive MIMO system,” in Proc. IEEE Int. Conf. Acoust., Speech, Signal Process. (ICASSP), Rhodes Island, Greece, Jun. 2023, pp. 1–5
work page 2023
-
[21]
N. Golyandina, V. Nekrutkin, and A. A. Zhigljavsky, Analysis of Time Series Structure: SSA and Related Techniques. Boca Raton, FL, USA: Chapman & Hall, Jun. 2001
work page 2001
-
[22]
Singular spectrum analysis based on L1-norm,
M. Kalantari, M. Yarmohammadi, and H. Hassani, “Singular spectrum analysis based on L1-norm,” Fluct. Noise Lett., vol. 15, no. 1, p. 1650009, Mar. 2016
work page 2016
-
[23]
Signal enhancement— Acomposite property mapping algorithm,
J. A. Cadzow, “Signal enhancement— Acomposite property mapping algorithm,” IEEE Trans. Acoust., Speech, Signal Pro- cess., vol. 36, no. 1, pp. 49–62, Jan. 1988
work page 1988
-
[24]
On low-rank Hankel matrix denois- ing,
M. Yin and R. S. Smith, “On low-rank Hankel matrix denois- ing,” IF AC-PapersOnLine, vol. 54, no. 7, pp. 198–203, 2021
work page 2021
-
[25]
H. Knirsch, M. Petz, and G. Plonka, “Optimal rank-1 Hankel approximation of matrices: Frobenius norm and spectral norm and Cadzow’s algorithm,” Linear Algebra and its Appl., vol. 629, pp. 1–39, Nov. 2021
work page 2021
-
[26]
Cadzow’s basic algorithm, alternating projections and singular spectrum analysis,
J. Gillard, “Cadzow’s basic algorithm, alternating projections and singular spectrum analysis,” Statist. and its Inter., vol. 3, no. 3, pp. 335–343, Jan. 2010
work page 2010
-
[27]
Sur le point pour lequel la somme des distances de n points donnés est minimum,
E. Weiszfeld, “Sur le point pour lequel la somme des distances de n points donnés est minimum,” Tohoku Math. J., vol. 43, pp. 355–386, 1937
work page 1937
-
[28]
Two decades of array signal processing research: The parametric approach,
H. Krim and M. Viberg, “Two decades of array signal processing research: The parametric approach,” IEEE Signal Process. Mag., vol. 13, no. 4, pp. 67–94, Jul. 1996
work page 1996
-
[29]
L. C. Godara, “Application of antenna arrays to mobile com- munications. II. beam-forming and direction-of-arrival consid- erations,” Proc. IEEE, vol. 85, no. 8, pp. 1195–1245, Aug. 1997
work page 1997
-
[30]
Direction- of-arrival analysis of airborne ice depth sounder data,
U. Nielsen, J.-B. Yan, S. Gogineni, and J. Dall, “Direction- of-arrival analysis of airborne ice depth sounder data,” IEEE Trans. Geosci. Remote Sens., vol. 55, no. 4, pp. 2239–2249, Apr. 2017
work page 2017
-
[31]
Coherent DOA estimation in sea surface observation with direction-finding HF radar,
J. Zhao, Y. Tian, B. Wen, and Z. Tian, “Coherent DOA estimation in sea surface observation with direction-finding HF radar,” IEEE Trans. Geosci. Remote Sens., vol. 59, no. 8, pp. 6651–6661, Aug. 2021
work page 2021
-
[32]
DOA estimation of underwater acoustic signals based on deep learning,
P. Li and Y. Tian, “DOA estimation of underwater acoustic signals based on deep learning,” in Proc. 2nd Int. Seminar Artif. Intell., Netw. Inf. Technol. (AINIT), Shanghai, China, Oct. 2021, pp. 221–225
work page 2021
-
[33]
S. K. Joshi, S. V. Baumgartner, A. B. C. da Silva, and G. Krieger, “Direction-of-arrival angle and position estimation for extended targets using multichannel airborne radar data,” IEEE Geosci. Remote Sens. Lett., vol. 19, pp. 1–5, Feb. 2022
work page 2022
-
[34]
Coverage analysis of UA Vs in millimeter wave networks: A stochastic geometry approach,
M. Boschiero, M. Giordani, M. Polese, and M. Zorzi, “Coverage analysis of UA Vs in millimeter wave networks: A stochastic geometry approach,” in Proc. Int. Wireless Commun. Mobile Comput. Conf. (IWCMC), Limassol, Cyprus, Jun. 2020, pp. 351–357
work page 2020
-
[35]
Millimeter wave wireless assisted robot navigation with link state classification,
M. Yin, A. K. Veldanda, A. Trivedi, J. Zhang, K. Pfeiffer, Y. Hu, S. Garg, E. Erkip, L. Righetti, and S. Rangan, “Millimeter wave wireless assisted robot navigation with link state classification,” IEEE Open J. Commun. Soc., vol. 3, pp. 493–507, Mar. 2022
work page 2022
-
[36]
Path planning under uncertainty to localize mmWave sources,
K. Pfeiffer, Y. Jia, M. Yin, A. K. Veldanda, Y. Hu, A. Trivedi, J. Zhang, S. Garg, E. Erkip, S. Rangan, and L. Righetti, “Path planning under uncertainty to localize mmWave sources,” in Proc. IEEE Int. Conf. Robot. Autom. (ICRA), London, U.K., Jun. 2023, pp. 3461–3467
work page 2023
-
[37]
Radio SLAM: A review on radio- based simultaneous localization and mapping,
B. Amjad, Q. Z. Ahmed, P. I. Lazaridis, M. Hafeez, F. A. Khan, and Z. D. Zaharis, “Radio SLAM: A review on radio- based simultaneous localization and mapping,” IEEE Access, vol. 11, pp. 9260–9278, 2023
work page 2023
-
[38]
Multiple emitter location and signal parameter estimation,
R. Schmidt, “Multiple emitter location and signal parameter estimation,” IEEE Trans. Antennas Propag., vol. 34, no. 3, pp. 276–280, Mar. 1986
work page 1986
-
[39]
ESPRIT—estimation of signal parame- ters via rotational invariance techniques,
R. Roy and T. Kailath, “ESPRIT—estimation of signal parame- ters via rotational invariance techniques,” IEEE Trans. Acoust., Speech, Signal Process., vol. 37, no. 7, pp. 984–995, Jul. 1989
work page 1989
-
[40]
Sensor array processing based on subspace fitting,
M. Viberg and B. Ottersten, “Sensor array processing based on subspace fitting,” IEEE Trans. Signal Process., vol. 39, no. 5, pp. 1110–1121, May 1991
work page 1991
-
[41]
Subspace direction finding with an auxiliary-vector basis,
R. Grover, D. A. Pados, and M. J. Medley, “Subspace direction finding with an auxiliary-vector basis,” IEEE Trans. Signal Process., vol. 55, no. 2, pp. 758–763, Jan. 2007
work page 2007
-
[42]
Tensor-based near-field localization using massive antenna arrays,
I. Podkurkov, G. Seidl, L. Khamidullina, A. Nadeev, and M. Haardt, “Tensor-based near-field localization using massive antenna arrays,” IEEE Trans. Signal Process., vol. 69, pp. 5830– 5845, Aug. 2021
work page 2021
-
[43]
Matrix pencil method for estimating parameters of exponentially damped/undamped sinusoids in noise,
Y. Hua and T. Sarkar, “Matrix pencil method for estimating parameters of exponentially damped/undamped sinusoids in noise,” IEEE Trans. Acoust., Speech, Signal Process., vol. 38, no. 5, pp. 814–824, May 1990
work page 1990
-
[44]
Elimination of the effects of mutual coupling in adaptive thin wire antennas,
R. S. Adve, “Elimination of the effects of mutual coupling in adaptive thin wire antennas,” Ph.D. dissertation, Syracuse University, Syracuse, NY, USA, 1996
work page 1996
-
[45]
MUSIC for single-snapshot spectral estimation: Stability and super-resolution,
W. Liao and A. Fannjiang, “MUSIC for single-snapshot spectral estimation: Stability and super-resolution,” Appl. Comput. Harmon. Anal., vol. 40, no. 1, pp. 33–67, Jan. 2016
work page 2016
-
[46]
On spatial smoothing for direction-of-arrival estimation of coherent signals,
T.-J. Shan, M. Wax, and T. Kailath, “On spatial smoothing for direction-of-arrival estimation of coherent signals,” IEEE Trans. Acoust., Speech, Signal Process., vol. 33, no. 4, pp. 806–811, Aug. 1985
work page 1985
-
[47]
Single snapshot DOA estimation based on spatial smoothing MUSIC and CNN,
C. Liu, W. Feng, H. Li, and H. Zhu, “Single snapshot DOA estimation based on spatial smoothing MUSIC and CNN,” in Proc. IEEE Int. Conf. Signal Process., Commun. Comput. (ICSPCC), Xi’an, China, Aug. 2021, pp. 1–5
work page 2021
-
[48]
On single snapshot direction-of-arrival estimation,
C. Degen, “On single snapshot direction-of-arrival estimation,” in Proc. IEEE Int. Conf. Wireless Space Extreme Environ. (WiSEE), Montréal, QC, Canada, Oct. 2017, pp. 92–97
work page 2017
-
[49]
S. L. M. Jr., Digital Spectral Analysis. Mineola, NY, USA: Dover Publications, Mar. 2019
work page 2019
-
[50]
Experimen- tal evaluation of AoA estimation for UA V to massive MIMO,
T. Rice, D. Pandey, D. Ramirez, and E. Knightly, “Experimen- tal evaluation of AoA estimation for UA V to massive MIMO,” in Proc. IEEE Mil. Commun. Conf. (MILCOM), Boston, MA, USA, Oct. 2023, pp. 839–844
work page 2023
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.